CN106021918B - 一种轨道交通网络化运行推演算法 - Google Patents
一种轨道交通网络化运行推演算法 Download PDFInfo
- Publication number
- CN106021918B CN106021918B CN201610329595.4A CN201610329595A CN106021918B CN 106021918 B CN106021918 B CN 106021918B CN 201610329595 A CN201610329595 A CN 201610329595A CN 106021918 B CN106021918 B CN 106021918B
- Authority
- CN
- China
- Prior art keywords
- now
- plat
- cell
- platform
- moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 101100206458 Mus musculus Them4 gene Proteins 0.000 claims description 40
- 238000011144 upstream manufacturing Methods 0.000 claims description 26
- QPFMBZIOSGYJDE-UHFFFAOYSA-N 1,1,2,2-tetrachloroethane Chemical compound ClC(Cl)C(Cl)Cl QPFMBZIOSGYJDE-UHFFFAOYSA-N 0.000 claims description 12
- 229920013660 Cellon Polymers 0.000 claims description 12
- 238000000034 method Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
本发明的技术方案是提供了一种轨道交通网络化运行推演算法,其特征在于,包括以下步骤:第一步、构建轨道网络拓扑结构;第二步、以一个周期时长为固定间隔循环计算轨道网络拓扑结构中每个站台及每个区间的状态。通过本发明提供的算法可以得到如下数据:实时区间列车数,实时区间乘客数,实时区间列车乘客拥挤度,实时站台停站列车数,实时站台停站乘客数,实时站台停站列车乘客拥挤度,实时站台滞留候车人数,实时站台候车区人数,从而为轨道交通的管理提供便利。
Description
技术领域
本发明涉及一种基于自动售票检票系统(AutoFareCollection,AFC)及自动列车监控系统(AutomaticTrainSupervision,ATS)的数据推算出轨道网络中实时客流状态信息的方法。
背景技术
城市轨道交通由于高密度运转、列车行车时间间隔短、行车速度高、列车编组辆数多而具有较大的运输能力。单向高峰每小时的运输能力最大可达到6万~8万人次(市郊铁道);地铁达到3万~6万人次,甚至达到8万人次;轻轨1万~3万人次,有轨电车能达到1万人次。城市轨道交通的运输能力远远超过公共汽车,因此成为了城市交通网络中的主干。
为了保证城市轨道交通安全、有序地运行,目前已建立了自动售票检票系统(AutoFareCollection,AFC)及自动列车监控系统(AutomaticTrainSupervision,ATS)。自动售票检票系统是集计算机技术、信息收集和处理技术、机械制造于一体的自动化售票、检票系统。自动列车监控系统是非安全设备系统,包括控制中心系统和集中控制站设备。上述两个系统的建立为城市轨道交通的运营提供了大量数据,但是这两个系统无法推算出轨道网络中实时客流状态信息。
发明内容
本发明的目的是提供一种推算出轨道网络中实时客流状态信息的方法。
为了达到上述目的,本发明的技术方案是提供了一种轨道交通网络化运行推演算法,其特征在于,包括以下步骤:
第一步、构建轨道网络拓扑结构,在轨道网络拓扑结构中每个轨道交通站点为一个集群,若当前轨道交通站点包含有N条轨道交通线路,N≥1,则当前轨道交通站点包含N个逻辑站点,每个逻辑站点为当前集群下的一个子集群,每个逻辑站点包含一个位于轨道交通上行线路上的站台及一个位于轨道交通下行线路上的站台,则两个站台分别为当前子集群下的两个节点,相邻两个轨道交通站点中通过同一根轨道可直接到达的两个站台在轨道网络拓扑结构中通过一个区间连接起来;
第二步、以一个周期时长为固定间隔循环计算轨道网络拓扑结构中每个站台及每个区间的状态,对于任意一个站台Plat_1及任意一个区间Cell_1而言,其状态的计算步骤包括:
步骤2.1、计算区间Cell_1流出乘客,计算得到如下变量:
Cell_1.CellCPOUT(now)=Cell_1.CellCTOUT(now)×Cell_1.TrainCAP,式中,Cell_1.CellCPOUT(now)表示当前周期时刻至前一周期时刻之间的时间段内流出区间Cell_1的列车所携带的人数,Cell_1.CellCTOUT(now)表示当前周期时刻至前一周期时刻之间的时间段内流出区间Cell_1的列车的列数,Cell_1.TrainCAP表示区间Cell_1所对应的轨道交通线路线路每列列车的定员人数;
Cell_1.CellARV(now)=Cell_1.CellCPOUT(now)×Cell_1.CellS(pre),式中,Cell_1.CellARV(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达区间Cell_1所对应的下游站台的总人数,Cell_1.CellS(pre)表示上一周期时刻区间Cell_1内的列车平均乘客饱和度;
Cell_1.CellU(now)=Cell_1.CellARV(now)×Cell_1.CellURATE(now),式中,Cell_1.CellU(now)表示当前周期时刻Cell_1.CellARV(now)中到站后不下车的乘客数,Cell_1.CellURATE(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达下游站台的乘客中到站后不下车的比例;
Cell_1.CellET(now)=Cell_1.CellARV(now)×Cell_1.CellETRATE(now),式中,Cell_1.CellET(now)表示Cell_1.CellARV(now)中当前周期时刻到站后下车出站点结束本次出行的乘客数,Cell_1.CellETRATE(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达下游站台的乘客中到站后出站点结束本次出行的比例;
Cell_1.CellTRs(p,now)
=p.PlatQTRs(Cell_1,now)
=Cell_1.CellARV(now)×Cell_1.CellTRRATES(p,now),式中:
Cell_1.CellTRs(p,now)表示当前周期时刻离开区间Cell_1换乘往站台p的乘客数,p.PlatQTRs(Cell_1,now)表示当前周期时刻换乘到站台p的乘客中来自区间Cell_1的乘客数,Cell_1.CellTRRATES(p,now)表示当前周期时刻离开区间Cell_1换乘往站台p的比例系数;
Cell_1.CellTR(now)=∑Cell_1.CellTRs(*,now),式中,Cell_1.CellTR(now)表示当前周期时刻Cell_1.CellARV(now)中到站后下车换乘的乘客数,∑Cell_1.CellTRs(*,now)表示当前周期时刻离开区间Cell_1换乘往各个不同站台的乘客数的和;
步骤2.2、计算站台Plat_1流出乘客,若站台Plat_1仅为一个区间Cell_11的上游站台,则计算得到如下变量:
Plat_1.PlatCPOUT(now)=Plat_1.PlatCTOUT(now)×Plat_1.TrainCAP,式中,Plat_1.PlatCPOUT(now)表示当前周期时刻至前一周期时刻之间的时间段内流出站台Plat_1轨道的列车所运载的人数,Plat_1.PlatCTOUT(now)当前周期时刻至前一周期时刻之间的时间段内流出站台Plat_1轨道的列车的列数,Plat_1.TrainCAP表示站台Plat_1所对应的轨道交通线路线路每列列车的定员人数;
Cell_11.CellCPIN(now)=Cell_11.CellCTIN(now)×Cell_11.TrainCAP,式中,Cell_11.CellCPIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_11的列车所运载的人数,Cell_11.CellCTIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_11的列车的列数,Cell_11.TrainCAP表示区间Cell_11所对应的轨道交通线路线路每列列车的定员人数;
cpout=Min{Plat_1.PlatCPOUT(now),Cell_11.CellCPIN(now)},式中,cpout表示站台Plat_1的载客出站列车所运载的人数的临时变量;
Plat_1.PlatCN(now)=Plat_1.PlatS(pre)×cpout,式中,Plat_1.PlatCN(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat_1出发进入下游各区间的总人数中原该站台区域轨道列车内的乘客数,Plat_1.PlatS(pre)表示上一周期时刻站台Plat_1区域轨道上的列车平均乘客饱和度,到站列车需下车乘客已下车而站台候车乘客还未上车时的饱和度;
Ctmp=cpout×Plat_1.TrainCD-Plat_1.PlatCN(now),式中,Ctmp表示站台Plat_1的出站列车剩余的能运载的人数的临时变量;
Plat_1.PlatQTR(now)=∑Plat_1.PlatQTRs(*,now),式中,Plat_1.PlatQTR(now)表示当前周期时刻至前一周期时刻之间的时间段内由各个方向通过通道、同站台、出站方式换乘到某站台候车区的总乘客数,∑Plat_1.PlatQTRs(*,now)表示当前周期时刻至前一周期时刻之间的时间段内各个方向通过通道、同站台、出站方式换乘到站台Plat_1候车区的乘客数的和;
Plat_1.PlatQIN(now)=Plat_1.PlatQTR(now)+Plat_1.PlatQGT(now),式中,Plat_1.PlatQIN(now)表示当前周期时刻至前一周期时刻之间的时间段内到达站台Plat_1候车区候车的乘客总数,Plat_1.PlatQGT(now)表示当前周期时刻至前一周期时刻之间的时间段内由闸机进站后到达站台Plat_1候车区的非换乘乘客;
若Ctmp<=0,则表示站台候车区乘客全部无法上车,则:
Plat_1.PlatON(now)=0,式中,Plat_1.PlatON(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat_1出发进入下游各区间的总人数中由该站台候车区上车的乘客数;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now),式中,Plat_1.PlatDPT(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat_1出发进入下游各区间的总人数;
若Ctmp>0且Ctmp>=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now),则表示站台候车区乘客可以全部上车,式中,Plat_1.PlatQ(pre)表示上一周期时刻站台Plat_1区域候车区内由于列车过饱和而无法及时上车的乘客数:
Plat_1.PlatON(now)=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now);
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
Ctmp其他取值情况时,表示站台候车区乘客无法全部上车:
Plat_1.PlatON(now)=Ctmp;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若站台Plat_1为区间Cell_11及区间Cell_21的上游站台,则计算得到如下变量:
Plat_1.PlatCPOUT(now)=Plat_1.PlatCTOUT(now)×Plat_1.TrainCAP;
Cell_11.CellCPIN(now)=Cell_11.CellCTIN(now)×Cell_11.TrainCAP;
Cell_21.CellCPIN(now)=Cell_21.CellCTIN(now)×Cell_21.TrainCAP,式中,Cell_21.CellCPIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_21的列车所运载的人数,Cell_21.CellCTIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_21的列车的列数,Cell_21.TrainCAP表示区间Cell_21所对应的轨道交通线路线路每列列车的定员人数;
cpout=Min{Plat_1.PlatCPOUT(now),Cell_1.CellCPIN(now)+Cell_1.CellCPIN(now)};
Plat_1.PlatCN(now)=Plat_1.PlatS(pre)×cpout;
Ctmp=cpout×Plat_1.TrainCD-Plat_1.PlatCN(now);
Plat_1.PlatQTR(now)=∑Plat_1.PlatQTRs(*,now);
Plat_1.PlatQIN(now)=Plat_1.PlatQTR(now)+Plat_1.PlatQGT(now);
若Ctmp<=0,则表示站台候车区乘客全部无法上车,
Plat_1.PlatON(now)=0;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若Ctmp>0且Ctmp>=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now),则表示站台候车区乘客可以全部上车:
Plat_1.PlatON(now)=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now);
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
Ctmp其他取值情况时,表示站台候车区乘客无法全部上车:
Plat_1.PlatON(now)=Ctmp;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若站台Plat_1的下游未与任何区间相连,则:
Plat_1.PlatCN(now)=0;Plat_1.PlatON(now)=0;Plat_1.PlatDPT(now)=0;
步骤2.3、计算区间Cell_1流入乘客,若区间Cell_1为位于其上游的站台Plat的唯一一条区间,则计算得到如下变量:
Cell_1.CellDPT(now)=Plat.PlatDPT(now),式中,Cell_1.CellDPT(now)表示当前周期时刻至前一周期时刻之间的时间段内由上游站台Plat出发进入区间Cell_1的总乘客数,Plat.PlatDPT(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat出发进入下游各区间的总人数;
Cell_1.CellCN(now)=Plat.PlatCN(now),式中,Cell_1.CellCN(now)表示Cell_1.CellDPT(now)中原上游站台区域轨道列车内的乘客数,Plat.PlatCN(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat出发进入下游各区间的总人数中原该站台区域轨道列车内的乘客数
Cell_1.CellON(now)=Plat.PlatON(now),式中,Cell_1.CellON(now)表示Cell_1.CellDPT(now)中由上游站台候车区上车的乘客数,Plat.PlatON(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat出发进入下游各区间的总人数中由该站台候车区上车的乘客数;
若位于区间Cell_1上游的站台Plat有两条相连的区间,分别为区间Cell_1及区间Cell_2,则若Cell_1.CellCPIN(now)=0且Cell_2.CellCPIN(now)=0,式中,Cell_1.CellCPIN(now)及Cell_2.CellCPIN(now)=0分别表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_1及区间Cell_2的列车所运载的人数,则:
Cell_1.CellDPT(now)=0;
Cell_1.CellCN(now)=0;
Cell_1.CellON(now)=0;
否则:
Cell_1.CellDPT(now)=
Plat.PlatDPT(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.Cell CPIN(now)};
Cell_1.CellCN(now)=
Plat.Pla_tCN(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.Cell CPIN(now)};
Cell_1.CellON(now)=
Plat.PlatON(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.Cell CPIN(now)};
步骤2.4、计算站台Plat_1的流入乘客,若站台Plat_1的上游仅有一条区间Cell_1,则:
Plat_1.PlatU(now)=Cell_1.CellU(now),式中,Plat_1.PlatU(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车由上游各区间到达站台Plat_1的总乘客数中不下车的乘客数;
若站台Plat_1的上游有两条区间,即区间Cell_1及区间Cell_2,则:
Plat_1.PlatU(now)=Cell_1.CellU(now)+Cell_2.CellU(now),式中,Cell_2.CellU(now)表示当前周期时刻Cell_2.CellARV(now)中到站后不下车的乘客数,Cell_2.CellARV(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_2到达区间Cell_2所对应的下游站台的总人数;
若站台Plat_1没有对应的上游区间,或没有对应的下游区间,则:
Plat_1.PlatU(now)=0;
步骤2.5、更新站台Plat_1的状态:
当站台Plat_1没有对应的上游区间,即为始发类型站台时:
Plat_1.PlatM(now)=0,式中,Plat_1.PlatM(now)表示当前周期时刻站台Plat_1区域轨道内的列车数;
Plat_1.PlatN(now)=0,式中,Plat_1.PlatN(now)表示当前周期时刻站台Plat_1区域轨道列车上的在途乘客总数;
Plat_1.PlatS(now)=0,式中,Plat_1.PlatS(now)表示当前周期时刻站台Plat_1区域轨道上的列车平均乘客饱和度,到站列车需下车乘客已下车而站台候车乘客还未上车时的饱和度;
Plat_1.PlatQ(now)=Max{0,Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now)-Plat_1.PlatON(now)},式中,Plat_1.PlatQ(now)表示当前周期时刻站台Plat_1区域候车区内由于列车过饱和而无法及时上车的乘客数;
Plat_1.PlatPQ(now)=Plat_1.PlatQ(now)+Plat_1.PlatQIN(now),式中,Plat_1.PlatPQ(now)表示当前周期时刻站台Plat_1区域候车区内滞留乘客及正常进站候车乘客的总和,即站台Plat_1候车区乘客总数;
当站台Plat_1没有对应的下游区间,即为终点类型站台时:
Plat_1.PlatM(now)=0;
Plat_1.PlatN(now)=0;
Plat_1.PlatS(now)=0;
Plat_1.PlatQ(now)=0;
Plat_1.PlatPQ(now)=0;
当站台Plat_1为除终点类型及始发类型外的其他类型的站台时:
Plat_1.PlatM(now)=Plat_1.PlatM(pre)+Plat_1.PlatCTIN(now)-Plat_1.PlatCTOU T(now);
Plat_1.PlatN(now)=Plat_1.PlatN(pre)+Plat_1.PlatU(now)-Plat_1.PlatCN(now);
Plat_1.PlatQ(now)=Max{0,Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now)-Plat_1.PlatO N(now)};
Plat_1.PlatPQ(now)=Plat_1.PlatQ(now)+Plat_1.PlatQIN(now);
若Plat_1.PlatM(now)>0,则:
Plat_1.PlatS(now)=Plat_1.PlatN(now)/{Plat_1.TrainCAP×Plat_1.PlatM(now)};
若Plat_1.PlatM(now)<=0,则:
Plat_1.PlatS(now)=0;
步骤2.5、更新区间Cell_1的状态:
Cell_1.CellM(now)=Cell_1.CellM(pre)+Cell_1.CellCTIN(now)-Cell_1.CellCTO UT(now),式中,Cell_1.CellM(now)表示当前周期时刻区间Cell_1内的列车数,Cell_1.CellM(pre)表示上一个周期时刻区间Cell_1内的列车数;
Cell_1.CellN(now)=Cell_1.CellN(pre)+Cell_1.CellDPT(now)-Cell_1.CellARV(now),式中,Cell_1.CellN(now)表示当前周期时刻区间Cell_1内列车上的在途乘客总数,Cell_1.CellN(pre)表示上一个周期时刻区间Cell_1内列车上的在途乘客总数;
若Cell_1.CellM(now)>0,则:
Cell_1.CellS(now)=Cell_1.CellN(now)/{Cell_1.TrainCAP×Cell_1.CellM(now)},式中,Cell_1.CellS(now)表示当前周期时刻区间Cell_1内的列车平均乘客饱和度;
若Cell_1.CellM(now)<=0,则:
Cell_1.CellS(now)=0。
优选地,在所述步骤2.1中,若区间Cell_1没有任何与其直接相连的下游站台,则将所述Cell_1.CellU(now)设为0。
本发明根据实时ATS运能数据推算轨道网络各站台及区间内的虚拟列车运行过程及状态,再根据列车推演结果以及AFC客流短期预测结果和客流选择行为特征参数推算轨道网络中的客流时空分布。
通过本发明提供的算法可以得到如下数据:实时区间列车数,实时区间乘客数,实时区间列车乘客拥挤度,实时站台停站列车数,实时站台停站乘客数,实时站台停站列车乘客拥挤度,实时站台滞留候车人数,实时站台候车区人数,从而为轨道交通的管理提供便利。
附图说明
图1为计算区间Cell_1流出乘客时的计算模型图;
图2为计算一种情况下的站台Plat_1流出乘客的计算模型图;
图3为计算另一种情况下的站台Plat_1流出乘客的计算模型图;
图4为计算一种情况下的区间Cell_1流入乘客的计算模型图;
图5为计算另一种情况下的区间Cell_1流入乘客的计算模型图;
图6为计算一种情况下的站台Plat_1流入乘客的计算模型图;
图7为计算另一种情况下的站台Plat_1流入乘客的计算模型图。
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下。
在进行本发明提供的一种轨道交通网络化运行推演算法前需要针对轨道交通构建轨道网络拓扑结构,在轨道网络拓扑结构中每个轨道交通站点为一个集群,若当前轨道交通站点包含有N条轨道交通线路,N≥1,则当前轨道交通站点包含N个逻辑站点,每个逻辑站点为当前集群下的一个子集群,每个逻辑站点包含一个位于轨道交通上行线路上的站台及一个位于轨道交通下行线路上的站台,则两个站台分别为当前子集群下的两个节点,相邻两个轨道交通站点中通过同一根轨道可直接到达的两个站台在轨道网络拓扑结构中通过一个区间连接起来。
在不考虑算法输入数据异常的情况下,本发明提供的一种轨道交通网络化运行推演算法可以简单概括为以下步骤:
步骤1、获取静态配置及离线参数数据;
步骤2、根据离线参数数据,初始化所有区间和站台的初始状态变量;
步骤3、以一个周期时长为固定间隔(30秒)循环计算轨道网络的状态(包括6个子过程);
步骤3-A、计算区间流出乘客;
步骤3-B、计算站台流出乘客;
步骤3-C、计算区间流入乘客;
步骤3-D、计算站台流入乘客;
步骤3-E、更新站台状态;
步骤3-F、更新区间状态。
其中,步骤1与步骤2在程序启动时首次执行,并且在程序持续运行过程中每天凌晨执行一次;通常在轨道网络中无列车运营时段可以以全零状态初始化所有区间和站台的状态,这样可以得到与真实情况完全一致的初始状态,对后续推演不会产生任何影响;然而在一些情况下无法选定这种“全零”的时刻来启动系统,如在早高峰时段出现程序异常关闭重启后,此时需要使用离线推演的轨道网络区间及站台状态(以计划ATS数据、历史客流数据作为输入)来作为程序启动时刻所有区间及站台的初始状态值;为了综合这两种初始化情况,统一采用离线状态来作为初始化的状态值,而在离线状态中非运营时段轨道网络状态可以设定为“全零”(通过剔除计划ATS数据中非运营时段调配列车产生的列车进出站记录,再进行离线推演得到离线状态)。另外,由于在某个正常运营时刻,轨道网络的真实状态与离线推演得到的状态几乎不可能是完全一致的,这种不一致会导致后续推演的结果出现额外的误差,更严重的后果是导致推演单元(区间或站台对象)在某些时刻的推演结果出现异常(如储车量M<0,需要使用局部离线模式或使推演单元的状态为0)。
由于步骤1需要每天被执行一次,步骤3则是每个周期30秒被执行一次,所以在程序实现中将会生成两个定时器任务执行线程,步骤2会在每天凌晨步骤1执行过后第一次执行步骤3时被执行一次;而步骤3包含6个子过程,为了提高实时推演的计算效率,同时也考虑到推演算法逻辑中并行计算的可行性,6个子过程会按序先后生成6个线程组进行计算,仅当前一个线程组中的所有线程全部结束后,后一个线程组才会被创建并开始计算,其中每个线程会负责完成一条线路的计算任务,如果轨道网络中共有k条线路,则通过这一多线程机制可以至多提高k倍推演计算效率。6个子过程的计算具有前后依赖性,计算先后顺序固定;而在每个子过程内部,所有计算对象(区间或站台)之间具有并行更新无序性。
步骤3的基本思想是:前一周期时刻的轨道网络状态+当前周期时刻至前一周期时刻间时段内的轨道网络状态转移量=得到当前周期时刻的轨道网络状态。
假设已知前一周期时刻轨道网络所有状态变量A(pre)和周期时刻至前一周期时刻间时段内的BP类轨道网络状态转移量BP(now),则计算当前周期时刻轨道网络状态A(now)的具体方法(步骤3的6个子过程)如下。
在对步骤3的6个子过程进行详细说明之前,对使用的变量归纳说明如下:
算法变量共包括两大类:状态变量(A类变量)和状态转移变量(B类变量)。其中状态变量表示轨道网络在某个周期时刻的瞬间的各种状态值;状态转移变量表示轨道网络在前后相邻的两个周期时刻之间的时间段内发生的各种状态值的变化量,状态转移变量又包括输入变量(BP类变量)和中间变量(BM类变量)。
前一周期时刻的状态变量,叠加上该时刻到之后一个周期时刻之间的时间段内的状态转移变量,即可得到后一周期时刻瞬间的状态变量。因此轨道网络在某一周期时刻的状态由前一周期时刻轨道网络的状态以及两个周期时刻之间的状态转移量来决定。
状态变量(A类变量)
CellM:某个周期时刻某区间内的列车数;
CellN:某个周期时刻某区间内列车上的在途乘客总数;
CellS:某个周期时刻某区间内的列车平均乘客饱和度(拥挤度或载客率);
PlatM:某个周期时刻某站台区域轨道内的列车数;
PlatN:某个周期时刻某站台区域轨道列车上的在途乘客总数;
PlatS:某个周期时刻某站台区域轨道上的列车平均乘客饱和度(拥挤度或载客率),到站列车需下车乘客已下车而站台候车乘客还未上车时的饱和度;
PlatQ:某个周期时刻某站台区域候车区内由于列车过饱和而无法及时上车的乘客数;
PlatPQ:某个周期时刻某站台区域候车区内滞留乘客及正常进站候车乘客的总和,即站台候车区乘客总数。
关于PlatPQ的解释:
站台候车区人数PQ通常会小于实际观测到的站台人数,因为真实的列车运营是站台攒够一个列车间隔时间(如5分钟)的进站台人数(如100人),然后列车一进站一下子把人全运走(假设列车很空),所以观测到的站台人数是从0~100缓慢均匀增加的。而推演模型中由于周期是30秒,所以每个时刻进站台人数是100/10=10人,并且在运能充足时每个周期的人都能被运走,相当于每次(30秒)赞够一小波人(10人)然后马上运走,所以就会出现模型中站台候车人数连续10个周期(共5分钟)一直保持站台候车区人数PQ=10,而现实中可能是在这5分钟里观测到的站台人数是从0~100缓慢均匀增加的!!!根本原因是推演模型的时间粒度统一到了30秒,而现实中列车间隔一定大于30秒,并且现实的列车间隔是不确定的。(从站台候车人数时变图形曲线上看就是,现实中的锯齿形曲线被滤波成了水平线)。另外,周期设30秒是因为站台轨道也要推演(保证列车运行推演的连续性),而列车停站时间较小,若周期取很大则违反了CTM模型一个周期内车辆不得穿越单元规定(违反会导致单元进出流量不守恒,模型计算出错),所以要将周期设为30秒,大多数停站时间都大于30秒,若个别小于30秒则处理成等于30秒。
状态转移变量(B类变量)---外部实时输入变量(BP类变量)
PlatCTIN(来源为实时ATS运能数据):某个周期时刻至前一周期时刻之间的时间段内流入某站台轨道的列车运能(列);
PlatCTOUT(来源为实时ATS运能数据):某个周期时刻至前一周期时刻之间的时间段内流出某站台轨道的列车运能(列);
CellCTIN(来源为实时ATS运能数据):某个周期时刻至前一周期时刻之间的时间段内流入某区间的列车运能(列);
CellCTOUT(来源为实时ATS运能数据):某个周期时刻至前一周期时刻之间的时间段内流出某区间的列车运能(列);
PlatQGT(来源为实时AFC客流数据):某个周期时刻至前一周期时刻之间的时间段内由闸机进站后到达某站台候车区的非换乘乘客(包括通道换乘、出站换乘、同站台换乘)。
状态转移变量(B类变量)---中间变量(BM类变量)
PlatCPIN:某个周期时刻至前一周期时刻之间的时间段内流入某站台轨道的列车运能(人);
PlatCPOUT:某个周期时刻至前一周期时刻之间的时间段内流出某站台轨道的列车运能(人);
CellCPIN:某个周期时刻至前一周期时刻之间的时间段内流入某区间的列车运能(人);
CellCPOUT:某个周期时刻至前一周期时刻之间的时间段内流出某区间的列车运能(人);
PlatQTR:某个周期时刻至前一周期时刻之间的时间段内由各个方向通过通道、同站台、出站方式换乘到某站台候车区的总乘客数;
PlatQTRs:某个周期时刻至前一周期时刻之间的时间段内各个方向通过通道、同站台、出站方式换乘到某站台候车区的乘客数的集合;
PlatQIN:某个周期时刻至前一周期时刻之间的时间段内到达某站台候车区候车的乘客总数,即PlatQGT与PlatQTR之和;
PlatDPT:某个周期时刻至前一周期时刻之间的时间段内由某站台出发进入下游各区间的总人数;
PlatCN:PlatDPT中原该站台区域轨道列车内的乘客数;
PlatON:PlatDPT中由该站台候车区上车的乘客数;
PlatARV:某个周期时刻至前一周期时刻之间的时间段内随列车由上游各区间到达某站台的总乘客数;
PlatU:PlatARV中不下车的乘客数;
PlatD:PlatARV中下车的乘客数;
CellARV:某个周期时刻至前一周期时刻之间的时间段内随列车离开某区间到达下游站台的总人数;
CellU:CellARV中到站后不下车的乘客数;
CellD:CellARV中到站后下车的乘客数;
CellTR:CellARV中到站后下车换乘的乘客数;
CellTRs:CellARV中到站后下车换乘到各个其他站台的乘客数的集合;
CellET:CellARV中到站后下车出站点结束本次出行的乘客数;
CellDPT:某个周期时刻至前一周期时刻之间的时间段内由上游站台出发进入某区间的总乘客数;
CellCN:CellDPT中原上游站台区域轨道列车内的乘客数;
CellON:CellDPT中由上游站台候车区上车的乘客数。
其他静态配置参数
Step:推演周期(步长),默认30秒;由于模型假设一个周期30秒内列车不会穿越任何区间和站台对象的完整范围,因此平滑处理成运能数据之前的原始列车进出站台及进出区间记录的进出时间差最小为一个周期30秒,若不足30秒则将该列车的“出”记录的时刻向后顺延若干秒以满足30秒时间差的要求,由于通常前后两列列车间隔不会小于30秒,因此这样的处理不会打乱前后班次列车的先后顺序;
Nowtime:当前周期时刻(完整日期格式);
Nowsec:当前周期时刻(近似到Step整数倍的当天绝对秒数);
Smooth:运能数据统计的平滑窗口大小(秒),默认20分钟即1200秒;
Moffset:CellM和PlatM的取值间隔,理论上CellM和PlatM的取值只可能是Moffset的整数倍的数值,Moffset=Step/Smooth;
CoachNUM:某条线路的列车编组车厢数;
CoachCAP:某条线路的车厢定员人数;
TrainCAP:某条线路每列列车的定员人数;
TrainCD:某条线路的列车最大饱和度参数,默认值140%;
CellURATE:某个周期时刻至前一周期时刻之间的时间段内随列车离开某区间到达下游站台的乘客中到站后不下车的比例;
CellETRATE:某个周期时刻至前一周期时刻之间的时间段内随列车离开某区间到达下游站台的乘客中到站后出站点结束本次出行的比例;
CellTRRATEs:某个周期时刻至前一周期时刻之间的时间段内随列车离开某区间到达下游站台的乘客中到站后换乘到各个其他站台的比例的集合;以上3组比例参数CellURATE、CellETRATE、CellTRRATEs归一化;
PlatINISTATE:某站台对象初始状态变量的集合;
CellINISTATE:某区间初始状态变量的集合。
对于上述任意一个变量而言,在该变量前增加前缀表示该变量所属的站台或区间,其中,站台用Plat表示,区间用Cell表示;增加后缀,表示该变量所属的周期时刻,其中,now表示当前周期时刻,pre表示上一个周期时刻。例如对于变量CellARV而言,Cell_1.CellARV(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达区间Cell_1所对应的下游站台的总人数,其他变量以此类推。
Step3-A、计算区间Cell_1流出乘客
如图1中一般情况下,计算Cell_1流出乘客:
Cell_1.CellCPOUT(now)=Cell_1.CellCTOUT(now)×Cell_1.TrainCAP;
Cell_1.CellARV(now)=Cell_1.CellCPOUT(now)×Cell_1.CellS(pre);
Cell_1.CellU(now)=Cell_1.CellARV(now)×Cell_1.CellURATE(now);
Cell_1.CellET(now)=Cell_J.CellARV(now)×Cell_1.CellETRATE(now);
令Cell_1.CellTRs(p,now)为当前周期时刻Cell_1换乘往站台p的乘客数,Cell_1.CellTRRATES(p,now)为当前周期时刻Cell_1换乘往站台p的比例系数,p.PlatQTRs(Cell_1,now)为当前周期时刻换乘到站台p的乘客中来自Cell_1的乘客数,则
Cell_1.CellTRs(p,now)
=p.PlatQTRs(Cell_1,now)
=Cell_1.CellARV(now)×Cell_1.CellTRRATES(p,now);
Cell_1.CellTR(now)=∑Cell_1.CellTRs(*,now),式中,∑Cell_1.CellTRs(*,now)表示当前周期时刻离开区间Cell_1换乘往各个不同站台的乘客数的和。
当Cell_1为终点类型区间的特殊情况时:
Cell_1.CellU(now)=0;
其他计算过程同一般情况下计算方法。
Step3-B、计算站台Plat_1流出乘客
如图2所示的第一种情况下,计算站台Plat_1流出乘客:
Plat_1.PlatCPOUT(now)=Plat_1.PlatCTOUT(now)×Plat_1.TrainCAP;
Cell_1.CellCPIN(now)=Cell_1.CellCTIN(now)×Cell_1.TrainCAP;
cpout为Plat_1的载客出站列车运能(人)临时变量,则有:
cpout=Min{Plat_1.PlatCPOUT(now),Cell_1.CellCPIN(now)};
Plat_1.PlatCN(now)=Plat_1.PlatS(pre)×cpout;
Ctmp为Plat_1的出站列车剩余运能(人)临时变量,则有:
Ctmp=cpout×Plat_1.TrainCD-Plat_1.PlatCN(now);
Plat_1.PlatQTR(now)=∑Plat_1.PlatQTRs(*,now);
Plat_1.PlatQIN(now)=Plat_1.PlatQTR(now)+Plat_1.PlatQGT(now)。
若Ctmp<=0,则表示站台候车区乘客全部无法上车:
Plat_1.PlatON(now)=0;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若Ctmp>0且Ctmp>=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now),则表示站台候车区乘客可以全部上车:
Plat_1.PlatON(now)=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now);
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
Ctmp其他取值情况时,表示站台候车区乘客无法全部上车:
Plat_1.PlatON(now)=Ctmp;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now)。
如图3所示的第二种情况下,计算站台Plat_1流出乘客:
Plat_1.PlatCPOUT(now)=Plat_1.PlatCTOUT(now)×Plat_1.TrainCAP;
Cell_1.CellCPIN(now)=Cell_1.CellCTIN(now)×Cell_1.TrainCAP;
Cell_2.CellCPIN(now)=Cell_2.CellCTIN(now)×Cell_2.TrainCAP;
cpout为Plat_1的载客出站列车运能(人)临时变量,则有:
cpout=Min{Plat_1.PlatCPOUT(now),Cell_1.CellCPIN(now)+Cell_1.CellCPIN(now)};
Plat_1.PlatCN(now)=Plat_1.PlatS(pre)×cpout;
Ctmp为Plat_1的出站列车剩余运能(人)临时变量,则有:
Ctmp=cpout×Plat_1.TrainCD-Plat_1.PlatCN(now);
Plat_1.PlatQTR(now)=∑Plat_1.PlatQTRs(*,now);
Plat_1.PlatQIN(now)=Plat_1.PlatQTR(now)+Plat_1.PlatQGT(now)。
若Ctmp<=0,则表示站台候车区乘客全部无法上车:
Plat_1.PlatON(now)=0;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若Ctmp>0且Ctmp>=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now),则表示站台候车区乘客可以全部上车:
Plat_1.PlatON(now)=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now);
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
Ctmp其他取值情况时,表示站台候车区乘客无法全部上车:
Plat_1.PlatON(now)=Ctmp;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now)。
当Plat_1为终点类型站台的特殊情况时:
Plat_1.PlatCN(now)=0;
Plat_1.PlatON(now)=0;
Plat_1.PlatDPT(now)=0。
Step3-C、计算区间Cell_1流入乘客
如图4所示的第一种情况下,计算Cell_1流入乘客:
Cell_1.CellDPT(now)=Plat_1.PlatDPT(now);
Cell_1.CellCN(now)=Plat_1.PlatCN(now);
Cell_1.CellON(now)=Plat_1.PlatON(now)。
如图5所示的第二种情况下,计算Cell_1流入乘客:
若Cell_1.CellCPIN(now)=0且Cell_2.CellCPIN(now)=0,则:
Cell_1.CellDPT(now)=0;
Cell_1.CellCN(now)=0;
Cell_1.CellON(now)=0;
否则:
Cell_1.CellDPT(now)=
Plat_1.PlatDPT(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.CellCPIN(now)};
Cell_1.CellCN(now)=
Plat_1.PlatCN(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.Ce llCPIN(now)};
Cell_1.CellON(now)=
Plat_1.PlatON(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.Ce llCPIN(now)}。
Step3-D、计算站台Plat_1流入乘客
如图6所示的第一种情况下,计算Plat_1流入乘客:
Plat_1.PlatU(now)=Cell_1.CellU(now)。
如图7所示的第二种情况下,计算Plat_1流入乘客:
Plat_1.PlatU(now)=Cell_1.CellU(now)+Cell_2.CellU(now);
当Plat_1为始发类型或终点类型站台的特殊情况时:
Plat_1.PlatU(now)=0。
Step3-E、更新站台Plat_1状态
当站台Plat_1为始发类型站台时:
Plat_1.PlatM(now)=0;
Plat__1.PlatN(now)=0;
Plat_1.PlatS(now)=0;
Plat_1.PlatQ(now)=Max{0,Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now)-Plat_1.PlatO N(now)};
Plat_1.PlatPQ(now)=Plat_1.PlatQ(now)+Plat_1.PlatQIN(now)。
当Plat_1为终点类型站台时:
Plat_1.PlatM(now)=0;
Plat_1.PlatN(now)=0;
Plat_1.PlatS(now)=0;
Plat_1.PlatQ(now)=0;
Plat_1.PlatPQ(now)=0。
当Plat_1为除始发类型及终点类型外的中间类型站台时:
Plat_1.PlatM(now)=Plat_1.PlatM(pre)+Plat_1.PlatCTIN(now)-Plat_1.PlatCTOU T(now);
Plat_1.PlatN(now)=Plat_1.PlatN(pre)+Plat_1.PlatU(now)-Plat_1.PlatCN(now)。
若Plat_1.PlatM(now)>0,则:
Plat_1.PlatS(now)=Plat_1.PlatN(now)/{Plat_1.TrainCAP×Plat_1.PlatM(now)};
若Plat_1.PlatM(now)<=0,则:
Plat_1.PlatS(now)=0。
Plat_1.PlatQ(now)=Max{0,Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now)-Plat_1.PlatO N(now)};
Plat_1.PlatPQ(now)=Plat_1.PlatQ(now)+Plat_1.PlatQIN(now)。
Step3-F、更新区间Cell_1状态
Cell_1.CellM(now)=Cell_1.CellM(pre)+Cell_1.CellCTIN(now)-Cell_1.CellCTO UT(now);
Cell_1.CellN(now)=Cell_1.CellN(pre)+Cell_1.CellDPT(now)-Cell_1.CellARV(now);
若Cell_1.CellM(now)>0,则:
Cell_1.CellS(now)=Cell_1.CellN(now)/{Cell_1.TrainCAP×Cell_1.CellM(now)};
若Cell_1.CellM(now)<=0,则:
Cell_1.CellS(now)=0。
Claims (2)
1.一种轨道交通网络化运行推演算法,其特征在于,包括以下步骤:
第一步、构建轨道网络拓扑结构,在轨道网络拓扑结构中每个轨道交通站点为一个集群,若当前轨道交通站点包含有N条轨道交通线路,N≥1,则当前轨道交通站点包含N个逻辑站点,每个逻辑站点为当前集群下的一个子集群,每个逻辑站点包含一个位于轨道交通上行线路上的站台及一个位于轨道交通下行线路上的站台,则两个站台分别为当前子集群下的两个节点,相邻两个轨道交通站点中通过同一根轨道可直接到达的两个站台在轨道网络拓扑结构中通过一个区间连接起来;
第二步、以一个周期时长为固定间隔循环计算轨道网络拓扑结构中每个站台及每个区间的状态,对于任意一个站台Plat_1及任意一个区间Cell_1而言,其状态的计算步骤包括:
步骤2.1、计算区间Cell_1流出乘客,计算得到如下变量:
Cell_1.CellCPOUT(now)=Cell_1.CellCTOUT(now)×Cell_1.TrainCAP,式中,Cell_1.CellCPOUT(now)表示当前周期时刻至前一周期时刻之间的时间段内流出区间Cell_1的列车所携带的人数,Cell_1.CellCTOUT(now)表示当前周期时刻至前一周期时刻之间的时间段内流出区间Cell_1的列车的列数,Cell_1.TrainCAP表示区间Cell_1所对应的轨道交通线路线路每列列车的定员人数;
Cell_1.CellARV(now)=Cell_1.CellCPOUT(now)×Cell_1.CellS(pre),式中,Cell_1.CellARV(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达区间Cell_1所对应的下游站台的总人数,Cell_1.CellS(pre)表示上一周期时刻区间Cell_1内的列车平均乘客饱和度;
Cell_1.CellU(now)=Cell_1.CellARV(now)×Cell_1.CellURATE(now),式中,Cell_1.CellU(now)表示当前周期时刻Cell_1.CellARV(now)中到站后不下车的乘客数,Cell_1.CellURATE(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达下游站台的乘客中到站后不下车的比例;
Cell_1.CellET(now)=Cell_1.CellARV(now)×Cell_1.CellETRATE(now),式中,Cell_1.CellET(now)表示Cell_1.CellARV(now)中当前周期时刻到站后下车出站点结束本次出行的乘客数,Cell_1.CellETRATE(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_1到达下游站台的乘客中到站后出站点结束本次出行的比例;
Cell_1.CellTRs(p,now)
=p.PlatQTRs(Cell_1,now)
=Cell_1.CellARV(now)×Cell_1.CellTRRATES(p,now),式中:
Cell_1.CellTRs(p,now)表示当前周期时刻离开区间Cell_1换乘往站台p的乘客数,p.PlatQTRs(Cell_1,now)表示当前周期时刻换乘到站台p的乘客中来自区间Cell_1的乘客数,Cell_1.CellTRRATES(p,now)表示当前周期时刻离开区间Cell_1换乘往站台p的比例系数;
Cell_1.CellTR(now)=∑Cell_1.CellTRs(*,now),式中,Cell_1.CellTR(now)表示当前周期时刻Cell_1.CellARV(now)中到站后下车换乘的乘客数,∑Cell_1.CellTRs(*,now)表示当前周期时刻离开区间Cell_1换乘往各个不同站台的乘客数的和;
步骤2.2、计算站台Plat_1流出乘客,若站台Plat_1仅为一个区间Cell_11的上游站台,则计算得到如下变量:
Plat_1.PlatCPOUT(now)=Plat_1.PlatCTOUT(now)×Plat_1.TrainCAP,式中,Plat_1.PlatCPOUT(now)表示当前周期时刻至前一周期时刻之间的时间段内流出站台Plat_1轨道的列车所运载的人数,Plat_1.PlatCTOUT(now)当前周期时刻至前一周期时刻之间的时间段内流出站台Plat_1轨道的列车的列数,Plat_1.TrainCAP表示站台Plat_1所对应的轨道交通线路线路每列列车的定员人数;
Cell_11.CellCPIN(now)=Cell_11.CellCTIN(now)×Cell_11.TrainCAP,式中,Cell_11.CellCPIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_11的列车所运载的人数,Cell_11.CellCTIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_11的列车的列数,Cell_11.TrainCAP表示区间Cell_11所对应的轨道交通线路线路每列列车的定员人数;
cpout=Min{Plat_1.PlatCPOUT(now),Cell_11.CellCPIN(now)},式中,cpout表示站台Plat_1的载客出站列车所运载的人数的临时变量;
Plat_1.PlatCN(now)=Plat_1.PlatS(pre)×cpout,式中,Plat_1.PlatCN(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat_1出发进入下游各区间的总人数中原该站台区域轨道列车内的乘客数,Plat_1.PlatS(pre)表示上一周期时刻站台Plat_1区域轨道上的列车平均乘客饱和度,到站列车需下车乘客已下车而站台候车乘客还未上车时的饱和度;
Ctmp=cpout×Plat_1.TrainCD-Plat_1.PlatCN(now),式中,Ctmp表示站台Plat_1的出站列车剩余的能运载的人数的临时变量;
Plat_1.PlatQTR(now)=∑Plat_1.PlatQTRs(*,now),式中,Plat_1.PlatQTR(now)表示当前周期时刻至前一周期时刻之间的时间段内由各个方向通过通道、同站台、出站方式换乘到某站台候车区的总乘客数,∑Plat_1.PlatQTRs(*,now)表示当前周期时刻至前一周期时刻之间的时间段内各个方向通过通道、同站台、出站方式换乘到站台Plat_1候车区的乘客数的和;
Plat_1.PlatQIN(now)=Plat_1.PlatQTR(now)+Plat_1.PlatQGT(now),式中,Plat_1.PlatQIN(now)表示当前周期时刻至前一周期时刻之间的时间段内到达站台Plat_1候车区候车的乘客总数,Plat_1.PlatQGT(now)表示当前周期时刻至前一周期时刻之间的时间段内由闸机进站后到达站台Plat_1候车区的非换乘乘客;
若Ctmp<=0,则表示站台候车区乘客全部无法上车,则:
Plat_1.PlatON(now)=0,式中,Plat_1.PlatON(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat_1出发进入下游各区间的总人数中由该站台候车区上车的乘客数;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now),式中,Plat_1.PlatDPT(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat_1出发进入下游各区间的总人数;
若Ctmp>0且Ctmp>=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now),则表示站台候车区乘客可以全部上车,式中,Plat_1.PlatQ(pre)表示上一周期时刻站台Plat_1区域候车区内由于列车过饱和而无法及时上车的乘客数:
Plat_1.PlatON(now)=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now);
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
Ctmp其他取值情况时,表示站台候车区乘客无法全部上车:
Plat_1.PlatON(now)=Ctmp;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若站台Plat_1为区间Cell_11及区间Cell_21的上游站台,则计算得到如下变量:
Plat_1.PlatCPOUT(now)=Plat_1.PlatCTOUT(now)×Plat_1.TrainCAP;
Cell_11.CellCPIN(now)=Cell_11.CellCTIN(now)×Cell_11.TrainCAP;
Cell_21.CellCPIN(now)=Cell_21.CellCTIN(now)×Cell_21.TrainCAP,式中,Cell_21.CellCPIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_21的列车所运载的人数,Cell_21.CellCTIN(now)表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_21的列车的列数,Cell_21.TrainCAP表示区间Cell_21所对应的轨道交通线路线路每列列车的定员人数;
cpout=Min{Plat_1.PlatCPOUT(now),Cell_1.CellCPIN(now)+Cell_1.CellCPIN(now)};
Plat_1.PlatCN(now)=Plat_1.PlatS(pre)×cpout;
Ctmp=cpout×Plat_1.TrainCD-Plat_1.PlatCN(now);
Plat_1.PlatQTR(now)=∑Plat_1.PlatQTRs(*,now);
Plat_1.PlatQIN(now)=Plat_1.PlatQTR(now)+Plat_1.PlatQGT(now);
若Ctmp<=0,则表示站台候车区乘客全部无法上车,
Plat_1.PlatON(now)=0;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若Ctmp>0且Ctmp>=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now),则表示站台候车区乘客可以全部上车:
Plat_1.PlatON(now)=Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now);
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
Ctmp其他取值情况时,表示站台候车区乘客无法全部上车:
Plat_1.PlatON(now)=Ctmp;
Plat_1.PlatDPT(now)=Plat_1.PlatON(now)+Plat_1.PlatCN(now);
若站台Plat_1的下游未与任何区间相连,则:
Plat_1.PlatCN(now)=0;Plat_1.PlatON(now)=0;Plat_1.PlatDPT(now)=0;
步骤2.3、计算区间Cell_1流入乘客,若区间Cell_1为位于其上游的站台Plat的唯一一条区间,则计算得到如下变量:
Cell_1.CellDPT(now)=Plat.PlatDPT(now),式中,Cell_1.CellDPT(now)表示当前周期时刻至前一周期时刻之间的时间段内由上游站台Plat出发进入区间Cell_1的总乘客数,Plat.PlatDPT(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat出发进入下游各区间的总人数;
Cell_1.CellCN(now)=Plat.PlatCN(now),式中,Cell_1.CellCN(now)表示Cell_1.CellDPT(now)中原上游站台区域轨道列车内的乘客数,Plat.PlatCN(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat出发进入下游各区间的总人数中原该站台区域轨道列车内的乘客数
Cell_1.CellON(now)=Plat.PlatON(now),式中,Cell_1.CellON(now)表示Cell_1.CellDPT(now)中由上游站台候车区上车的乘客数,Plat.PlatON(now)表示当前周期时刻至前一周期时刻之间的时间段内由站台Plat出发进入下游各区间的总人数中由该站台候车区上车的乘客数;
若位于区间Cell_1上游的站台Plat有两条相连的区间,分别为区间Cell_1及区间Cell_2,则若Cell_1.CellCPIN(now)=0且Cell_2.CellCPIN(now)=0,式中,Cell_1.CellCPIN(now)及Cell_2.CellCPIN(now)=0分别表示当前周期时刻至前一周期时刻之间的时间段内流入区间Cell_1及区间Cell_2的列车所运载的人数,则:
Cell_1.CellDPT(now)=0;
Cell_1.CellCN(now)=0;
Cell_1.CellON(now)=0;
否则:
Cell_1.CellDPT(now)=
Plat.PlatDPT(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.CellCPIN(now)};
Cell_1.CellCN(now)=
Plat.PlatCN(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.CellCPIN(now)};
Cell_1.CellON(now)=
Plat.PlatON(now)×Cell_1.CellCPIN(now)/{Cell_1.CellCPIN(now)+Cell_2.CellCPIN(now)};
步骤2.4、计算站台Plat_1的流入乘客,若站台Plat_1的上游仅有一条区间Cell_1,则:
Plat_1.PlatU(now)=Cell_1.CellU(now),式中,Plat_1.PlatU(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车由上游各区间到达站台Plat_1的总乘客数中不下车的乘客数;
若站台Plat_1的上游有两条区间,即区间Cell_1及区间Cell_2,则:
Plat_1.PlatU(now)=Cell_1.CellU(now)+Cell_2.CellU(now),式中,Cell_2.CellU(now)表示当前周期时刻Cell_2.CellARV(now)中到站后不下车的乘客数,Cell_2.CellARV(now)表示当前周期时刻至前一周期时刻之间的时间段内随列车离开区间Cell_2到达区间Cell_2所对应的下游站台的总人数;
若站台Plat_1没有对应的上游区间,或没有对应的下游区间,则:
Plat_1.PlatU(now)=0;
步骤2.5、更新站台Plat_1的状态:
当站台Plat_1没有对应的上游区间,即为始发类型站台时:
Plat_1.PlatM(now)=0,式中,Plat_1.PlatM(now)表示当前周期时刻站台Plat_1区域轨道内的列车数;
Plat_1.PlatN(now)=0,式中,Plat_1.PlatN(now)表示当前周期时刻站台Plat_1区域轨道列车上的在途乘客总数;
Plat_1.PlatS(now)=0,式中,Plat_1.PlatS(now)表示当前周期时刻站台Plat_1区域轨道上的列车平均乘客饱和度,到站列车需下车乘客已下车而站台候车乘客还未上车时的饱和度;
Plat_1.PlatQ(now)=Max{0,Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now)-Plat_1.PlatON(now)},式中,Plat_1.PlatQ(now)表示当前周期时刻站台Plat_1区域候车区内由于列车过饱和而无法及时上车的乘客数;
Plat_1.PlatPQ(now)=Plat_1.PlatQ(now)+Plat_1.PlatQIN(now),式中,Plat_1.PlatPQ(now)表示当前周期时刻站台Plat_1区域候车区内滞留乘客及正常进站候车乘客的总和,即站台Plat_1候车区乘客总数;
当站台Plat_1没有对应的下游区间,即为终点类型站台时:
Plat_1.PlatM(now)=0;
Plat_1.PlatN(now)=0;
Plat_1.PlatS(now)=0;
Plat_1.PlatQ(now)=0;
Plat_1.PlatPQ(now)=0;
当站台Plat_1为除终点类型及始发类型外的其他类型的站台时:
Plat_1.PlatM(now)=Plat_1.PlatM(pre)+Plat_1.PlatCTIN(now)-Plat_1.PlatCTOUT(now);
Plat_1.PlatN(now)=Plat_1.PlatN(pre)+Plat_1.PlatU(now)-Plat_1.PlatCN(now);
Plat_1.PlatQ(now)=Max{0,Plat_1.PlatQ(pre)+Plat_1.PlatQIN(now)-Plat_1.PlatON(now)};
Plat_1.PlatPQ(now)=Plat_1.PlatQ(now)+Plat_1.PlatQIN(now);
若Plat_1.PlatM(now)>0,则:
Plat_1.PlatS(now)=Plat_1.PlatN(now)/{Plat_1.TrainCAP×Plat_1.PlatM(now)};
若Plat_1.PlatM(now)<=0,则:
Plat_1.PlatS(now)=0;
步骤2.6、更新区间Cell_1的状态:
Cell_1.CellM(now)=Cell_1.CellM(pre)+Cell_1.CellCTIN(now)-Cell_1.CellCTOUT(now),式中,Cell_1.CellM(now)表示当前周期时刻区间Cell_1内的列车数,Cell_1.CellM(pre)表示上一个周期时刻区间Cell_1内的列车数;
Cell_1.CellN(now)=Cell_1.CellN(pre)+Cell_1.CellDPT(now)-Cell_1.CellARV(now),式中,Cell_1.CellN(now)表示当前周期时刻区间Cell_1内列车上的在途乘客总数,Cell_1.CellN(pre)表示上一个周期时刻区间Cell_1内列车上的在途乘客总数;
若Cell_1.CellM(now)>0,则:
Cell_1.CellS(now)=Cell_1.CellN(now)/{Cell_1.TrainCAP×Cell_1.CellM(now)},式中,Cell_1.CellS(now)表示当前周期时刻区间Cell_1内的列车平均乘客饱和度;
若Cell_1.CellM(now)<=0,则:
Cell_1.CellS(now)=0。
2.如权利要求1所述的一种轨道交通网络化运行推演算法,其特征在于,在所述步骤2.1中,若区间Cell_1没有任何与其直接相连的下游站台,则将所述Cell_1.CellU(now)设为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610329595.4A CN106021918B (zh) | 2016-05-18 | 2016-05-18 | 一种轨道交通网络化运行推演算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610329595.4A CN106021918B (zh) | 2016-05-18 | 2016-05-18 | 一种轨道交通网络化运行推演算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021918A CN106021918A (zh) | 2016-10-12 |
CN106021918B true CN106021918B (zh) | 2019-03-08 |
Family
ID=57097785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610329595.4A Active CN106021918B (zh) | 2016-05-18 | 2016-05-18 | 一种轨道交通网络化运行推演算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021918B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945355A (zh) * | 2017-11-29 | 2018-04-20 | 中铁程科技有限责任公司 | 信息处理方法及装置、计算机可读存储介质 |
CN108132913B (zh) * | 2017-12-11 | 2021-10-29 | 深圳先进技术研究院 | 一种轨道交通客流移动估计方法、系统及电子设备 |
CN113184024B (zh) * | 2021-05-21 | 2023-02-24 | 北京全路通信信号研究设计院集团有限公司 | 一种面向有轨电车的到站时间静态预测方法及其系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984993A (zh) * | 2014-05-13 | 2014-08-13 | 东南大学 | 一种轨道交通客流od分布实时推测方法 |
CN105224999A (zh) * | 2015-09-10 | 2016-01-06 | 北京市交通信息中心 | 基于afc数据的城市轨道交通实时客流预测方法及系统 |
-
2016
- 2016-05-18 CN CN201610329595.4A patent/CN106021918B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984993A (zh) * | 2014-05-13 | 2014-08-13 | 东南大学 | 一种轨道交通客流od分布实时推测方法 |
CN105224999A (zh) * | 2015-09-10 | 2016-01-06 | 北京市交通信息中心 | 基于afc数据的城市轨道交通实时客流预测方法及系统 |
Non-Patent Citations (3)
Title |
---|
城市轨道交通线路实时客流推算方法研究;孙海瑞;《中国优秀硕士学位论文全文数据库(电子期刊)》;20121015;C033-142 |
城市轨道交通网络化乘客信息服务系统研究及应用;沈峰等;《第九届中国智能交通年会》;20150528;1096-1103 |
基于票务信息的城市轨道交通客流实时监测系统的研究;刘海涛;《中国优秀硕士学位论文全文数据库(电子期刊)工程科技Ⅱ辑》;20130115;C033-56 |
Also Published As
Publication number | Publication date |
---|---|
CN106021918A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110866649B (zh) | 一种短期地铁客流预测方法、系统及电子设备 | |
Wang et al. | Passenger-demands-oriented train scheduling for an urban rail transit network | |
Zaki et al. | Online bus arrival time prediction using hybrid neural network and Kalman filter techniques | |
CN104183119B (zh) | 基于路段od反推的实时交通流分布预测方法 | |
Treethidtaphat et al. | Bus arrival time prediction at any distance of bus route using deep neural network model | |
CN106021918B (zh) | 一种轨道交通网络化运行推演算法 | |
CN109872535A (zh) | 一种智慧交通通行预测方法、装置及服务器 | |
Gong et al. | Sequential framework for short-term passenger flow prediction at bus stop | |
D'Acierno et al. | The definition of a model framework for managing rail systems in the case of breakdowns | |
CN107025805A (zh) | 基于候补机制的协同决策系统时刻计算方法及对应系统 | |
Zhang et al. | Coupling analysis of passenger and train flows for a large-scale urban rail transit system | |
CN103927433A (zh) | 基于流控制的多源点多出口人群应急疏散调度方法 | |
Ma et al. | Dynamic system optimal routing in multimodal transit network | |
Sarango et al. | Mathematical model for predicting the vehicle occupancy for urban passenger transport on the route | |
Fanti et al. | A Petri net model for fleet sizing of electric car sharing systems | |
CN110415518A (zh) | 一种客流与运力的监控预警方法、装置、设备及介质 | |
Vateekul et al. | Recurrent Neural‐Based Vehicle Demand Forecasting and Relocation Optimization for Car‐Sharing System: A Real Use Case in Thailand | |
CN108734345A (zh) | Ota机票客服话务量的预测方法及系统 | |
Di Febbraro et al. | INTRANET: A new simulation tool for intermodal transportation systems | |
CN107886190B (zh) | 一种突发事件下乘客出行选择行为仿真方法及装置 | |
Jin et al. | Route planning of rescue vehicles in the process of dynamic change of traffic volume under emergency conditions | |
Zhang et al. | Dynamic multi-bus dispatching strategy with boarding and holding control for passenger delay alleviation and schedule reliability: A combined dispatching-operation system | |
CN112182838A (zh) | 一种基于afc数据的短时客流预测动态模型 | |
SHARIAT et al. | Creating bus timetables under stochastic demand | |
Nguyen et al. | Scheduling problem for bus rapid transit routes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |