CN108876064B - 基于客流量的城市轨道交通网络关键车站识别方法 - Google Patents
基于客流量的城市轨道交通网络关键车站识别方法 Download PDFInfo
- Publication number
- CN108876064B CN108876064B CN201810967392.7A CN201810967392A CN108876064B CN 108876064 B CN108876064 B CN 108876064B CN 201810967392 A CN201810967392 A CN 201810967392A CN 108876064 B CN108876064 B CN 108876064B
- Authority
- CN
- China
- Prior art keywords
- passenger
- station
- passengers
- time
- delta
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012546 transfer Methods 0.000 claims abstract description 49
- 238000007781 pre-processing Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000006399 behavior Effects 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims 1
- 238000005070 sampling Methods 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 44
- 238000012423 maintenance Methods 0.000 abstract description 4
- 238000010606 normalization Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000000630 rising effect Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
本发明提供了一种基于客流量的城市轨道交通网络关键车站识别方法,属于城市轨道交通运行控制技术领域。首先收集AFC设备中的乘客数据并进行预处理,获取OD对数据、进站乘客量和出站乘客量;然后结合交通网络拓扑结构,根据OD对数据对乘客行为进行仿真,获取站点的换乘乘客量和运送乘客量;最后根据进站乘客量、出站乘客量、换乘乘客量和运送乘客量获取站点的重要性程度。本发明考虑每个站点的进入乘客数量、离开乘客数量、转乘乘客数量和运送乘客数量这四项指标,使用基于多智能体的仿真和最小—最大归一化的分析方法得到每个站点的重要性程度,以此识别城市轨道交通网络中的关键站点,对城市轨道交通管理人员的维护工作提供理论支持。
Description
技术领域
本发明涉及轨道交通列车运行控制技术领域,具体涉及一种基于客流量的城市轨道交通网络关键车站识别方法。
背景技术
城市轨道交通运输由于其可靠性高、容量大、污染小等优点,在全世界范围内得到了飞速的发展。大城市的轨道交通运输已成为一个大型复杂的网络,在城市交通系统中发挥着主要作用。城市轨道交通网络中存在着一些对整个网络正常运行起着重要作用的关键站点,识别这些站点并对其投入更多的人力物力来管理和维护,对保证整个城市轨道交通网络的正常运行有着重要的意义。
识别关键站点的传统研究方法是基于交通网络的拓扑结构,未能考虑到各个站点的实际客流数据,不能根据实际客流量来确定实际的关键站点。
发明内容
本发明的目的在于提供一种能够根据每一站点的实际客流量,确定城轨交通系统中的关键站点,为优化轨道交通提供依据的基于客流量的城市轨道交通网络关键车站识别方法,以解决上述背景技术中存在的技术问题。
为了实现上述目的,本发明采取了如下技术方案:
本发明提供的一种基于客流数据的城市轨道交通网络关键站点识别方法,包括如下步骤:
步骤S110:收集AFC设备中的乘客数据并进行预处理,获取OD对数据、进站乘客量和出站乘客量;
步骤S120:结合交通网络拓扑结构,根据所述OD对数据对乘客行为进行仿真,获取站点的换乘乘客量和运送乘客量;
步骤S130:根据所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量获取站点的重要性程度;
进一步的,所述步骤S120具体包括:
根据所述OD对数据确定一个乘客集合,根据所述乘客集合中的乘客状态仿真模拟所述乘客的路径信息,根据所述路径信息确定所述换乘乘客量和所述运送乘客量。
进一步的,所述乘客集合中的乘客状态包括:
Sδ=0,表示乘客δ刚刚进入城市轨道交通系统;
Sδ=1,表示乘客δ刚刚从进站口进入一个站点,正走向站台准备等车;
Sδ=2,表示乘客δ刚到站台,开始等车;
Sδ=3,表示乘客δ刚乘上车;
Sδ=4,表示乘客δ刚在换乘站下车,开始换乘另一条城市轨道交通线路;
Sδ=5,表示乘客δ刚到目的站点下车,准备检票出站;
Sδ=6,表示乘客δ已经完成检票出站,并已退出城市轨道交通系统。
进一步的,所述根据所述乘客集合中的乘客状态仿真模拟所述乘客的出行路径包括:
步骤1.1:初始化参数tb、te、Δt和Td,设置时钟t=tb,第k次仿真周期的时间区间为[tb+(k-1)×Td,tb+k×Td],k=1,2,3...K-1,
其中,tb表示仿真开始的时间,te表示仿真结束的时间,Δt表示仿真步长,Td表示每个仿真的仿真周期,K=(te-tb)÷Td表示仿真周期的个数;
步骤1.2:第k个仿真中的乘客集合为{δ(O,D)k},对其中的任意一个乘客δ(O,D)k,设置初始状态Sδ=0,将{δ(O,D)k}添加至乘客总体£;其中,
定义(O,D)表示一个从站点O到站点N的OD对,且所述第k个仿真满足 t=ts+(k-1)×Td,t<ts+k×Td,t=t+Δt;
步骤1.3:当Sδ=0时,设置tδ=aδn,n站点表示起始站点O,将站点n 增加至路径信息pδ中,判断tδ∈[t,t+Δt]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=1,每重新计算一次,k增加1;其中,tδ表示第k个仿真采样时的乘客的当前时间,aδn表示乘客进入站点n的时间;
当Sδ=1时,判断tδ∈[t,t+Δt]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=2,否则,设置tδ=tδ+t(0,1)n,每重新计算一次,k增加1;其中, t(0,1)n表示乘客从站点n的进站口到站台所需的步行时间;
当Sδ=2时,根据列车时刻表得到即将驶入站点n的列车为i,判断 tδ∈[ain,din]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=3,每重新计算一次,k增加1;其中,ain表示列车i到达站点n的时间,din表示列车i离开站点 n的时间;
当Sδ=3时,识别列车i与乘客路径相同的最后一站m1,设置n=m1、 tδ=am1i,将站点m1作为换乘站点添加至路径信息pδ中,判断tδ∈[ain,din]且 n=D是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=5,否则,乘客δ(O,D)k进入状态Sδ=4,每重新计算一次,k增加1;其中,am1i表示列车i到达m1站点的时刻,D表示最终站点;
当Sδ=4时,寻找乘客δ(O,D)k的目标站点m2,设置n=m2、 tδ=tδ+t(4,3)n,将站点m2作为终点站添加至路径信息pδ中,判断tδ∈[t,t+Δt] 是否成立,若成立,则乘客δ(O,D)k的状态重新设置为Sδ=2;其中,t(4,3)n表示乘客在站点m1的换乘时间;
当Sδ=5时,设置tδ=tδ+t(5,6)n,判断tδ∈[t,t+Δt]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=6,每重新计算一次,k增加1;其中,t(5,6)n表示乘客从站点m2的站台到出站口所需的步行时间;
当Sδ=6时,计算结束,将乘客δ(O,D)k从乘客总体£中移除;
步骤1.4:根据步骤1.3的计算方法,确定集合{δ(O,D)k}中所有乘客的出行路径。
进一步的,所述步骤S130具体包括:
使用最小—最大值方程对每个站点的所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量进行归一化:
其中,xn表示进入站点n的日均乘客数量,Xmax表示进入一个站点的最大日均乘客数量;Xmin表示进入一个站点的最小日均乘客数量;yn表示离开站点n的日均乘客数量;Ymax表示离开一个站点的最大日均乘客数量;Ymin表示离开一个站点的最小日均乘客数量;un表示站点n日均转乘乘客数量;Umax表示一个站点的最大日均转乘乘客数量;Umin表示一个站点的最小日均转乘乘客数量;wn表示站点n的日均运送乘客数量;Wmax表示一个站点的最大日均运送乘客数量;Wmin表示一个站点的最小日均运送乘客数量。
进一步的,根据所述归一化的加权求和,确定站点的重要性程度:
其中,c1、c2、c3和c4分别表示所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量的加权系数。
本发明有益效果:从实际客流量的角度出发,考虑每个站点的进入乘客数量、每个站点的离开乘客数量、每个站点的转乘乘客数量、每个站点的运送乘客数量这四项指标,使用基于多智能体的仿真和最小—最大归一化的分析方法得到每个站点的重要性程度,以此识别城市轨道交通网络中的关键站点,对城市轨道交通管理人员的维护工作提供理论支持。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一所述的基于客流数据的城市轨道交通网络关键车站识别方法流程图。
图2为本发明实施例二所述的基于客流数据的城市轨道交通网络关键车站识别方法流程图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的模块。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或模块,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、模块和/或它们的组。
需要说明的是,在本发明所述的实施例中,除非另有明确的规定和限定,术语“连接”、“固定”等应做广义理解,可以是固定连接,也可以是可拆卸连接,或成一体,可以是机械连接,也可以是电连接,可以是直接连接,也可以是通过中间媒介间接连接,可以是两个元件内部的连通,或两个元件的相互作用关系,除非具有明确的限定。对于本领域技术人员而言,可以根据具体情况理解上述术语在本发明实施例中的具体含义。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语 (包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
为便于对本发明实施例的理解,下面将结合附图以具体实施例为例做进一步的解释说明,且实施例并不构成对本发明实施例的限定。
本领域普通技术人员应当理解的是,附图只是一个实施例的示意图,附图中的部件或装置并不一定是实施本发明所必须的。
实施例一
如图1所示,本发明实施例一提供一种基于客流量的城市轨道交通网络关键站点识别方法,该方法包括如下步骤:
步骤S110:收集AFC设备中的乘客数据并进行预处理,获取OD对数据、进站乘客量和出站乘客量;
步骤S120:结合交通网络拓扑结构,根据所述OD对数据对乘客行为进行仿真,获取站点的换乘乘客量和运送乘客量;
步骤S130:根据所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量获取站点的重要性程度。
在具体的实施例一中,所述步骤S120具体包括:
根据所述OD对数据确定一乘客集合,根据所述乘客集合中的乘客的乘客状态仿真模拟所述乘客的路径信息,根据所述路径信息确定所述换乘乘客量和所述运送乘客量。
在具体实施例一中,所述乘客集合中的乘客的乘客状态包括:
Sδ=0,表示乘客δ刚刚进入城市轨道交通系统;
Sδ=1,表示乘客δ刚刚从进站口进入一个站点,正走向站台准备等车;
Sδ=2,表示乘客δ刚到站台,开始等车;
Sδ=3,表示乘客δ刚乘上车;
Sδ=4,表示乘客δ刚在换乘站下车,开始换乘另一条城市轨道交通线路;
Sδ=5,表示乘客δ刚到目的站点下车,准备检票出站;
Sδ=6,表示乘客δ已经完成检票出站,并已退出城市轨道交通系统。
在本发明具体实施例一中,所述根据所述乘客集合中的乘客的乘客状态仿真模拟所述乘客的出行路径包括:
步骤1.1:初始化仿真参数tb、te、Δt和Td,设置仿真时钟t=tb,第k次仿真周期的时间区间为[tb+(k-1)×Td,tb+k×Td],k=1,2,3...K-1,
其中,tb表示仿真开始的时间,te表示仿真结束的时间,Δt表示仿真步长,Td表示每个仿真的仿真周期,K=(te-tb)÷Td表示仿真周期的个数;
步骤1.2:第k个仿真中的乘客集合为{δ(O,D)k},对其中的任意一个乘客δ(O,D)k,设置初始状态Sδ=0,将{δ(O,D)k}添加至乘客总体£;其中,
定义(O,D)表示一个从站点O到站点N的OD对,且所述第k个仿真满足 t=ts+(k-1)×Td,t<ts+k×Td,t=t+Δt;
步骤1.3:当Sδ=0时,设置tδ=aδn,n站点表示起始站点O,将站点n 增加至路径信息pδ中,判断tδ∈[t,t+Δt]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=1,每重新计算一次,k增加1;其中,tδ表示第k个仿真采样时的乘客的当前时间,aδn表示乘客进入站点n的时间;
当Sδ=1时,判断tδ∈[t,t+Δt]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=2,否则,设置tδ=tδ+t(0,1)n,每重新计算一次,k增加1;其中, t(0,1)n表示乘客从站点n的进站口到站台所需的步行时间;
当Sδ=2时,根据列车时刻表得到即将驶入站点n的列车为i,判断 tδ∈[ain,din]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=3,每重新计算一次,k增加1;其中,ain表示列车i到达站点n的时间,din表示列车i离开站点 n的时间;
当Sδ=3时,识别列车i与乘客路径相同的最后一站m1,设置n=m1、 tδ=am1i,将站点m1作为换乘站点添加至路径信息pδ中,判断tδ∈[ain,din]且 n=D是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=5,否则,乘客δ(O,D)k进入状态Sδ=4,每重新计算一次,k增加1;其中,am1i表示列车i到达m1站点的时刻,D表示最终站点;
当Sδ=4时,寻找乘客δ(O,D)k的目标站点m2,设置n=m2、 tδ=tδ+t(4,3)n,将站点m2作为终点站添加至路径信息pδ中,判断tδ∈[t,t+Δt] 是否成立,若成立,则乘客δ(O,D)k的状态重新设置为Sδ=2;其中,t(4,3)n表示乘客在站点m1的换乘时间;
当Sδ=5时,设置tδ=tδ+t(5,6)n,判断tδ∈[t,t+Δt]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=6,每重新计算一次,k增加1;其中,t(5,6)n表示乘客从站点m2的站台到出站口所需的步行时间;
当Sδ=6时,计算结束,将乘客δ(O,D)k从乘客总体£中移除;
步骤1.4:根据步骤1.3的计算方法,确定集合{δ(O,D)k}中所有乘客的出行路径。
在本发明的具体实施例一中,所述步骤S130具体包括:
使用最小—最大值方程对每个站点的所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量进行归一化:
其中,xn表示进入站点n的日均乘客数量,Xmax表示进入一个站点的最大日均乘客数量;Xmin表示进入一个站点的最小日均乘客数量;yn表示离开站点n的日均乘客数量;Ymax表示离开一个站点的最大日均乘客数量;Ymin表示离开一个站点的最小日均乘客数量;un表示站点n日均转乘乘客数量;Umax表示一个站点的最大日均转乘乘客数量;Umin表示一个站点的最小日均转乘乘客数量;wn表示站点n的日均运送乘客数量;Wmax表示一个站点的最大日均运送乘客数量;Wmin表示一个站点的最小日均运送乘客数量。
在本发明的具体实施例一中,根据所述归一化的加权求和,确定站点的重要性程度:
其中,c1、c2、c3和c4分别表示所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量的加权系数。
实施例二
本发明实施例二提供的一种基于客流量的城市轨道交通网络关键车站识别方法,该方法包括如下流程步骤:
步骤1,对收集得到的智能卡数据进行预处理;
位于每个站点入口和出口的交通智能卡自动收费(AFC)设备可以收集全部的乘客数据,包括乘客进站时间、出站时间和每对起止站点的OD对。对智能卡数据进行预处理可以直接得到OD对数据、每个站点的进站乘客数量和每个站点的出站乘客数量。
步骤2,得到每个站点的换乘乘客数量和运送乘客数量;
因为换乘过程的不确定性,所以不能从智能卡数据中直接获得每个站点的换乘乘客数量和运送乘客数量。因此,在本发明实施例二中,提出了一个基于多智能体的仿真过程,综合了OD数据、时间表图和交通网络拓扑结构,准确地计算得到每个站点的换乘乘客数量和运送乘客数量;此仿真可以近似模拟乘客在路线选择中的行为,并通过实际的人工测量数据和大量的历史数据对所提的仿真参数进行验证。
在本发明实施例二中,上述基于多智能体的仿真过程使用如下参数进行计算:
定义N来表示一组有限的节点(即站点),并定义(O,D)来表示一个从站点O到站点N的OD对,其中O,D∈N,定义一个乘客δ,用sδ来表示其状态,该乘客δ有七种不同的状态,即:
sδ=0表示乘客δ刚刚进入城市轨道交通系统;
sδ=1表示乘客δ刚刚从进站口进入一个站点,正走向站台准备等车;
sδ=2表示乘客δ刚到站台,开始等车;
sδ=3表示乘客δ刚乘上车;
sδ=4表示乘客δ刚在换乘站下车,开始换乘另一条城市轨道交通线路;
sδ=5表示乘客δ刚到目的站点下车,准备检票出站;
sδ=6表示乘客δ已经完成检票出站,并已退出城市轨道交通系统。
在实际应用中,并非所有的乘客旅程都包含这7个状态,也就是说,在一个乘客旅途中可能没有发生过换乘或进行过不止一次换乘。
此外,在仿真过程中也使用了以下参数和变量,即:
t表示仿真时钟;
Δt表示仿真步长;
tb表示仿真开始的时间;
te表示仿真结束的时间;
Td表示每个仿真周期的时间长短;
K表示仿真周期的数量,即K=(te-tb)/Td;
k表示仿真指标,对于每个k=1,2,…K-1,第k次仿真周期的时间区间为 [ts+(k-1)Td,ts+kTd);
δ(O,D)k表示某一个在第k次仿真周期内进入城市轨道交通系统的乘客,起点为站点O,终点为站点D,其中O,D∈N,k=1,2,…,K-1;
{δ(O,D)k}表示在第k次仿真周期内进入城市轨道交通系统的所有乘客δ(O,D)k的集合;
£表示整个仿真过程中所有乘客的总体;
tδ表示乘客δ的当前时间;
n表示乘客δ当前所在站点;
aδn表示乘客δ进入站点n的时间;
t(0,1)n表示站点n从进站口到站台所需的平均步行时间;
t(4,3)n表示换乘站n的平均换乘时间;
t(5,6)n表示站点n从站台到出站口所需的平均步行时间;
i表示列车的编号;
ain表示列车i到达站点n的时间;
din表示列车i驶出站点n的时间;
pδ表示乘客δ的路径信息,即乘客δ经过的站点。
具体的,在本发明的实施例二中,上述城市轨道交通系统中乘客基于多智能体仿真过程计算程序如下:
初始化:初始化仿真参数tb,te,Δt和Td。设置仿真时钟t=tb,设置仿真指数 k=1。
过程:
步骤3,对每个站点的四项指标进行归一化;
因为这四项指标的影响程度是不同的,以最小—最大值方程对每个站点的四项指标进行归一化,具体的参数和符号定义如下:
xn表示进入站点n的日均乘客数量;
Xmax表示进入一个站点的最大日均乘客数量;
Xmin表示进入一个站点的最小日均乘客数量;
yn表示离开站点n的日均乘客数量;
Ymax表示离开一个站点的最大日均乘客数量;
Ymin表示离开一个站点的最小日均乘客数量;
un表示站点n日均转乘乘客数量;
Umax表示一个站点的最大日均转乘乘客数量;
Umin表示一个站点的最小日均转乘乘客数量;
wn表示站点n的日均运送乘客数量;
Wmax表示一个站点的最大日均运送乘客数量;
Wmin表示一个站点的最小日均运送乘客数量;
使用最小—最大值方程对每个站点的四项指标进行归一化:
步骤4,计算每个站点的重要性程度;
根据四归一化指标的加权求和来求解各个站点的重要性程度:
实施例三
本发明实施例三提供一种分析计算北京市地铁交通网络中279个站点的重要性程度的方法,所用的交通网络拓扑和客运量数据从2016年1月1日至2016 年12月31日的数据。
步骤1,对收集得到的智能卡数据进行预处理;
这里列举了北京市10个站点的智能卡数据,数据来源于位于每个站点入口和出口的交通智能卡自动收费(AFC)设备,对数据进行预处理,显示这10 个站点的进站人数、每个站点的出站人数,如表1所示。
表1数据预处理
站点 | NDAPE1 | NDAPE2 | 站点 | NDAPE1 | NDAPE2 |
西直门 | 64474 | 64728 | 呼家楼 | 29577 | 29924 |
国贸 | 67319 | 70628 | 宋家庄 | 46540 | 45143 |
西单 | 65601 | 68677 | 西二旗 | 61166 | 67215 |
北京南站 | 95703 | 70292 | 北京西站 | 95779 | 67900 |
大望路 | 63397 | 64844 | 朝阳门 | 52349 | 58398 |
缩写:进入某站点的日均乘客数量(NDAPE1),离开某站点的日均乘客数量(NDAPE2)。
步骤2,得到每个站点的换乘乘客数量和运输乘客数量;
通过使用本发明提出城市轨道交通系统中乘客基于多智能体仿真过程计算每个站点的换乘乘客数量和运输乘客数量,计算结果如表2所示。
表2每个站点的换乘乘客数量和运输乘客数量
站点 | NDATP | NDAPP | 站点 | NDATP | NDAPP |
西直门 | 234153 | 845385 | 呼家楼 | 172362 | 1230322 |
国贸 | 130271 | 1235133 | 宋家庄 | 195147 | 466742 |
西单 | 87421 | 1165278 | 西二旗 | 131688 | 354186 |
北京南站 | 49640 | 628385 | 北京西站 | 43705 | 350461 |
大望路 | 98923 | 872879 | 朝阳门 | 83591 | 947488 |
缩写:某站点的日均转乘乘客数量(NDATP),某站点的日均运送乘客数量(NDAPP)。
步骤3,以最小—最大值方程对每个站点的四项指标进行归一化:
步骤4,计算每个站点的重要性程度;
其中加权系数c1=0.3,c2=0.2,c3=0.3,c4=0.2
计算结果如表3所示。
表3北京市10个站点的重要性程度和排名
站点 | 重要性程度(排名) | 站点 | 重要性程度(排名) |
西直门 | 0.7903(1) | 呼家楼 | 0.5809(6) |
国贸 | 0.7432(2) | 宋家庄 | 0.5760(7) |
西单 | 0.6669(3) | 西二旗 | 0.5749(8) |
北京南站 | 0.6303(4) | 北京西站 | 0.5723(9) |
大望路 | 0.6183(5) | 朝阳门 | 0.5606(10) |
计算结果表明,关键站点位于北京市城市中心枢纽(如西直门,国贸,西单)和连接市中心和郊区大型住宅区的咽喉区域(如大望路,宋家庄,西二旗),城市轨道交通管理人员应优先维护这些站点。
综上所述,本发明实施例所述的基于客流量的城市轨道交通网络关键车站识别方法,从实际客流量的角度出发,考虑每个站点的进入乘客数量、每个站点的离开乘客数量、每个站点的转乘乘客数量、每个站点的运送乘客数量这四项指标,使用基于多智能体的仿真和最小—最大归一化的分析方法得到每个站点的重要性程度,以此识别城市轨道交通网络中的关键站点,对城市轨道交通管理人员的维护工作提供理论支持。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种基于客流量的城市轨道交通网络关键站点识别方法,其特征在于,包括如下步骤:
步骤S110:收集交通智能卡自动收费AFC设备中的乘客数据并进行预处理,获取起始站点至最终站点OD对数据、进站乘客量和出站乘客量;
步骤S120:结合交通网络拓扑结构,根据所述OD对数据对乘客行为进行仿真,获取站点的换乘乘客量和运送乘客量;
步骤S130:根据所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量获取站点的重要性程度;
所述步骤S120具体包括:
根据所述OD对数据确定一乘客集合,根据所述乘客集合中的乘客状态仿真模拟所述乘客的路径信息,根据所述路径信息确定所述换乘乘客量和所述运送乘客量;
所述乘客集合中的乘客状态包括:
Sδ=0,表示乘客δ刚刚进入城市轨道交通系统;
Sδ=1,表示乘客δ刚刚从进站口进入一个站点,正走向站台准备等车;
Sδ=2,表示乘客δ刚到站台,开始等车;
Sδ=3,表示乘客δ刚乘上车;
Sδ=4,表示乘客δ刚在换乘站下车,开始换乘另一条城市轨道交通线路;
Sδ=5,表示乘客δ刚到目的站点下车,准备检票出站;
Sδ=6,表示乘客δ已经完成检票出站,并已退出城市轨道交通系统;
所述根据所述乘客集合中的乘客状态仿真模拟所述乘客的出行路径包括:
步骤1.1:初始化参数tb、te、△t和Td,设置时钟t=tb,第k次周期的时间区间为[tb+(k-1)×Td,tb+k×Td],k=1,2,3...K-1,其中,tb表示开始的时间,te表示结束的时间,△t表示步长,Td表示每个周期的时间长度,K=(te-tb)÷Td表示周期的个数;
步骤1.2:第k个周期中的乘客集合为{δ(O,D)k},对其中的任意一个乘客δ(O,D)k,设置初始状态Sδ=0,将{δ(O,D)k}添加至乘客总体其中,(O,D)表示一个从站点O到站点N的OD对,且所述第k个周期满足t=ts+(k-1)×Td,t<ts+k×Td,t=t+△t;
步骤1.3:当Sδ=0时,设置tδ=aδn,n站点表示起始站点O,将站点n增加至路径信息pδ中,判断tδ∈[t,t+△t]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=1,每重新计算一次,k增加1;其中,tδ表示第k个周期采样时的乘客的当前时间,aδn表示乘客进入站点n的时间;
当Sδ=1时,判断tδ∈[t,t+△t]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=2,否则,设置tδ=tδ+t(0,1)n,每重新计算一次,k增加1;其中,t(0,1)n表示乘客从站点n的进站口到站台所需的步行时间;
当Sδ=2时,根据列车时刻表得到即将驶入站点n的列车为i,判断tδ∈[ain,din]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=3,每重新计算一次,k增加1;其中,ain表示列车i到达站点n的时间,din表示列车i离开站点n的时间;
当Sδ=3时,识别列车i与乘客路径相同的最后一站m1,设置n=m1、tδ=am1i,将站点m1作为换乘站点添加至路径信息pδ中,判断tδ∈[ain,din]且n=D是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=5,否则,乘客δ(O,D)k进入状态Sδ=4,每重新计算一次,k增加1;其中,am1i表示列车i到达m1站点的时刻,D表示最终站点;
当Sδ=4时,寻找乘客δ(O,D)k的目标站点m2,设置n=m2、tδ=tδ+t(4,3)n,将站点m2作为终点站添加至路径信息pδ中,判断tδ∈[t,t+△t]是否成立,若成立,则乘客δ(O,D)k的状态重新设置为Sδ=2;其中,t(4,3)n表示乘客在站点m1的换乘时间;
当Sδ=5时,设置tδ=tδ+t(5,6)n,判断tδ∈[t,t+△t]是否成立,若成立,则乘客δ(O,D)k进入状态Sδ=6,每重新计算一次,k增加1;其中,t(5,6)n表示乘客从站点m2的站台到出站口所需的步行时间;
步骤1.4:根据步骤1.3的计算方法,确定集合{δ(O,D)k}中所有乘客的出行路径。
2.根据权利要求1所述的基于客流量的城市轨道交通网络关键站点识别方法,其特征在于,所述步骤S130具体包括:
使用最小—最大值方程对每个站点的所述进站乘客量、所述出站乘客量、所述换乘乘客量和所述运送乘客量进行归一化:
其中,xn表示进入站点n的日均乘客数量,Xmax表示进入一个站点的最大日均乘客数量;Xmin表示进入一个站点的最小日均乘客数量;yn表示离开站点n的日均乘客数量;Ymax表示离开一个站点的最大日均乘客数量;Ymin表示离开一个站点的最小日均乘客数量;un表示站点n日均转乘乘客数量;Umax表示一个站点的最大日均转乘乘客数量;Umin表示一个站点的最小日均转乘乘客数量;wn表示站点n的日均运送乘客数量;Wmax表示一个站点的最大日均运送乘客数量;Wmin表示一个站点的最小日均运送乘客数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810967392.7A CN108876064B (zh) | 2018-08-23 | 2018-08-23 | 基于客流量的城市轨道交通网络关键车站识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810967392.7A CN108876064B (zh) | 2018-08-23 | 2018-08-23 | 基于客流量的城市轨道交通网络关键车站识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108876064A CN108876064A (zh) | 2018-11-23 |
CN108876064B true CN108876064B (zh) | 2022-03-18 |
Family
ID=64321842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810967392.7A Active CN108876064B (zh) | 2018-08-23 | 2018-08-23 | 基于客流量的城市轨道交通网络关键车站识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108876064B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110188974A (zh) * | 2019-04-01 | 2019-08-30 | 浙江交通职业技术学院 | 一种地铁交通网络脆弱性评估方法 |
CN110111140A (zh) * | 2019-04-25 | 2019-08-09 | 华南理工大学 | 一种基于ic卡数据的热点换乘枢纽识别方法 |
CN111367900B (zh) * | 2020-02-26 | 2023-08-08 | 北京交通大学 | 基于afc数据的城市轨道交通网络常态限流强度的计算方法 |
CN116822762A (zh) * | 2023-03-29 | 2023-09-29 | 北京交通大学 | 地铁网络系统稳定性监测控制方法、系统和设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529815A (zh) * | 2016-11-15 | 2017-03-22 | 同济大学 | 城市轨道交通网络乘客出行时空轨迹估计方法及其应用 |
CN106779190A (zh) * | 2016-12-02 | 2017-05-31 | 东南大学 | 一种城市轨道交通乘客出行路径建议方法及系统 |
CN108198002A (zh) * | 2018-01-30 | 2018-06-22 | 武汉地铁运营有限公司 | 一种轨道交通票价测算、对比方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009049923A1 (de) * | 2009-10-19 | 2011-05-05 | Eads Deutschland Gmbh | Passagierbewegungsvorhersage- und Optimierungssystem |
-
2018
- 2018-08-23 CN CN201810967392.7A patent/CN108876064B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529815A (zh) * | 2016-11-15 | 2017-03-22 | 同济大学 | 城市轨道交通网络乘客出行时空轨迹估计方法及其应用 |
CN106779190A (zh) * | 2016-12-02 | 2017-05-31 | 东南大学 | 一种城市轨道交通乘客出行路径建议方法及系统 |
CN108198002A (zh) * | 2018-01-30 | 2018-06-22 | 武汉地铁运营有限公司 | 一种轨道交通票价测算、对比方法 |
Non-Patent Citations (1)
Title |
---|
基于网络客流传播的轨道交通关键站点识别;胡映月等;《西南交通大学学报》;20171022;第52卷(第6期);第1193-1200页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108876064A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108876064B (zh) | 基于客流量的城市轨道交通网络关键车站识别方法 | |
Chen et al. | Uncertainty in urban mobility: Predicting waiting times for shared bicycles and parking lots | |
CN111582750B (zh) | 一种轨道交通反向乘客识别及乘车班次确定方法及系统 | |
CN108009972B (zh) | 一种基于多源数据校核的多方式出行o-d需求估计方法 | |
CN107194491A (zh) | 一种基于公交客流和站间行程时间预测的动态调度方法 | |
CN109949098B (zh) | 一种高速公路收费站车道数的计算方法 | |
Wilson et al. | The potential impact of automated data collection systems on urban public transport planning. | |
CN108806302A (zh) | 一种车辆调度方法及装置 | |
CN113177657A (zh) | 轨道交通客流预测方法及装置 | |
Cheon et al. | Data-driven stochastic transit assignment modeling using an automatic fare collection system | |
CN114399099A (zh) | 一种城市轨道交通断面客流的预测方法及装置 | |
CN106846214A (zh) | 分析交通枢纽可达性对区域公共运输方式竞争影响的方法 | |
Zhang et al. | Short‐Term Passenger Flow Forecast of Rail Transit Station Based on MIC Feature Selection and ST‐LightGBM considering Transfer Passenger Flow | |
CN107704972A (zh) | 一种航线旅客需求预测方法 | |
Takeuchi et al. | Evaluation method of robustness for train schedules | |
Zheng et al. | Improved iterative prediction for multiple stop arrival time using a support vector machine | |
CN110657817A (zh) | 行程路线的推荐方法及装置 | |
CN109447352B (zh) | 一种公交出行od矩阵反推的组合方法 | |
Kim et al. | Predicting express train choice of metro passengers from smart card data | |
CN111367900A (zh) | 基于afc数据的城市轨道交通网络常态限流强度的计算方法 | |
CN116451867A (zh) | 一种基于时空出行路径标定的地铁短时客流预测方法 | |
Ivanchev et al. | Stochastic bus traffic modelling and validation using smart card fare collection data | |
CN110796301B (zh) | 一种基于ic卡数据的乘客流量预测方法及装置 | |
CN112182838A (zh) | 一种基于afc数据的短时客流预测动态模型 | |
Liu et al. | A Short‐Turn Dispatching Strategy to Improve the Reliability of Bus Operation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |