Managing Recurrent Congestion of Subway Network in Peak Hours with Station Inflow Control

Station inflow control (SIC) is an important and effective method for reducing recurrent congestion during peak hours in the Beijing, Shanghai, and Guangzhou subway systems. This work proposes a practical and efficient method for establishing a static SIC scheme in normal weekdays for large-scale subway networks. First, a traffic assignment model without capacity constraint is utilized to determine passenger flow distributions on the network. An internal relationship between station inflows and section flows is then constructed. Second, capacity bottlenecks are identified by considering the transport capacity of each section. Then, a feedback-based bottleneck elimination strategy is established to search target control stations and determine their control time and control strength. To validate the effectiveness of the proposed approach, a decision support system coded in the C# programming language was developed, and the Beijing subway was used as a case study. The results indicate that the proposed method and tool are capable of practical applications, and the generated SIC plan has better performance over the existing SIC plan. This study provides a practical and useful method for operation agencies to construct SIC schemes in the subway system.


Introduction
1.1.Motivation.Station inflow control (SIC) is to relieve congestion and ensure operational safety by controlling the number of passengers entering a station [1].It is usually applied during excessive congestion such as peak hours and other expected high passenger flows.Though SIC causes travel delays and inconvenience to passengers, it plays an important role in maintaining safety when travel demand greatly exceeds capacity.In April 2011, the Beijing subway first employed SIC at 17 stations during morning peak hours [2].Since then, SIC has been gradually implemented in other subway systems in China, such as the Guangzhou and Shanghai subways.With the continuous growth of travel demand, the number of controlled stations and the length of control time have increased greatly.At the end of 2016, the number of controlled stations in the Beijing subway reached 74, covering almost 25% of all stations [3].
Up to now, a practical and efficient method to generate SIC schemes for large-scale transit networks is still lacking.Current SIC scheme is established almost depending on operator's experience.According to the regulations on operation safety management for urban rail transit, which was released by the Beijing municipal commission of transport, a station should be controlled when passenger inflows reach the warning state of 70% of the load capacity of the station [1].This regulation is the only official guidance for carrying out SIC actions.However, it does not provide a clear and system-wide optimization approach.For instance, congestion of a certain station may not be caused by large inflows but by limited train capacity.Hence, a system-wide optimization approach that considers not only large inflows but also train capacities should be constructed.
To address the SIC issue in subway systems, previous works have proposed mathematical optimization models to establish a cooperative SIC plan within multiple stations [4,5].Nevertheless, these methods have low computational efficiency and lack capability to address the problem in largescale networks, because of complex solution and oversimplification of the practical problem in formulating model.It is well known that capacity bottlenecks are the critical reason for congestion.Hence, we start from this basic point and construct a new feedback-based bottleneck elimination strategy to determine the SIC plan.Compared to the existed mathematical models, the proposed method has direct meaning and high computational efficiency and can be easily expanded by considering practical factors in different situations.The major characteristics of the model include the following.(a) A feedback-based bottleneck elimination algorithm is used to generate the control scheme, which can provide quantitative support for operation managers, including control stations, control times, and control strength.(b) The method has high computational efficiency and capability for applications in large-scale rail transit networks.(c) The proposed approach has good extensibility with regard to practical factors, such as platform load capacity and traffic environment outside the station.Subway.The Beijing subway developed rapidly in recent years, consisted of 18 lines with a total length of 574 km by the end of 2016.The ridership reached 3.66 billion trips in 2016 (based on the count of unlinked trips, as every transfer in the system is counted as an additional trip), which means more than 10 million trips per day [6].To relieve congestion, train headways in some lines have been reduced to the minimum time (less than two minutes).However, travel demand is much greater than the transport capacity, and the congestion is becoming serious.

SIC Applications in Beijing
The Beijing Subway firstly employed SIC measures in April of 2011, which controlled seventeen stations in Line 1 and Line BT [2].The control time starts from 7:00 am to 9:00 am on weekdays.Then, the number of controlled stations continues to grow as travel demand increases.Figure 1 presents the number of controlled stations and annual ridership in the last five years.
In practice, three control levels are generally employed based on the crowding extent in a station.The first level limits passengers entering platform when platform capacity is insufficient.The second level limits passengers going to purchasing area, such as the area close to ticket gates.If the congestion cannot be relieved by these two control actions, then a third level is adopted, which limits passengers entering station by setting fences to slow down flow speed or closing station gates.Control level increases with the extent of congestion growth.Representative SIC measures applied in the Beijing subway are shown in Figure 2.
According to a survey of the Beijing subway, SIC schemes will be updated every three months in accordance with the variation of travel demand.If a station meets with heavy congestion, station masters will apply for a SIC action to the Beijing Subway Operation Company.Then, operation agencies decide whether to authorize the application depending on flow status and previous experience.However, there exist drawbacks in the redesign process for establishing a new SIC plan.First, station masters only make decisions based on a single station's crowdedness, not from the system wide.Second, there is no scientific method to support operation managers to decide which station should be controlled.Moreover, a SIC scheme only determines which station and what time should be controlled, not including the control strength.

Objective and
Organization.This work aims to provide a practical method for generating static SIC schemes for largescale subway networks only in peak hours and guides SIC actions through a quantization reference of control strength.
A proper SIC scheme could ensure operational safety and improve the service level for passengers.The remainder of this paper is organized as follows.Section 2 provides related works on SIC actions and congestion remission measures for rail transit system.The methodology to determine a control plan based on a feedback-based bottleneck elimination approach is given in Section 3. Afterwards, a discussion of its implementation in the Beijing subway is provided in Section 4. Finally, conclusions and future research are discussed.

Literature Review
Congestion is an intractable problem for worldwide traffic and transit systems.Two important strategies are usually used to manage congestion, which are capacity enhancement and travel demand management.For subway systems, it is difficult to improve capacity because physical facilities have upper limits on capacity, and a long time is required to build new lines.Measures of demand management become effective and positive tools for relieving congestion.In practice, two types of congestion-relief measures are widely used.The first one is time-varying price approach which has been used in the London [7], Melbourne subway systems [8,9], and so on.The second is station inflow control which is broadly implemented in the subway systems of China.Next, we will summarize the related works on these two aspects.2.1.Time-Dependent Pricing.The differential price approach to shift demand from peak periods is not new to either research or practice, which has been implemented in many large cities in the US from the 1970s [10].There are two basic strategies of differential pricing measure, peak surcharges and off-peak discounts.In the field of subway system, Santiago underground system firstly employed time-varying pricing strategy to reduce peak congestion in 1986 [11].Thereafter, many subway systems have implemented the varying price strategy to manage peak demand.A recent review by LEK finds that around 40% of urban rail networks worldwide provide some form of peak surcharge and/or off-peak discount [12].The most famous differential pricing strategy is "early bird tickets" in Melbourne metro, which provides free tickets if trips are completed before 7:00 am.After a half year of the implementation of the policy, it is found that 23% of commuters transfer from peak hours to off-peak periods, and the departure time of travelers moves forward about 42 minutes in average.The pricing strategy reduces peak demand and relieves congestion significantly [9].More applications of differential price policy in subway systems can refer to the reviews of Hale and Charles [13] and Liu and Charles [14].The mechanism of differential pricing for congestion relief is to influence passenger's travel behaviors and even the demand in temporal-spatial domains.Pricing strategy can intervene with the choices of passengers on departure time, destination, travel mode, and travel route.Usually the effect of a new price policy needs a long time to appear.Compared with the elastic pricing strategy, inflow control manages peak demand in a direct and forced way, by removing the mountain peak of the demand.For a mature subway system, pricing strategy is more preferable than control actions, while in developing subway systems where travel demand grows and fluctuates greatly, inflow control is an efficient measure for reliving congestion and maintaining operational safety.

