# The Aggregation Mechanism Mining of Passengers’ Flow with Period Distribution Based on Suburban Rail Lines

Authors: Xiaobing Ding; Hua Hu; Zhigang Liu; Haiyan Zhu
Journal: Discrete Dynamics in Nature and Society
 (2017)

---

## Abstract

Stranding too many passengers at the stations will reduce the service level; if measures are not taken, it may lead to serious security problems. Deeply mining the time distribution mechanism of passenger flow will guide the operation enterprises to make the operation plans, emergency evacuation plans, and so on. Firstly, the big data theory is introduced to construct the mining model of temporal aggregation mechanism with supplement and correction function, then, the clustering algorithmTime_clusterkm,n is used to mine the peak time interval of passenger flow, and the passenger flow time aggregation rule is studied from the angle of traffic dispatching command. Secondly, according to the rule of mining traffic aggregation, passenger flow calculation can be determined by the time of train lines in the suburbs of vehicle speed ratio, to match the time period of the uneven distribution of passenger flow. Finally, an example is used to prove the superiority of model in determining train ratios with the experience method. Saving energy consumption improves the service level of rail transit. The research can play a positive role in the operation of energy consumption and can improve the service level of urban rail transit.


---

## Body

## 1. Introduction

The success of urbanization in China has also brought some problems as follows: urban traffic congestion in megacities is becoming more and more serious, and the focus of urban life is shifting to suburbs, and the passenger flow between suburbs increases. The passenger flow of the rail transit has a great change that is difficult to predict, resulting in increased difficulty of passenger organization; large passenger flow stranded will lead to major safety issues; for example, the serious trampling accident in Chen Yi Square in Shanghai Bund on New Year’s Day of 2015 was a security alarm. By means of big data theory analysis and advanced data mining technology, the study on the aggregation rule of rail passenger flow at the station can obtain the characteristics of passenger flow at different stations and different means in time and help the operators adjust the driving schedule in time and provide decision-making basis for daily transportation dispatching.In many passenger travel distribution models, the gravity model and the growth coefficient method are widely used. In 1940, Stouffer first proposed the interventional opportunity model. In 1955, Casey proposed a gravity model, which was subsequently interpreted NT in terms of maximum entropy and maximum likelihood. In 1965, Fumess proposed the well-known growth coefficient method. Subsequently, the linear regression model was brought into use in traffic planning to predict travel occurrence and attractiveness in the United States. At the end of 1960, the British scholars put forward a cross-classification method; on this basis, Gordon proposed a nonfamily travel model [1]. Yam et al. developed a special model to predict official travel [2]. Bowman and Ben-Akiva established a model to predict the travel occurrence from the traveler’s activities as the starting point [3]. These foreign scholars made a more in-depth study in the passenger flow prediction and distribution models, which laid a theoretical basis. In recent years, many studies have been conducted in China: Xiaoguang proposed a polynomial distribution lagging model for short-term traffic flow forecasting, which gave overall considerations to factors impacting on the traffic state time series except its own lagging term, and the model calculation is more accurate [4]. Enjian analyzed the theoretical flaws in the traditional disaggregated models and constructed an OD distribution forecasting model based on aggregated data [5]. On the joint model of OD distribution and path selection, Professor Peng used the extrapolation method on parameter calibration because of the lack of survey results on destination selection, which resulted in the separation of two phases and affected the prediction accuracy of the model [6]. A short-term passenger flow OD estimation model based on the state space method was proposed, and the accuracy of the prediction results is improved [7]. In the OD matrix, the objective function was estimated by using the generalized least square model as the OD matrix, and the corresponding solving process using the Lagrangian algorithm was given. However, the accuracy of back stepping results was barely satisfactory [8, 9]. A multiobjective mathematical model was proposed to address a situation in which several projects are candidate to be invested completely or partially. The computational results were compared with solutions obtained by NSGA-II algorithm [10]. A multiobjective model for relief resources distribution facilities under an uncertain condition is investigated in two ways of demand satisfaction by considering the relief resources accessibility and demand satisfaction in a fuzzy logic. To solve the problem and to do the sensitivity analysis, the NSGA-II algorithm was presented. A constraint method was also proposed to evaluate the performance of the proposed algorithm [11].The above research results have not been studied deep in the mass data preprocessing, and there have still been no supplementary mechanisms for coping with lack of data field value and so on. The application of passenger data mainly focuses on the calibration of model parameters. If we can carry on the fine pretreatment and the missing supplement to the passenger flow data, and on this basis, we explore deeper to probe the station time aggregation mechanism of passengers, which will provide the data support for the optimization of the scheme. Based on the passenger flow data of AFC inbound and outbound traffic from the operators, this paper studies the information tracking based on passengers’ transportation card numbers and travel time and calculates the travel time of passengers by combining relevant time parameter data, so as to construct the model of mining passenger aggregation mechanism. It is of great help on train operating plan making and temporary adjustment and can provide theoretical and methodic references for supplementation of passenger flow forecasting and emergency relief model.


## 2. Urban Rail Transit Passenger Flow Data Cleaning and Processing

After many years of operation, the rail transportation enterprise has accumulated a great deal of operation-related data, and the Oracle database has a large capacity, and because of not-so-reasonable design of fields, stored procedures, data rules, and paradigms in the early stage of database construction, the stored operational data contains redundancy or error that not only occupies a lot of storage space and causes a waste of storage resources but also greatly reduces the effectiveness of data analysis and stability. Therefore, for abnormal data identification, missing data complement, error data correction, redundant data cleaning, and so on [12], big data analysis and processing methods are of great significance.Data redundancy, missing data, and nonstandard quality problems will have a fatal impact on big data application, so there is a need for cleaning big data having quality problems. You can use the parallel technology to achieve high scalability of big data cleaning. However, due to the lack of effective design, common data cleaning methods based on MapReduce have redundant computation, which leads to the performance degradation. To optimize the data cleaning to improve the efficiency of cleaning operations may start from the process shown in Figure 1 to design the optimization process.Figure 1
The flowchart of data cleaning and preprocessing.In the data generated by the rail transit operation, some are filled in by hand, some are automatically generated by the machine, and some are inevitably missed. For this reason, a missing value filling mechanism based on naive Bayesian classification is proposed. See Figure2 for the basic process:(1)PCi∣X=PX∣Ci×PCiPX.Figure 2
The flowchart of missing value filling.P ( C i ∣ X ) means that when incident X occurs, the probability of incident Ci occurs, P(X∣Ci) means that when incident Ci occurs, the probability of incident X occurs, PCi means the probability of incident Ci, and P(X) means the probability of incident X.(1  ) Parameter Estimation. The entire missing value is populated with the value of the maximum probability as the filling value calculated by Bayesian classification [13]. The task of the parameter estimation module is to calculate all the probabilities using (1), where all values of P(X) are constant, so only the calculation of P(X∣Ci)×P(Ci) is needed. If the prior probability of missing values is unknown, then it can be assumed to be equal probability, and P(X∣Ci) can be calculated using (2)PX∣Ci=PX1∣Ci×PX2∣Ci×⋯×PXn∣Ci.P ( X ∣ C i ) means that when incident Ci occurs, the probability of incident X occurs, P(X1∣Ci) means that when incident Ci occurs, the probability of incident X1occurs,…,P(Xn∣Ci) means that when incident Ci occurs, the probability of incident Xn occurs.So the entire parameter estimation module is used to calculate the probabilityP(X∣Ci) of each value of all attributes.From the perspective of probability theory, when the sample space is large enough, the probability can be approximately equivalent to that of the frequency, and the system uses frequency of occurrence of each attribute value in tuples without missing values to calculate probabilityP(Xn∣Ci).(2  ) Connection Module. According to (2), the probability of various values to be filled in the tuples containing missing values that are determined by its dependency attribute value can be calculated. However, sinceMapReduce can only process one record at a time in theMap phase andReduce phase, the system must make the dependency attribute value and its conditional probability correlate; this is existence necessity and the problem to be solved of the connection module. The input data of the connection module is the output data of the parameter estimation module and the original data to be filled. The output data is the file correlating the dependent attribute values in the tuples containing the missing data with the conditional probability of the values.(3  ) Filling Module. The filling module is implemented with MapReduce. First, the output result and the original input data of the connection module are subject to connective operation with offset as the key. TheMap phase is similar to the connection module. TheReduce phase uses (2) to calculate the conditional probability corresponding to each Ci (a possible value of an attribute to be filled) and select Ci with the maximum P(Ci∣X) probability for filling the missing value.


## 3. The Time_clusterkm,n Model of Aggregation Mechanism Mining

When the database is changed, the fuzzy frequent attribute set may not be frequent in the updated database. In this paper, the definition of negative bounds of the fuzzy frequent attribute sets is proposed for the fuzzy sets of the fuzzy frequent attribute set, the negative boundary, and the fuzzy attribute set in the original database and used for mining the time passenger flow.


### 3.1. Mining Association Rules of Passenger Flow Time Series in Rail Transit

In fuzzy discretization of time series, set the time seriesS=(x1,x2,x3,…,xn), the time window of the width w acts on the subsequence Si=(x1,x2,…,xi+w-1) formed by S, and a series of subsequences x1,x2,x3,…,xn-w+1 of the width w is formed by single-step slip(3)Ws,w=Si∣i=1,2,…,n-w+1.W ( s , w ) is regarded as n-w+1 points in w Euclidean space, randomly assigned to the class k, to calculate each type of center as a representative of each category and the elements Sii=1,2,…,n-w+1 of the set W(s,w), and the membership attribute function uj(Si) belonging to the class-j representatives is(4)ujSi=1/si-xj21/b-1∑c=1k1/si-xc2j=1,2,…,k;b>1,where b>1 denotes a constant that controls the fuzzy degree of the clustering results. si-xj2 represents the square of the distance from each point to the representative point of the class j.In the calculation of the support, ifA occurs, B occurs in time T; that is, we need to determine the frequency of occurrence of A:(5)FA=∑i=1n-w+1uASi,where uA(Si) is the membership of point Si belonging to the Ath representative attribute.The supportc of the correlation rule A→TB is c(A→TB)=F(A,B,T)/F(A), where F(A,B,T) denotes frequency of occurrence of B in time T after A happens, with reference to the C-means method of the rule A→TB(6)CBT:A=pA×pBT∣Alog⁡pBT∣ApBT+1-pBT∣Alog⁡1-pBT∣A1-pBT,where p(A) denotes probability of occurrence of the rule in A, p(BT∣A) represents occurrence forms of B in time T after A occurs, and the right side of (5) represents the information acquisition process from a priori probability p(BT) to posteriori probability p(BT∣A).


### 3.2. Mining of Temporal Passenger Flow OD Pairs Based on the Support of Correlation Rules

By preprocessing the AFC passenger flow data of the rail transit stations, the obtained data have been more accurate and in accord with the data format and redundancy requirement of the application data mining tools. The main algorithm is used to explore temporal aggregation correlation rules for passenger data for the purpose of mining the temporal passenger flow and the passenger flow OD pairs with traffic card numbers as the major cord and the designed algorithm is shown as follows.Step 1.
CalculateuA(Si) first and then c(A→TB) according to the support of passenger flow of the rail transit, and set minsupport=c(A→TB).Step 2.
Initialize AFC passenger flow database according to the mining condition with the cleaned and preprocessed data as initial data for correlation rule mining, scanning the transaction databaseTID, and finding all the item sets whose length is k=1, forming candidate 1-item set (C1), substituting C1 into (5) and (6), and calculating the support of each item, comparing with the minimum support one by one, and forming the frequent 1-item set if the support is greater than the minsupport.Step 3.
Generate candidate(k+1) item set Ck+1 according to the frequent k-item set, if Ck+1≠Φ, go to the next step; otherwise the loop ends.Step 4.
Scan the database to calculate and determine the support for each candidate item set.Step 5.
Delete candidate items with the support smaller thanminsupport, forming a frequent item set Lk+1 of the (k+1) item.Step 6.
k = k + 1, go to Step 3.Step 7.
This results in the frequent item set correlation rules of temporal passenger flow aggregation.Step 8.
Output 2-hour particle size of the temporal passenger flow aggregation and the results of OD pairs to guide the preparation and adjustment of rail transit operation plans.


## 3.1. Mining Association Rules of Passenger Flow Time Series in Rail Transit

In fuzzy discretization of time series, set the time seriesS=(x1,x2,x3,…,xn), the time window of the width w acts on the subsequence Si=(x1,x2,…,xi+w-1) formed by S, and a series of subsequences x1,x2,x3,…,xn-w+1 of the width w is formed by single-step slip(3)Ws,w=Si∣i=1,2,…,n-w+1.W ( s , w ) is regarded as n-w+1 points in w Euclidean space, randomly assigned to the class k, to calculate each type of center as a representative of each category and the elements Sii=1,2,…,n-w+1 of the set W(s,w), and the membership attribute function uj(Si) belonging to the class-j representatives is(4)ujSi=1/si-xj21/b-1∑c=1k1/si-xc2j=1,2,…,k;b>1,where b>1 denotes a constant that controls the fuzzy degree of the clustering results. si-xj2 represents the square of the distance from each point to the representative point of the class j.In the calculation of the support, ifA occurs, B occurs in time T; that is, we need to determine the frequency of occurrence of A:(5)FA=∑i=1n-w+1uASi,where uA(Si) is the membership of point Si belonging to the Ath representative attribute.The supportc of the correlation rule A→TB is c(A→TB)=F(A,B,T)/F(A), where F(A,B,T) denotes frequency of occurrence of B in time T after A happens, with reference to the C-means method of the rule A→TB(6)CBT:A=pA×pBT∣Alog⁡pBT∣ApBT+1-pBT∣Alog⁡1-pBT∣A1-pBT,where p(A) denotes probability of occurrence of the rule in A, p(BT∣A) represents occurrence forms of B in time T after A occurs, and the right side of (5) represents the information acquisition process from a priori probability p(BT) to posteriori probability p(BT∣A).


## 3.2. Mining of Temporal Passenger Flow OD Pairs Based on the Support of Correlation Rules

By preprocessing the AFC passenger flow data of the rail transit stations, the obtained data have been more accurate and in accord with the data format and redundancy requirement of the application data mining tools. The main algorithm is used to explore temporal aggregation correlation rules for passenger data for the purpose of mining the temporal passenger flow and the passenger flow OD pairs with traffic card numbers as the major cord and the designed algorithm is shown as follows.Step 1.
CalculateuA(Si) first and then c(A→TB) according to the support of passenger flow of the rail transit, and set minsupport=c(A→TB).Step 2.
Initialize AFC passenger flow database according to the mining condition with the cleaned and preprocessed data as initial data for correlation rule mining, scanning the transaction databaseTID, and finding all the item sets whose length is k=1, forming candidate 1-item set (C1), substituting C1 into (5) and (6), and calculating the support of each item, comparing with the minimum support one by one, and forming the frequent 1-item set if the support is greater than the minsupport.Step 3.
Generate candidate(k+1) item set Ck+1 according to the frequent k-item set, if Ck+1≠Φ, go to the next step; otherwise the loop ends.Step 4.
Scan the database to calculate and determine the support for each candidate item set.Step 5.
Delete candidate items with the support smaller thanminsupport, forming a frequent item set Lk+1 of the (k+1) item.Step 6.
k = k + 1, go to Step 3.Step 7.
This results in the frequent item set correlation rules of temporal passenger flow aggregation.Step 8.
Output 2-hour particle size of the temporal passenger flow aggregation and the results of OD pairs to guide the preparation and adjustment of rail transit operation plans.


## 4. Case Study and Analysis

### 4.1. Cleaning and Missing Filling of Passenger Flow Data

