Adjust train schedules to minimize passenger waiting times during disruptions

Huijuan Zhou (Beijing Key Laboratory of Intelligent Control Technology for Urban Road Traffic, North China University of Technology, Beijing, China)
Rui Wang (Beijing Key Laboratory of Intelligent Control Technology for Urban Road Traffic, North China University of Technology, Beijing, China)
Dongyang Weng (Beijing Key Laboratory of Intelligent Control Technology for Urban Road Traffic, North China University of Technology, Beijing, China)
Ruoyu Wang (Beijing Key Laboratory of Intelligent Control Technology for Urban Road Traffic, North China University of Technology, Beijing, China)
Yaoqin Qiao (Beijing Key Laboratory of Intelligent Control Technology for Urban Road Traffic, North China University of Technology, Beijing, China)

Smart and Resilient Transportation

ISSN: 2632-0487

Article publication date: 14 December 2023

Issue publication date: 11 March 2024

201

Abstract

Purpose

The interruption event will seriously affect the normal operation of urban rail transit lines,causing a large number of passengers to be stranded in the station and even making the train stranded in the interval between stations. This study aims to reduce the impact of interrupt events and improve service levels.

Design/methodology/approach

To address this issue, this paper considers the constraints of train operation safety, capacity and dynamic passenger flow demand. It proposes a method for adjusting small loops during interruption events and constructs a train operation adjustment model with the objective of minimizing the total passenger waiting time. This model enables the rapid development of train operation plans in interruption scenarios, coordinating train scheduling and line resources to minimize passenger travel time and mitigate the impact of interruptions. Regarding the proposed train operation adjustment model, an improved genetic algorithm (GA) is designed to solve it.

Findings

The model and algorithm are applied to a case study of interruption events on Beijing Subway Line 5. The results indicate that after solving the constructed model, the train departure intervals can be maintained between 1.5 min and 3 min. This ensures both the safety of train operations on the line and a good match with passengers’ travel demands, effectively reducing the total passenger waiting time and improving the service level of the urban rail transit system during interruptions. Compared to the GA algorithm, the algorithm proposed in this paper demonstrates faster convergence speed and better computational results.

Originality/value

This study explicitly outlines the adjustment method of using short-turn operation during operational interruptions, with train departure times and station stop times as decision variables. It takes into full consideration safety constraints on train operations, train capacity constraints and dynamic passenger demand. It has constructed a train schedule optimization model with the goal of minimizing the total waiting time for all passengers in the system.

Keywords

Citation

Zhou, H., Wang, R., Weng, D., Wang, R. and Qiao, Y. (2024), "Adjust train schedules to minimize passenger waiting times during disruptions", Smart and Resilient Transportation, Vol. 6 No. 1, pp. 2-20. https://doi.org/10.1108/SRT-10-2023-0009

Publisher

:

Emerald Publishing Limited

Copyright © 2023, Huijuan Zhou, Rui Wang, Dongyang Weng, Ruoyu Wang and Yaoqin Qiao.

License

Published in Smart and Resilient Transportation. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence maybe seen at http://creativecommons.org/licences/by/4.0/legalcode


1. Introduction

Urban rail transit is essential for providing residents with fast, efficient and environmentally friendly transportation. However, interruptions in its operation can lead to significant delays and lower service quality. To address this, it is crucial to develop practical and reliable train operation adjustment strategies using scientific decision models. These strategies are vital for effectively managing interruptions and minimizing train delays, ensuring smooth and efficient rail transit operations.

Train operation adjustment studies under delay conditions primarily focus on constructing optimization models based on real fault scenarios. These models are used to optimize train schedules by incorporating relevant train scheduling strategies. Canca et al. (2016) considered the impact of transfer connections and developed a network-level train timetable optimization model based on short-turning operations. Zhou et al. (2022) considered the impact of transfer connections, a network-level train timetable optimization model is constructed based on interlacing route operations. Jia et al. (2022) considered the stochastic nature of passenger flow and proposed a coordinated optimization model. It integrates short-turning plans, passenger flow control strategies and train timetable adjustments based on the dynamic relationship between passenger flow evolution and train operations. Zhu et al. (2020a) addressed the uncertainty in the duration of interruptions by formulating the train timetable rescheduling problem as a deterministic equivalent form of a rolling time-space two-stage stochastic programming problem. They calculate the optimal timetable solution for interruption scenarios based on the probabilities associated with the interruption duration. Binder et al. (2021) designed a train timetable adjustment algorithm considering train capacity. It aims to achieve a user equilibrium between timetables and passengers by optimizing objectives, including operating costs, deviation from the original timetable and passenger convenience. Zhang et al. (2021) considered cancellation, delay and speed reduction strategies and developed a mixed integer nonlinear model for adjusting train speeds and timetables. Ghaemi et al. (2016) created a macroscopic train scheduling model using mixed integer linear programming to obtain the optimal turnaround time for trains in the event of complete line disruptions. Zhu et al. (2020b) developed a model using mixed integer linear programming to minimize generalized travel time for passengers during interruptions. The model considers passenger reactions, decision-making behavior and station capacity, aiming to minimize onboard time, waiting time at origin and transfer stations and the number of transfers. Lin et al. (2018) established an optimization model for train operations, considering constraints such as the number of train units and passenger travel time. The objective of the model is to optimize the train schedule, prioritizing the rapid evacuation of stranded passengers. Yang et al. (2021) proposed a coordinated optimization method for train adjustment and passenger flow control to improve train operation order and passenger travel experience. Their approach involves a two-level linear programming model that includes skip-stop scheduling and multi-station passenger flow control.

Most of the above literature only makes quantitative analysis of train operation adjustment models in urban rail transit during disruption scenarios; the considered passenger flow data is typically fixed and does not account for dynamic changes. Therefore, based on the discussion of the above literature, this paper adopts a train operation adjustment method using short-turning operations. The objective function is to minimize the total waiting time for all passengers. It involves reconfiguring the train timetable during disruption scenarios and dynamically adjusting train departure and stopping times. Meanwhile, an improved genetic algorithm (IGA) is proposed and applied to an actual disruption event on Beijing Metro Line 5 to supplement the research on train operation adjustment under disruption scenarios.

2. Adjustment objectives and problem analysis

2.1 Adjustment objectives

When the urban rail transit system experiences a disruption event, the train operation adjustment needs to quickly formulate a train operation organization plan under the disruption scenario according to the passenger demand, taking into account a series of condition constraints such as line conditions, train numbers, headways and dwell times. The purpose is to reduce the travel time of passengers on the line by coordinating train and line resources, and to achieve that the train can resume the planned timetable operation as soon as possible after the disruption ends. Based on existing results, the objectives of train operation adjustment can generally be constructed from various optimization objectives from the perspectives of passengers and operators.

2.1.1 Passenger demand.

From the perspective of passenger demand, the objectives of train operation adjustment can be divided into three aspects: first, to reduce the total travel time of all passengers; second, to reduce the total waiting time of all passengers; and third, to reduce the in-vehicle congestion cost. The reduction of total travel time is the time difference from when passengers enter the platform to wait until they get off at the destination station, including platform waiting time, section running time and train dwell time; the reduction of total waiting time includes platform waiting time and in-vehicle passenger waiting time; the in-vehicle congestion cost is to make the train load factor as close as possible to the range that passengers can tolerate. The three are interrelated and are important factors that affect passenger travel choice, which are related to passenger satisfaction with the service level of the urban rail transit system.

2.1.2 Operator demand.

From the perspective of operator demand, the objectives of train operation adjustment can be divided into three aspects: first, to reduce the total delay time of trains; second, to ensure the balance of train arrivals and departures; and third, to reduce the number of people gathering on the station platforms. The reduction of total delay time of trains is to minimize the deviation between the actual timetable and the planned timetable of all affected trains; the balance of train arrivals and departures is to avoid the situation that the arrival interval of two adjacent trains is too small, resulting in a large difference in train load factor between the two trains; the reduction of the number of people gathering on the station platforms is to ensure that the platform waiting population is within the safe capacity range, and to avoid the occurrence of large passenger flow events at the station. The three are closely related and are important factors that operators consider for disruption scenarios.

The above-mentioned adjustment objectives are interrelated and contradictory; therefore, different degrees of adjustment need to be made according to the actual scenario in the process of formulating the optimization objectives of train operation adjustment.

2.2 Problem analysis

Large passenger flow is the decisive factor for subway failures. The morning and evening peak periods increase the frequency of vehicle use, vehicle wear and tear and failure frequency. According to the official Weibo announcement of Beijing Subway, signal failure is the most frequent and common failure event in subway operations, which often causes bidirectional delay and interruption of the line. Therefore, this paper takes a subway line that suffers from a signal failure event causing bidirectional interruption as the research object and uses the line passenger flow demand, train operation data, etc., as inputs to construct an optimization model of train timetable under the interruption event, so as to quickly and accurately replan the train timetable under the interruption event, and thus guarantee the capacity and passenger travel of the urban rail transit system to a certain extent.

In the abstracted line operation scenario depicted in Figure 1, when there is an interruption event in the N-3 interval of a bidirectional line, both directions are affected and normal operations are disrupted. If the interruption persists without resolution, timely emergency adjustment methods are required to ensure the transportation capacity of the line. This paper adopts the short-turn operation adjustment method, which allows for flexible modifications to train departure times and station dwell times.

Solution approach: Assuming Station N-3 or nearby stations can accommodate short-turning operations, trains can be organized to perform early turnarounds at Station N-3. This maintains segment operations between Station 1 and Station N-3, aiming to preserve service levels, minimize passenger waiting time, reduce station crowding and ensure passenger travel and safety.

3. Train schedule optimization model based on passenger waiting time

3.1 Model assumptions and symbol meanings

To facilitate research and modeling, the following assumptions are made in this study:

  • The stations where train reversal operations can be conducted are not used for train storage.

  • The railway line does not have the conditions for overtaking.

  • The location of the interruption is known, and the duration of the interruption can be predicted.

  • The passenger arrival rate and alighting rate at each station are calculated based on actual card-swiping data.

  • During the boarding and alighting process, the “alight before boarding” and “first-come, first-serve” principles are followed, and the incoming passenger flow at each station is considered to be uniformly distributed.

  • At any given time, only one train is allowed to stop at the platform of each station in both directions. All trains follow a station-to-station stopping pattern during operation.

Model parameters and symbol definitions as shown in Table 1.

3.2 Objective function

This paper focuses on passenger demand and operational considerations, aiming to optimize train operation adjustments based on the smoothness of passenger flow. The objective is to minimize the total waiting time for passengers. The model takes into account onboard time, waiting time for passengers between consecutive train arrivals, and secondary waiting time for stranded passengers. To ensure timely emergency decision-making during disruption events, the calculation of waiting time between consecutive train arrivals is simplified for improved computational speed. The objective function is calculated as follows:

(1) minT=WT1+WT2+WT3

3.2.1 Total waiting time of passengers on board.

The in-vehicle time of passengers is the product of the number of remaining passengers in the vehicle after the train arrives at the station and the actual dwell time. The number of remaining passengers in the vehicle after the train arrives at the station is obtained by calculating the difference between the number of passengers in the vehicle when the train arrives at the station and the number of passengers getting off when the train arrives at the station. Then, the calculation is performed by adding up all the stations and trains on the line. The calculation formula is as follows:

(2) WT1=k=1|K|u=2|U1|(Qk,u1inDk,ual)(tk,udeptk,uarr)

3.2.2 Total waiting time of passengers arriving within the interval between adjacent service trains.

The waiting time of passengers arriving at the station within the departure interval of two adjacent trains is the integral of the product of the passenger arrival rate at the station and the passenger waiting time for departure within the departure interval of two adjacent trains, and then the calculation is performed by adding up all the stations and trains on the line. The calculation formula is as follows:

(3) WT2=k=1|K|u=1|U1|tk1,udeptk,udepλu(t)(tk,udept)dt

Based on the above analysis, to improve the speed of model iteration calculation, the above formula is simplified. The calculation formula is as follows:

(4) WT2=k=1|K|u=1|U1|[λun(tk,udeptk1,udep)22]

3.2.3 Total waiting time of passengers stranded on the platform.

The secondary waiting time of passengers who cannot leave the station by the current train and are stranded at the station is the product of the number of platform stranded passengers when the train leaves the station and the time difference between the departures of two adjacent trains, and then the calculation is performed by adding up all the stations and trains on the line. The calculation formula is as follows:

(5) WT3=k=1|K|u=1|U1|(Wk,uhall(tk+1,udeptk,udep))

3.3 Constraints related to train operation safety

  • Interval running time constraint:

    (6) tk,u+1arrtk,udep=tk,urun

The model assumes that the interval running time for trains is a fixed value:
  • The constraints on train dwell time at stations:

    (7) tk,udeptk,uarr=tk,udwell
    (8) tmindwelltk,udwelltmaxdwell

The dwell time in this paper is a variable value. Equation (8) ensures that the train dwell time at stations satisfies the actual operational requirements, with the dwell time being a variable value:

  • Train tracking interval constraint:

    (9) tk+1,udeptk,udepH1
    (10) tk+1,uarrtk,uarrH2
    (11) tk+1,uarrtk,udepH3

In the model assumption, it is assumed that no overtaking will occur during the train operation process. Then, equations (9)–(11) are to ensure that the difference between the arrival and departure times of two adjacent trains at station u is within the range of the minimum safe driving interval so as to ensure the safety of the line operation.

  • Train turnaround operation constraint:

    (12) tk,u+1arrtk,udep=δz

This paper proposes to adopt the adjustment method of short-turn operation to ensure the service level under system operation interruption. To ensure the safety of train operation, equation (12) is the difference between the arrival time of the down train after the up train completes the turnback operation at the turnback station and the departure time of the up train at the turnback station, which is the turnback operation time of the train:

3.4 Constraints related to passenger flow control

  • Actual number of passengers boarding the train:

    (13) Ak,u,von={min{Pk,u,vwait,Cap},if u=1min{Pk,u,vwait,CapQk,u1in+Dk,ual},if 1<u|U|
    (14) Ak,uon=v=u+1|U|1Ak,u,von

  • Number of passengers alighting from the train:

    (15) Dk,ual={0,   if u=1v=1u1Ak,von·ρu(t),  if  1<u|U|

  • Number of passengers waiting to board the train:

    (16) Pk,u,vwait=Wk1,u,vhall+Ek,u,v
    (17) Pk,uwait=v=u+1|U|1Pk,u,vwait=v=u+1|U|1(Wk1,u,vhall+Ek,u,v)

  • Number of passengers stranded at the station:

    (18) Wk,u,vhall=Pk,u,vwaitAk,u,von
    (19) Wk,uhall=Pk,uwaitAk,uon=v=u+1|U|1(Pk,u,vwaitAk,u,von)

  • Number of passengers on board the train:

    (20) Qk,uin={Ak,uon,  if  u=1Qk1,uinDk,ual+Ak,uon,  if  1<u<|U|

  • Train capacity constraint:

    (21) Qk,uinCap

3.5 Constraint linearization

In the above-mentioned constraints related to passenger flow control, equation (20) represents a nonlinear expression. According to linearization theory, the constraints can be linearized by introducing two binary variables, α,ε ∈ [0,1], to transform the nonlinear relationships into linear inequalities. The specific linearization process is as follows: if Pk,u,vwait>Cap is true, then α = 0; otherwise, α = 1. Similarly, the comparison between Pk,u,vwait and CapQk,u1in+Dk,ual is linearized in the same manner:

(22) if  u=1:{Ak,u,vonPk,u,vwaitAk,u,vonPk,u,vwaitM(1α)Ak,u,vonCapAk,u,vonCapMα
(23) if  1<u|U|:{Ak,u,vonPk,u,vwaitAk,u,vonPk,u,vwaitM(1ε)Ak,u,vonCapQk,u1in+Dk,ualAk,u,vonCapQk,u1in+Dk,ualMε

4. Solving algorithm

4.1 Algorithm principle

This paper presents a train operation adjustment model for urban rail transit based on fault interruption events, with the objective of minimizing passenger waiting time. The model involves a large number of stations and train services, and the iterative process of solving the model is complex. In emergency scenarios with unexpected interruptions, the response time is crucial. However, precise solution algorithms can be time-consuming for large-scale problems. Therefore, this article adopts an IGA to solve the model. The genetic algorithm (GA) is modified and designed to solve the model, building upon traditional GAs.

This paper proposes an IGA based on the GA, as shown in Figure 2. The improvements are as follows:

  • To adjust the crossover rate and mutation rate adaptively according to the number of iterations and to realize the nonlinearity of the adaptive adjustment curve of the crossover rate and mutation rate so as to improve the global search ability and convergence speed of the algorithm. At the same time, in the adjustment process, the sigmoid function, which is a smoother neuron activation function in neural networks, is used as the adaptive adjustment curve so as to retain the excellent individuals as much as possible.

  • To adopt the elitist strategy, there is a problem in the IGA that it can adjust the crossover rate and mutation rate adaptively, but their values may gradually increase, resulting in the elimination of good individuals. Therefore, the elitist strategy is adopted to avoid eliminating the individual with the highest fitness value in the previous generation and to ensure that the algorithm can converge better.

4.2 Algorithm flow

Step 1: Input parameter initialization. Load passenger flow data, timetable data, interval running time, planned stop time and other initial data, set train capacity, population size, number of iterations, crossover rate, mutation rate and calculation times and other related parameters.

Step 2: Encoding. The decision variables of the model in this paper are the departure time and dwell time of the trains at each station on the line, which is a continuous parameter optimization problem. Therefore, real number encoding is used to encode the departure time and dwell time. Because the running time between stations is a fixed value, it is only necessary to determine the departure time of the train at the origin station, and then determine the dwell time at each station, and then the departure time of the train at the subsequent stations can be calculated.

Step 3: Generate the initial population. Considering the safety of train operation, use a random method to generate an initial array X within the range of [0, 1] according to the number of iterations. Under the condition of meeting the departure interval constraint and dwell time constraint, determine the departure time t of the first train within the range of [1, hmin], determine the dwell time t of the first train within the range of [tmindwell,tmaxdwell], determine the departure time of subsequent trains according to T = X (hmaxhmin) + hmin and determine the dwell time of subsequent trains in the same way. This generates an initial population of N chromosomes.

Step 4: Calculate the fitness of each individual in the population. Individual fitness is the evaluation criterion of the individual’s superiority or inferiority in the population, so it is also called the evaluation function, which mainly judges the individual’s fitness by the characteristics of the population individuals. The objective function of the model in this paper is to minimize the total waiting time of passengers, and the smaller the objective function is, the better the individual solution is. Therefore, the algorithm selects the objective function as the calculation formula of the fitness function:

(24) Fitness=min(T)
Step 5: Selection operation. Selection operation is also called replication operation, which is to use a method to select a certain number of excellent individuals from the current population as parents, and start breeding the next generation. It should be noted that the probability of being selected is determined by the individual fitness. This paper adopts the most commonly used roulette wheel selection method. The detailed calculation steps are as follows:
  • First, calculate the fitness function value fitness (vi) of each individual vi = (i = 1, 2 … popsize) in the population.

  • Then, calculate the total fitness value of the population: Fitness=i=1popsizefitness(vi).

  • Then calculate the probability pi = fitness(vi)/fitness and cumulative probability qi = Σpi of each individual vi = (i = 1, 2…popsize) in the population being selected.

  • Generate a random number z within the range of [0,1]; if z < qi, then select the first chromosome; otherwise, select the ith chromosome vi (i=2…popsize) that makes qi−1< z < qi true.

Step 6: Crossover operation. Crossover is the process of producing new individuals by selecting and exchanging genes from both chromosomes of two individuals. This paper adopts the simulated binary crossover method, which divides the parent individuals into two parts, performs crossover N/2 times, selects two individuals for crossover, uses the minimum departure interval hmin as the judgment condition and considers the number of iterations to use the neural network activation function to improve and adjust the crossover probability Pc, thus obtaining a new population:

(25) Pc=k1+logsig[k2×(iter/iterationsk3)]×k4
(26) logsig(x)=11+ex

Step 7: Mutation operation. Mutation is the process of randomly changing the individual gene characteristics in the population with a very small mutation probability. Generally, the mutation probability is set to a very small value. This paper considers the number of iterations to use the neural network activation function to improve and adjust the mutation probability Pm and adopts the polynomial mutation method to generate a new population:

(27) Pm=h1+logsig[h2×(iter/iterationsh3)]×h4

The mutation form xk = ck + δ ×(popmax-popmin) of the polynomial mutation method, then the process of generating new individuals by mutating the parent individuals is: first, randomly generate a random number ψ, ψ ∈ [0,1], select the chromosome position where the random number is less than the mutation probability, and then calculated according to the following formula:

(28) δ={[2u+(12u)(1δ1)ηm+1]1/(ηm+1)1,               u0.51[2(1u)+2(u0.5)(1δ2)ηm+1]1/(ηm+1),  u>0.5

Where, δ1=ckpopminpopmaxpopmin, δ2=popmaxckpopmaxpopmin, u ∈ [0,1].

Step 8: Elite preservation strategy. The algorithm is prone to eliminate excellent individuals in the iterative calculation process of crossover and mutation, so the elite preservation strategy is adopted to prevent the elimination of excellent individuals. The main idea of the elite preservation strategy is that when the calculation iterates to the nth generation, the optimal individual in the population is c (n), and if there is no individual better than c (n) in the next generation population, then c (n) individual is taken as the optimal individual of the next generation population.

Step 9: Algorithm termination judgment. Judge whether the maximum number of iterations of the algorithm has met the preset number of iterations. If it is greater than the set number of iterations, the calculation stops and outputs the optimal solution in the iterative calculation process. Otherwise, return to Step 5 and perform the calculation again until the number of iterations is equal to the set number of iterations, and then output the optimal solution in the iterative calculation process.

5. Case studies

5.1 Basic data and parameters

The study focuses on Beijing Subway Line 5 as the research object. The line has a total of 23 stations, and all transfer stations on the line are equipped with crossover tracks. The direction of travel is from Songjiazhuang to Tiantongyuan North. The research period is set from 17:45 to 18:30 on February 19, 2019. During this time, a switch failure occurred at Datunlu East Station, resulting in a disruption of train operations in the section between Datunlu East Station and Tiantongyuan North Station.

Based on the above interruption event in the section between Songjiazhuang and Datunlu East, a short-turn operation is conducted, and the revised train timetable is developed using the train operation adjustment model proposed in this study. For ease of statistical analysis, the stations in the direction from Songjiazhuang to Tiantongyuan North are numbered as Station 1 to Station 23, as shown in Table 2. The dwelling time of trains during the research period is shown in Table 3, and the running time of intervals is shown in Table 4. The section numbers are numbered in sequence according to the station numbers. Table 5 lists the parameters in the model and algorithm.

5.2 Results analysis

This article uses two algorithms (GA and IGA) to solve the model. Figure 3 shows the change in fitness values with the number of iterations, and Table 6 presents the results of the objective function calculation for both algorithms. It can be observed that IGA can achieve the solution in fewer iterations compared to the GA algorithm, and the calculation results show a time saving of 1,830,741 s for passengers within the system. IGA exhibits faster convergence speed and better computational results, which can enhance the operational service quality of urban rail transit systems under interruption conditions to some extent.

The initial and adjusted train timetables are shown in Figures 4 and 5, respectively. The initial train timetable data was obtained from the official website of Beijing Subway for the study period. From Figure 5, it can be observed that there are 20 train services operating in the up direction at Station 17 as part of the short-turn operation. The dashed lines in the figure represent the canceled train services in the down direction toward Station 23, totaling 22 train services. Table 6 shows that the replanned train timetable significantly reduces the waiting time for passengers within the system. Once the interruption event is over, all subsequent trains will resume their planned timetable operations.

Figures 68 show that the number of passengers staying at the station is significantly reduced after optimization, reducing passenger aggregation and ensuring the travel safety of passengers. Figure 7 shows that the number of passengers in all trains leaving the station is significantly reduced after optimization, which can not only match the passenger flow and traffic well but also improve the travel experience of passengers. Figure 8 shows that the number of passengers getting off at each station when all trains arrive is reduced after optimization, indicating that the optimized train operation plan can further ensure the travel of passengers to a greater extent, thereby reducing passenger aggregation on the platform.

These three evaluation indicators all show that the train timetable after optimization can further demonstrate the effectiveness of the model in this paper, which can reduce the degree of passenger aggregation and improve the travel experience to a certain extent.

In the model, both departure time and station dwell time are decision variables. After solving the model using the algorithm and processing the results, a comparative analysis was conducted on the optimized departure intervals (Table 7) and the optimized station dwell times (Table 8) compared to their initial values. From Table 7, it can be observed that there are significant differences between the optimized departure intervals and the initial ones. The planned departure intervals were fixed at 2 min or 3 min, which provided a balanced distribution of departure times but did not align well with passenger travel demands. However, the model’s solution allowed for departure intervals ranging from 1.5 min to 3 min, ensuring both safe train operations and better alignment with passengers’ travel demands. This effectively reduced the total waiting time for passengers. Table 8 shows that the optimized station dwell times differ significantly from the initial fixed values. Station 1 experienced the largest increase in dwell time, from 30 s to 56 s, whereas Station 12 experienced the largest decrease, from 60 s to 46 s. Other stations, such as Stations 2, 3 and 4, had some increases in dwell time, while Stations 6, 7 and 16 had some decreases. Overall, the total dwell time increased by 119 s, allowing more passengers to board the train and alleviating congestion on the platforms.

In conclusion, after the occurrence of an interruption event, the use of the small-loop operation adjustment method and the optimization model of train timetable based on passenger waiting can effectively reduce passenger waiting time and crowd congestion within the system, thereby improving the service level of the system under interruption events.

6. Conclusions

This study comprehensively considers the constraints of train operation safety, capacity and dynamic passenger demand. It adopts the small-loop operation adjustment method to handle disruption events and constructs a train timetable adjustment model with the objective of minimizing the total passenger waiting time. An IGA is proposed as the solving algorithm. The model and algorithm are applied to a real disruption event on Beijing Metro Line 5 for validation. The results demonstrate the following:

  • After a disruption event occurs, implementing a small-loop operation plan allows for the continued operation of some trains, thereby reducing the impact of the disruption on the rail transit system. By limiting the scope of the disruption, its effects can be contained within a smaller area, allowing other stations to continue normal operations and avoiding a complete shutdown of the entire line. This approach contributes to improved overall operational reliability and passenger satisfaction by minimizing the extent of the disruption and its consequences.

  • The proposed train operation adjustment model, as validated through practical examples, demonstrates that the initial departure intervals exhibit good balance but fail to align well with passenger travel demands. However, after optimization, the departure intervals range between 1.5 min and 3 min, ensuring both the safety of train operations and a better alignment with passengers’ travel needs. This optimization effectively reduces the overall waiting time for passengers and enhances the service level of the urban rail transit system during disruptions.

Figures

Illustration of train operation organization during interruption

Figure 1.

Illustration of train operation organization during interruption

Steps for solving with improved genetic algorithm

Figure 2.

Steps for solving with improved genetic algorithm

Variation curve of fitness values with iteration count

Figure 3.

Variation curve of fitness values with iteration count

Train schedule before model optimization adjustment

Figure 4.

Train schedule before model optimization adjustment

Train schedule after model optimization adjustment

Figure 5.

Train schedule after model optimization adjustment

Comparison of the number of stranded passengers at different stations

Figure 6.

Comparison of the number of stranded passengers at different stations

The number of passengers on board all trains when departing from the station

Figure 7.

The number of passengers on board all trains when departing from the station

The number of passengers disembarking at each station when all trains arrive

Figure 8.

The number of passengers disembarking at each station when all trains arrive

Model-related parameters and symbol definitions

Parameters Symbol definitions Parameters Symbol definitions
|K| Service train set, m, k are train indices WT1 Total waiting time of passengers on board
|U| Set of stations on the line, u, v are train indices WT2 Total waiting time of passengers arriving within the interval between adjacent service trains
Cap Train capacity WT3 Total waiting time of passengers stranded on the platform
H1 Minimum time interval between the departures of two consecutive trains from the same station λu (t) The arrival rate of passengers at station u at time t
H2 The minimum time interval between the arrival of two consecutive trains at the same station λun The arrival rate of passengers at station u during the nth time period
H3 The minimum time interval between the arrival and departure of two consecutive trains at the same station ρu (t) The alighting rate of passengers at station u at time t
tk,udep The departure time of train k at station u Dk,ual The number of passengers alighting from train k at station u
tk,udwell The station dwell time of train k at station u T The minimum total waiting time for all passengers
tmindwell Minimum station dwell time of the train tk,uarr The arrival time of train k at station u
tmaxdwell Maximum station dwell time of the train δz Train reversal time
tk,urun The travel time of train k between stations u and u + 1 Ak,u,von The actual number of passengers boarding train k at station u when train k departs from station u to station v
Pk,u,vwait The number of passengers waiting to board train k at station u, planning to travel to station v, when train k arrives at station u Ak,uon The actual number of passengers boarding train k at station u when train k departs from station u
Pk,uwait The number of passengers waiting to board train k at station u when train k arrives at station u Wk,u,vhall The number of passengers who were unable to board train k at station u before its departure to station u when train k departs from station u
Ek,u,v The number of passengers traveling from station u to station v, who arrived at station u within the departure interval of adjacent service trains k and k − 1 Wk,uhall The number of passengers who were unable to board train k at station u when train k departs from station u
Qk,uin The number of passengers on board train k when it departs from station u

Source: Authors’ work

Station number

Station name Station no. Station name Station no.
Songjiazhuang 1 Hepingli Beijie 13
Liujiayao 2 Heping Xiqiao 14
Puhuangyu 3 Huixin Xijie Nankou 15
Temple of Heaven East Gate 4 Huixin Xijie Beikou 16
Ciqikou 5 Datunludong 17
Chongwenmen 6 Beiyuanlubei 18
Dongdan 7 Lishuiqiaonan 19
Dengshikou 8 Lishuiqiao 20
Dongsi 9 Tiantongyuannan 21
Zhangzizhong Lu 10 Tiantongyuan 22
Beixinqiao 11 Tiantongyuanbei 23
Yonghegong Lama Temple 12

Source: Authors’ work

Dwell time

Station Dwell time (s) Station Dwell time (s) Station Dwell time (s)
1 30 9 40 17 30
2 40 10 35 18 40
3 40 11 35 19 45
4 35 12 60 20 40
5 45 13 30 21 45
6 50 14 30 22 40
7 50 15 40 23 40
8 35 16 60

Source: Authors’ work

Interval running time

Interval Interval running time (s) Interval Interval running time (s) Interval Interval running time (s)
1 122 9 86 17 253
2 80 10 76 18 101
3 128 11 79 19 99
4 95 12 86 20 116
5 78 13 87 21 85
6 76 14 85 22 84
7 84 15 76
8 78 16 122

Source: Authors’ work

Other related parameters of the model

Parameters Value Parameters Value
Train capacity 1,440 people Minimum departure interval 1.5 min
Reversal time 1.5 min Maximum departure interval 3 min
Population size 80 Minimum station dwell time 45 s
Number of iterations 200 Maximum station dwell time 60 s
Crossover probability 0.5 h1 0.1
Mutation probability 0.1 h2 5
k1 0.5 h3 0.5
k2 5 h4 0.1
k3 0.5 Popmax 1
k4 0.3 Popmin 0
η 1 m 1

Source: Authors’ work

Comparison of calculation results

Objective function (s) WT1/(s) WT2/(s) WT3/(s)
Initial 16,157,070.64 2,781,076.80 3,931,921.70 9,444,072.14
GA 15,383,045.94 2,469,551.18 3,530,028.24 9,383,466.52
IGA 13,552,304.17 2,102,633.45 2,263,137.96 9,186,532.76

Source: Authors’ work

Optimization results of departure intervals

Train Before/(s) After/(s) Train Before/(s) After/(s) Train Before/(s) After/(s)
1–2 120 93 8–9 120 168 15–16 120 114
2–3 180 108 9–10 120 139 16–17 120 153
3–4 120 97 10–11 120 158 17–18 180 158
4–5 120 153 11–12 120 127 18–19 120 136
5–6 120 129 12–13 180 131 19–20 120 115
6–7 120 180 13–14 120 137
7–8 180 118 14–15 120 107

Source: Authors’ work

Comparison of optimization results for dwell time

Station Before/(s) After/(s) Station Before/(s) After/(s) Station Before/(s) After/(s)
1 30 56 7 50 47 12 30 45
2 40 47 8 35 46 14 30 45
3 40 47 9 40 47 15 40 48
4 35 46 10 35 48 16 60 47
5 45 46 11 35 45 17 30 52
6 50 46 12 60 46

Source: Authors’ work

References

Binder, S., Maknoon, M.Y., Azadeh, S.S., et al. (2021), “Passenger-centric timetable rescheduling: a user equilibrium approach”, Transportation Research Part C: Emerging Technologies, Vol. 132, p. 103368.

Canca, D., Barrena, E., Laporte, G. and Ortega, F.A. (2016), “A short-turning policy for the management of demand disruptions in rapid transit systems”, Annals of Operations Research, Vol. 246 Nos 1/2, pp. 145-166.

Ghaemi, N., Goverde, R.M.P. and Cats, O. (2016), “Railway disruption timetable: short-turnings in case of complete blockage”, 2016 IEEE International Conference on Intelligent Rail Transportation (ICIRT), IEEE, pp. 210-218.

Jia, B., Zhu, L., Li, S.K. and Liu, J.L. (2022), “Collaborative optimization strategy of short-turning plan and passenger flow control in urban rail transit”, Journal of Transportation Systems Engineering and Information Technology, Vol. 22 No. 1, pp. 124-132.

Lin, L., Meng, X.L. and Song, Z.Z. (2018), “A research on train operation plan under the condition of high-speed railway lines interruption”, Railway Transport and Economy, Vol. 40 No. 7, pp. 7-12. (in Chinese).

Yang, T.Y., Zhao, P., Yao, X.M. and Zhang, P. (2021), “Coordinated optimization model of rail transit delayed train stop-skip pattern and passenger flow control”, Journal of Transportation Systems Engineering and Information Technology, Vol. 21 No. 2, pp. 105-110, 132, (in Chinese).

Zhang, J.J, Hu, W.F., Peng, T. and Yang, C.H. (2021), “Moving block principle-based multi-strategy optimal scheduling method for trains in case of segment blockages”, Scientia Sinica Informationis, Vol. 51 No. 3, pp. 413-429, (in Chinese).

Zhou, W.T., Qian, L., Han, B.M. and Ran, J.L. (2022), “Train rescheduling of urban rail transit under bi-directional interruption in operation section”, Journal of Southeast University (Natural Science Edition), Vol. 52 No. 4, pp. 770-779, (in Chinese), doi: 10.3969/j.issn.1001-0505.2022.04.019.

Zhu, Y.Q. and Goverde, R.M.P. (2020a), “Dynamic and robust timetable rescheduling for uncertain railway disruptions”, Journal of Rail Transport Planning and Management, Vol. 15, p. 100196, doi: 10.1016/j.jrtpm.2020.100196.

Zhu, Y.Q. and Goverde, R.M.P. (2020b), “Integrated timetable rescheduling and passenger reassignment during railway disruptions”, Transportation Research Part B: Methodological, Vol. 140, pp. 282-314, doi: 10.1016/j.trb.2020.09.001.

Further reading

Author, F. and Author, S. (2016), “Title of a proceedings paper”, in Editor, F., Editor, S. (Eds), Conference 2016, LNCS, Springer, Heidelberg, Vol. 9999, pp. 1-13.

Zhow, W.T., Qian, L., Han, B.M., et al. (2022), “Train rescheduling of urban rail transit under bi-directional interruption in operation section”, Journal of Southeast University, Vol. 52 No. 4, pp. 770-779, (Natural Science Edition).

Corresponding author

Rui Wang can be contacted at: 17671610731@163.com

Related articles