Optimized Skip-Stop Metro Line Operation Using Smart Card Data

Skip-stop operation is a low cost approach to improving the efficiency of metro operation and passenger travel experience. This paper proposes a novel method to optimize the skip-stop scheme for bidirectional metro lines so that the average passenger travel time can be minimized. Different from the conventional “A/B” scheme, the proposed Flexible Skip-Stop Scheme (FSSS) can better accommodate spatially and temporally varied passenger demand. A genetic algorithm (GA) based approach is then developed to efficiently search for the optimal solution. A case study is conducted based on a real world bidirectional metro line in Shenzhen, China, using the time-dependent passenger demand extracted from smart card data. It is found that the optimized skip-stop operation is able to reduce the average passenger travel time and transit agencies may benefit from this scheme due to energy and operational cost savings. Analyses are made to evaluate the effects of that fact that certain number of passengers fail to board the right train (due to skip operation). Results show that FSSS always outperforms the all-stop scheme even when most passengers of the skipped OD pairs are confused and cannot get on the right train.


Introduction
The time required for stops at stations is usually a major travel time component for urban transit systems.Cortés et al. [1] discussed different approaches and operation strategies that can reduce total passenger travel time.With the goal of reducing travel time and increasing overall transit operation efficiency, skip-stop operation is a scheme that can have fleets not making all designated stops along a route.Skip-stop operation has been widely applied in bus transit based on the A/B mode.Under the A/B mode, operating fleets are labeled as either type A or type B; stations are also predetermined as type A, type B, or type AB stations.Type A fleets can only stop at type A and type AB stations and type B fleets can only stop at type B and type AB stations [2,3].Eberlein et al. [4] proposed a real-time deadheading strategy to determine the dispatching time.Deadheading vehicles and skip operations were optimized to minimize the total passenger cost.Sun and Hickman [5] formulated the real-time skip-stop problem for bus operation and presented an enumeration method with good computation speed.Yu et al. [6] studied the service reliability of a bus route in the city of Dalian and proposed a heuristic algorithm-based deadheading optimization strategy.Fu et al. [7] presented a new strategy to have all-stop and skip-stop trains operate in pairs.The interests of operators and passengers were balanced in the model by using an exhaustive search procedure.Skip-stop operation is applied in practice by the Transmilenio system in Bogota and Metro Rapid system in Los Angeles for bus operation and has been proved to be highly effective according to the US National Research Council [8] and Leiva et al. [9].
For urban metro systems, the A/B mode-based skip-stop strategy was first developed for the metro systems in Chicago city since 1948 [10,11].It was later implemented in Philadelphia and the City of New York [12].However, due to budget constraints and large train headways caused by operating different types of trains, it switched back to the all-stop mode after several years.A/B mode has also been implemented in the metro system of Santiago, Chile, since 2007, which has demonstrated a significant benefit to passengers and transit operators [10].Lee et al. [13] and Jamili and Aghaee [14] reported that, under skip-stop operation, passengers experienced shorter travel time and operators benefited from the scheme by having reduced fleet size and operational costs.
In addition, skip-stop operation has been commonly used as a recovery strategy in case of special events such as demand surges and major disruptions according to Cao et al. [15].
With respect to the use of skip-stop metro operation, there are some challenges that have not been fully addressed in the existing literature.The first is that since overtaking is prohibited on single-track metro lines, a larger departure headway is needed for the case that the target train travels faster (due to skip operation) than the previous train.This may also reduce the line capacity and cause some difficulty in timetable design.The second challenge is brought by the infrastructure constraints at the terminal stations, where trains need to wait and/or be stored at the turnaround terminal before serving the backward direction.Therefore the safe and efficient operation at the turnaround terminals needs to be properly considered.The third challenge is the precise estimation of dynamic passenger demand pattern, which significantly affects the design of skip-stop scheme.Last but not least, skip-stop operation usually leads to passenger confusion.As a result, certain number of passengers may fail to board the right train.
To better address the first challenge, recent research interest has shifted from the conventional A/B skip-stop mode to Flexible Skip-Stop Schemes (FSSS).A summary and comparison of the existing skip-stop schemes for metro operation is provided in Table 1.Sun and Hickman [5], Zheng et al. [16], Freyss et al. [10], Lee et al. [13], and Cao et al. [15] used different approaches to optimize A/B skip-stop scheme for a single direction urban transit line based on assumed origin-destination (OD) information.Sogin et al. [17] developed a mixed integer program that could be solved by commercial solvers for small scale problems.A genetic algorithm (GA) based solution approach was developed to solve larger scale problems with more stations and trains considered.Niu et al. [19] optimized a skip-stop scheme based on predetermined train stop plans; a mixed integer nonlinear programming model was constructed and solved by a GAMS solver in their paper.Wang et al. [18] developed optimization models that offer real-time operations with Flexible Skip-Stop Schemes; that is, trains are allowed to adjust their skip-stop patterns in real-time to minimize passenger travel time and energy consumption.Jamili and Aghaee [14] built a model to determine the "uncertain skip patterns" which means that trains can flexibly skip or stop at stations along the line so that the train operation speed is optimized.
Another often overlooked challenge is the operation at turnaround terminals.The only work that considers turnaround operation is Wang et al. [18], which treated the turnaround terminal as an intermediate station.The authors believe it is more realistic to consider turnaround operation with departure time optimization and storage capacity.This enables more flexible bidirectional operation that can lead to more efficient line utilization.For the third challenge, recent research has shown that transit smart card data can be utilized to better interpret the temporal and spatial dynamics of passenger demand [20,21].Niu et al. [19] used time-varying passenger demand obtained from the high-speed rail system.For metro operation, researchers (e.g., Sogin et al. [17]; Wang et al. [18]; Jamili and Aghaee [14]) usually assume uniform passenger arrival rates.Little work has been done using dynamic passenger demand from real operation data.With respect to the last challenge, although it is widely recognized as a major drawback of skip-stop operation, the effect of "fail to board the right train" has often been overlooked and little relevant analysis has been reported in literature.
In this paper, we develop a FSSS for bidirectional metro operation during off-peak periods.The problem is formulated as a mixed integer nonlinear programming (MINLP) model that minimizes the average passenger travel time.The proposed scheme is able to consider the train operation at the turnaround terminal with departure time optimization and storage capacity.The departure times at the two terminal stations are constrained by three types of departure intervals, that is, the minimum safety headway, headway that meets the passenger demand, and the turnaround headway.A GAbased solution approach is then developed to efficiently solve the problem.The proposed scheme is tested using a real metro line in Shenzhen, China.The case study considers timedependent demand information inferred from the smart card data of the metro system.Compared to the all-stop scheme, we find that the proposed FSSS can reduce the average passenger travel time.Transit agencies may also benefit from this scheme due to energy and operational cost savings.For the cases that certain percentage of passengers fail to board the right train, results show that FSSS always outperforms the all-stop scheme even when most passengers of the skipped OD pairs are confused and cannot get on the right train.The remainder of this paper is organized as follows.Section 2 presents the FSSS model formulation.Section 3 proposes a GA-based solution approach to efficiently solve the model.The experiment and numerical results are summarized in Section 4, followed by the concluding remarks in Section 5.