Passenger flow data of this paper comes from the AFC gate of Shanghai Rail Transit Line 9, and passenger flow period covers 6:00–22:00. The main data field includes card number, arrival time, inbound station, departure time, outbound station, P + R (the indication if the P + R parking lots is used), and total time (travel time), many of which are redundant data such as fare evasion, transportation card lost, and other reasons that resulted in the invalid data, and some raw data are shown as Table1.Table 1
The AFC passenger flow data of Shanghai Metro Line 16 (original).
Card number Arrival time Inbound station Departure time Outbound station P + R Total time (minute) U12452172127 09:12 Jiuting 09:43 Xujiahui  31 U07452172874 05:58 University Town 06:49 Zhaojiabang Positive U31842254712 06:58 University Town 07:28 Yishan Rd. Positive U12452172127 08:31 Jiuting 09:03 Xujiahui U07452172874 06:21 University Town 07:23 Zhaojiabang Positive U31842254712 06:24 University Town 07:28 Yishan Rd. Positive 30 U07452172874 16:46 Zhaojiabang Rd. 17:51 University Town Positive U07452172874 17:02 Xujiahui    30 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮From the original data we can see that some data fields are missing, especially the fields that have to be filled; the data cleaning and filling method designed in Section2 is used to deal with the original data to facilitate data mining.In order to fill the loss values of the original data, substituting (2) into (1) yields(7)PCi∣X=PX∣Ci=PX1∣Ci×PX2∣Ci×⋯×PXn∣Ci×PCiPX.Based on the card number field value of AFC data, travel times and the OD percentage of each passenger are collected and retrieved by SQL statement and substituted into (7), and through parameter estimation, connection, filling, and other modules, supplementing the original data yields results of Table 2.Table 2
The AFC passenger flow data of Shanghai Metro Line 16 (supplementary and correction).
Card number Arrival time Inbound station Departure time Outbound station P + R Total time (minute) U12452172127 09:12 Jiuting 09:43 Xujiahui Negative 31 U07452172874 05:58 University Town 06:49 Zhaojiabang Rd. Positive 51 U31842254712 06:58 University Town 07:28 Yishan Rd. Negative 30 U12452172127 08:31 Jiuting 09:03 Xujiahui Negative 32 U07452172874 06:21 University Town 07:10 Zhaojiabang Rd. Positive 41 U31842254712 06:24 University Town 07:28 Yishan Rd. Positive 30 U07452172874 16:46 Zhaojiabang Rd. 17:32 University Town Positive 46 U07452172874 17:02 Xujiahui 17:48 University Town Positive 30 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮


### 4.2. Mining of Passenger Flow Period Aggregation