Station Inflow Control. Station inflow control has been
widely implemented in the subway systems of China, including Beijing, Tianjin, Shanghai, Guangzhou, and Chengdu subway.With the peak congestion being serious, how to perform SIC actions becomes an absorbing topic for researchers and practitioners.Previous works can be divided into three categories: macroscopic principles for SIC strategy, mesoscopic methods for constructing a SIC plan, and microscopic measures for carrying out SIC actions in a station.This work falls into the second category on how to construct SIC schemes.
Macroscopic principles provide overall guidance for SIC.Liu and Jiang pointed out earlier that SIC actions should be executed at three levels of "station-line-network" [15].In practice, congestion of a station may not be caused by large inflows but for limited train capacity, because the stations ahead of the congestion station occupy most of the train capacity.Hence, collaborative control within multiple stations should be considered.Based on a balance principle, Liu et al. proposed a cooperative control strategy between two stations [16].Huang   a single line [17].Cooperative inflow control within multiple stations has been a consensus agreement for carrying out SIC actions.
On the mesoscopic level of method construction for SIC plans, Zhao et al. firstly proposed a mathematical optimization model to establish SIC schemes for a single line, with the objective of minimizing travel delay and maximizing the turnover of passenger flows [4].Lu et al. also proposed a linear integer programming model on the basis of network topology and definition of passenger demand; however, the model is only for a line [18].Then, Yao et al. established a coordinated passenger inflow control model for networks.In order to improve the efficiency and capability of the model for a network, a static assignment model is used to establish the internal flow relationship [5].With the view of reducing travel delay, Guo et al. built a cooperative model for a network with constraints on the capacity of stations [19], and Wang et al. established an integer programming model based on an analysis of passenger delay and alight and board processes [20].In conclusion, these models try to use optimization method to solve the SIC problem.Though cooperative control strategy can be considered, large numbers of parameters and complex solution limit their capability in real large-scale networks.The critical point of these models is establishing internal flow relationships, whose complexity increases greatly with network scale grows.In practice, congestion usually takes place in a regional wide.When we construct a SIC plan, there is no need to consider the stations far away from the target congestion station.Different from optimization methods, we establish a heuristic inflow control algorithm to establish SIC schemes.
Moreover, some works focus on SIC actions implemented in stations, such as ticket gates optimization and physical fences setting for inflow control [21].Dou et al. provide a cloud model to determine what time to start control actions for a station [22].Xie studied the detail control actions for transfer stations [23].In addition to normal control actions, Coulson et al. propose a new strategy of monetary incentives to reduce congestion by redirecting of passengers to less crowded exits [24].However, it is difficult to use a normalized method to determine control actions, because each station has specific physical structure and organization rules.How to establish detail SIC actions is out of range of this work.

Methodology
This section discusses a feedback-based bottleneck elimination algorithm for generating SIC schemes.Firstly, the framework to generate SIC scheme is provided.Then, the internal relationship between station inflows and section flows is established using a capacity-unconstrained assignment model, which provides important parameters for constructing SIC plans.Thirdly, capacity bottlenecks are identified by considering section flows and corresponding transport capacities.Fourthly, the elimination processes for a single bottleneck and multiple bottlenecks are presented separately.Lastly, the flowchart for generating a time-dependent SIC plan for a network is given.

Framework for Generating an SIC Scheme.
It is well known that bottleneck caused by the imbalance between supply and demand is the origin of congestion.We proceed from this point to establish a feedback-based elimination algorithm to remove bottlenecks and then search the corresponding target control stations.Usually, a SIC scheme contains three parts, namely, control stations, time, and strength.The framework for SIC scheme generation is shown in Figure 3.