FSSS Timetabling
In this section, we first introduce the notations, followed by a brief description of bidirectional train operation problem.We then introduce the proposed FSSS model formulation. bpg : total number of the passengers whose in-vehicle time is reduced by skip-stop operation.

Bidirectional Metro Operation
Problem.Consider a bidirectional metro line.The system is comprised of  tr + 1 trains and   stations, indexed from 0 to  tr and 1 to   , respectively.Train 0 is an all-stop train that departs at time 0; this all-stop train is operated to ensure the system is in operation at the beginning of the optimization period.Station 1 is the departure terminal and station   is the turnaround terminal.The line has only one track in each direction as shown in Figure 1.To avoid confusion, stations are denoted as 1, 2, . . .,   ,   +1, . . ., 2  , in which stations 1 (2  ) indicate departure terminal; stations   and   + 1 represent the turnaround terminal; {1, 2  }, {2, 2  − 1}, . . ., {  ,   + 1} are used to represent the same stations, which share the same dwell time.The section travel times between two adjacent stations are assumed to be the same for both directions.In real practices, there are usually more than one turnaround line at the turnaround terminal where trains are allowed to be stored on the lines.As illustrated in Figure 1, train  − 1 up to train  −  − 1 are allowed to wait at turnaround terminal before train  −  departs from station   + 1.A proof is given in Appendix to show that more turnaround lines (e.g.,  ≥ 1) can lead to better train utilization and therefore increase the line capacity.The solution of the metro operation problem is to make skip/stop decisions and determine the departure time at each station.Thereby a complete operation timetable can be obtained.
The authors notice that there is a debate regarding whether skip-stop is more suitable for off-peak or peak period operation.For example, Niu et al. [19] built mathematical models to solve the rail optimization problem during peak period.As reported in Suh et al. [12], the benefits received in off-peak period (due to skip-stop) are more significant compared to the peak period.The results are largely due to the fact that there are larger passenger demands at most stations during the peak period.The skip operation may lead to excessive waiting time at the skipped stations and it may also cause station congestion.In this work, the authors consider the off-peak operation.The train capacity is not considered as a constraint since our early study [22] found out that trains are far from crowded during the off-peak period.