According to the time series discretization design in Section3.1, the passenger flow period is divided according to the 2-hour particle size, and single-step slip forms a series of time series S=(x1,x2,x3,…,xn) with width of 2.CalculatingW(s,w)={Si∣i=1,2,…,n-w+1} and the membership function ujSi=1/si-xj21/b-1/∑c=1k1/si-xc2j=1,2,…,k;b>1 obtains the membership of station passenger flow aggregation of every 2-hour particle size.Passenger flow period covers 6:00 to 20:00 and is divided as 2-hour particle size into 6:00–8:00, 8:00–10:00, 10:00–12:00, 12:00–14:00, 14:00–16:00, 16:00–18:00, 18:00–20:00, and 20:00–22:00, and membershipujSi, calculation of support C(A-B), and passenger flow of each period are shown in Table 3.Table 3
Characteristic values of 2-hour particle size of Da Xue Cheng Station.
Time series Class Membershipuj(Si) Calculation of supportC(A-B) Passenger flow 6:00–8:00 0.85 0.90 47,804 8:01–10:00 0.74 0.86 31,012 10:01-12:00 0.91 0.84 11,430 12:01–14:00 0.90 0.80 9,010 14:01–16:00 0.89 0.84 13,014 16:01–18:00 0.82 0.81 30,254 18:01–20:00 0.69 0.83 30,147 20:01–22:00 0.92 0.84 29,874The above-mentioned mining algorithm can be used to calculate 2-hour particle size period passenger flow of the Shanghai Metro Line 9: Guilin Road, Caohejing Development Zone, Hechuan Road, Xingzhong Road, Qibao, Zhongchun Road, Jiuting, Sijing, Sheshan, Dongjing, Da Xue Cheng Station. The results are shown in Table4.Table 4
The aggregation passenger flow based on distribution time.
Time range Station St1 St2 St3 St4 St5 St6 St7 St8 St9 St10 St11 6:01–8:00 7,031 6,031 9,214 15,274 22,215 14721 16,014 8,721 10,721 30,264 28,214 8:01–10:00 12,031 13,031 8,610 11,201 20,203 6,321 11,782 6,321 9,843 19,721 18,734 10:01–12:00 10,574 7,574 7,141 9,845 6,201 5014 9,487 7,034 7,123 10,351 16,321 12:01–14:00 9,414 3,414 6,573 9,216 5,845 4782 5,014 6,791 8,787 9,014 10,024 14:01–16:00 9,018 10,034 5,362 8,142 9,216 3587 6,571 8,482 5,591 7,782 9,789 16:01–18:00 8,124 8,124 5,217 9,345 13,201 8721 20,362 6,217 10,258 9,587 8,137 18:01–20:00 10,250 22,230 6,057 9,714 9,231 6321 18,217 9,354 9,239 8,214 10,354 20:01–22:00 9,217 13,217 5,325 8,521 7,227 5014 9,057 8,242 8,427 9,241 9,792Figure3 shows the change trend of distribution time data. It can show the passenger flow distribution of the stations in the 2-hour particle size range visually. The distribution of the passenger flow can guide the operation directly.Figure 3
The aggregation passenger flow based on distribution time.The mining results in Table4 can provide the data support for the travel rush hour of the passengers and traffic adjustment scheme of the rail transit station, so that the train capacity can be matched with the passenger flow, and the passengers are quickly distributed to the road network. It can also provide data support for the determination of the speed ratio of express/slow trains of the subway lines of the rail transit. The ratio of the express/slow trains can be determined according to the temporal passenger flow, thus saving energy consumption and passenger travel time.


### 4.3. Mining of Passenger Flow OD Pairs Based on Correlation Rules

According to Section3.2, c(A→TB) of the OD pairs of each traffic card number and compared with minsupport=c(A→TB), the pairs with values higher than it are determined as the OD pairs. The mining results are shown in Table 5. The mining results can provide accurate decision-making support for the determination of the arrival scheme of the express trains in the express/slow mode.Table 5
The mining of passengers’ travel OD pair rule (upward direction).
Card number O D P + R U12452172127 Jiuting Xujiahui Negative U07452172874 Songjiang University Town Zhaojiabang Rd. Positive U31842254712 Songjiang University Town Yishan Rd. Positive U01552172876 Songjiang New Town Qibao Positive U07452172874 Songjiang New Town Yishan Rd. Positive U12452376572 Songjiang Sports Center Century Ave. Negative U06571172546 Dongjing Century Ave. Negative U318745684324 Sheshan Xujiahui Negative U01555476158 Qibao Century Ave. Negative U15478475786 Jiuting Century Ave. Negative ⋮ ⋮ ⋮ ⋮It can be seen from Table5 that if we can combine the travel time of passengers to mine the OD pairs, it can provide more elaborated decision support for the optimization of the express/slow train program. To a certain extent, the accessibility of passenger travel is increased; especially in the peak commuting, accessibility is particularly important to reduce the unnecessary transfer, and transfer passenger flow in the road network is relatively reduced.


## 4.1. Cleaning and Missing Filling of Passenger Flow Data

Passenger flow data of this paper comes from the AFC gate of Shanghai Rail Transit Line 9, and passenger flow period covers 6:00–22:00. The main data field includes card number, arrival time, inbound station, departure time, outbound station, P + R (the indication if the P + R parking lots is used), and total time (travel time), many of which are redundant data such as fare evasion, transportation card lost, and other reasons that resulted in the invalid data, and some raw data are shown as Table1.Table 1
The AFC passenger flow data of Shanghai Metro Line 16 (original).
Card number Arrival time Inbound station Departure time Outbound station P + R Total time (minute) U12452172127 09:12 Jiuting 09:43 Xujiahui  31 U07452172874 05:58 University Town 06:49 Zhaojiabang Positive U31842254712 06:58 University Town 07:28 Yishan Rd. Positive U12452172127 08:31 Jiuting 09:03 Xujiahui U07452172874 06:21 University Town 07:23 Zhaojiabang Positive U31842254712 06:24 University Town 07:28 Yishan Rd. Positive 30 U07452172874 16:46 Zhaojiabang Rd. 17:51 University Town Positive U07452172874 17:02 Xujiahui    30 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮From the original data we can see that some data fields are missing, especially the fields that have to be filled; the data cleaning and filling method designed in Section2 is used to deal with the original data to facilitate data mining.In order to fill the loss values of the original data, substituting (2) into (1) yields(7)PCi∣X=PX∣Ci=PX1∣Ci×PX2∣Ci×⋯×PXn∣Ci×PCiPX.Based on the card number field value of AFC data, travel times and the OD percentage of each passenger are collected and retrieved by SQL statement and substituted into (7), and through parameter estimation, connection, filling, and other modules, supplementing the original data yields results of Table 2.Table 2
The AFC passenger flow data of Shanghai Metro Line 16 (supplementary and correction).
Card number Arrival time Inbound station Departure time Outbound station P + R Total time (minute) U12452172127 09:12 Jiuting 09:43 Xujiahui Negative 31 U07452172874 05:58 University Town 06:49 Zhaojiabang Rd. Positive 51 U31842254712 06:58 University Town 07:28 Yishan Rd. Negative 30 U12452172127 08:31 Jiuting 09:03 Xujiahui Negative 32 U07452172874 06:21 University Town 07:10 Zhaojiabang Rd. Positive 41 U31842254712 06:24 University Town 07:28 Yishan Rd. Positive 30 U07452172874 16:46 Zhaojiabang Rd. 17:32 University Town Positive 46 U07452172874 17:02 Xujiahui 17:48 University Town Positive 30 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮


## 4.2. Mining of Passenger Flow Period Aggregation

According to the time series discretization design in Section3.1, the passenger flow period is divided according to the 2-hour particle size, and single-step slip forms a series of time series S=(x1,x2,x3,…,xn) with width of 2.CalculatingW(s,w)={Si∣i=1,2,…,n-w+1} and the membership function ujSi=1/si-xj21/b-1/∑c=1k1/si-xc2j=1,2,…,k;b>1 obtains the membership of station passenger flow aggregation of every 2-hour particle size.Passenger flow period covers 6:00 to 20:00 and is divided as 2-hour particle size into 6:00–8:00, 8:00–10:00, 10:00–12:00, 12:00–14:00, 14:00–16:00, 16:00–18:00, 18:00–20:00, and 20:00–22:00, and membershipujSi, calculation of support C(A-B), and passenger flow of each period are shown in Table 3.Table 3
Characteristic values of 2-hour particle size of Da Xue Cheng Station.
Time series Class Membershipuj(Si) Calculation of supportC(A-B) Passenger flow 6:00–8:00 0.85 0.90 47,804 8:01–10:00 0.74 0.86 31,012 10:01-12:00 0.91 0.84 11,430 12:01–14:00 0.90 0.80 9,010 14:01–16:00 0.89 0.84 13,014 16:01–18:00 0.82 0.81 30,254 18:01–20:00 0.69 0.83 30,147 20:01–22:00 0.92 0.84 29,874The above-mentioned mining algorithm can be used to calculate 2-hour particle size period passenger flow of the Shanghai Metro Line 9: Guilin Road, Caohejing Development Zone, Hechuan Road, Xingzhong Road, Qibao, Zhongchun Road, Jiuting, Sijing, Sheshan, Dongjing, Da Xue Cheng Station. The results are shown in Table4.Table 4
The aggregation passenger flow based on distribution time.
Time range Station St1 St2 St3 St4 St5 St6 St7 St8 St9 St10 St11 6:01–8:00 7,031 6,031 9,214 15,274 22,215 14721 16,014 8,721 10,721 30,264 28,214 8:01–10:00 12,031 13,031 8,610 11,201 20,203 6,321 11,782 6,321 9,843 19,721 18,734 10:01–12:00 10,574 7,574 7,141 9,845 6,201 5014 9,487 7,034 7,123 10,351 16,321 12:01–14:00 9,414 3,414 6,573 9,216 5,845 4782 5,014 6,791 8,787 9,014 10,024 14:01–16:00 9,018 10,034 5,362 8,142 9,216 3587 6,571 8,482 5,591 7,782 9,789 16:01–18:00 8,124 8,124 5,217 9,345 13,201 8721 20,362 6,217 10,258 9,587 8,137 18:01–20:00 10,250 22,230 6,057 9,714 9,231 6321 18,217 9,354 9,239 8,214 10,354 20:01–22:00 9,217 13,217 5,325 8,521 7,227 5014 9,057 8,242 8,427 9,241 9,792Figure3 shows the change trend of distribution time data. It can show the passenger flow distribution of the stations in the 2-hour particle size range visually. The distribution of the passenger flow can guide the operation directly.Figure 3
The aggregation passenger flow based on distribution time.The mining results in Table4 can provide the data support for the travel rush hour of the passengers and traffic adjustment scheme of the rail transit station, so that the train capacity can be matched with the passenger flow, and the passengers are quickly distributed to the road network. It can also provide data support for the determination of the speed ratio of express/slow trains of the subway lines of the rail transit. The ratio of the express/slow trains can be determined according to the temporal passenger flow, thus saving energy consumption and passenger travel time.


## 4.3. Mining of Passenger Flow OD Pairs Based on Correlation Rules

According to Section3.2, c(A→TB) of the OD pairs of each traffic card number and compared with minsupport=c(A→TB), the pairs with values higher than it are determined as the OD pairs. The mining results are shown in Table 5. The mining results can provide accurate decision-making support for the determination of the arrival scheme of the express trains in the express/slow mode.Table 5
The mining of passengers’ travel OD pair rule (upward direction).
Card number O D P + R U12452172127 Jiuting Xujiahui Negative U07452172874 Songjiang University Town Zhaojiabang Rd. Positive U31842254712 Songjiang University Town Yishan Rd. Positive U01552172876 Songjiang New Town Qibao Positive U07452172874 Songjiang New Town Yishan Rd. Positive U12452376572 Songjiang Sports Center Century Ave. Negative U06571172546 Dongjing Century Ave. Negative U318745684324 Sheshan Xujiahui Negative U01555476158 Qibao Century Ave. Negative U15478475786 Jiuting Century Ave. Negative ⋮ ⋮ ⋮ ⋮It can be seen from Table5 that if we can combine the travel time of passengers to mine the OD pairs, it can provide more elaborated decision support for the optimization of the express/slow train program. To a certain extent, the accessibility of passenger travel is increased; especially in the peak commuting, accessibility is particularly important to reduce the unnecessary transfer, and transfer passenger flow in the road network is relatively reduced.


## 5. Conclusions

Suburban rail transport passenger transport organization problems have troubled operation and management enterprises over years. Passenger transport organizations are involved in key indicators such as the temporal passenger flow, passenger OD, and speed car ratio of the express/slow trains. Passenger flow data of the rail transit are complicated in structure, in which redundancy and missing coexist. In this paper, we use the big data processing method to purify and supplement the original data and then explore the time distribution rule and the OD pairs of the passenger flow in depth. The passenger flow AFC data of Shanghai Rail Transit Line 9 are taken as inputs and are revised, supplemented, and mined to give out the time distribution characteristics of the passenger flow of Line 9 in the end and verify the feasibility of the method. In the future, the time distribution rule and OD pairs of passenger flow under the network condition will be mined in order to guide the establishment of the operation plan of the whole road network.



---

*Source: 3819304-2017-06-22.xml*