Flow Relationship Construction.
In a rail transit system, transport capacity has an upper limit, and passengers should wait for the next train if the capacity of the incoming train is insufficient.In order to obtain the original status of flow distributions on a network, the capacity limitation of section is not considered in our traffic assignment model but will be taken into consideration later.A Logit-based stochastic user equilibrium (SUE) assignment model is applied to establish the relationship between station inflows and section (link) flows [25].In the model, there are mainly two types of relationships: (1) the inflows of a station that travel through given sections and their corresponding travel through rates and (2) the link flows of a section that comes from particular stations and their corresponding capacity occupation rates.
A subway network is represented by a directed graph  = (, ), where  = {1, 2, . . ., , } is the set of stations (nodes) and  = {1, 2, . . ., , } is the set of sections (links).The analysis time (such as 7:00 am to 9:00 am) is divided into short time spans with the same length, which are denoted by  = {1, 2, . . ., }.Let   () be the travel demand from station  to station  in time interval , where   () and   () are the section flow and transport capacity of section  at time , respectively.Define   as the set of feasible routes for the OD (Origin-Destination) pair (, ), and   as the th route,   ∈   .Let    be the generalized travel cost of the th route.
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 can be calculated by the distance corresponding to the walk link and the average walking speed of travelers.Normally, the on-board time for passengers travelling through sections is constant, which can be fixed from train schedules (timetable).The waiting time includes two components: (a) waiting time at the origin station; (b) waiting time at the transfer station if a transfer is required.The average waiting time at an origin/transfer station is equal to half of the train headway time.The transfer time is similar to the access or egress walk time.Then the generalized cost of a feasible route can be formulated as follows: where are waiting time at origin and transfer station, respectively.The weight factors for different travel processes can be estimated by the maximum likelihood method through travel surveys.Referring to the previous work on traffic assignment for the Beijing subway [26],   = 0.21,   = 0.14,   = 0.28, and   = 0.37.The value of factors may be variant for different subway systems.Suppose the OD matrix at time  is given, which can be accurately determined from automatic fare collection (AFC) records.Take an OD pair (, ) as an example to show the process of establishing the internal flow relationship.According to the Logit-based SUE model, the probability for passengers in OD pair (, ) choosing the th route is given by ( 2) and the flow volume of the th route can be computed using (3).
where    is the choice probability and  is the average cost of all feasible routes in the OD pair (, ).
Assign the whole demand of the OD pair (, ) to all feasible routes; then passenger flows of section  that comes from this OD can be determined, as shown in the following: where    () refers to the passenger flow of section , which comes from OD pair (, ), and     is a binary variable: if     = 1, section  belongs to the route   ; otherwise     = 0.
Next, assign all demands of all OD pairs to the network; the total section flow of section  can be calculated by In this work, two important parameters are defined to describe the internal relationship between station inflows and section flows.These are the station-section travel through rate    and the section-station capacity occupation rate    .The two parameters describe the correlation from opposite perspectives, where    starts from station view to explain that the inflow of a station will go through specific sections, and    aims to explain that section flows come from particular stations.An example of these two parameters is shown in Figure 4.In Figure 4(a), assume the inflow volume of station A is 1000 passengers and the flows traveling through the subsequent sections are 1000, 800, 600, and 500.Therefore, the station-section travel rates    are 100%, 80%, 60%, and 50% respectively.The station-section travel rate can be computed by where    () represents the percentage of passengers who depart from station  traveling through section  at time .
Suppose the section flows between stations C and D are 500 passengers, as shown in Figure 4(b), in which the numbers of passengers from stations A, B, and C are 100, 150, and 250, respectively.Then, the section-station capacity occupation rates are 20%, 30%, and 50%.The formulation of the section-station capacity occupation rate is given by , ∀ ∈ ,  ∈ ,  ∈ , where    () represents the percentage of passenger flows from station  in all flows traveling through section  at time .
The section-station capacity occupation rate    should satisfy

Capacity Bottlenecks Identification.
A bottleneck emerges when section flows exceed its transport capacity.In practice, passengers should queen on the platform for unable boarding.Note that there is no capacity limitation in our traffic assignment model; then the bottleneck is exhibited as the surplus flow rather than section capacity.Define Δ  as the surplus flow of section m, as shown in the following: The transport capacity can be calculated from the train timetable directly or from train headways, as shown in (10).Because of the randomness of passengers' arrival and imbalance in the distribution of passengers on different vehicles, the available capacity of a train is usually less than the theoretical capacity.Here, the peak hour factor (PHF) is used to describe the reduction of transport capacity [27] and is determined using (11).
where   is the number of trains running through section  during time t, D means the nominal load capacity for a train (usually 1440 passengers per train),  PHF is the PHF (0.25 ≤  PHF ≤ 1.0), and  is the permissible overload rate ( ≤ 130%).
where  60 is the inflows of a certain line in an hour and  15 is the maximum peak flow during 15 minutes in this hour.

