CN113744526B - 一种基于lstm和bf的高速公路风险预测方法 - Google Patents
一种基于lstm和bf的高速公路风险预测方法 Download PDFInfo
- Publication number
- CN113744526B CN113744526B CN202110979482.XA CN202110979482A CN113744526B CN 113744526 B CN113744526 B CN 113744526B CN 202110979482 A CN202110979482 A CN 202110979482A CN 113744526 B CN113744526 B CN 113744526B
- Authority
- CN
- China
- Prior art keywords
- time
- risk
- lstm
- accident
- traffic
- 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
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
- G08G1/0129—Traffic data processing for creating historical data or processing based on historical data
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/017—Detecting movement of traffic to be counted or controlled identifying vehicles
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Probability & Statistics with Applications (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供了一种基于LSTM和BF的高速公路风险预测方法,包括线下风险预警预测模型训练和线上风险预警模型实时预测;线下风险预警预测模型训练包括构建区分事故与安全状态的LSTM瞬时风险判别模型,以及建立基于递减系数、先验概率和阈值的BF序列风险预测模型;线上风险预警模型实时预测时,以规定周期内各观测点上LSTM瞬时风险判别模型的输出序列为输入向量,以模式内的观测频数占比为先验概率,输入BF序列风险预测模型,最终得到未来路段风险状态预测结果。本发明提高了高速公路实时风险预测的精度。
Description
技术领域
本发明涉及交通安全评价和智能交通技术领域,特别是一种基于LSTM(LongShort-Term Memory,长短期记忆神经网络)和BF(Bayesian Filtering,贝叶斯滤波)的高速公路风险预测方法。
背景技术
汽车保有量前六的国家中,我国交通事故致死率远超美国、德国、日本等,且重大交通安全事故伤亡人数与财物损失一直维持在高位,交通安全已成为影响我国社会经济发展的突出问题。研究高速公路交通事故风险状态预测方法,为高速公路交通安全早期预警提供可靠依据,有助于消除安全隐患,预防和减少高速公路交通事故发生。
高速公路交通事故风险预测一般基于海量交通流数据,首先对交通事故风险特征进行挖掘和提取,然后根据提取特征建立预测模型,实现对某一路段/区域在未来一段时间内交通事故风险的预测。传统交通事故风险预测通常采用参数模型方法,如多项逻辑回归模型和贝叶斯网络模型等。近年来,也有一些研究引入了非参数的机器学习模型,如K近邻法、支持向量机、随机森林法等。随着深度学习技术的不断发展,许多研究者也将深度网络模型应用在了风险评估领域,如递归神经网络、图卷积网络等。然而目前大部分模型通常需要详细的交通流量、占有率和车速等参数信息,在实际应用方面,特别是对于普遍仅拥有车辆车牌识别数据类型的高速公路卡口***,仍存在一定局限性。
发明内容
针对现有技术中存在不足,本发明提供了一种基于LSTM和BF的高速公路风险预测方法,通过LSTM瞬时风险判别模型学习历史风险数据中存在的时间依赖关系,并通过BF序列风险预测模型融合LSTM瞬时风险判别模型预测结果提高实时风险预测效果,解决通过高速公路卡口数据高效预测事故风险状态的问题。
本发明是通过以下技术手段实现上述技术目的的。
一种基于LSTM和BF的高速公路风险预测方法,包括步骤:
S1,线下风险预警预测模型训练
基于历年高速公路历史交通事故数据,得到K类事故交通时空模式库,统计得到各事故交通时空模式的发生频数;基于事故发生前指定时间内上下游卡口间的时空范围,通过卡口车牌识别历史数据提取多步高频时变变量,构建区分事故与安全状态的LSTM瞬时风险判别模型;以规定周期内各观测点上LSTM瞬时风险判别模型判别结果的序列为输入向量,利用低频时变变量和常变量匹配的交通时空模式,建立基于先验概率的BF序列风险预测模型;
S2,线上风险预警模型实时预测
实时获取预测时刻的低频时变变量和预测路段常变量,实时匹配预测时空的交通时空模式;实时采集每一辆经过卡口车辆的牌照、车辆类型、经过时刻和所在车道,提取时间窗内上下游卡口间的多步高频时变变量,输入LSTM瞬时风险判别模型,得到t时刻的LSTM瞬时风险判别模型判别结果;以规定周期内各观测点上LSTM瞬时风险判别模型的输出序列为输入向量,以匹配交通时空模式内的观测频数占比为先验概率,输入BF序列风险预测模型,最终得到未来路段风险状态预测结果。
进一步地,所述LSTM瞬时风险判别模型的构建过程,具体为:
1)基于历史交通事故,以预测时刻t0前tw内上下游卡口M1M2之间的时空范围为输入时间窗,提取T步输入高频时变特征向量X={x1,x2,…,xt,…,xT},输出从预测时刻起未来tp时长内有事故发生,得到事故样本;
2)对每一个事故样本同一时段、同一路段的无事故时空范围进行随机下采样,以预测时刻t′0前tw内上下游卡口M1M2之间的时空范围为输入时间窗,提取T步输入高频时变特征向量X′={x′1,x′2,…,x′t,…,x′T},输出从预测时刻起未来tp时长内无事故发生,得到非事故样本;
3)融合事故样本和非事故样本,训练LSTM瞬时风险判别模型,使得损失函数最小,得到最终标定的LSTM瞬时风险判别模型。
更进一步地,所述预测时刻t0是事故发生时刻tc前tp时间,tp为预测时长。
进一步地,所述LSTM瞬时风险判别模型的输出序列值视为二分类随机变量y∈{0,1},随机变量y∈{0,1}的取值由参数θ决定,参数θ的期望值表示为:
其中:y为LSTM瞬时风险判别模型的输出序列,m和q分别代表LSTM输出序列事故类和非事故类变量的个数,a和b是beta分布的超参数,分别代表事故类和非事故类的先验概率。
更进一步地,所述a和b的值随每个样本所处的交通时空模式变化,参数a的值为该样本所处交通时空模式TMc内事故记录的观测频数fc占所有模式总频数的占比,即相应参数其中fk为各交通时空模式内事故记录的观测频数。
更进一步地,所述每个样本所处的交通时空模式是通过以下方法确定的:
1)基于高速公路历史交通事故数据,以事故样本的低频时变变量和路段常变量为特征变量[z1,z2,...,zs,...,zS],通过K-modes聚类方法得到K类事故模式,每一类事故模式的聚类中心为TMk=[z1k,z2k,...,zsk,...,zSk],k=1,2,...,K;
2)以待确定样本的低频时变变量和路段常变量为该样本的交通时空模式特征变量 TMt=[z1t,z2t,...,zst,...,zSt],在TMk=[z1k,z2k,...,zsk,...,zSk]中基于汉明距离进行模式匹配,与 TMt汉明距离最小的模式TMc即为该样本所处的交通时空模式。
更进一步地,所述未来路段风险状态预测是根据所述E(θ|y)与阈值的大小进行判定的,当E(θ|y)<τ1:预测路段处于低风险状态,当τ1≤E(θ|y)<τ2:预测路段处于中风险状态,当τ2≤E(θ|y):预测路段处于高风险状态,其中τ1、τ2分别代表中、高风险状态的阈值。
进一步地,所述LSTM瞬时风险判别模型的隐藏层和Dense全连接输出层中间设有随机失活Dropout,且LSTM瞬时风险判别模型配合使用输入连接权重L2正则化和学习率衰减。
进一步地,所述高频时变变量包括上下游卡口流量、车道间流量差、车道平均流量、路段流量密度和大型/小型车流量比。
进一步地,所述低频时变变量包括季节、周天数、日时段和天气类型。
进一步地,所述常变量包括路段车道数、道路线形、路段进出口匝道间距和路段进出口匝道个数,并将所述路段进出口匝道间距连续变量按取值区间离散化。
本发明的有益效果为:
(1)本发明无需传感器采集详细的交通流量、占有率和车速等参数信息,适用于普遍仅拥有车辆车牌识别数据类型的高速公路卡口***,具有较高的实用性。
(2)本发明通过构建LSTM-BF预测框架,其中LSTM瞬时风险判别模型学习历史风险数据中存在的时间依赖关系,BF序列风险预测模型融合规定周期内LSTM瞬时风险判别模型的预测结果,提高了高速公路实时风险预测精度。
(3)本发明区分高频时变变量、低频时变变量和道路常变量共三种类型变量,通过高频时变变量构造LSTM瞬时风险判别模型的多步特征输入,通过低频时变变量和道路常变量匹配交通时空模式得到BF序列风险预测模型风险先验概率,充分考虑了不同类型变量对风险预测的差异化影响,提高了高速公路实时风险预测模型的准确性。
附图说明
图1为本发明的基于LSTM和BF的高速公路风险预测方法流程图;
图2为本发明的LSTM模型输入时空范围的概念示意图;
图3为本发明的LSTM模型多步特征输入(步数T=4)的概念示意图;
图4为本发明的BF模型输入序列的概念示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图对本发明的具体实施方式做详细说明。
如图1所示,一种基于LSTM和BF的高速公路风险预警方法,具体包括如下步骤:
步骤一,线下风险预警预测模型训练
根据交通事故信息库,获得历年高速公路历史交通事故数据,基于历年高速公路历史交通事故数据,提取事故样本的低频时变变量和路段常变量,通过聚类分析方法得到K类事故模式库,并统计得到各事故模式的发生频数。利用卡口实时信息采集***采集的车牌识别数据,基于事故发生前指定时间内上下游卡口间的时空范围,采用移动时间窗形式提取多步高频时变特征变量,构建区分事故与安全状态的LSTM瞬时风险判别模型(即深度学习神经网络)。以规定周期内各观测点上LSTM瞬时风险判别模型判别结果的序列为输入向量,建立基于先验概率的BF序列风险预测模型。
步骤一具体为:
步骤(1),根据高速公路沿线卡口位置,将高速公路主线划分为L个路段,根据交通流方向确定路段上下游卡口关联关系,获取路段常变量,包括路段车道数、道路线形、路段进出口匝道间距、路段进出口匝道个数,并将其中路段进出口匝道间距连续变量按取值区间离散化。
步骤(2),基于高速公路历史交通事故数据,获取步骤(1)划分的各路段历史交通事故记录数据,提取事故发生路段、发生时间(即低频时变变量,具体包括所在季节、周天数(一周内星期几)、日时段(一天内小时数))、天气类型(属于低频时变变量)和步骤(1)所述路段常变量为特征变量[z1,z2,...,zs,...,zS]。基于所述特征变量,通过K-modes聚类方法得到K 类事故模式,每一类事故模式的聚类中心为TMk=[z1k,z2k,...,zsk,...,zSk],k=1,2,...,K,并统计得到各事故模式内事故记录的观测频数f1,f2,…,fK。
步骤(3),利用卡口实时信息采集***采集的车牌识别数据,基于事故发生前指定时间内上下游卡口间的时空范围,采用移动时间窗形式提取多步高频时变特征变量,构建区分事故与安全状态的LSTM瞬时风险判别模型。所述步骤(3)具体通过以下子步骤来实现:
步骤3-1),基于卡口实时信息采集***历史采集的经过卡口的车辆牌照、车辆经过时刻、车辆所在车道和车辆类型(小型车/大型车)信息,统计单位时间内通过各卡口的总流量、分车道流量和分车型流量。
步骤3-2),采用移动时间窗形式构建LSTM瞬时风险判别模型的事故样本输入和输出,所述步骤3-2)具体通过以下子步骤来实现:
步骤3-2-1),记高速公路历史交通事故数据为D个,根据交通流方向确定路段上下游卡口关联关系确定事故发生路段的上游卡口M1、下游卡口M2和发生时刻tc。
步骤3-2-2),如图2所示,以事故发生时刻tc前tp时间为预测时刻t0,即t0=tc-tp,tp为预测时长。以t0时刻前tw(单位为分钟min)内上下游卡口M1M2之间的时空范围为输入时间窗,以5min为单位步长,提取T=tw/5步输入特征向量X={x1,x2,…,xt,…,xT},每步输入xt包括的高频时变变量有:上下游卡口流量(总流量)、车道间流量差(分车道流量的差值)、车道平均流量(分车道流量的均值)、路段流量密度(总流量除以上下游卡口之间的间距)、大型/小型车流量比(分车型流量的比值),其中各时变变量均采用z-score方法进行标准化处理;输入时间窗的移动步长为1min;输出标签记为1,即从预测时刻起未来tp时长内有事故发生,如图3所示。
步骤3-2-3),对所有D个记录事故按步骤3-2-2)所述方法构建LSTM瞬时风险判别模型正样本的输入和输出,得到共D个LSTM瞬时风险判别模型正样本。按照下采样比例为正样本(事故样本)数:负样本(非事故样本)数=1:3,对每一个事故样本同一时段、同一路段的无事故时空范围进行随机下采样,并按步骤3-2-2)所述方法构建LSTM瞬时风险判别模型负样本的输入,输出标签记为0,即从预测时刻起未来tp时长内无事故发生,得到共3D个LSTM瞬时风险判别模型负样本。
步骤3-2-4),为了避免模型深度过高致使数据过拟合,设LSTM瞬时风险判别模型的隐藏层层数为1,隐藏层节点个数由试验获取。同时,为了防止模型过拟合,在所述隐藏层和 Dense全连接输出层中间加入随机失活Dropout,并配合使用输入连接权重L2正则化和学习率衰减技术来改善过拟合情况。
步骤3-2-5),融合步骤3-2-3)中的正、负样本为总样本(共D+3D=4D个),按训练集、验证集和测试集分别占60%、20%和20%随机抽取样本构成相应集合,以二分类交叉熵损失函数为损失函数,采用Adam优化器,训练步骤3-2-4)所述LSTM瞬时风险判别模型使得损失函数最小,得到最终标定的LSTM瞬时风险判别模型。
步骤(4),以规定周期内LSTM瞬时风险判别模型判别结果的序列为输入向量,建立基于先验概率的BF序列风险预测模型。所述步骤(4)具体通过以下子步骤来实现:
步骤4-1),基于步骤(3)获取的4D个总样本,对于每个样本的预测时刻t0,按时间先后顺序将t0前间隔为1min的N-1个观测点和t0本身作为该样本的观测周期(即t0-N+1,t0-N+2,…,t0共N个观测点)。
步骤4-2),如图4所示,以步骤4-1)得到的某样本观测周期内的每一个观测点为预测时刻t0,按步骤3-2-2)所述方法得到对应N个观测点的N组LSTM瞬时风险判别模型输入特征,分别输入步骤3-2-5)得到的LSTM瞬时风险判别模型,得到观测周期内N个观测点上LSTM瞬时风险判别模型判别结果构成的序列y=[y1,...,yN],并以该序列为BF模型的输入向量。BF输出标签与该样本的正负标签一致,若原样本为正样本,则BF输出标签标记为 1,否则标记为0。
步骤4-3),将LSTM瞬时风险判别模型的输出视为二分类随机变量y∈{0,1}的观测样本,建立基于递减系数、先验概率和阈值的BF序列风险预测模型,所述步骤4-3)具体通过以下子步骤来实现:
步骤4-3-1),随机变量y∈{0,1}的取值由参数θ决定:
p(y=1|θ)=θ
其中参数θ代表随机变量y=1(事故类)的概率。采用beta分布作为参数θ的先验概率分布形式:
其中Γ(x)是gamma函数,a和b是beta分布的超参数,分别代表两种类别的先验知识,可由两种类别的初始观测频率表示。
步骤4-3-2),当给定LSTM输出序列y=[y1,...,yN],根据贝叶斯概率公式,参数θ的后验概率p(θ|y)为先验概率beta(θ|a,b)与二项分布概率bin(m|N,θ)的乘积并进行标准化,即:
其中m和q分别代表LSTM输出序列中y=1(事故类)和y=0(非事故类)变量的个数,并满足m+q=N。则参数θ的期望值可简单表示为:
步骤4-3-3),考虑到越靠近当前观测时刻t的LSTM瞬时风险判别模型输出应予以更大的权重,加入随观测时间变化的递减函数rn,以提高最终预测结果的准确性:
其中n=1,2,...,N代表LSTM序列的第n个输出值的序号(N为序列最后一个输出值的序号);C∈(0,1]为递减常数,其大小将显著影响BF的预测效果,当C取1时,表示LSTM序列中每个输出结果的权重相同;mn和qn代表LSTM序列第n个输出结果的相应标签值,并满足mn+qn=1,如若输出为事故类(yn=1),则标记mn=1,qn=0;反之非事故类(yn=0) 为mn=0,qn=1。
步骤4-3-4),通过设定阈值τ为BF序列风险预测模型判定最终预测类别,即若E(θ|y)>τ,则判定为事故类,反之则判定为非事故类(安全类)。
步骤4-4),基于获取的BF样本训练步骤4-3)所述BF序列风险预测模型,得到最终标定的BF序列风险预测模型。所述步骤4-4)具体通过以下子步骤来实现:
步骤4-4-1),对步骤4-1)中所有的样本按步骤4-2)所述方法构建BF序列风险预测模型样本的输入和输出,最终得到共4D个BF序列风险预测模型样本。
步骤4-4-2),根据高速公路历史交通事故数据,提取每个样本对应的季节、周天数、日时段和天气类型低频时变变量。以低频时变变量和步骤(1)获取的样本对应路段常变量为该样本的交通时空模式特征变量TMt=[z1t,z2t,...,zst,...,zSt],基于汉明距离在步骤(2)获得的K 类事故模式TMk=[z1k,z2k,...,zsk,...,zSk],k=1,2,...,K中进行模式匹配。汉明距离计算方法为:初始化距离值为0;对每一个特征变量zs,若zst=zsk,则距离值加1,否则距离值保持不变;循环该过程,至所有特征变量[z1,z2,...,zs,...,zS]遍历结束。选取与TMt汉明距离最小的模式TMc为该样本所处的交通时空模式。
步骤4-4-4),按训练集和测试集分别占75%和25%随机抽取步骤4-4-1)所获样本构成相应集合,以预测F1值为目标函数,通过网格搜索法获取BF序列风险预测模型其他参数C、 N、τ的最优值C*、N*、τ*,其中各参数备选值:C∈[0.5,1]、N∈[5,15]、τ∈[0.1,1)。最终得到基于{C=C*,N=N*,τ=τ*}最优参数组合的BF序列风险预测模型。
步骤二,线上风险预警模型实时预测
卡口实时信息采集***实时采集每一辆经过卡口车辆的牌照、车辆类型、经过时刻和所在车道,通过步骤一所述特征提取方法实时提取时间窗内上下游卡口间的多步高频时变特征变量,输入LSTM瞬时风险判别模型,得到t时刻的LSTM瞬时风险判别模型判别结果。分别通过卡口实时信息采集***和交通地理信息库实时获取步骤一所述低频时变变量和路段常变量,基于汉明距离在步骤一获得的事故模式库中实时匹配当前交通时空模式,得到当前模式频数与其他模式频数的占比。以规定周期内各观测点上LSTM瞬时风险判别模型判别结果的序列为输入向量,以模式频数占比为风险先验概率,输入BF序列风险预测模型,最终得到未来路段风险状态预测结果。
步骤二具体为:
步骤(1),对于预测时刻t,通过卡口实时信息采集***实时获取步骤一所述低频时变变量,通过交通地理信息库实时获取步骤一所述预测路段常变量。以低频时变变量和预测路段常变量为特征变量TMt=[z1t,z2t,...,zst,...,zSt],基于汉明距离在步骤一获得的K类事故模式 TMk=[z1k,z2k,...,zsk,...,zSk],k=1,2,...,K中进行实时模式匹配。汉明距离计算方法与步骤一中所述一致。选取与TMt汉明距离最小的模式TMc为预测时刻预测路段所处的交通时空模式。
步骤(2),对于预测路段l,利用该路段上游和下游卡口信息***分别对每一辆经过上游和下游卡口的车辆拍照并自动识别牌照,具体采集的车辆信息包括车辆牌照、车辆经过时刻、车辆所在车道和车辆类型(小型车/大型车)。基于所述采集信息,按步骤一所述BF序列风险预测模型输入向量构建方法实时确定预测时刻t的观测周期(即t-N+1,t-N+2,…,t共N个观测点),并按步骤一所述LSTM模型多步高频时变特征变量构建方法得到对应N个观测点的N 组LSTM瞬时风险判别模型输入特征,分别输入步骤一训练得到的LSTM瞬时风险判别模型,得到周期内N个观测点上LSTM判别结果构成的序列y=[y1,...,yN]。
步骤(3),以步骤(1)得到的模式TMc内事故记录的观测频数fc占所有模式总频数的占比为步骤一所述BF序列风险预测模型中参数a的值,即参数以步骤(2)得到的LSTM判别结果序列y=[y1,...,yN]为BF序列风险预测模型的输入向量,输入步骤一所述基于参数组合的BF序列风险预测模型,得到预测时刻预测路段的事故风险概率E(θ|y)。
步骤(4),根据步骤(3)得到的事故风险概率制定相应的高速公路风险预警策略:当步骤(3)所述BF序列风险预测模型输出风险概率E(θ|y)<τ1:预测路段处于低风险状态,高速公路管理部门无需采取任何处理措施;当τ1≤E(θ|y)<τ2:预测路段处于中风险状态,高速公路管理部门应加强对该路段的监控,做好预警准备;当τ2≤E(θ|y):预测路段处于高风险状态,高速公路管理部门应对进出该路段交通流采取平稳措施,对驾驶员予以预警提醒,并做好事故救援准备。其中τ1,τ2∈[0,1],且τ1<τ2,分别代表中、高风险状态的阈值。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (6)
1.一种基于LSTM和BF的高速公路风险预测方法,其特征在于,包括步骤:
S1,线下风险预警预测模型训练
基于历年高速公路历史交通事故数据,得到K类事故交通时空模式库,统计得到各事故交通时空模式的发生频数;基于事故发生前指定时间内上下游卡口间的时空范围,通过卡口车牌识别历史数据提取多步高频时变变量,构建区分事故与安全状态的LSTM瞬时风险判别模型;以规定周期内各观测点上LSTM瞬时风险判别模型判别结果的序列为输入向量,利用低频时变变量和常变量匹配的交通时空模式,建立基于先验概率的BF序列风险预测模型;
S2,线上风险预警模型实时预测
实时获取预测时刻的低频时变变量和预测路段常变量,实时匹配预测时空的交通时空模式;实时采集每一辆经过卡口车辆的牌照、车辆类型、经过时刻和所在车道,提取时间窗内上下游卡口间的多步高频时变变量,输入LSTM瞬时风险判别模型,得到t时刻的LSTM瞬时风险判别模型判别结果;以规定周期内各观测点上LSTM瞬时风险判别模型的输出序列为输入向量,以匹配交通时空模式内的观测频数占比为先验概率,输入BF序列风险预测模型,最终得到未来路段风险状态预测结果;
所述LSTM瞬时风险判别模型的构建过程,具体为:
1)基于历史交通事故,以预测时刻t0前tw内上下游卡口M1M2之间的时空范围为输入时间窗,提取T步输入特征向量X={x1,x2,…,xt,…,xT},输出从预测时刻起未来tp时长内有事故发生,得到事故样本;所述预测时刻t0是事故发生时刻tc前tp时间,tp为预测时长;
2)对每一个事故样本同一时段、同一路段的无事故时空范围进行随机下采样,以预测时刻t0′前tw内上下游卡口M1M2之间的时空范围为输入时间窗,提取T步输入特征向量X′={x1′,x2′,…,xt′,…,x′T},输出从预测时刻起未来tp时长内无事故发生,得到非事故样本;
3)融合事故样本和非事故样本,训练LSTM瞬时风险判别模型,使得损失函数最小,得到最终标定的LSTM瞬时风险判别模型;
所述LSTM瞬时风险判别模型的输出序列值视为二分类随机变量y∈{0,1},随机变量y∈{0,1}的取值由参数θ决定,参数θ的期望值表示为:
其中:y为LSTM瞬时风险判别模型的输出序列,m和q分别代表LSTM输出序列事故类和非事故类变量的个数,a和b是beta分布的超参数,分别代表事故类和非事故类的先验概率;
将LSTM瞬时风险判别模型的输出视为二分类随机变量y∈{0,1}的观测样本,建立基于递减系数、先验概率和阈值的BF序列风险预测模型:
加入随观测时间变化的递减函数rn,参数θ的期望值表示为:其中:n=1,2,...,N代表LSTM序列的第n个输出值的序号,N为序列最后一个输出值的序号,C∈(0,1]为递减常数,mn和qn代表LSTM序列第n个输出结果的相应标签值;
通过设定阈值τ为BF序列风险预测模型判定最终预测类别;
每个样本所处的交通时空模式是通过以下方法确定的:
1)基于高速公路历史交通事故数据,以事故样本的低频时变变量和路段常变量为特征变量[z1,z2,...,zs,...,zS],通过K-modes聚类方法得到K类事故模式,每一类事故模式的聚类中心为TMk=[z1k,z2k,...,zsk,...,zSk],k=1,2,...,K;
2)以待确定样本的低频时变变量和路段常变量为该样本的交通时空模式特征变量TMt=[z1t,z2t,...,zst,...,zSt],在TMk=[z1k,z2k,...,zsk,...,zSk]中基于汉明距离进行模式匹配,与TMt汉明距离最小的模式TMc即为该样本所处的交通时空模式;
所述未来路段风险状态预测是根据E(θ|y)与阈值的大小进行判定的,当E(θ|y)<τ1:预测路段处于低风险状态,当τ1≤E(θ|y)<τ2:预测路段处于中风险状态,当τ2≤E(θ|y):预测路段处于高风险状态,其中τ1、τ2分别代表中、高风险状态的阈值,E(θ|y)为参数θ的期望值。
3.根据权利要求1所述的基于LSTM和BF的高速公路风险预测方法,其特征在于,所述LSTM瞬时风险判别模型的隐藏层和Dense全连接输出层中间设有随机失活Dropout,且LSTM瞬时风险判别模型配合使用输入连接权重L2正则化和学习率衰减。
4.根据权利要求1所述的基于LSTM和BF的高速公路风险预测方法,其特征在于,所述高频时变变量包括上下游卡口流量、车道间流量差、车道平均流量、路段流量密度和大型/小型车流量比。
5.根据权利要求1所述的基于LSTM和BF的高速公路风险预测方法,其特征在于,所述低频时变变量包括季节、周天数、日时段和天气类型。
6.根据权利要求1所述的基于LSTM和BF的高速公路风险预测方法,其特征在于,所述常变量包括路段车道数、道路线形、路段进出口匝道间距和路段进出口匝道个数,并将所述路段进出口匝道间距连续变量按取值区间离散化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110979482.XA CN113744526B (zh) | 2021-08-25 | 2021-08-25 | 一种基于lstm和bf的高速公路风险预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110979482.XA CN113744526B (zh) | 2021-08-25 | 2021-08-25 | 一种基于lstm和bf的高速公路风险预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113744526A CN113744526A (zh) | 2021-12-03 |
CN113744526B true CN113744526B (zh) | 2022-12-23 |
Family
ID=78732681
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110979482.XA Active CN113744526B (zh) | 2021-08-25 | 2021-08-25 | 一种基于lstm和bf的高速公路风险预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113744526B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114582131B (zh) * | 2022-03-17 | 2023-08-29 | 中远海运科技股份有限公司 | 一种基于匝道智能控流算法的监控方法及*** |
CN114944055B (zh) * | 2022-03-29 | 2023-04-18 | 浙江省交通投资集团有限公司智慧交通研究分公司 | 基于电子收费门架的高速公路碰撞风险动态预测方法 |
CN114863680B (zh) * | 2022-04-27 | 2023-04-18 | 腾讯科技(深圳)有限公司 | 预测处理方法、装置、计算机设备及存储介质 |
CN115376308B (zh) * | 2022-05-26 | 2024-06-04 | 南京工程学院 | 一种汽车行驶时间的预测方法 |
CN115565373B (zh) * | 2022-09-22 | 2024-04-05 | 中南大学 | 高速公路隧道事故实时风险预测方法、装置、设备及介质 |
CN116030627B (zh) * | 2022-12-31 | 2024-04-30 | 东南大学 | 一种融合预测交通风险变量的道路交通事故分析方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646534A (zh) * | 2013-11-22 | 2014-03-19 | 江苏大学 | 一种道路实时交通事故风险控制方法 |
CN106355883A (zh) * | 2016-10-20 | 2017-01-25 | 同济大学 | 基于风险评估模型的交通事故发生概率获取方法及*** |
CN107742193A (zh) * | 2017-11-28 | 2018-02-27 | 江苏大学 | 一种基于时变状态转移概率马尔可夫链的行车风险预测方法 |
WO2018103313A1 (zh) * | 2016-12-06 | 2018-06-14 | 杭州海康威视数字技术股份有限公司 | 发生交通事故的风险预测方法、装置及*** |
CN108198415A (zh) * | 2017-12-28 | 2018-06-22 | 同济大学 | 一种基于深度学习的城市快速路事故预测方法 |
CN112949999A (zh) * | 2021-02-04 | 2021-06-11 | 浙江工业大学 | 一种基于贝叶斯深度学习的高速交通事故风险预警方法 |
CN112990545A (zh) * | 2021-02-08 | 2021-06-18 | 东南大学 | 一种高速公路交织区交通安全状态预测方法 |
-
2021
- 2021-08-25 CN CN202110979482.XA patent/CN113744526B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646534A (zh) * | 2013-11-22 | 2014-03-19 | 江苏大学 | 一种道路实时交通事故风险控制方法 |
CN106355883A (zh) * | 2016-10-20 | 2017-01-25 | 同济大学 | 基于风险评估模型的交通事故发生概率获取方法及*** |
WO2018103313A1 (zh) * | 2016-12-06 | 2018-06-14 | 杭州海康威视数字技术股份有限公司 | 发生交通事故的风险预测方法、装置及*** |
CN107742193A (zh) * | 2017-11-28 | 2018-02-27 | 江苏大学 | 一种基于时变状态转移概率马尔可夫链的行车风险预测方法 |
CN108198415A (zh) * | 2017-12-28 | 2018-06-22 | 同济大学 | 一种基于深度学习的城市快速路事故预测方法 |
CN112949999A (zh) * | 2021-02-04 | 2021-06-11 | 浙江工业大学 | 一种基于贝叶斯深度学习的高速交通事故风险预警方法 |
CN112990545A (zh) * | 2021-02-08 | 2021-06-18 | 东南大学 | 一种高速公路交织区交通安全状态预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113744526A (zh) | 2021-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113744526B (zh) | 一种基于lstm和bf的高速公路风险预测方法 | |
Park et al. | Real-time prediction and avoidance of secondary crashes under unexpected traffic congestion | |
CN104834977B (zh) | 基于距离度量学习的交通警情等级预测方法 | |
CN102750824B (zh) | 基于多种神经网络分类器投票的城市道路交通状态检测方法 | |
CN112215487B (zh) | 一种基于神经网络模型的车辆行驶风险预测方法 | |
Florio et al. | Neural-network models for classification and forecasting of freeway traffic flow stability | |
Yan et al. | Spatial-temporal chebyshev graph neural network for traffic flow prediction in iot-based its | |
CN112085947A (zh) | 一种基于深度学习和模糊聚类的交通拥堵预测方法 | |
CN109345832B (zh) | 一种基于深度递归神经网络的城市道路超车预测方法 | |
CN101216998A (zh) | 基于模糊粗糙集的证据理论城市交通流信息融合方法 | |
CN112990545B (zh) | 一种高速公路交织区交通安全状态预测方法 | |
Wang et al. | Crash prediction for freeway work zones in real time: A comparison between Convolutional Neural Network and Binary Logistic Regression model | |
Pawar et al. | Classification of gaps at uncontrolled intersections and midblock crossings using support vector machines | |
CN113449905A (zh) | 一种基于门控循环单元神经网络的交通拥堵预警方法 | |
CN117238126A (zh) | 一种连续流道路场景下的交通事故风险评估方法 | |
Zhu et al. | Early identification of recurrent congestion in heterogeneous urban traffic | |
Yang et al. | A traffic dynamic operation risk assessment method using driving behaviors and traffic flow Data: An empirical analysis | |
Florio et al. | Neural network models for classification and forecasting of freeway traffic flow stability | |
Sun et al. | Vision-based traffic conflict detection using trajectory learning and prediction | |
CN116798223A (zh) | 基于宏观基本图/fcm聚类的子区划分及状态识别方法 | |
Huang | SVM‐Based Real‐Time Identification Model of Dangerous Traffic Stream State | |
Obereigner et al. | Methods for Traffic Data Classification with regard to Potential Safety Hazards | |
Shi et al. | CPT‐DF: Congestion Prediction on Toll‐Gates Using Deep Learning and Fuzzy Evaluation for Freeway Network in China | |
Abeyratne et al. | Applying big data analytics on motor vehicle collision predictions in New York City | |
Srinivas et al. | Integrating Smart Technologies for Enhanced Traffic Management and Road Safety: A Data-Driven Approach |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20221208 Address after: 557100 Yunguan Office Building, Yongle Road, Danjiang Town, Leishan County, Qiandongnan Miao and Dong Autonomous Prefecture, Guizhou Province Applicant after: Guizhou Qiantong Zhilian Technology Co.,Ltd. Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301 Applicant before: JIANGSU University |
|
GR01 | Patent grant | ||
GR01 | Patent grant |