Simulation-Based Dynamic Passenger Flow Assignment Modelling for a Schedule-Based Transit Network

The online operationmanagement and the offline policy evaluation in complex transit networks require an effective dynamic traffic assignment (DTA) method that can capture the temporal-spatial nature of traffic flows. The objective of this work is to propose a simulation-based dynamic passenger assignment framework and models for such applications in the context of schedule-based rail transit systems. In the simulation framework, travellers are regarded as individual agents who are able to obtain complete information on the current traffic conditions. A combined route selection model integrated with pretrip route selection and entrip route switch is established for achieving the dynamic network flow equilibrium status. The train agent is operated strictly with the timetable and its capacity limitation is considered. A continuous time-driven simulator based on the proposed framework and models is developed, whose performance is illustrated through a large-scale network of Beijing subway. The results indicate that more than 0.8million individual passengers and thousands of trains can be simulated simultaneously at a speed ten times faster than real time. This study provides an efficient approach to analyze the dynamic demand-supply relationship for large schedule-based transit networks.


Introduction
Urban rail transit has developed rapidly in China during the last ten years.By the end of 2015, there were 26 cities operating the rail transit with a total length of 3,618 km [1].With the travel demand growing radically, some intractable issues emerge for the operation management, such as the recurrent congestion in peak hours and the train schedule construction for the travel demand under extremely unbalanced conditions.It has been conjectured that the development of Dynamic Traffic Management Systems (DTMS) can be a feasible approach for addressing these prominent problems.However, assessing the benefits and influences of such systems are difficult as these strategies are highly dependent on the travel behavior of individual passenger in response to travel information and control actions.Hence, it is necessary to develop useful methods and tools for assessing the DTM policies and control strategies.
Dynamic traffic assignment (DTA) models are regarded as a valuable tool to evaluate the DTM performance for its capabilities in capturing the dynamic nature of traffic flows and describing the formation and propagation of traffic congestions.DTA models can be divided into two categories: the mathematic-based DTA models and simulationbased DTA models.Mathematic-based DTA models, such as optimization programmes [2], variational inequalities [3], and optimal control [4], have solid theory foundations and perfect analyzability.However, few of them are implemented in real-world traffic networks because of their limitations in computation efficiency and complex parameter requirements.Simulation-based DTA models are considered more suitable for real-world applications [5,6].There have been several DTA simulators developed and applied to road networks, such as DYNAMIT (https://its.mit.edu/software/dynamit) and DYNASMART (http://mctrans.ce.ufl.edu/featured/dynasmart/). Unfortunately, the previous models and tools are mainly developed for road traffic networks.The traffic flow characteristics in transit system are much different within road traffic, where the passenger flows are highly nonlinear.Hence, DTA models in road traffic systems are not suitable for transit system.The DTA tools for the specific schedule-based rail transit networks are underdeveloped.
In the sphere of transit assignment, the DTA models also can be divided into two categories, the frequencybased models and the schedule-or timetable-based models [7].Frequency-based approaches consider services in terms of sets of lines, where the run scheduled times are not considered explicitly, while the schedule-based approach refers to services in terms of runs, using the real vehicle arrival/departure time to obtain attributes that can be explicitly considered in the run choice.With the development of transit modelling, the DTA approaches have gradually changed from frequency-based models to schedule-based models.The capacity limitation of vehicles is a critical point that should be considered in transit systems, because passengers can only travel through sections by trains that have strict capacity constraints.However, previous approaches, whether frequency-or schedule-based, have mainly focused on line capacities not vehicle capacities [8].The "fail-to-board probabilities" is usually used to describe the approximate congestion conditions of vehicles [9,10].In this work, the train is seen as an individual agent and its capacity constraints can be strictly described by limiting the number of boarding passengers.Therefore, queuing processes for passengers waiting on platforms can be captured.For the schedule-based transit assignment models, the space-time or diachronic graph is usually used to represent the transit services network, which contains the service, demand, and access/egress subgraphs [11][12][13][14].However, this kind of spacetime network is extremely complex when transit lines and transit services (train runs) are large.It is difficult to find a shortest path quickly in the networks.In order to enhance the computation efficiency of the DTA models, a two-stage path choice process is constructed in this work.First, search the  shortest paths for each OD pair without considering the passenger flows and train services.This process seeks to establish the interrelationship between links and ODs.Second, update the cost of each route with the simulation clock advancing by considering the real-time passenger flows and services.When passengers enter the transit system, they can choose the "best" path immediately through the current traffic conditions.The two-stage path choice process has high efficiency because the speed for updating route cost is much quicker than searching a new route.
This work aims to propose an efficient and practical simulation-based DTA framework and models for the schedule-based rail transit network.The characteristics of the proposed approaches include (1) describing the travel behavior at the level of individual passengers and presenting the detail travel processes for passengers travelling through the network, such as walking within the station, waiting, transferring, boarding, and alighting; (2) integrating the pretrip and enroute path choice behavior together, which not only allows describing the route selection processes in normal conditions but also enables to capture the enroute switch behavior under exceptional conditions; (3) the capacity which is considered by the accurate train load capacity instead of the line capacity and the train runs which can adapt to different types of timetable; (4) the simulation models having high computational efficiency for large scale (large number of passengers and trains) rail transit networks.
The remainder of the paper is organized as follows.The simulation-based modelling framework for DTA in the schedule-based transit network is given first.Following are the major models under the simulation framework, mainly including passenger generation, route selection, and network-loading models.Then, model implementation in the Beijing subway networks are presented in Section 4. Finally, Section 5 presents the conclusion and future works.

Modelling Framework
Generally, any simulation-based DTA systems consist of two critical components: a route selection module and a networkloading module [16].The route selection module describes the principle for travellers choosing their route and then determines the macroscopical flow status on the network.The network-loading module attempts to represent the temporalspatial evolution of traffic flows on the network after the routes of passengers are determined.Figure 1 presents the simulation-based DTA framework and model structure for a schedule-based transit network, where the major parts are highlighted.
In the simulation framework, travellers are regarded as individual agents who can obtain the current traffic conditions on the network and always select the best path.Trains are operated strictly with the preset timetable, which allows considering the strict arrival/departure time.The simulation processes for obtaining the macroscopical network flow status are to describe the whole travel procedures of all individual passengers in microscopic view, such as egress/access, waiting for the train, and transferring through passageways.The detailed introduction of these two critical components is provided in the next section.Additionally, the network structure and the travel demand generation that produces the input for DTA models are important and will also be introduced in Section 3.

Route Selection Module.
There are several types of route choice models for different applications, such as equilibriumbased or information-based, pretrip or enroute, and withinday or day-to-day.Detail descriptions of passenger choice models can infer to the review of Szeto and Wong [17].In this work, the within-day dynamic user equilibrium (DUE) is used to describe the principle for traveller route choice, which can be seen as follows: when travelling, travellers should always select their "best" path at each decision point, and the cost of used routes is no larger than the unused routes.In a rail transit network, passengers usually have no secondary route when they are at the enroute decision points (transfer stations) because of the low connectivity of the network.Hence, the pretrip user equilibrium principle can be suitable in normal conditions.However, passengers should change their routes in some special conditions, such as switching their routes when they cannot go to the end for incidents.In order to enhance the capability of models in describing the passengers routing behavior under different situations, a combined

Route selection
The k-shortest path route selection method which integrated the pretrip and enroute models is used in this work.The pretrip route selection model which follows the dynamic user equilibrium is established under normal conditions.Moreover, the enroute path switch model based on a bounded rational rule is applied for special situations.This combined model is consistent well with the path selection behaviors of rail transit passengers.
Though the rail transit network has low connectivity, there may also exist alternative route in an OD pair.The -shortest paths are used to describe the discrepancy in route choice behavior of different passengers.How to find the "best" route immediately for each traveller is of great influence on the computation efficiency of a simulation system.For a large-scale network, it is almost impossible to search the routes online.A two-stage path choice process is established in this work, which contains (1) initial -shortest paths search and (2) the real-time route cost update.The process of initial -shortest paths search is to establish the interconnections between links and OD paths, and the route cost updating procedure aims to describe the real-time traffic conditions on the network.The route cost will be refreshed in a small preset time interval (such as two minutes).When a passenger enters the station, he/she can choose the "best" path quickly by considering the approximate real-time traffic conditions.Hence, the flow state of the network is not a rigorous but an approximate user equilibrium status.

Network-Loading Module.
Once the traveller's route is determined, the network-loading module is to simulate the movement of travellers through the network, with the output of dynamic flow distributions of nodes (stations/platforms) and sections (trains).In this work, the continues-time-driven simulation approach is utilized to simulate the processes of all passengers moving on the network.The networkloading module contains three major components: (1) passenger walking within the station, such as access, egress, and transfer; (2) train movement operation; and (3) interactions between passengers and trains on platforms.The detailed description of these components will be presented in Section 3.

Network and Train Schedule.
The rail transit network is represented by a directed graph  = (, , ), where  is the node set,  is the link set, and  is the set of transit lines.The nodes consist of two types: source and sink nodes (entrance and exit gates) for passenger arriving and leaving, and transit nodes (platforms) for train stopping.A station is composed of several interconnected nodes including the source/sink nodes and the transit nodes.According to the connected node type, the links connecting two neighbouring nodes are divided into two types.The first one is the walk link for passenger walking and transferring and the second one is the transit link for a train running between two stations.A transit line is a fixed path along which vehicles periodically run, corresponding to the train route.Note that more than one transit line may exist on a physical line, such as for long and short train routes.Figure 2 presents the topology structure for rail transit networks.
Let  in  be a transit line where   is a set of transit nodes and   is a set of links,  = (  ,   ).Denote  , and  , as the th transit node and th transit link in the transit line ,   , is the arrival time for vehicle  at th transit node (station), and   , is the departure time for vehicle  at th transit node.The timetable  , for vehicle  running on transit line  can be represented by groups of arrival time   , and departure time   , ,  , = {(  , ,   , )}.The dwell time   , for vehicle  is given as (1).The relationship among arrival time, departure time, and link running time is shown in (2).Table 1 provides an example of a train schedule (timetable).where ℎ , is the time for a vehicle running through link  and link  connects nodes  , and  ,+1 .

Individual Passenger Generation.
There are two types of traffic demand information which can be used as input for generating individual passengers in the DTA simulation system.The first one is the time-dependent origin-destination (OD) matrixes that are usually for online traffic management applications.The method of dynamic OD estimation and prediction can be inferred to the works of Yao et al. [18,19].
The second one is the automatic fare collection (AFC) data that is primarily for offline policy evaluations.The generation processes in the simulation system should be as consistent as possible with the actual passenger flow characteristics in arrival time and spatial distributions.When a passenger is generated, the origin station, destination station, and departure time (arrival station time) should be produced.Suppose that passengers arrive uniformly if the time span  OD is very short, for example, no more than five minutes.Then the passenger arriving rate    of source node  (belonging to station ) in time interval  can be calculated as

Generation from
The number of passengers arriving at node  during a simulation time step is where Δ represents the time step of the simulation clock and [] represents the down rounding computation.However, while the simulation time step is very small (such as 1 second), there may not be a single passenger arriving at the station in a simulation step.Referring to other related simulation works [20], we assume that the time intervals between two consecutive passengers follow the negative exponential distributions.Hence, the probability of a passenger arriving in time   is as follows: where   represents the arrival time of the previous passenger.When a passenger is generated at station , his/her destination can be determined from the trip distribution fractions calculated from the O-D tables.The probability of selecting node  as the destination station is 3.2.2.Generation from Automatic Fare Transaction Records.
The AFC system has been used widely in public transport and records accurate trip information for each passenger, including origin and destination stations, check-in and check-out times, and card type.The AFC data provides a new and massive data source for traffic planning and management.When the AFC records are used as the input for the simulation system, the passenger agents can be generated accurately whose arrival patterns follow the nature of the actual travel demand.Unfortunately, it can only be used to analyze the historical traffic status or evaluate the traffic policies corresponding to offline applications.Table 2 presents the data structure of the AFC record information used for generating passengers.

Route Selection.
For the particular structure of the rail transit network, passengers can only change their route at transfer stations, and usually they are unable to choose a secondary route once they start their trip.Hence, in this work, a combined route selection procedure is used to achieve the dynamic network equilibrium conditions: passengers follow the predetermined equilibrium principles under normal circumstances, and a bounded rational rule is adopted for passenger route changes under emergencies.The advantages of the combined route selection model are (1) high computational efficiency, (2) conformance to the route choice behavior for rail transit travellers, and (3) ability to capture the behavior change under special conditions.

Equilibrium-Based Pretrip Path Choice.
Assume that all travellers can obtain all of the information on the current traffic conditions and each traveller will select the best (usually the lowest cost) path.The travel pattern of pretrip dynamic equilibrium is defined as follows: for all travellers who leave their origin at any time in any O-D pair, the costs of any used routes are equal and minimal and smaller or equal to the cost of unused routes.The route cost is represented by a generalized cost function which is composed of four parts: (1) access and egress walk time; (2) on-board time; (3) waiting time; and (4) transfer time.
The walk time contains access and egress walk time, which is calculated by the distance corresponding to the walk link and the average walking speed of travellers.The total walk time  walk can be represented as where  in walk represents the access walk time and  out walk is the egress walk time.
Normally, the on-board time for passengers travelling through sections is constant, which can be fixed with the train schedules.Define  in as the total on-board time of passengers travelling through all sections within a route.The waiting time includes three components: (1) waiting time at the origin station; (2) waiting time at the transfer station if a transfer is required; and (3) the delay waiting time due to overloaded trains.The average waiting time at an origin/transfer station is equal to half of the train headway time corresponding to travel line.However, the rail transit system has strict capacity restraints; some passengers cannot board the first train that arrives because it is overloaded.This delayed waiting time should be considered if there is no sufficient capacity.Assume that passengers who fail to board the first arriving train will wait for the next train.Thus, the total waiting time  wait can be represented as where  ori wait represents the waiting time at the origin station;  trans wait is the waiting time at the transfer station;  = 1 if there is transfer, otherwise,  = 0; and  is the weighted factor for delay waiting time  delay wait .When a passenger is unable to board the first train, he/she needs to wait the headway time for the next train.The delay time changes greatly for different stations and should be frequently updated with the simulation clock.The delay waiting time is formulated as where   is the number of passengers who are delayed;   represents the numbers of passengers who are able to board the first train; and ℎ  is the headway time of the transit line .
Define   as a set of the  shortest routes in the O-D pair (, ),   represents a route, and   ∈   ;    is the generalized cost of route   .The time-dependent generalized cost    of route   during time step  can be formulated as follows: where   ,   ,   , and   are the weighted factors for the cost corresponding to different travel processes and    transfer represents the total transfer walk time in route   .These weighted factors can be estimated by the maximum likelihood method through travel surveys.There is no doubt that passengers are unable to know the exact cost of each route but merely estimate the cost based on their experiences and make trade-off decisions.The deviation between the practical travel time and the perceived time makes a difference in route choice.In this work, the Multinomial Logit (MNL) model which has been a popular method for the probabilistic choice is used to describe the uncertainty in the route selection.Suppose the deviations follow the Gumbel distributions.Then the probability for a passenger selecting the route   during simulation time  can be described as where  min   () represents the minimum cost of all alternative paths in the O-D pair (, ).
To improve the simulation efficiency, the initial  shortest paths of all O-D pairs will be searched and stored in the initialization procedure without considering passenger flows.The route cost will be updated in a preset short time interval (such as two minutes) with the simulation clock running.Therefore, passengers could choose the best path when they start at the origin station, according to the approximate realtime traffic conditions.

Rule-Based Enroute Path Switch.
Because of unexpected events such as incidents and disturbances, the cost of the originally selected route could significantly change after the passenger leaves the origin station.The influenced passengers must switch their routes or cancel their trips, as shown in Figure 3. Hence, the behavior change of passengers in response to special conditions should be added to the simulation framework.Based on the bounded rational model proposed by Mahmassani [21], a modified model is constructed to describe the switch rule.If the gains from a route change satisfy and suffice for a specific value, passengers should switch their routes; if the cost of an alternative route exceeds the highest expected value, passengers will cancel their trip and choose other travel modes; otherwise, passengers should keep on waiting for the traffic recover.The switching rule is formulated as where  represents an indicator variable equal to 1 when a passenger switches from the current route to the alternate, −1 when the passenger cancels the rail transit travel, and 0 if the current path is maintained.TTC() is the current travel cost of the original route from decision node  to the destination, and TTC  () represents the historical travel cost of the original route from decision node  to the destination; TTA() is the travel cost of the alternative route from the decision node to the destination;  and  are the relative difference threshold values; in this work,  = 0.3 and  = 2.0; and  min is the absolute minimum travel cost advantage needed for a route switch.The values of  and  can be estimated from traffic survey and disparate for different transit systems.
From the example shown in Figure 3, it can be found that the event locations are important for a passenger behavior change.In case 1, passengers at decision node  can switch to the alternative route if the travel cost satisfies the condition in (12).However, passengers have to cancel their travel in case 2. Note that although the train movement in our framework is operated strictly by the timetable, the influence of special events can also be evaluated by inputting different train schedules under the same travel demand.

Network Loading.
The function of network loading is to move the individual passengers through the network, producing the output of the time-varying node and link flows.Compared to the network-loading processes for an urban road network, the differences for rail transit include the following: (1) travellers can only move through sections by train, (2) the travel time for transit links (sections) does not change for congestion, and (3) the passenger boarding and alighting interaction is complex; travellers may fail to get on a train that has reached its capacity.For a rail transit system, the network-loading processes contain three major components: passenger walking within the station, train movement control, and the interaction of passenger boarding and alighting.

Queue-Based Model for Passenger Walking and Waiting.
There are several microscopic models for simulating the pedestrian walk behavior in a station, such as the cellular automata model [22] and the social force model [23].However, these models require high computational resources and are not suitable for large-scale network simulation.To offer a computationally efficient traffic simulator, this work uses a simplified queue-based model to describe the walking and waiting behavior for passengers within stations.
Define a walk link (access or egress) as a route for passengers walking from the gate to the platform or back and transfer link as a route from one platform to another.Then, a station can be represented by interconnected walk links, transfer links, platform nodes, and gate nodes.To capture the passenger flow congestion, virtual stack queues related to the walk links were constructed and named entrance queues and exit queues, as shown in Figure 4.The stack queue has no storage capacity limitation.Here, the pass capacity of the link is used to capture the queuing processes of passengers.The pass capacity of a walk link is set to be the minimum pass capacity of all key points within the route, such as gates and downstairs, which can be calculated by referring to the station design handbook.A more realistic queue with spatial storage capacity is constructed for the platform node, which attempts to capture the passenger waiting process for trains.The storage capacity can be fixed by the width and length of the platforms and the congestion passenger density.The connection between physical spatial queues and virtual stack queues can be further reviewed in the work of Hurdle and Son [24].
Newell's simplified model based on cumulative inflow and outflow counts is used to capture the passenger congestion formation and dissipation procedures in queues [15], which is illustrated in Figure 5. Entering the station process can be used as an example to describe the passenger queuing procedure.For a passenger  added to the entrance queue, his/her arrival time is   (), and the departure time will be   (),   () =   () + FFWT under free flow conditions.The free flow walk time (FFWT) can be computed by the link distance and average walking speed.When the simulation clock advances to   (), if the outflow capacity (platform storage capacity) is available and passenger  is at the end of the entrance queue, this passenger can move to the platform waiting queue for trains; otherwise, this passenger must stay in the entrance queue and wait for available outflow capacity.The congestion delay time (CDT) for passenger is   () −   ().
Define  as the number of entrance queues in a station;  as the index of entrance queues,  ∈ ;  pass  as the unit pass capacity of queue  during a simulation interval;  store as the storage capacity of the platform;  , as the number as the transfer flow from the entrance queue to the wait queue.
Based on these definitions, the calculation procedure for a passenger entering the station is shown in Algorithm 1.
The egressing and transferring processes are similar to the entering procedures.Additionally, the number of boarding passengers determines the outflows of the wait queue on the platform.Then, the status of the wait queue changes in a discrete nonlinear way, which will be described in next section.

Train Movement Control. According to the report from Beijing Mass Transit Railway Operation Corporation
Limited, the average punctuality rate of trains is more than 99.7%.It is reasonable to assume that trains run strictly with the scheduled timetables.The train is generated at the first station of the transit line.Suppose that a train runs on transit line .The section  , of transit line  connects the nodes  , and  ,+1 ,  , = (   ,    ) and  ,+1 = ( +1  ,  +1  ), where (, ) represents the node location.Figure 6 shows the schematic diagram for updating the train position.The train location can be calculated as follows: where (  ,   ) represents the train position in time  and Δ and ℎ , have same meanings as before.
The train position within a transit link can be computed by the coordinate transition from the corresponding start and end node locations.If a train stops at a platform for passengers to board and alight, the train position will not be updated until it departs.The train movement may be disturbed for emergencies, such as incidents and congestion, which means the train operating processes should be controlled dynamically.However, controlling the train operating processes from the microscopic view is complicated and challenging.Future works are in progress to incorporate the dynamic train control process based on the moving block train control system into our framework.

Boarding and Alighting
Interaction. Figure 7 presents the interaction flowchart of passenger boarding and alighting when a train arrives.The interaction follows the principle of "first-alight-then-board." Assume that the preset dwell time is sufficient for passenger boarding and alighting and the passengers who fail to board the first train will wait for the next train.
When a train arrives at a station, all of the passengers on the train should decide whether to continue riding or alight, according to their O-D and preselected route.For each passenger, if the stop is (a) the destination station or (b) the transfer station, the passenger should alight from the train.After passenger alighting has finished, the passengers waiting on the platform should determine whether or not to board.Because different transit lines may share the same station, it is important for passengers to board the correct train.Passengers who can board the train should satisfy the following: (i) The next stop of the train is the same as the next stop for the passenger to make sure the train is going the correct direction.(ii) The train will stop at the target stations for passengers, which represents the destination station or enroute transfer stations.(iii) The train has sufficient capacity for the boarding passengers.
Suppose a train numbered  arrives at the station .Define   as the train loading capacity,  −1  as the number of passengers on the train when the train departs from the previous station  −1 ,    as the number of passengers waiting on the platform, and    as the number of passengers needing to alight.Using these definitions, the number of passengers who can board is where    represents the number of passengers who can board.If    <    , some passengers should wait for the next train, and the delayed passengers    are

Simulator Description.
A simulator based on the proposed models was developed in the C# programming language using object-oriented and continuous time-driven approaches.The tool was executed on an Intel PC under Windows 7 with a 3.8 GHz CPU and 4 G RAM. Figure 8 shows the main interface of the simulator.The Beijing rail transit, one of the busiest transit systems in the world, is applied to test the models and simulator.It consists of 344 stations and 18 lines with a total length of 631 kilometres in 2016.The ridership is approximately 10 million trips (each transfer is counted as an additional trip) and 7,539 trains operate daily.The minimum train headway time in peak hours reaches 103 seconds [25].The headway and dwell times are measured in seconds in the operation of rail transit.Hence, the time step of the simulation clock is set to be 1 second and the triggering time for iteration is 100 milliseconds, which is ten times faster than the real time.Test results show that (1) the traffic flow status for a large-scale transit network can be simulated with high efficiency and only 120 minutes is needed to simulate a whole day (from 5:00 to 24:00) and (2) more than 1,000 trains and 0.8 million individual passengers can be simulated simultaneously.The main capabilities of the tool include the following:   3.
There are two indicators used to quantify the model accuracy.The first one is the average relative deviation ARD  between the simulated travel time and the actual travel time for each passenger, as shown in (16).The other one is the average relative deviation ARD  between the simulated section flows and the flows from the official revenue clearance centre, which is calculated by (17).
where   means the simulated travel time for passenger ; where   () represents the simulated section flow for section  during time interval ;    () is the assignment section flow from the revenue clearance centre; and  is the number of sections.
Figure 9 presents the ARD  distributions from 80,000 randomly selected travellers.We can observe that the deviation for most of the passengers is about 5%∼25%.The main reasons for deviation are the random variable influence in the MNL-based model that causes the wrong route choice and the different walking speed for passengers when they go through the passageways.However, the weighted average relative deviation is about 17%, which are still acceptable for a real-world application.
Figure 10 provides the time-dependent ARD  deviations from 05:00 a.m. to 12:00 a.m.We can find that most of the  deviations are less than 6%.Therefore, it can be concluded that the proposed models have an acceptable accuracy and are suitable for DTA application in the schedule-based rail transit networks.

Simulator Applications.
Two case studies are presented to show the capabilities of the simulator.The first one shows the descriptive capability for describing the spatial-temporal passenger flow patterns, and the second one provides the evaluation capability for the event influence.4.
In Beijing rail transit, the direction for lines is defined as the up direction is from east to west and from south to north for straight lines and counterclockwise for the circle lines; the down direction represents the opposite directions.Figure 11 shows the dynamic network passenger flow status during morning peak hours (7:00-9:00), where  means the capacity utilization rates which corresponds to different  colours.We can see that the up direction of Line 1 and the down direction of Lines 5 and 6 and BT are the most congested lines in the Beijing subway during the morning peak hours.To maintain safe transit operations, a forced inflow control measure is widely used in the Beijing rail transit, which limits the number of passengers entering the stations.Most of the controlled stations are located on these congested lines.
A significant amount of passenger flow information for different levels can be output by the simulator, such as platform, train, line, and network.Figure 12 shows the passenger density on the down direction platform of the Tiantongyuan    (TTY) station every 30 seconds.The TTY station is one of the most congested stations on Line 5. We can find that the maximum passenger density is as high as 3.0 passengers per square metres, which almost reaches the maximum safe density for platforms.The platform passenger density can be used to support train schedule regulation decisions and transport capacity allocations.
Figure 13 presents the time-varying section flows in the up direction of Line 1.The maximum section flow of the lines is the key parameter for calculating the train headway time.We can see that the pivotal sections of Line 1 are sections 6 to 11; nearly 120,000 passengers travel across these sections in the morning.
Figure 14 presents the network inflows, outflows, and passengers loading distribution every 30 minutes.There is no doubt that the reason for these two sharp peaks in the morning and evening is commuter travelling.The imbalanced travel demand makes it difficult for transport organizations and capacity utilization.We can find that the peak of the network load is approximately 500,000 passengers, which is less than the maximum capacity of the simulator.Hence, the  developed simulator is efficient enough for large-scale transit network applications.

Special Events Evaluation.
The traffic management for special events such as concerts, large conferences, and sports events is a difficult task for transport corporations.Generally, the transport schedules should be adjusted for the mass traffic flows caused by these events.However, precise quantitative empirical data are usually insufficient due to the lack of observation data and analysis tools.In rail transit systems, the detailed trip records of each traveller can be stored in the AFC system.However, it is also impossible to determine which passenger is part of the event audience.Considering the high regularity of rail transit passengers, a comparison method between the event day and the corresponding weekday of the previous week is used to evaluate the traffic influence of special events.A football game held in Worker's Stadium on 23 April 2016, starting at 19:35 and ending at 21:40, is used as our case study.The Worker's Stadium is one the most important centres for entertainment and sports in Beijing, and it can accommodate approximately 150,000 people.The location of the Stadium is shown in Figure 15.The Dongsishitiao (DSST) Station is the nearest railway station to the Worker's Stadium, so it is the most used station for passengers going to the stadium.In this study, the comparison day of the previous week is 16 April 2016.
Figure 16 shows the time-varying inflow and outflow distributions of the DSST station.It can be clearly observed that (1) there is a high arrival peak flow (station outflow) before the game, with the peak time at 18:30; (2) there is also a high departure peak flow (station inflow) after the game, with the peak at 22:00; and (3) the number of departure passengers entering the DSST station is as much as four times than the normal day, and the departure peak is much higher than the arrival flow.If there is no additional transport capacity added to support these mass flows, congestion and safety issues may emerge for the event stations.The number of extra section flows is used as the evaluation indicator for the event influence, which equals the section flows in an event day minus the flows in a normal day.After simulating the traffic status of two comparison days, the time-dependent extra section flow distribution after the event ends is shown in Figure 17. means the number of flows, which also corresponds to different colours.The red colour means the extra section flow is over 800 passengers, and the other colours represent smaller numbers of passengers.It can be found that (1) an outward evacuation process exists after the sports game ends, with a peak 20 minutes later; (2) the extra section flow (event flow) mainly travels through Lines 2, 6, and 1; hence, the train schedule adjustment should be made for these lines; and (3) fifty minutes later the event flows almost disappear from the network.
The maximum section flow in a line determines the train headway time.Table 5 gives the maximum section flows for the most influenced lines, as shown in Figure 17.We can determine that (1) the most influenced line is Line 2 and the influenced time is from 21:30-22:30 and (2) the most influenced time for the network is 22:00-22:30.Hence, some regulatory measures should be made for transporting these additional event flows, such as reducing headway time or running a short route train for these influenced sections.

Conclusions
A simulation-based DTA framework and models for the schedule-based rail transit network were proposed in this work.The dynamic equilibrium status of network flows is achieved by integrating the pretrip equilibrium and enroute path switch models, which is consistent with the route choice

Figure 1 :
Figure 1: Simulation-based DTA framework and model structure for a schedule-based transit system.

Figure 3 :
Figure 3: An example of a route switch procedure.

Figure 4 :
Figure 4: Passenger queues for passenger walking within a station.
(i) Modelling the complex routing behavior integrated with the pretrip route selection and enroute switch processes for rail transit travellers.(ii) Modelling the train movement for different types of train schedules, such as long and short train routes and a variable number of train units.(iii) Representing dynamics of network flows at multiple levels (network, line, station, section, train, and platform) and for any time span.(iv) Tracking the target passengers for behavior changes under special conditions.

Figure 10 :
Figure 10: Time-dependent deviations of section flows.

4. 3 . 1 .
Dynamics of Passenger Flows.The descriptive capability for dynamic flow status is the fundamental and core function of DTA models and simulators.The same traffic demand input and train schedules used in the model verification are applied in this case study.The simulation time is from 5:00 to 24:00.Part of time-dependent transport capacities for each line are given in Table

2 )Figure 12 :
Figure 12: Passenger flow density on the down direction platform of the TTY station.

Figure 13 :
Figure 13: Time-varying section flows in the up direction of Line 1.

Figure 16 :
Figure 16: Time-dependent inflows and outflows of the DSST Station.

Table 1 :
Figure 2: Topology structure for a rail transit network.An example of a train schedule.

Table 2 :
Data structure of AFC record information.
Calculate  , by checking the time  with the expected departure time   for each passenger in the entrance queue Step 2. Calculate transfer flow  trans , ,  trans , = min{ , ,  pass  ,  store −  ,−1 } Step 3. Transfer passengers from the entrance queue to the wait queue  , =  ,−1 +  trans

Table 3 :
Information of parameters.

Table 5 :
Comparison results between the event day and a normal day.