Single Bottleneck Elimination Algorithm.
The elimination algorithm for a single bottleneck is the basis for removing multiple bottlenecks.We will describe how to remove a single bottleneck first.There are two types of control strategies for eliminating a bottleneck: single-station control and cooperative multistations control, as shown in Figure 5. Consider the bottleneck in Figure 5 as an example to show the elimination process.Suppose section C-D is a capacity bottleneck and Δ C-D is the overload flow.While Δ C-D is low, the bottleneck can be removed by controlling station B (assume the flows of section C-D mainly come from station B).When the congestion is heavy, multiple stations, such as stations A, B, and C, should be controlled.Two important questions should be answered when we remove the bottleneck, that is, (1) how to select the target control station/stations and (2) how to determine the control strength of each station that ensures the bottleneck can be removed.

Target Control Stations Selection.
The purpose of target control station selection is to determine which station/stations should be controlled.There is no doubt that it is more effective to control the stations with larger flows traveling through the bottleneck section.We can choose the target control stations according to the section-station capacity occupation rate    .Define overload rate  as the congestion pressure of the bottleneck, as shown in The higher the congestion pressure is, the more the stations should be controlled.Moreover, the more stations controlled, (1) Flow relationship: the section-station capacity occupation rate represents the capacity utilization of sections.It is the most important parameter used to determine the control station weight.The larger the rate, the more important the station in removing the bottleneck.
(2) Response time: response time describes the temporal connection between control stations and bottlenecks, which is expressed by the running time of the train between the control station and the bottleneck.When the control station is near the bottleneck, the effect of the control action is more obvious.
(3) Traffic conditions outside the station: traffic condition is assistant factor when we carry out SIC actions.Two representative factors are considered in this work, which are the area of the square outside the station, and the bus operation condition.First, there should be sufficient space for passengers queuing outside the station.The smaller the station square, the lower the control strength.Second, if there are enough buses around the station, travelers have high possibility to transfer to a bus instead of metro.In fact, the frequency, service area, and level of buses influent passengers' transfer behavior.However, it is a hard work to consider these influences in depth, because each bus line has its specific service area and characteristics.For available data, a simple index of the number of bus lines outside the station is used to describe the bus operation condition in this work.Usually, the more the bus lines, the higher the control strength.
(4) Platform load capacity: platform load capacity is another major factor that determines the control strength, which can be measured based on the platform area and the designed flow density.It is safer and robust for a platform with a larger capacity.Hence, if the platform capacity of a station is insufficient, greater control strength should be used to maintain safety.
The station weight can be expressed as a function of the above factors, as shown in where    is the set of target control stations for bottleneck section ; Δ   is the train running time between station  and bottleneck ;    is the area of the station square;    is the area of platform;   is the number of bus lines around the station;  indicates the importance of these factors.In this work,  1 = 0.4,  2 = 0.1,  3 = 0.1,  4 = 0.1,  5 = 0.3.
For a given bottleneck, the weights of target control stations should satisfy the constraint of ( 14).The station weights can be normalized by (15).
where    () is the final station weight and γ  () is the initial control weight from (13).

Control Strength Calculation.
After the target control stations and the corresponding control weights are obtained, the control strength needs to be determined.The inflow control rate is defined to quantify the control strength, which is the percentage of passengers who are limited to entering the station relative to the travel demand, as shown in the following: where   () is the control rate,   () is the whole travel demand of station  at time , and    () is the number of passengers who are limited to entering the station.To remove bottleneck  with overload flow Δ  (), the effective inflows that need to be controlled can be calculated by Δ  ()   ().Note that this is the effective flow, not the real flow, because not all passengers who enter station  will travel through bottleneck .Here, we can use the station-section travel through rate to revise the control inflows, as shown in The control rate of the stations is given by