Assumptions.
The key assumptions used in this paper are summarized and discussed as below.
Assumption 1.The proposed method targets offline timetabling design and the passenger information (e.g., OD) is assumed to be known and can be accurately calibrated using smartcard data collected from their day-to-day travel activities.
Assumption 2. The number of the supplied trains and the storage capacity at the departure terminal are unlimited.Assumption 3.All stations can only accommodate one train at a time for one direction; overtaking is prohibited at any point of the line.This assumption is consistent with the common infrastructure setup at most metro systems in China.Assumption 4. The travel time between two adjacent stations and the dwell time at each station varies along the line but they are assumed to be fixed for all trains.The time loss associated with a stop (due to braking and acceleration) is assumed to be a constant for all trains and all stops.Assumption 5.It is assumed that if train  stops at station  and is about to skip station , certain percentage of passengers from OD pairs  to  will mistakenly board the train (this percentage is referred to as the confusion rate).It is assumed that these passengers shall notice their mistake (e.g., via train radio broadcast) once they get on the train, and they will alight at the next stopping station.These passengers will then wait at the station until being transported by the next available (correct) train to their destination .

Objective Function.
The objective of FSSS is to determine the operation scheme and departure times to minimize the average passenger travel time including the in-vehicle travel time and waiting time.The operation for train  at station  is defined as a binary variable  , .
The total number of onboard passengers  pg , the total passenger waiting time  wait , and the total passenger in-vehicle time  veh can be calculated using (2), in which  , pg ,  , wait , and  , veh are the number of onboard passengers, the passenger waiting time, and the passenger in-vehicle time for train  at station , respectively.
We use  , to denote the departure time of train  at station .The dwell time at station  is   dw , which is determined based on the real metro operation data.The arrival time at station  can be easily acquired by subtracting   dw from  , .Variable    is used to denote the train that stopped at station  prior to the arrival of train , which can be referred to as    = max{  |   < ,  , =    , = 1} according to Niu et al. [19].
The time-dependent passenger OD rate between stations  and  during time period  is denoted as  ,  and the length of each period is .veh can be calculated as in ( 6)-( 8).The objective function is then formulated to minimize the average passenger travel time; see formulation (9).Period Period "g" "g + 1" safety and the capacity of a metro line.Three types of train headways are considered in this study, which are the minimum safety headway  min , the headway that meets the passenger demand   , and the minimum turnaround headway   .Different from the work in Wang et al. [18] who took the turnaround terminal as a normal station, we specifically consider the departure time optimization and storage capacity at the turnaround station.In this subsection, the headway-related constraints are discussed below, followed by the operational constraints.10) and (11) represent that train 0 departures at time 0; after arriving at the turnaround terminal, the new departure time of train 0 is calculated by adding the section travel time, station dwell time, and the turnaround time.For safe operation, the departure interval of two consecutive trains should satisfy formulation ( 12) at all stations.

