CN112232415A - 一种高铁全网络站点延误传播的识别方法、设备及介质 - Google Patents
一种高铁全网络站点延误传播的识别方法、设备及介质 Download PDFInfo
- Publication number
- CN112232415A CN112232415A CN202011109315.1A CN202011109315A CN112232415A CN 112232415 A CN112232415 A CN 112232415A CN 202011109315 A CN202011109315 A CN 202011109315A CN 112232415 A CN112232415 A CN 112232415A
- Authority
- CN
- China
- Prior art keywords
- delay
- speed rail
- network
- station
- propagation
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
- G06F18/24155—Bayesian classification
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/10—Operations, e.g. scheduling or time tables
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Mechanical Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Train Traffic Observation, Control, And Security (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
本发明公开了一种高铁全网络站点延误传播的识别方法、设备及介质,方法为:将运行时段划分为若干时间窗;根据高铁列车历史运行情况,计算高铁全网络每个站点每个时间窗内所有延误列车的平均延误时长;根据高铁全网列车时刻表中每两个站点之间的运行时长确定候选边;基于平均延误时长和候选边,使用NOTEARS算法构建初始的高铁全网络站点延误传播贝叶斯网络结构,删除其中权重小于预设权重阈值的边,得到最终的高铁全网络站点延误传播贝叶斯网络结构;当某个站点发生延误时,根据最终的高铁全网络站点延误传播贝叶斯网络结构,计算其余站点受到延误影响的强度。本发明可有效识别发生站点延误在高铁全网络中传播情况,具有较高的精度和广度。
Description
技术领域
本发明属于交通技术领域,具体涉及一种基于贝叶斯网络的高铁全网络站点延误传播的识别方法、设备及介质。
背景技术
高速列车作为一种中长距离的出行方式,具有速度快、安全性高、正点率高、能耗小、运输能力大等优点。但是高速列车在实际运营的过程中,受到多种因素(极端天气、设备故障、交通事故、运营管理不当等)的影响,经常发生列车到站延误事件。高速列车发生延误不仅会打乱列车运行计划,影响高速列车到站的运营管理;而且会使乘客的出行受到影响,降低乘客的乘车舒适度。研究理解高速列车延误传播机理,识别高速列车在站点间的延误传播,对管理者制定列车运行计划,提高车站的运营管理水平具有十分重要的现实意义。但是,现有识别高速列车在站点间的延误传播的方法存在以下问题:
(1)通过专家的领域知识手动建立高速列车在站点间的延误传播关系不具备客观性,且此方法只能识别两个站点之间是否具有传播关系,难以获得站点间延误传播的强弱关系。
(2)现在只具有识别某条高铁单线的延误传播关系的方法,但是识别高铁单线站点之间的延误传播关系的方法并不能识别延误在整个网络中的传播情况,因为整个高铁网络的延误传播较某条高铁单线的延误传播更为复杂。
综上所述,现在识别高铁全网络站点延误传播的方法仍有不足,因为高铁网络规模巨大,结构复杂,如果采用上述方法,很难有效识别高铁全网络的站点延误传播。为了更好的理解高速列车在延误传播机理,制定列车到站计划,提高车站的运营管理水平,迫切需要一种基于贝叶斯网络的高铁全网络站点延误传播的识别方法。
发明内容
本发明解决的技术问题是,针对现有技术的不足,提供了一种基于贝叶斯网络的高铁全网络站点延误传播的识别方法、设备及介质,通过采取数据驱动的方式,有效地识别所发生的站点延误在高铁全网络中传播情况,具有较高的精度和广度,适用于列车延误传播和预测研究,有利于后续列车到站计划的制定与调整。
为了实现上述技术目的,本发明提供的技术方案为:
一种基于贝叶斯网络的高铁全网络站点延误传播的识别方法,包括:
步骤1,将高铁历史每天的运行时段平均划分为若干时间窗;
步骤2,根据高铁列车的历史运行情况,计算高铁全网络的每个站点在每个时间窗内所有延误列车的平均延误时长;
步骤3,根据高铁全网列车时刻表中每两个站点之间的运行时长,确定贝叶斯网络的候选边;
步骤4,基于步骤2得到的平均延误时长和步骤3得到的候选边,使用NOTEARS算法构建初始的高铁全网络站点延误传播贝叶斯网络结构;
步骤5,对初始的高铁全网络站点延误传播贝叶斯网络结构,删除其中边权重小于预设权重阈值的边,得到最终的高铁全网络站点延误传播贝叶斯网络结构;
步骤6,当高铁全网络的某个站点发生延误时,根据最终的高铁全网络站点延误传播贝叶斯网络结构,判断高铁全网络中受到该延误站点影响的站点以及影响强度。。
进一步的,所述步骤2的具体步骤为:
步骤2.1,将每个延误列车分配至对应的站点和时间窗内;
步骤2.2,通过以下公式计算每个站点分别在每个时间窗内所有延误列车的平均延误时长:
进一步的,所述步骤3的具体步骤为:
步骤3.1,在高铁全网列车时刻表中,查找任意两个站点之间是否有列车运行:如果有,按照步骤3.2进行判断;如果没有,则该两个站点相连接的边不能成为候选边;
步骤3.2,选择该两个站点之间列车运行时间最短的时长:如果该时长小于时间窗的时长,则该两个站点相连接的边可以成为候选边;否则,该两个站点相连接的边不能成为候选边。
进一步的,所述步骤4的具体步骤为:
步骤4.1,初始化贝叶斯网络的网络结构:将每一个站点作为贝叶斯网络的节点,以步骤3得到的候选边作为贝叶斯网络的候选边;
步骤4.2,将步骤2得到的平均延误时长,整理为m行n列的站点平均延误时长矩阵,其中每行对应一个时间窗,每列对应一个站点;
步骤4.3,根据站点平均延误时长矩阵和初始化得到的贝叶斯网络的网络结构,采用NOTEARS算法进行贝叶斯网络结构学习,得到初始的高铁全网络站点延误传播贝叶斯网络结构。
进一步的,步骤5中的预设权重阈值,设置为初始的高铁全网络站点延误传播贝叶斯网络结构中渗流相变点对应的权重值。
进一步的,每个时间窗的时长为1h。
本发明还提供一种电子设备,包括存储器及处理器,所述存储器中存储有计算机程序,其特征在于,所述计算机程序被所述处理器执行时,使得所述处理器实现上述任一项所述的方法。
本发明还提供一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现上述任一项所述的方法。
有益效果
本发明通过构建高铁全网络站点延误传播贝叶斯网络,当某一站点的高速列车发生延误时,即可得知因此站点延误导致其他发生延误的站点,从而提前通知相关高铁站点的乘客和工作人员做好准备,通过车站站内即将到站的列车进行调度调整,可以对未发生延误的列车降低列车发生延误的概率,对已经发生延误的列车降低列车延误时长,同时可以为高速列车运行图的制定与调整提供参考。根据高铁全网络站点延误传播贝叶斯网络,可以判断某个站点的高速列车延误之后,通过网络中边的权重判断网络中其他站点受到延误影响的强度,进而关注网络中影响其他站点的关键车站以及受到较大延误影响的站点。这对于整个高速列车网络的站点拓扑关系的改进有着重大的指引意义。
本发明和现有技术相比,有以下优点:
(1)传统通过专家的领域知识手动建立高速列车在站点间的延误传播关系不具备客观性,且传统方法只能识别两个站点之间是否具有传播关系,难以获得站点间延误传播的强弱关系。本发明方法通过数据驱动的方式,更能精准反应网络中延误传播的实际情况,且可以识别受到延误传播重大影响的车站。
(2)本方法通过NOTEARS算法构建得到高铁全网络站点延误传播贝叶斯网络,解决了以往只能识别高铁单线站点之间的延误传播关系而不能识别大型综合网络中站点之间的延误传播关系的问题,为研究全网络的延误传播以及延误预测提供指导意义。
附图说明
图1为本发明实施例所述方法的流程图;
图2为本发明实施例所运用的高铁站点地理位置示意图;
图3为本发明实施例采用渗流相变点预设权重阈值时,最大团簇和第二大团簇大小随权重阈值的变化图;
图4为本发明实施例最终的高铁全网络站点延误传播贝叶斯网络。
具体实施方式
本实施例以本发明的技术方案为依据开展,给出了详细的实施方式和具体的操作过程,对本发明的技术方案作进一步解释说明。
实施例1:
本实施例提供一种基于贝叶斯网络的高铁全网络站点延误传播的识别方法,如图1所示,包括:
步骤0,获取高速列车运行基本信息和高速列车站点基本信息;
本实施例采集的高速列车运行基本信息包括列车时刻表信息和列车运行实绩数据信息,采集的时间段为2018年1月20日至3月20日。高速列车站点基本信息包括高速列车站点名信息和高速列车站点经纬度信息。其中,高铁列车运行时刻表信息在中国铁路总公司12306官方网站上获得,列车时刻表中记录了4079辆开行的高速列车;列车运行实绩数据信息来源于我们使用python开发的爬虫工具对中国铁路总公司12306官方网站公布的列车实时到站情况,包含了高速列车车次、到达站、图定到站时间、实际到站时间、延误时长、日期等信息;高速列车站点基本信息数据来源于腾讯地图,其中包括站点名和站点经纬度信息,除去有数据缺失的高铁站点,一共有562个高铁站点,如图2所示。
步骤1,将高铁一天的运行时段平均划分为若干时间窗;
本实施例中,基于中国高速列车一天主要运行时段6:00-24:00,将18个小时分割为18个时间间隔为1h的时间窗。本实施例的历史运行时段包括2018年1月20日至3月20日共60天,共有1080个时间窗。
步骤2,根据高铁列车的历史运行情况,计算高铁全网络的每个站点分别在每个时间窗内所有延误列车的平均延误时长;具体步骤为:
步骤2.1,将每个延误列车分配至对应的站点和时间窗内;
步骤2.2,通过以下公式计算每个站点分别在每个时间窗内所有延误列车的平均延误时长:
步骤3,根据高铁全网列车时刻表中每两个站点之间的运行时长,确定高铁全网络站点延误传播贝叶斯网络的候选边;具体步骤为:
步骤3.1,在高铁全网列车时刻表中,查找任意两个站点之间是否有列车运行:如果有,按照步骤3.2进行判断;如果没有,则该两个站点相连接的边不能成为候选边;
步骤3.2,选择该两个站点之间列车运行时间最短的时长:如果该时长小于时间窗的时长,则该两个站点相连接的边可以成为候选边;否则,该两个站点相连接的边不能成为候选边。
步骤4,基于步骤2得到的平均延误时长和步骤3得到的候选边,使用NOTEARS算法构建初始的高铁全网络站点延误传播贝叶斯网络结构;具体步骤为:
步骤4.1,初始化贝叶斯网络的网络结构:将每一个站点作为贝叶斯网络的节点,以步骤3得到的候选边作为贝叶斯网络的候选边;
步骤4.2,将步骤2得到的平均延误时长,整理为m行n列的站点平均延误时长矩阵,其中每行对应一个时间窗,每列对应一个站点;
步骤4.3,根据站点平均延误时长矩阵和初始化得到的贝叶斯网络的网络结构,采用NOTEARS算法进行贝叶斯网络结构学习,得到初始的高铁全网络站点延误传播贝叶斯网络结构。本实施例中,初始高铁全网络站点延误传播贝叶斯网络中共有548个节点,1168条边。
步骤5,对初始的高铁全网络站点延误传播贝叶斯网络结构,删除其中边权重小于预设权重阈值的边,得到最终的高铁全网络站点延误传播贝叶斯网络结构;
其中的预设权重阈值,设置为初始的高铁全网络站点延误传播贝叶斯网络结构中渗流相变点对应的权重值(相变点,是指初始的高铁全网络站点延误传播贝叶斯网络发生相变的临界点,代表着网络前后不同状态的临界点,且相变点位置的确定是通过绘制图3所示的最大团簇和第二大团簇随权重的变化规律而得到的,相变点对应的权重值为第二大团簇有最大的大小时(纵坐标的值最大)所对应的权重)。本实施例最后得到具有472个节点,463条边的网络结构,如图4所示。初始的高铁全网络站点延误传播贝叶斯网络结构中,每条边都有对应的权重大小,权重越大的边对网络来说就越重要。如果不断调整权重阈值,在网络中留下大于权重阈值的边,网络中存在的边会越来越少,同时网络也因此而瓦解,瓜分为了多个团簇。本实施例中渗流概念中相变点的确定,是根据图3中第二大团簇中团簇最大时进行确定。
步骤6,当高铁全网络的某个站点发生延误时,根据最终的高铁全网络站点延误传播贝叶斯网络结构,判断高铁全网络中受到该延误站点影响的站点以及影响强度。
实施例2:
本实施例提供一种电子设备,包括存储器及处理器,所述存储器中存储有计算机程序,其特征在于,所述计算机程序被所述处理器执行时,使得所述处理器实现实施例1所述的方法。
实施例3:
本实施例提供一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现实施例1所述的方法。
以上实施例为本申请的优选实施例,本领域的普通技术人员还可以在此基础上进行各种变换或改进,在不脱离本申请总的构思的前提下,这些变换或改进都应当属于本申请要求保护的范围之内。
Claims (8)
1.一种基于贝叶斯网络的高铁全网络站点延误传播的识别方法,其特征在于,包括:
步骤1,将高铁历史每天的运行时段平均划分为若干时间窗;
步骤2,根据高铁列车的历史运行情况,计算高铁全网络的每个站点在每个时间窗内所有延误列车的平均延误时长;
步骤3,根据高铁全网列车时刻表中每两个站点之间的运行时长,确定贝叶斯网络的候选边;
步骤4,基于步骤2得到的平均延误时长和步骤3得到的候选边,使用NOTEARS算法构建初始的高铁全网络站点延误传播贝叶斯网络结构;
步骤5,对初始的高铁全网络站点延误传播贝叶斯网络结构,删除其中边权重小于预设权重阈值的边,得到最终的高铁全网络站点延误传播贝叶斯网络结构;
步骤6,当高铁全网络的某个站点发生延误时,根据最终的高铁全网络站点延误传播贝叶斯网络结构,判断高铁全网络中受到该延误站点影响的站点以及影响强度。
3.根据权利要求1所述的方法,其特征在于,所述步骤3的具体步骤为:
步骤3.1,在高铁全网列车时刻表中,查找任意两个站点之间是否有列车运行:如果有,按照步骤3.2进行判断;如果没有,则该两个站点相连接的边不能成为候选边;
步骤3.2,选择该两个站点之间列车运行时间最短的时长:如果该时长小于时间窗的时长,则该两个站点相连接的边可以成为候选边;否则,该两个站点相连接的边不能成为候选边。
4.根据权利要求1所述的方法,其特征在于,所述步骤4的具体步骤为:
步骤4.1,初始化贝叶斯网络的网络结构:将每一个站点作为贝叶斯网络的节点,以步骤3得到的候选边作为贝叶斯网络的候选边;
步骤4.2,将步骤2得到的平均延误时长,整理为m行n列的站点平均延误时长矩阵,其中每行对应一个时间窗,每列对应一个站点;
步骤4.3,根据站点平均延误时长矩阵和初始化得到的贝叶斯网络的网络结构,采用NOTEARS算法进行贝叶斯网络结构学习,得到初始的高铁全网络站点延误传播贝叶斯网络结构。
5.根据权利要求1所述的方法,其特征在于,步骤5中的预设权重阈值,设置为初始的高铁全网络站点延误传播贝叶斯网络结构中渗流相变点对应的权重值。
6.根据权利要求1所述的方法,其特征在于,每个时间窗的时长为1h。
7.一种电子设备,包括存储器及处理器,所述存储器中存储有计算机程序,其特征在于,所述计算机程序被所述处理器执行时,使得所述处理器实现如权利要求1~6中任一项所述的方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1~6中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011109315.1A CN112232415B (zh) | 2020-10-16 | 2020-10-16 | 一种高铁全网络站点延误传播的识别方法、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011109315.1A CN112232415B (zh) | 2020-10-16 | 2020-10-16 | 一种高铁全网络站点延误传播的识别方法、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112232415A true CN112232415A (zh) | 2021-01-15 |
CN112232415B CN112232415B (zh) | 2021-07-20 |
Family
ID=74118447
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011109315.1A Active CN112232415B (zh) | 2020-10-16 | 2020-10-16 | 一种高铁全网络站点延误传播的识别方法、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112232415B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927872A (zh) * | 2014-04-28 | 2014-07-16 | 武汉大学 | 一种基于浮动车数据预测多时段行程时间分布的方法 |
CN108764526A (zh) * | 2018-04-20 | 2018-11-06 | 西南交通大学 | 一种基于列车运行秩序的列车晚点人工神经网络识别方法 |
CN109948204A (zh) * | 2019-03-05 | 2019-06-28 | 清华大学 | 一种基于贝叶斯网的高铁接触网动态风险分析方法 |
CN110782070A (zh) * | 2019-09-25 | 2020-02-11 | 北京市交通信息中心 | 一种城市轨道交通突发事件客流时空分布预测方法 |
CN110969295A (zh) * | 2019-11-25 | 2020-04-07 | 北京交通大学 | 列车区段延误预测误差控制方法 |
CN111091225A (zh) * | 2019-11-08 | 2020-05-01 | 北京交通大学 | 一种基于深度学习的城市轨道交通短期客流预测方法 |
-
2020
- 2020-10-16 CN CN202011109315.1A patent/CN112232415B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927872A (zh) * | 2014-04-28 | 2014-07-16 | 武汉大学 | 一种基于浮动车数据预测多时段行程时间分布的方法 |
CN108764526A (zh) * | 2018-04-20 | 2018-11-06 | 西南交通大学 | 一种基于列车运行秩序的列车晚点人工神经网络识别方法 |
CN109948204A (zh) * | 2019-03-05 | 2019-06-28 | 清华大学 | 一种基于贝叶斯网的高铁接触网动态风险分析方法 |
CN110782070A (zh) * | 2019-09-25 | 2020-02-11 | 北京市交通信息中心 | 一种城市轨道交通突发事件客流时空分布预测方法 |
CN111091225A (zh) * | 2019-11-08 | 2020-05-01 | 北京交通大学 | 一种基于深度学习的城市轨道交通短期客流预测方法 |
CN110969295A (zh) * | 2019-11-25 | 2020-04-07 | 北京交通大学 | 列车区段延误预测误差控制方法 |
Non-Patent Citations (2)
Title |
---|
XUNZHENG,BRYON ARAGAM,PRADEEP RAVIKUMAR,ERIC P.XING: "《DAGs with NO TEARS: Continuous Optimization for Structure Learning》", 《ARXIV》 * |
孙宝丹; 王培超; 周鋆: "《NO TEARS算法在XSS攻击检测中的应用研究》", 《小型微型计算机系统》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112232415B (zh) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108109382B (zh) | 基于复合网络的拥堵点、拥堵线、拥堵区域的发现方法 | |
CN102087788B (zh) | 基于浮动车车速置信度的交通状态参数估计方法 | |
US20210407284A1 (en) | Vehicle traffic flow prediction method with missing data | |
CN108596202B (zh) | 基于移动终端gps定位数据计算个人通勤时间的方法 | |
CN109558988B (zh) | 一种基于大数据融合的电动汽车能耗预测方法及系统 | |
CN111008223B (zh) | 一种基于时空关联规则的区域交通拥堵相关性计算方法 | |
CN103984994B (zh) | 一种城市轨道交通客流高峰持续时间预测方法 | |
CN104408203A (zh) | 一种预测移动对象轨迹终点的方法 | |
CN104123833A (zh) | 一种道路状况的规划方法和装置 | |
CN107195177B (zh) | 基于分布式内存计算框架对城市交通路况的预测方法 | |
CN105740643A (zh) | 一种基于城市区域网格自适应的pm2.5浓度推测方法 | |
CN112232722B (zh) | 一种基于建筑信息模型技术的建筑设备管理系统 | |
CN112215416B (zh) | 智能规划巡检航线系统及方法 | |
CN111739299A (zh) | 稀疏轨迹的车辆排队长度确定方法、装置、设备及介质 | |
CN109410576A (zh) | 多源数据融合的路况分析方法、装置、存储介质及系统 | |
CN114021883A (zh) | 一种高峰时期地铁接驳共享单车的调度方法 | |
CN112036757A (zh) | 基于手机信令和浮动车数据的停车换乘停车场的选址方法 | |
CN111199247A (zh) | 一种公交运行仿真方法 | |
CN115496316A (zh) | 一种基于多源数据的公交线路客流动态分配方法 | |
CN105243841A (zh) | 基于网页地图的实时路况采集与预测方法 | |
CN112232415B (zh) | 一种高铁全网络站点延误传播的识别方法、设备及介质 | |
CN106600030A (zh) | 一种基于路网行程时间可靠性的出行预留时间计算方法 | |
CN117273538A (zh) | 一种交通碳排放空间网络关联结构分析方法及系统 | |
Tan et al. | Statistical analysis and prediction of regional bus passenger flows | |
CN112630836B (zh) | 基于电网微地形覆冰威胁路径延伸分析的监测布点方法及系统 |
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 |