Multibottlenecks Elimination Algorithm.
There are always many bottlenecks on the network, and these bottlenecks are interconnected through passenger flows.For example, if we remove a bottleneck with heavy congestion, another light bottleneck close to this bottleneck may disappear.Consider the example in Figure 6, which shows the elimination process for multiple bottlenecks.Suppose both section C-D and section D-E are bottlenecks, and their target control station sets are {A, B, C} and {B, C, D}, respectively.When we control the inflows of stations A, B, and C to remove bottleneck C-D, the link flows through section D-E may reduce, or the bottleneck may even disappear.Hence, the connections among bottlenecks should be considered when eliminating multiple bottlenecks.The internal flow relationship among these bottlenecks is illustrated in Figure 6(b).The reduced flows of the related links can be computed by the stationsection travel through rate, as shown in where   () is the updated section flow of section  after station  is under control.
To rationally remove all bottlenecks, an iterative elimination algorithm is employed in this work: (1) Sort all bottlenecks according to their overload rate in descending order and select the first one as the target bottleneck.(2) Establish the inflow control plan for the target bottleneck by using the single bottleneck elimination algorithm; see Section 3.4.
(3) Update the section flows of the related sections using (19).
(4) Identify the bottlenecks on the network after updating the section flows.
(5) If all the bottlenecks are removed, the SIC plan is finished; otherwise, go to step (1).
3.6.SIC Scheme Generation for a Network.Based on the above approach for generating a SIC plan for a single time interval, the algorithm for construction of a SIC scheme throughout the peak hours (7:00 am to 9:00 am) will be established in this section.A SIC scheme usually contains three elements, namely, control stations, control time, and strength.Discretize the continuous time into short time periods of equal length (such as 30 minutes).For each time span, the proposed algorithm can be used to generate a portion of the SIC plan.
It should be noted that the passengers who are limited to entering the station in time interval  will influence the travel demand in next time period  + 1.Moreover, travel behaviors of passengers can be influenced by SIC actions, such as travel mode change and departure time reschedule.Then the temporal-distribution of demand can be changed.However, it is very hard to consider the interactive influence when we construct the SIC plan.In this work, it is supposed that travel demand is constant, and passengers will wait outside the station and not transfer to other travel modes.
Then, the travel demand at time  + 1 should be updated by considering the passengers delayed at time , as shown in (20).Moreover, the section flows related to the control station should also be updated, as shown in (21).
where   () is the updated inflow of station  in time  + 1, if the station is controlled in time .where   () is the updated section flow of section  in time  + 1.
The iterative algorithm based on the bottleneck elimination approach to generate the whole SIC scheme during peak hours is presented in Figure 7.

Tool Description.
A decision support tool based on the proposed method was developed in the C# programming language and is used to validate the efficiency and performance of the approach.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 interfaces of the tool.Thanks to the Beijing Subway Operation Co. Ltd. for providing data, the tool has been tested in the Beijing subway system.The results demonstrate that it has high efficiency for applications in large transit networks, and reasonable SIC plans can be obtained.
The main capabilities of the tool include the following: (a) constructing time-dependent internal-relationship between station inflows and section flows and representing them in graphical visualization; (b) generating SIC schemes based on the iterative bottleneck elimination algorithm and automatically exporting the plan in the form of a report; (c) evaluating the performance of the generated SIC plan; (d) tracking the detailed elimination process for a particular bottleneck.