Minimum Safety Headway. Equations (
Note that skip operation leads to shorter travel time and it may cause some safety concerns.An illustrative example is provided in Figure 3. Suppose train 1 stops at all stations and train 2 skips station 2 and station 3.If train 2 only keeps the minimum safety headway  min at the departure terminal, the minimum safety headway will not suffice at station 3.This situation poses safety issues and requires further intervention control.To avoid this, an extra amount of time Δ  is added to the departure time (of train ) at the departure terminal; see the trajectory of train 2.
To calculate the value of Δ  , set   and set    are defined as the skip-stop plans of train  for the forward direction and backward direction, respectively.Each set is comprised of a series of the "skip arcs"; each "skip arc" connects two stations where train  stops, between which at least one station has been skipped.For example, as shown in Figure 4  − The next step is to calculate the maximum travel time difference between train  − 1 and train  at any point along the metro line, which can be expressed as . ., }.In the cases that this maximum travel time difference is larger than zero (i.e., train  travels faster than train  − 1), an extra amount of time is added to the departure time of train  to avoid collision.In the case that the maximum travel time difference is smaller than zero, no adjustment is needed; see formulation (15).Similarly, we use Δ   to represent the maximum travel time difference for the backward direction.
To ensure minimum safety headway between two consecutive trains at any point of the line, the departure time of train  is constrained at terminal stations 1 and   + 1 by formulation (16a) and formulation (16b), respectively.
This constraint, however, is too strict for skip-stop operation.In the cases that the safety operation headway is larger than the demand headway (due to skip operation), we cannot guarantee that the demand headway is satisfied at all times.Therefore we consider constraints (18a) and (18b) as a compromise to constraints (17a) and (17b).
In practice, some stations such as terminal stations and transfer stations should not be skipped.For this, we further define a set SS that contains the stations that can be skipped.For the stations that do not belong to this set, they should not be skipped; see constraint (21).Here SS ∪ SS = SN, where SN is the set that contains all metro stations. , = 0, ∀ ∈ SS. (21)

Solution Approach
In the previous section, the FSSS model is formulated as a mixed integer nonlinear programming (MINLP) problem, which is comprised of the objective function (9), constraint (1), constraint (10)-constraint (11), constraints (16a) and (16b), and constraints (18a) and (18b)-constraint (21).Note that the model requires departure time optimization at each station, which imposes high computational complexity.To simplify  the problem, we apply a rule-based simulation (( 22)-( 24)) that generates departure times to satisfy constraints (16a) and (16b) and constraints (18a) and (18b)-constraint (19).A similar rule-based computation of train departure times was implemented by Salzborn [23].First, we define the terminal departure time as (22); the departure time at the terminal departure takes the maximum allowed value that satisfies the safety and passenger demand constraints; therefore it should result in the smallest required number of trains.Similarly, the departure time at the turnaround terminal should take (23).Equation ( 24) is then used to update the departure times at the intermediate stations.
It is particularly important to note that if the turnaround terminal is considered as an intermediate station (as in most existing literature), an extra amount of the maximum travel time difference, which takes the form of max{Δ  , Δ   }, needs to be added to the departure terminal.In the cases that Δ   is larger than Δ  , the departure time at the departure terminal can be expressed by (25).And this leads to delayed departure time compared to (22), which results in line capacity loss.
Furthermore, it is widely known that NP-hard problems suffer from the "Curse of Dimensionality"; separating the forward and backward operation scheme can significantly reduce the dimension of the problem; that is, the computational complexity reduces from (2 2   tr ) to (2    tr +1 ).
Using the simulation scheme, the departure times can be determined, and the problem is thus reduced to the optimization of the binary variables  , .The problem is NP-hard and an exhaustive search of the skip-stop schemes has an exponential time complexity (2    tr ).To avoid enumerating all possible outcomes, we develop a genetic algorithm (GA) approach to guide the searching process.GA is a metaheuristic method that imitates the process of natural evolution.
GA is motivated by the principles of natural selection and survival-of-the-fittest individuals [24].It is proved to be efficient in solving NP-hard optimization problems in the transportation field [25][26][27][28].This method was also used by Lee et al. [13] to find a solution to the A/B mode-based skipstop operation problem.In their work, gene is defined as the station type and it can take the choices of A, B, or AB.The time complexity of their problem is therefore 3   .Compared to Lee et al. [13], the proposed FSSS model offers flexible schemes and it can further consider the operation at the turnaround terminal.Realistic operational and infrastructure constraints are also considered.Therefore the proposed method is different from the one in Lee et al. [13].The solution procedure of the developed GA algorithm is shown in Figure 5.
In our solution approach, a gene is defined as the operation choice at the specific station (0 or 1 which denotes skip or stop) and therefore a chromosome is defined as a binary string that embodies { , |  ∈ SS}.And this ensures that constraint (21) is satisfied throughout the searching process.In the initialization step, a number of   chromosomes are randomly generated.Next, we check the feasibility of the generated chromosomes, using constraint (20).An arbitrarily large objective (fitness) value is assigned to the infeasible chromosomes to penalize the infeasible predecessors; for the feasible predecessors, their fitness values need to be evaluated using the objective function.We then apply the GA operators Equ. ( 14)-( 16a) and (16b) to the list of chromosomes to generate a list of new solutions.The GA operators include (i) a selection operator, which is used to select  ps fitter solutions based on the Roulette Wheel Selection Algorithm [29]; (ii) a crossover operator, which is used to generate the child solutions from the selected solutions through a recombination process (with crossover rate   ; see Back [30]); (iii) a mutation operator, which uses a small probability (i.e.,   ) to model the genetic mutation process.We further interpret the binary strings as Gray-coded integers, which is a proven practice that can be used to avoid the hamming cliff problem (see Back [30]).We then check the generated new solutions to see if they satisfy the stopping criteria; that is, the number of iterations exceeds a predefined value () or convergence is reached.If so, the final solution is recorded; otherwise, return to check the feasibility.Readers can refer to Jong [24] for more detailed information of GA.

Data Processing and Experiment
Setup.Smart card data collected from urban transit systems have been widely used in transportation applications.Pelletier et al. [31] reviewed the application of smart card data in public transit field and categorized the data usage into three levels: the strategic level for long-term planning, the tactical level for service adjustment, and the operational level for performance measurement.Munizaga and Palma [32] proposed a method to estimate metro passenger OD using GPS data and smart card data in Santiago.Hong et al. [33] developed an approach that can match passenger trips to trains.Zhang et al. [22] formulated a model to optimize the urban rail operation based on the detailed passenger transaction data during weekdays.It is revealed by the literature that smart card data can better reflect the passenger OD pattern.The proposed FSSS model was tested using the smart card data and operation parameters from a real world bidirectional metro line in Shenzhen, China.The metro system is shown in Figure 6.We selected Line 1 (with 28 intermediate stations and two terminals; highlighted in green) to test the proposed FSSS model.
The dynamic passenger OD rates were obtained from the smart card data using the following steps: (i) a transfer rule was defined based on the least number of stations; that is, for a trip that starts and ends at different lines, the trip will use the transfer station that yields the least number of stations  during the trip; (ii) the transfer rule was applied to the records that either start or terminate at some station of Line 1; (iii) the origins and destinations of all the Line 1 trips were then identified; (iv) by aggregating the OD data for the weekday off-peak period (13:30-17:00) during one month, the dynamic passenger OD rates ( ,  ) were calculated for each time period , which was set to be 15 minutes in the experiment.
This metro line currently operates under an all-stop scheme.The round trip travel time is about 138.8 minutes. tr was set to be 10.The minimum safety headway, the demand headway, and the minimum turnaround headway were set to be 2 minutes, 6 minutes, and 3 minutes, respectively.According to the real world practice of Chicago Metro [11], the maximum acceptable waiting time   was 12 minutes during the peak hours.As our study targets off-peak operation,   was set to be 15 minutes.The time loss  was calculated based on the reachable maximum speed  max using formulation (26).Here  is the maximum acceleration rate and  is the maximum deceleration rate.Parameters  max , , and  were set to be 80 km/h, 0.8 m/s 2 , and 1 m/s 2 , respectively; the calculated  was 25 seconds, using (26).Further discussions of this formulation can be found in Zheng et al. [16].
The off-peak boarding and alighting counts at each station are shown in Figure 7(a).It is found that the number of passengers is relatively low at station 21, station 22, and station 23.An example of the time-dependent passenger demand from station 8 to station 15 is illustrated in Figure 7(b).The section travel time between two consecutive stations and the dwell time at each station are shown in Table 2.The dwell time at transfer stations {3, 4, 8, 9, 15, 24} is 45 seconds, which is higher compared to the 35-second dwell time at other intermediate stations.The minimum turnaround time was set to be 120 seconds.The section travel times and dwell times were set according to the real world operating conditions.Parameter  was set to be 2 which means that one train is allowed to be stored at the turnaround line.
The all-stop operation scheme was simulated by setting parameter  , to 1 for all the trains and stations.We then use formulation (1)-formulation (11) and formulation ( 22)-formulation (24) to calculate the performance of the system under the all-stop scheme.It is found that the average waiting time  pg wait and the average in-vehicle time  pg tv are 3.0 minutes and 13.5 minutes, respectively.The average travel time under the all-stop scheme is 16.5 minutes.The total number of onboard passengers  pg is 26146.

