Abstract
A station disruption is an abnormal operational situation that the entrance or exit gates of a metro station have to be closed for a certain of time due to an unexpected incident. The passengers’ travel behavioral responses to the alternative station disruption scenarios and the corresponding controlling strategies are complex and hard to capture. This can lead to the hardness of estimating the changes of the networkwide passenger demand, which is the basis of carrying out a response plan. This paper will establish a model to solve the metro station disruption problem by providing optimal additional busbridging services. Two main contributions are made: a threelayer discrete choice behavior model is developed to analyze the dynamic passenger flow demand under station disruption; and an integrated algorithm is designed to manage and control the station disruption crisis by providing additional busbridging services with the objective of minimizing the total travel time of affected passengers and the operating cost of bridgingbuses. Besides, the multimodal transport modes, including metro, bridgingbus, sharedbike, and taxi, are considered as passengers’ alternative choices in face of the station disruption. A numerical study based on the Beijing metro network shows that additional busbridging services can significantly eliminate the negative impact of the station disruption.
1. Introduction
Station disruption, also called station closure, is an abnormal and usually unplanned operational situation in which operators must close the entrances or exits of a metro station for various reasons, such as unexpected incidents or taking steps to avoid overcrowding [1]. In these situations, passengers cannot use the closed station as their departure station or destination stations during the closure time. Generally, disruption means a serious deviation from the planned operations in the rail transit context [2–6]. Station disruption is a typical case of metro disruption, and it can strongly affect both the service and demand of the metro system. At the service level, a station disruption means that the trains, which are planning to stop at the closed station, have to adopt other alternative routes to detour the closed station, e.g., passing the closed station without stopping. At the demand level, the station disruption may cause significant changes in passenger flow demand at the closed stations and nearby stations.
To our knowledge, limited existing studies are available related to the topic of station disruption. Silva et al. [7] proposed a datadriven statistical method to determine the effect of station disruptions on the macroscopic passenger flow demand based on the smart card data from the day of the station disruption and the historical data. This approach was proved to be effective for identifying the impact of a station disruption, but it is not applicable for predicting the impact of a station disruption that does not actually occur or is about to occur. Another effective mathematical model that could capture the behavior of passengers affected by a station disruption was proposed for analyzing alternative disruption scenarios and their likely outcomes in our previous study [1]. But this model ignored the randomness of the passengers’ choice behaviors and the influence of the controlling strategies on passenger behaviors. Some other researchers studied the impacts of traffic incidents on road traffic, which could provide references for this study in metro system [8, 9].
Although we can capture the dynamic influence of the station disruption on passenger flow demand to some extent, how to minimize the severity of potential negative impacts of the station disruption by carrying the affected passengers to the destination with high service quality and low operation cost is still unsolved. Rescheduling timetables plays a limited role in solving the station disruption problem [10–18], because the station disruption defined in this paper does not interrupt the trains running process on the metro lines. But the busbridging service has been recognized as one of the most critical components of effective disruption responses in metro networks. Kepaptsoglou and Karlaftis [19] firstly proposed a methodological framework for planning and designing an efficient busbridging network. Wang et al. [20] explored the demand modeling problem for busbridging disruption in metro services based on the theory of compound Poisson processes and formulated it as a bulk queuing problem involving balking and reneging. Jin et al. [21] introduced a localized metrobus integration approach that aimed at enhancing the urban transit networks resilience to disruption. Based on the previous researches, Jin et al. [22] presented an optimizationbased approach of designing an optimal busbridging service for network disruptions by considering the commuter travel demand at the time of the disruption. But the travel demand was estimated from the historical data. In other words, the dynamic responses of passengers to the disruption were not considered.
In this paper, we attempt to manage and control the station disruption by providing additional busbridging services for passengers. Our contributions mainly consist of the following two parts: a threelayer discrete choice behavior model is developed to analyze the dynamic passenger flow demand under station disruption; and an integrated algorithm is designed to manage and control the station disruption crisis by providing additional busbridging services with the objective of minimizing the total travel time of affected passengers and the operating cost of bridgingbuses. Besides, the multimodal transport modes, including metro, bridgingbus, sharedbike, and taxi, are considered as passengers’ alternative choices in face of the station disruption.
The paper is organized as follows: Section 2 presents the detailed description of station disruption problem and the integrated optimization approach based on agentbased simulation. The case study results are reported in Section 3. Finally, conclusions are given in Section 4.
2. BusBridging Service for the Station Disruption Problem
In this section, we firstly introduce the metro station disruption problem in the metro network. Then, we describe this problem with an optimization approach integrated with a three discrete choice behavior model. In this approach, an optimal busbridging service plan is generated by a binary integer programming model. Besides, the passengers’ response to the busbridging information and its randomness are also formulated in the behavior model. Finally, a solution algorithm is developed to solve the problem.
2.1. Problem Description
Consider a scenario that an emergent event causes a temporary closure of a station in the metro network. We also assume that the metro operator can accurately estimate the closure duration time and then inform the networkwide passengers immediately for simplifying the problem. Not only passengers whose departure station or destination station is the closed station are affected, but also the passengers around the closed station are potentially affected. As illustrated in Figure 1, station is temporarily closed. The affected passengers planning to travel by metro have to alter to other alternative transport modes, such as taxi, sharedbike, bridgingbus, and walking.
Consequently, for passengers, they cannot board or alight the trains passing through without stopping at the closed station . For the passengers with their departure stations as station B, they may choose other alternative stations near the closed station as their new departure stations, wait at the closed station for recovery, or take another transit mode to finish their travels. For the passengers that station is their planned destinations, they have to decide whether to choose an alternative destination station or take another transit mode. The passengers make travel decisions according to their travel experiences and the released guidance information from public transit operators.
For metro operators, they must quickly carry out a responsive plan to replenish the services for the affected passengers. A common strategy is to restore the connectivity between the closed station and the rest of the metro network by running busbridging services. And, in station disruption problem, we consider designing the optimal direct busbridging routes between the closed station and the rest of the rail stations with an objective of minimizing total travel time of affected passengers and the operation cost.
In planning for responses to station disruption, challenges are to estimate the dynamic passenger demand interacting with the operational information and then to decide which stations should be connected with the closed station by busbridging services. Notice that when bridgingbuses are operated between the closed station and nearby stations, the passengers will also decide whether to choose it to finish their travels or not. Therefore, it is a complicated decision problem: which stations among the networkwide stations should be selected to connect with the closed station with the objective of minimizing total travel time of affected passengers and the operating cost of bridgingbuses, considering the dynamics of passenger behavior and demand. To describe and solve the station disruption problem, a hybrid optimization approach based on binary integer programming and passenger behavior simulation is constructed. And, in order to simplify the problem, we constrain our research to the scenario where only a single station is closed.
The basic framework of the model is shown in Figure 2. The objective function of the optimization model is minimizing the total travel time of affected passengers and the operating cost of bridgingbuses. The simulation procedure is used to describe the feedback of the passenger to the bridgingbuses and then to output the dynamic traffic flow for the optimization part. While the urban rail transit system involves largescale passengers and their dynamic travel decisions under station disruption, the passenger behavior cannot be described by the precise mathematical method. Therefore, an agentbased passenger behavior simulation model is considered to describe the behavioral response of the passengers to the station disruption scenario and its corresponding busbridging schemes.
2.2. Model Notations
Main notations used in the paper are listed as follows. All variables are assumed to be integer numbers to satisfy the engineering requirements.
Parameters and Sets G(): directed graph to denote the metro network. S indicates the station set, and A indicates the set of directed arcs. i, j, s: index of a station, i, j, . : index of an directed arc, . M: the vector of transport modes, , and the bus stands for the bridgingbus. : the vector of transport modes except for metro, . m: index of transport modes, m∈M. : the carrying capacity of the transport mode departing from the station i. : distance of the shortest path from station to station by the transport mode m. : travel time of the shortest path from station to station by the transport mode m. b: index of the closed station, . N: the total number of the available bridgingbuses.
Decision Variables P: set of passengers in the current metro system. : set of passengers from station to station at time t, . : index of the passenger with the origin station i, destination station j, and departure time of time t, . : the value of the travel time of urban transport mode m, which can be found in local studies in China. : the comfort penalty parameter of urban transport mode m, and there exits according to our common sense. : the fee for taking transport mode with the shortest path from station to station j. : the generalized cost of taking transport mode from station to station at time considering the travel time and price. : the total number of passengers taking transport mode at station at time t.
Intermediate Variables : indicates that whether operating bridgingbuses from the closed station to station i at time and then back to closed station at time . 1 denotes “YES”; 0 stands for “NO.”
Station Disruption Parameters : the main decision variable denotes whether station should be chosen as the alternative origin or destination station for the passenger p, where : the decision variable denotes whether the transport mode should be the alternative transport mode from station i to station j at time for the passenger p, where : the constant parameter denoting whether station is closed for the passenger p. equals 1 when station is out of service and is 0 otherwise. : the beginning time of the disruption of station i. : the actual ending time of the disruption of station i. : the extra journey time cost of metro transit from station to station at time due to the station disruption:
2.3. Model Formulations
The proposed integrated optimization model includes two parts. The optimization part determines which stations should be set as a bridge with the closed station while minimizing the operating cost of bridgingbuses and the total travel time of affected passengers. However, the behavior model part is used to capture the passenger choice behavior response to station disruption with a threelayer discrete choice model, determining which alternative travel strategy, which alternative origin or destination, and which replaceable transport mode will be chosen by the affected passengers.
2.3.1. Optimization Model
The optimization model is used to decide which stations should be connected to the closed station around the network with a goal of minimizing the operating cost of bridgingbuses and the total travel time of affected passengers. The model can be written as follows:
As indicated in (4), the objective function contains two parts: The first part is the operating cost of bridgingbuses, running between the disrupted station and other stations. The second part is the total travel time of all affected passengers travelling with various transport modes.
Equation (5) defines the function to calculate the total travel time of all affected passengers. In the function, the affected passengers are defined as two types: one is the passenger whose planned origin station is closed, and the other is the passenger whose planned destination is closed. The travel time of two types of passengers is calculated differently.
Equations (6)(7) indicate the constraints of the passenger behavioral choices. Equation (6) shows that an alternative origin or destination station will be chosen for a particular passenger . Constraint (7) indicates that one new transport mode will bridge the invalid origin station and the new replaceable station, or the new destination station to the invalid destination station.
Equations (8)(9) indicate the capacity constraints of multimodal transport resources. Constraint (8) denotes the capacity constraints of the transport mode at the closed station . Constraint (9) denotes the capacity constraints of the transport mode at the closed station s.
Equation (10) shows the decision variables are binary.
2.3.2. Passenger Behavior Model
The lower model is to capture the behavioral responses to the station disruption and its corresponding busbridging strategies and considering the randomness of the choice behavior on alternative stations or transport modes, e.g., bridgingbuses, taxis, metro, sharedbicycles, and walking.
Generate the Optimal Travel Decision. For a passenger , the objective of passenger behavior model is to minimize the generalized travel cost of making a new travel decision for the station disruption. The decision variables are and .
Subject to
Objective function (11) minimizes the generalized travel cost. Equations (12)–(14) indicate the basic constraints of the station disruption situation for a particular passenger, as illustrated in Figure 3. Constraint (12) indicates the generalized travel cost of making a new travel decision if the passenger’s origin station is closed. Constraint (13) denotes the generalized travel cost of making a new travel decision if the passenger’s destination station is closed. Constraint (14) indicates that only one case occurs for a particular passenger, which means that either the origin station () or the destination station () is closed.
(a)
(b)
Equations (15)(16) denote the limits of the decision variables. Equation (15) denotes that another transport mode will bridge the closed origin station and the new departure station, or the new destination station and the closed destination station except the metro. Equation (16) denotes that an alternative station will be chosen to replace the closed one.
Equation (17) indicates the generalized cost of taking transport mode from station to station at time considering the travel time, price, and comfort penalty. is defined as a very large constant. When , it means that the generalized travel cost will be very large if the passenger still chooses the bridgingbus as his alternative transport mode. It also indicates that only when bridgingbuses are operated between stations and the closed station b, the bridgingbuses can be available for passengers (). It can be seen that the passenger’s travel behavior is affected by the upper decision variables of controlling measures.
Equation (18) denotes that the time interval for a particular passenger to make a decision is only within the period of the station disruption. In other words, the model is applicable when the affected passenger has just been faced with the station disruption.
ThreeLayer Discrete Choice Behavior Model for Alternative Travel Strategies. In a previous section, a behavior optimization model is proposed to generate an optimal travel strategy with a minimum generalized travel cost for the case of a station disruption. But, in an actual disruption scenario, the passengers’ travel experience partially or completely turns into failure, and the necessary information for decisionmaking is usually lacking. Consequently, passengers cannot make a rational or optimal travel decision. Therefore, the randomness of passenger travel behavior should be considered in the newly proposed model.
Based on the stochastic utility theory, a threelayer nestedlogit model for passengers’ choice behavior under station disruption is constructed, which is illustrated in Figure 4.
In the proposed threelayer nestedlogit model, the first layer is to make a trip decision among quitting the metro journey (labeled as the event ), waiting for the recovery of station service (labeled as event ) and taking an alternative metro station (labeled as event ); the second layer is to make travel decisions among the alternative stations set for the passengers taking an alternative metro station, or choices among various transport modes to continue the journey; the third layer is to make choices among available transportation modes to bridge the new alternative metro station and the planned origin or destination station. The multimodal transport resources contain bridgingbus, taxi, shared bicycle, and walking, which are labeled as m1, m2, m3, and m4, respectively. Accordingly, the general cost of each behavior is , , and , respectively, for events , , and . Take a passenger as an example, and we can have these remarks.
(a) For the passenger who gives up the metro journey, the choice behavior is that the passenger will take the planned destination as his new departure station in the situation of and , or take the planned origin station as his new destination in the situation of and . Then, the probability of choosing the model can be calculated as
So can be calculated based on (12) ~ (14):
(b) For the passenger who waits for the recovery of station service, the choice behavior is that the passenger will take the planned origin station as his new departure station. Once the destination is closed, it is hard to decide which station should be chosen to wait at. So waiting for the recovery of station service can be one option only in the case that one’s planned departure station is closed. So can be calculated based on (12) ~ (14):
(c) For the passenger who takes an alternative metro station, the choice behavior is that the passenger will take an alternative station as his new departure station or destination station. If the station is selected, the generalized travel cost of the passenger choosing the transport mode can be calculated asFirstly, in this situation, the probability of choosing the model can be calculated as
We label as the generalized travel cost of the passenger who chooses the station as his new alternative departure or destination station:
Then, the probability of choosing the station can be calculated as
The probability of each passenger behavior decision is
The joint probability of passenger choosing a station as an alternative station can be calculated:
Furthermore, the joint probability of choosing transport mode can be
2.4. An Integrated Solution Algorithm
2.4.1. Genetic Algorithm
The decision variable in the upper model is , indicating whether the bridgingbus should be operated between the closed station and the station i∈S at the time interval t. 1 denotes “YES” and 0 stands for “NO.” Therefore, in the genetic algorithm, the traditional binary coding method can be used to design a feasible solution as the corresponding gene sequence or chromosome. A chromosome is represented by a vector as indicated in (31). Besides, the population is represented by in (32).
Generally, the genetic algorithm consists of three important operators: the crossover operator, the mutation operator, and the selection operator. We assume that the chromosome intersects with a certain probability and randomly selects the intersection with a certain crossover probability . In the algorithm, the crossover method is given as shown in Figure 5. Since the chromosome is composed of a specific gene fragment which cannot be divided only the same position of the gene fragment will be crossed.
The mutation operator can help avoid local optimal solutions. In this paper, we assume that each chromosome will be a certain probability of mutation operation and randomly select the mutation position, as shown in Figure 6.
The selection operator is to ensure the best individuals can survive in the next generation. And it is designed based on the fitness function, defined as , where is defined by formula (4). The more the fitness is, the more likely the best individual can survive. We adopt the Roulette method to determine whether each individual can survive in the next generation or not.
2.4.2. Solution Algorithm of Passenger Behavior Model
The general simulation procedure can be attributed to Yin et al. [1]. A new behavior model solver in the simulation algorithm is developed for the threelayer discrete choice behavior model, which is summarized as shown in Algorithm 1:

2.4.3. The General Solution Algorithm
In this paper, a genetic algorithm and a simulation based method are integrated to solve the proposed model. The genetic algorithm is used to solve the binary integer programing model and the simulation procedure is to generate different types of the passenger flow demand. And the basic simulation procedure is attributed to Yin et al. [1]. The general algorithm framework for solving the problem of station disruption controlling optimization is described as follows.
Step 1 (initialize). Initialize the parameters and set the maximum number of iterations to G, and the population size is . Initialize the disruption’s time duration as [1,n]. Set the current time interval t = 1.
Step 2. Set an initial set of busbridging scheme or initial population. Set the current generation .
Step 3. Solve the lower model. For a particular busbridging scheme, the simulation model and algorithm are called to capture the passenger travel behavior and each passengers’ travel time. And calculate Z.
Step 4. Solve the upper model. According to each passengers’ travel time, the fitness of each chromosome in the initial population can be calculated and through the following operations, to generate the upper optimization of the busbridging scheme strategy .
Step 4.1. Mutation operation: using the mutation operator to create a new chromosome.
Step 4.2. Cross operation: cross the operation on the chromosome.
Step 4.3. Selection operation: using the selection operator to ensure that the optimal individual genetic to the next generation.
Step 5. Verify the end condition 1. If the number of iterations is greater than , then the optimal solution of the model will be obtained; otherwise, , and the iteration will return to Step 3.
Step 6. Verify the end condition 2. If , end; otherwise, make , and return to Step 2.
3. Numerical Example
In this section, we present the application of suggested model and the solution algorithm to an example urban metro network based on real and virtual data. Figure 7 illustrates the example station disruption scenario, which contains 2 directional lines and 2 transfer stations extracted from Beijing Subway network, which is one of the busiest metro systems [23]. Figure 8 illustrates the practical train timetables. We assume that the Xizhimen station is closed. The cost of taking alternative transport modes is given in Table 1.
In order to compute and test the model more easily, we assume that each OD volume shares the same distribution, as shown in Figure 9. Different values of can be investigated for different scales of passenger demand. And other values of main parameters are defined in Table 2 for simplification.where t denotes the tth time interval of 15 minutes and indicates the maximum passenger demand of each two different stations per time interval, indicated in Figure 9.
3.1. Passenger Flow Demand Analysis under a Station Disruption
3.1.1. Analysis of Individual Behavioral Choices
Passengers’ Behavior Choices under Different the Closure Duration. For a particular OD from the Xizhimen station to the Pingguoyuan station, Figure 10 shows that, with increasing closure duration, the passengers will more likely alter their planned origin stations or give up their metro journeys, and fewer passengers will wait for its recovery of station service.
Alternatives of Transport Modes. Figure 11 shows the passengers’ choices in multimodal transports. With increasing the bridging distance up to 6 km, more passengers prefer to choose the taxi as their alternative transport mode.
(a)
(b)
Alternatives of Origin Stations. Figure 12 shows the top stations that passengers prefer to choose as their alternative origins when their planned origin station Xizhimen station is closed for different durations. The Fuxingmen, which is near the Xizhimen station, is always chosen by most passengers. The result also shows that fewer passengers will wait at the closed station and more passengers will give up their planned origin stations to nearby stations if the closure duration becomes larger.
(a) 15min
(b) 30min
(c) 45min
(d) 60min
3.1.2. Passenger Flow Demand at Affected Stations
Passenger Flow Dynamics at the Closed Station. In a 60 minutes’ station disruption, the number of entries at the Xizhimen station is changing over time, as is shown in Figure 13. It also shows that when the disruption comes to an end, the passengers waiting at closed station will gradually enter the closed station, which leads to more entries than normal. And, consequently, it takes longer than the closure duration itself for the entries to recover to a normal level, as shown in Figure 13.
(a)
(b)
Passenger Flow Dynamics at the Nearby Stations. The disruption of Xizhimen station results in a significant increase in the number of entries at the nearby stations, especially at the Fuxingmen station, as is shown in Figure 14.
(a)
(b)
3.2. Optimization Results
Optimization Results under Different Disruption Durations. As is indicated in Figure 15, a set of disruption scenarios are numerically computed for a determined OD volume of V=240, and the duration of station disruption ranges from 0 to 120 minutes. The results show that when the duration of station disruption becomes larger, the total travel time will be more significantly reduced if the bridgingbuses are operated. The improvement can reach up to 49.38% under a station disruption of 120 minutes. That is because when the station is closed for a longer time, more passengers will choose the bridgingbuses or other transport modes instead of waiting for recovery of the station service, which leads to a better controlling effect. It also indicates that, in the minor disruption scenario, there is no need to operate the bridgingbuses.
The blue line in Figure 15 also shows that our developed solution algorithm is timeexpensive. It costs up to 305 minutes in an only 120 minutes’ station disruption scenario in the small numerical study, which seems not to be efficient enough for largescale network cases or online emergent situations. That is because our proposed model must calculate all behavioral responses to the station disruption for every affected passenger in each iteration in the solution algorithm.
Optimization Results under Different Passenger Flow. As is indicated in Figure 16, another set of disruption scenarios are numerically computed for a disruption duration of 60 minutes, and the OD volume of ranges from 60 to 420. The results show that, with increasing the passenger flow demand, more total travel time will be significantly reduced if the bridgingbuses are operated.
4. Conclusion
An integrated optimization approach based on behavior model is proposed to manage and control the disruption crisis by providing additional busbridging services with the objective of minimizing the total travel time of affected passengers and the operating cost of bridgingbuses. Among them, a threelayer discrete choice behavior model is proposed to analyze the passenger flow demand under station disruption. Furthermore, the multimodal transport modes, including metro, bridgingbus, sharedbike, and taxi, are considered as passengers’ alternative choices in face of the station disruption. A numerical study based on Beijing metro network is conducted finally. Some conclusions can be drawn as follows.
(a) The effects of a station disruption with different durations on passenger flow demand are analyzed, and the results show that, with increasing closure duration, more passengers will alter their planned origin stations or quit their metro journey, and fewer passengers will wait at the closed station for its recovery. And, with increasing the bridging distance up to 6 km, more passengers prefer to choose the taxi as their alternative transport mode.
(b) The optimization results show that additional busbridging services can eliminate the station disruption impact significantly, especially with the increase of closure duration and the passenger flow demand. The improvement can reach up to 49.38% under a station disruption of 120 minutes and a determined OD volume of V=240 in our numerical study.
(c) Our developed solution algorithm is timeexpensive and seems not to be efficient and applicable enough for largescale network cases or realtime emergent situations. A more efficient algorithm should be developed in our future studies.
Data Availability
All data generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by NSFC (nos. 71771018, 71621001, 71701013, and 71601017), the Fundamental Research Funds for the Central Universities (no. 2018JBM074), Research Foundation of State Key Laboratory of Rail Traffic Control and Safety (no. RCS2018ZT004), and Technology Foundation of Beijing Jiaotong University (no. 2018RC010). The authors thank these agencies.