Timeframe for Updating SIC Plan.
Though the proposed method has no limitation on the iteration step length for updating SIC plan, other factors from applicability and robustness should be considered in practice.Firstly, we will analyze the stability of the model inputs, which are timedependent OD tables.The relative deviation shown in (22) is used to measure the stability of OD flows.Figure 9 shows the deviations of OD flows between two related weekdays (the day and the same weekday of last week) in different time intervals.It can be clearly found that the deviations increase greatly with the time interval shortening.If the time step is five minutes, the average deviation is over 50%.There is no doubt that SIC plans are not credible if a short time span is chosen to update the plan.where RE() is the average relative deviation of OD flows in time interval ;  is the total number of OD pairs; and   () and   () mean the volume of flows in a certain day and the contrastive day of last week, respectively.Secondly, the feasibility of implementation of control actions should also be considered when we make the plan.Because rail transit stations lack flexible control equipment, it is difficult to change control actions frequently, such as fences.Consider the variation of travel demand and the feasibility of implementation in practice; we suggest that the timeframe (time length) for updating SIC plan should not be less than 30 minutes.Beijing Subway 4.3.1.Data.The Beijing subway, one of the largest and most congested transit systems in China, was analyzed as a case study.To avoid the fluctuating influence of travel demand, the average travel demand over five weekdays from Monday April 10th to Friday April 14th in 2017 is used in the model.The demand was obtained from the automatic fare collection (AFC) system.Time-dependent inflows of each line are provided in Figure 10.The SIC action is usually carried out during the morning and evening peak hours.However, only morning peak hours (7:00 am to 9:00 am) are discussed in the study.The time span for updating the SIC scheme is set to be 30 minutes.Hence, there are four small time intervals during the morning peak hours.Transport capacity is an important parameter for identifying bottlenecks.In this work, the train timetable is employed to calculate the transport capacity for each line or section, which is presented in Table 2.

Case Study of
Moreover, traffic conditions around the station (such as the bus and square area) are obtained from a traffic survey, and platform load capacity is determined from the design map of stations.These parameters are crucial to determining the station weight when removing a particular bottleneck.Considering the copious amount of information in these parameters, they are not represented here in detail.
Note that Line 4, Line 14, and Line 16, which belong to Beijing Mass Transit Rail (MTR) Corporation, are excluded from the SIC scheme.The main reason is the different operational concepts between the Beijing subway and MTR Corporation.Regular SIC actions are never applied to these lines.

Results
. In April of 2017, the Beijing subway controlled 61 stations during morning peak hours.The generated SIC plan from the proposed method includes 63 stations.Comparing with the actual SIC plan, 52 stations are the same as the actual controlled stations, 9 stations should be cancelled, and 11 new stations should be controlled.Figure 11 shows the distribution of the actual and generated SIC plans for the network, in which the "canceled" means the station is controlled in practice but not included in the generated plan, and "newly added" has an opposite meaning.Table 3 presents the detailed plan of each line.Consider too many items in the SIC scheme, only an example of the detailed plan for Line 1 and BT is shown in Table 4.The control rate represents the control strength, which equals the percentage of limited passengers and all inflows.It can also be regarded as the reduction of speed of passengers entering the station.In Table 4, the station names are represented by their acronyms.
Table 5 provides the detailed elimination process for key bottlenecks (the heaviest congestion) in Line 1 and BT between 7:30 am and 8:00 am.For a certain bottleneck, the target control stations, station weights, and control volumes can be tracked.Through tracking the bottleneck elimination process, we can determine the function of each control station for a certain bottleneck.

Performance Evaluation.
Measuring the performance of SIC schemes is important for carrying out a new plan in practice.However, it is difficult because the actual plan of the Beijing subway is not clear.We can only know the control stations and a rough range of control times (such as 7:00 am to 9:00 am), and we do not know pivotal information about the control strength of each controlled station.Therefore, the number of delayed passengers cannot be determined.In this work, the section flows are used to evaluate the performance in an indirect way.First, time-varying link flows of each section under the actual SIC plan can be obtained from the Revenue Clearance Center (RCC) of the Beijing subway, which represents the performance of the practical plan.Then, section flows under the generated SIC scheme can be gained from the flow assignment model.The relative deviation of section flows between the practical and generated plan is used to measure the SIC performance, as shown in where () is the average relative deviation of all section flows,    () is the section flow of section  at time  using the generated SIC scheme,    () is the actual flow from the RCC, and   is the set of sections with  elements.
Table 6 shows the average relative deviation of section flows in each line.Note that the travel demand used to generate SIC plan is the average OD tables from five consecutive weekdays in April of 2017.Therefore, the actual section flows are the corresponding mean value.The following can be seen: (1) In examining the results over all the peak hours, the change of relative deviations is very small.It is easy to understand that SIC action only influences the departure time of travelers rather than the travel demand.Hence, the total flows through the section will not change greatly.
(2) In view of the small time span, the section flows in the controlled lines generally increase during 7:00 am∼ 8:00 am and decrease during 8:00 am∼9:00 am to a certain degree.This indicates that passengers enter the station earlier and that the new plan has a positive impact on reducing passengers' travel delays.
(3) The larger the section flow values change, the more controlled stations there are in these lines, such as Line 1 and Line 5.These lines are the most congested lines on the network.Put another way, it is necessary to adjust the actual SIC plan in these lines to improve the flow management.

4.3.4.
Guidance for SIC Actions.After the SIC scheme for a network is obtained, there is another important work that should be considered, which is how to determine the detail control actions for a specific station.In practice, station masters decide how to control flows according to their experience and the specific operational environment of the station.This is a difficult work because each station may have different physical structure, organization rules for passenger flows, flow characteristics of different entrances and exits, and so on.This work cannot provide a normalized method to determine the detail station control actions.However, the control rate/strength of the scheme could provide a guidance for setting control actions.The control rate represents the percentage of passengers who are limited to entering the station relative to the travel demand, which can also be regarded as the flow speed reduction rate.Then, what we need to do is to use some proper measures to slow down the inflow speed.We think the microsimulation tool, such as LEGION and VISSIM, which is widely used for evacuation evaluation for stations, is a useful method to determine the detail actions.

Conclusion
In this study, an algorithm based on a bottleneck elimination strategy for generating SIC schemes was proposed, which has high computational efficiency for application to large-scale subway networks.The approach has a clear meaning and can be expanded easily by considering other factors.A tool based on the proposed method has been developed in the C# programming language.A case study of the Beijing subway validates that there are good performance and practical value of the model.The study can support managers in constructing scientifically SIC plans and help station masters adopt proper actions in reference to the control rate.
A static SIC plan is suitable for recurrent congestion during peak hours because commute travel is very stable in the subway system.However, there are also sudden changes

Figure 1 :
Figure 1: The number of controlled stations and annual ridership of the Beijing subway.

Figure 2 :
Figure 2: Representative SIC measures in the Beijing subway.

Figure 3 :
Figure 3: Framework for generating a SIC scheme.

Figure 4 :
Figure 4: Schematic of relationship between station inflows and section flows.

Figure 5 :
Figure 5: Schematic of the single bottleneck elimination process.

Figure 6 :
Figure 6: Schematic for the internal relationship between multiple bottlenecks.

1 Figure 7 :
Figure 7: Algorithm for constructing SIC plan in whole peak hours.

Figure 8 :Figure 9 :
Figure 8: Main interfaces of the decision support tool.

4 Figure 10 :Figure 11 :
Figure 10: Time-dependent inflows of each line during morning peak hours.

Table 2 : 3 :
Transport capacity for each line during peak hours.The number of controlled stations in each line.
et al. studied propagation pattern of congestion and proposed a collaborative control strategy in is the generalized cost of the route   .,   ,   ,   are weight factors for the cost corresponding to different travel processes;  = 1 if there is a transfer; otherwise,  = 0; wait

Table 1 :
The proposed values for determining the number of target control stations.

Table 1 .
Once the number of control stations is determined, the top  stations according to the sectionstation capacity occupation rate can be selected as the target control stations.
3.4.2.Station Weight Calculation.If we only consider the internal relationship between station inflows and section flows, the control strength can be determined directly.However, in practice, other factors should be taken into consideration.Define a station weight    to describe the contribution of the control station  to eliminate the bottleneck .Four typical factors are considered in this work.