Skip-Stop Optimization Results
. The proposed FSSS model is then solved using the GA-based solution approach.There are five major parameters used in GA, including the size of the chromosome list (  ), population's size ( ps ), maximum number of iterations (), crossover rate (  ), and mutation rate (  ).The first three parameters were set to be 200, 500, and 2000, respectively.The crossover rate and the mutation rate are particularly important and are often found to affect the modeling performance (Back [30]).In the experiment, these parameters were set to be 0.8 and 0.005.The selection of parameters is based on a standard sensitivity analysis approach; details are not presented here.
For the FSSS, we manually set SS to include 20 low demand stations, that is, SS = {10, 13,14,18,21,22,23,27,28,29,32,33,34,38,39,40,43,47, 48, 51}.The corresponding chromosome length is 200 (20 binary variables multiple 10 trains).The problem is solved on a computer with 2.0 GB RAM on a 2.0 GHz Intel Core i7 processor.The GA runtime is around 2 hours.However, as the FSSS model is proposed for offline operations, the computation time is not a big concern.By solving the MINLP, the optimized FSSS is listed in Table 3, only for the stations that belong to set SS.For the stations that belong to SS, they cannot be skipped.It is found that the stations that have the lower passenger demand (e.g., station    agencies.The performance of FSSS is then compared to the all-stop operation, as shown in Table 3. The results present in Table 4 show that the objective value of FSSS decreases by 5%, from 990 seconds to 940 seconds, which indicates a travel time saving of about 50 seconds per passenger.It is found that the average in-vehicle time is reduced by 6.5% with slightly increased (less than 4 seconds) passenger waiting time.As indicated by the value of  bpg , about 26.9% of total passengers (6604 out of 24576) benefited from FSSS.Compared to the all-stop operation, the number of onboard passengers under FSSS slightly decreases by 6%.This is mainly because the last train skips some stations and some passengers are not served by this train.These passengers are assumed to be picked up by the subsequent trains [18].
The timetables of the two operation schemes are illustrated in Figure 8.The light line denotes the all-stop operation, while the dark line represents the FSSS operation.It is found that FSSS has a slightly higher travel speed (the dark lines travel faster than the light lines).In specific, the average round trip travel time under FSSS is 135.7 minutes, which is more than 3 minutes shorter compared to the all-stop operation.The average departure interval is about 6 minutes which is similar to the all-stop operation.
Figure 9 indicates the seat occupancy and maximum section passenger load under the two schemes.Seat occupancy is defined as the passenger load over the seat numbers.It is found that, compared with all-stop scheme, passengers who board the skip-stop trains have a better chance to find empty seats, meaning FSSS is more favorable to the passengers as it enhances their level of service.By analyzing the maximum section passenger load, we find that the highest passenger load under FSSS is 558, which is far below the train capacity (1860).Therefore FSSS does not cause congestion on the trains.

Analysis of the Passenger Confusion Rate.
It is widely recognized that the major problem in applying skip-stop operation is to keep the passengers informed regarding the skip-stop plans.The solution to it relies on information dissemination such as radio broadcast in stations/trains.Nonetheless, the authors believe the effect of certain percentage of passengers being confused by the skip operation and failing to find the right train is worth studying.The percentage is referred to as the passenger confusion rate hereafter in the paper.As noted in Assumption 5, we assume that these confused passengers shall notice their mistake (e.g., via train radio broadcast) once they board the train, and they will get off at the next stopping station.They will then wait at the station until being transported by the next available (correct) train to their destination.Variable   is introduced here to denote the confusion rate, and a sensitivity analysis was conducted to see the effects of   on the modeling performance.The variable   was tested by taking values from 0 to 100%, with an increment of 20%.
The calculation of the boarded passengers (on train  which stops at station , with a destination ) can be divided 0 tw For Case 2, the number of confused passengers can be calculated using (30).They will get off at the nearest stop denoted as   = min{  |   > ,  ,  = 1} after departure from station .If   <  (destination has not been passed), these passengers shall wait until   = min{  |   > ,    ,     , = 1} arrives.If   > , these passengers need to take a train in the opposite direction, noted as , where   is the total number of stations.The associated waiting time can be obtained using (31), and the in-vehicle time can be calculated using (32) By taking the summation of the above variables for all OD pairs of all trains, the total boarded passenger number, total waiting time, and in-vehicle time could be calculated, and the objective value could then be obtained using (9).The results are provided in Table 5.
It is observed that as the confusion rate increases (e.g., poorer information dissemination), the number of confused  passengers increases proportionally.However, this number only represents a small proportion of the some 25,000 passengers in total.This is because only the OD pairs with skipped destinations are affected, and the skipped stations have comparably low volume.For the confused passengers, it is found that their in-vehicle time and especially their waiting time increase.The increased in-vehicle time is because some passengers missed their destination and have to take the trains to the inverse direction (as in Case 2,   > ).And the increased waiting time is mainly because the correct train has skipped station   and the confused passengers are not able to get on the correct train at that station.Per the specific experiment, we found that FSSS always outperforms the skipstop scheme even when all the passengers of the skipped OD pairs get confused (i.e.,   = 1).Therefore the system performance is stable under FSSS.

Conclusion
In this study, we developed a FSSS approach for bidirectional metro line operation during off-peak periods to better accommodate the spatially and temporally varied passenger demand.The problem was formulated as a mixed integer nonlinear programming (MINLP) model that minimizes the average passenger travel time.The proposed scheme is able to optimize the skip-stop scheme and it also models the train departure time at the departure and turnaround terminals.The departure intervals are constrained by three types of headways, that is, the minimum safety headway, headway that meets the passenger demand, and the minimum turnaround headway.A GA-based solution approach was then developed to efficiently solve the problem.The proposed scheme was tested using the smart card data collected at a real metro line in Shenzhen, China.Time-dependent passenger demands were considered in the experiment and the passenger confusion rate was also analyzed in the case study part.
The model performance was assessed with respect to the passenger benefits and the operational benefits.Overall, the FSSS is effective in timetable design.It was found that the FSSS is able to reduce the average passenger travel time by 0.82 minutes, which corresponds to 5% travel time saving.About 26.9% of total passengers benefit from the skip-stop scheme.was also observed that the average round trip travel time under the skip-stop scheme is 135.7 minutes, which is more than 3 minutes shorter compared to the all-stop operation.The average train departure interval is 6 minutes which is the same compared to the all-stop scheme.From the timetable, it was found that 35 stations are skipped.Transit agencies may benefit from this scheme due to the savings of energy and operational costs.Through the sensitivity analysis of passenger confusion rate, we found that the proposed scheme always outperforms the all-stop scheme even when most passengers of the skipped OD pairs are confused and fail to get on the right train.
This paper mainly considers the skip-stop operation during off-peak scenarios in which the demands for some stations are quite low.As shown by the results, the skip operation does not lead to excessive waiting time at these stations.For peak-hour operations, since the passenger demands at most stations are relatively large, the skip operation may lead to longer waiting time at the skipped stations.In extreme cases, the number of waiting passengers may exceed the design capacity of the station.In future work, the tradeoff between in-vehicle travel time and passenger waiting time needs to be properly leveraged.The FSSS model currently requires a rulebased simulation to calculate the train departure time.Under the rule-based simulation, the terminal departure times need to be adjusted by the maximum travel time difference between two consecutive trains.This adjustment, however, may correspond to suboptimal solutions since departure choice can be made more wisely at each station including the intermediate stations.The GA-based solution approach requires a predetermined set (SS) that specifies the stations that can be skipped.This set needs to be determined using a more rigorous approach.The authors are also interested in extending the current methodology to optimize the skip-stop operation for a network-wide metro system.More results and findings will be provided in our subsequent studies.

Figure 1 :
Figure 1: Illustration of a bidirectional metro system.

Figure 2 :
Figure 2: The departure of two consecutive trains; (a) trains  and    depart in the same time period; (b) trains  and    depart in different time periods.

Figure 3 :
Figure 3: The illustration of the method to obtain the train safety headway.
Time-dependent passenger OD rates from station 8 to station 15

Table 1 :
Existing skip-stop schemes for metro operation.

Table 1 :
Continued. , pg : total number of onboard passengers for train  at station ;  , wait : total waiting time for passengers who board train  at station ;  , veh : total in-vehicle time for passengers who board train  at station ;  pg : number of onboard passengers;  wait : total passenger waiting time, in seconds;  veh : total passenger in-vehicle time, in seconds; SS: the candidate skipping station set that contains the stations that can be potentially skipped;   : the set that records the skip-stop schemes of train  for the forward direction; : the set that records the skip-stop schemes of train  for the backward direction;    : starting station of arc  of train ;    : ending station of arc  of train ; arc    ,   : skip arc that connects two stations (   and    ), between which at least one station is skipped; Δ    ,   : maximum travel time difference between    and    ; Δ  : maximum travel time difference between train − 1 and train  for the forward direction; Δ   : maximum travel time difference between train − 1 and train  for the backward direction.ss : number of skipped stations; To calculate the values of  , In this case, the passenger arrival rate is uniform.The resulting  , pg ,  , wait , and  , veh can be calculated as in (3)-(5).Here  , −     , represents the departure interval between trains  and    at station ; therefore  ,  ( , −     , ) represents the amount of passengers who arrive at station  and plan to travel to station  during this time interval. , −  , −   dw refers to the in-vehicle travel time for passengers travel from  to .The second case is when trains  and    depart in different time periods (i.e.,     , ∈ [( − 1), ) ∧  , ∈ [, ( + )) ∀ ≥ 1).The corresponding  , pg ,  , wait , and  , pg ,  , wait , and  , veh , two cases are considered; see Figure 2. The first case is when trains  and    depart in the same time period (i.e.,  , ,     , ∈ [( − 1), )).
, arc   1 ,  1 denotes the fact that train  stops at station   1 = 1 and stops again at station   1 = 4, between which station 2 and station 3 are skipped.Here    refers to the starting station of arc ;   , )(  dw +), where   dw is the dwell time at station  and  is the time loss due to the acceleration and deceleration associated with a stop at a station.Therefore, the maximum travel time difference between train  − 1 and train  from stations refers to the ending station of arc .The last "skip arc" in set   is denoted as arc    ,   . and    , denoted as Δ    ,   , can be represented as Δ    ,   =  (1 −  , ) (  dw + ) Here   is the headway that meets the passenger demand.,1≤, ∀1 ≤  ≤  tr ,(17a) ,  +1 ≤  0,  +1 +   , ∀1 ≤  ≤  tr .
(20)3.Turnaround Headway.In our model, we specifically consider  ( ≥ 1) turnaround lines at the turnaround terminal.Under the infrastructure geometry, train  −  should depart from station   + 1 before train  departs from station   .That is, formulations (18a) and (18b) should be satisfied.Here   refers to the minimum turnaround headway at the turnaround station.The formulation can better utilize the storage capacity at the turnaround station.,  +   −  −,  +1 ≥ 0, ∀ ≤  ≤  tr .(19)2.5.4.Operational Constraint.To avoid intolerable waiting time at the trip origins, we assume that each viable OD pair should be served in an acceptable time period which is marked as   ; that is, constraint(20)needs to be satisfied.Here  = max{  |   < ,  ,  ,  ,  > 0∧ ,  ,

Table 4 :
(27)ormance of all-stop operation and FSSS.Further define set    to represent the set of trains between trains    and  that stop at station  and skip station ; that is,    = { |  , = 1∧ , = 0,    ≤  < }; the cardinality of the set is denoted as ;  ℎ is used to denote the value of the ℎ-th element in set    .For Case 1, the number of onboard passengers from station  to station  via train  can be calculated using(27).
into two cases.Case 1 is that train  will stop at station ; therefore passengers (from OD  to ) board the right train.In Case 2, train  will skip station ; people who board the train are the confused passengers.

Table 5 :
Sensitive analysis of the passenger confusion rate.