CN106571034B - 基于融合数据的城市快速路交通状态滚动预测方法 - Google Patents

Abstract

本发明涉及一种基于多源异构数据融合的城市快速路交通状态滚动预测方法，包括步骤如下：针对待预侧路段获取检测点的卡口、微波数据；通过卡口的匹配车辆获取其真实的旅行时间，通过微波检测器获取检测点的真实速度，按照其空间相对位置，进行多源异构数据的融合；通过非线性插值算法获得任意非检测点速度，并构建虚拟轨迹线获取其估计旅行时间；通过优化目标函数(速度和旅行时间的均方根误差)，对其参数进行标定；利用标定好的参数对下一时间段的速度和旅行时间进行预测。本发明方法简便有效，便于操作，避免的传统方法在参数标定的过程中，一经标定，即不发生改变的弊端。

Description

基于融合数据的城市快速路交通状态滚动预测方法

技术领域

本发明涉及交通状态的预测领域，具体是一种基于多源异构数据融合的交通状态滚动预测方法。

背景技术

在城市道路交通中，交通流是由成千上万的人、车参与的，每个交通的参与者都会基于自己的行为目标与行为准则，依据不同的交通运行环境做出符合自身利益的决策。这种决策主要来源于对未来交通流运行状况的估计，只有准确的把握交通流运行过程中的行为变化，才能使制定出的决策更富有针对性。因此将交通流的动态性预测作为日常监测和管理的重要内容，是交通管理决策科学的基础，也是主动、全方位交通管理的关键。

对于出行者来说，需要实时交通状态信息和对未来交通状态的准确预测，从而做出有效的决策，例如决定出发时间，进行路径选择等。对于实时交通状态建立合理的模型，设计简单高效的算法来准确地预测未来的交通状态。交通流模型描述了车辆和驾驶者的动态行为，传统的交通状态预测模型参数依赖于对历史数据的校准，而且模型参数一经标定，便不再改变，在此情况下进行标定的参数就存在了一些不足的地方，不能够适应实时的交通状态预测，存在着应用的局限性。

发明内容

本发明为克服上述不足之处，提供一种基于多源异构数据融合的交通状态滚动预测方法，本发明通过获取卡口、微波数据，并在此基础上进行数据分析和参数的标定，实现对下一时刻的路况进行预测。本发明方法简便有效，便于操作，克服了传统的静态参数标定方法不能应用于实时交通状态预测的不足，解决了现有交通参数标定方法适用范围的局限性问题。

一种基于多源异构数据融合的城市快速路交通状态滚动预测方法，包括步骤如下：

(1)针对交通状态待预测的路段，获取检测点的位置，所述检测点包括布设在该路段的卡口和微波检测器。

(2)获得路段的真实旅行时间T(t_j)和检测点(x_i,t_i)处的真实速度v(x_i,t_i)，其中，t_j为在卡口处获取车辆驶入上匝道的时间，结合该车辆驶出下匝道时间，从而得到车辆的真实旅行时间T(t_j)，根据实测道路长度计算车辆在该路段的空间平均速度，并以该平均速度作为车辆在卡口处的行驶速度，即获得卡口处的真实行驶速度；通过微波探测器获取车辆的时间平均速度，从而得到该路段上微波检测器位置处的行驶速度，即得到微波检测器位置点的真实速度。

(3)将该路段上检测点的行驶速度，按照其空间相对位置进行多源异构数据融合。

(4)在步骤(3)的融合数据基础上，重建时空速度场。根据已知检测点的真实速度v(x_i,t_i)，求该路段上除检测点以外的非检测点的估计速度，具体采用如下非线性插值算法：

V(x,t)＝w(x,t)V_cong(x,t)+[1-w(x,t)]V_free(x,t)

式中，V_thr临界速度，ΔV为自由流和阻塞流转换之间的带宽，c_free为自由流速度，c_cong消散波速度。t_i，x_i，v_i分别为已知检测点所在的时刻、位置、速度。t，x，V分别为任意非检测点所在的时刻、位置、速度。σ为空间平滑宽度，τ是时间平滑宽度。

(5)计算检测点的估计速度，具体为：针对某一检测点，通过其余检测点的真实速度，结合步骤(4)所述的非线性插值算法，计算出检测点的估计速度

(6)根据步骤(2)中检测点的真实速度v(x_i,t_i)和步骤(4)中非检测点的估计速度，构建虚拟轨迹线，对旅行时间进行估计，得到估计旅行时间

(7)微波数据估计速度的标准均方根误差为

卡口数据估计旅行时间的标准均方根误差为

其中，v(x_i,t_i)为真实速度，为估计速度，T(t_j)为真实旅行时间，为估计旅行时间，N为微波检测数据量，M为卡口旅行时间数据量。

(8)采用交通状态预测滚动优化算法，通过计算其目标函数，对空间平滑宽度σ、时间平滑宽度τ两个参数进行标定：

式中：D_microwave(t₀-T_h→t₀)为从t₀-T_h时刻到t₀时刻的微波数据，D_LPR(t₀-T_h→t₀)为同一时间段的卡口数据，t₀为当前时刻，T_h为历史时间窗口，g为目标函数，即两种数据的标准均方根误差之和。

(9)根据步骤(8)标定后的空间平滑宽度σ、时间平滑宽度τ，采用步骤(4)所述的非线性的插值算法，结合当前时刻采集的真实速度，对该路段上任意位置的下一时刻行驶速度进行预测，进一步根据预测后的速度，采用步骤(6)所述的虚拟轨迹线算法，对下一时刻的旅行时间进行预测。

本发明的有益效果是采用滚动优化思想，在进行交通预测的过程中，依据实时数据，不断调整和更新模型参数，使得预测结果更加贴合实际，从而避免传统方法在参数标定过程中，一经标定，即不发生改变，造成交通状态预测结果与实际出现很大偏差的弊端。

附图说明

图1路段拓扑结构和检测器布设位置；

图2望江路上匝道，登云路下匝道所匹配的车辆为例，构建其虚拟轨迹线；

图3是本发明交通状态预测滚动方法示意图。

具体实施方式

本发明是基于国家自然科学基金青年基金项目(51508505)和浙江省自然科学基金杰出青年项目(LR17E080002)的研究，提出一种基于多源异构融合数据的城市快速路交通状态滚动预测方法。下面结合具体实施例对本发明进行进一步描述，但本发明的保护范围并不仅限于此。

以杭州市上塘-中河高架快速路为例，选取由南向北行驶的一段长度为9.7km的路段。其中共包含9个上下匝道(4个上匝道和5个下匝道)，以及6个微波检测器。路段拓扑结构和检测器布设位置如图1所示。通过以下方法对该路段的交通状态进行短时预测。

(1)针对交通状态待预测的路段，获取检测点的位置，所述检测点包括布设在该路段的卡口和微波检测器。

(2)获得路段的真实旅行时间T(t_j)和检测点的真实速度v(x_i,t_i)：在卡口处获取车辆驶入上匝道和驶出下匝道时间，从而得到车辆的真实旅行时间T(t_j)，根据实测道路长度计算车辆在该路段的空间平均速度，并以该平均速度作为车辆在卡口处的行驶速度，即获得卡口处的真实行驶速度，其中上下匝道匹配车辆数和平均旅行时间如表1所示；通过微波探测器获取车辆的时间平均速度，从而得到该路段上微波检测器位置处的行驶速度，即得到微波检测器位置点的真实速度。

表1

注：表中数字为匹配车辆数、括号中数字为匹配车辆的平均旅行时间(单位：秒)。

(3)将该路段上检测点的行驶速度，按照其空间相对位置进行多源异构数据融合。

(4)在步骤(3)的融合数据基础上，重建时空速度场。根据已知检测点的真实速度v(x_i,t_i)，求该路段上除检测点以外的非检测点的估计速度，具体采用如下非线性插值算法：

V(x,t)＝w(x,t)V_cong(x,t)+[1-w(x,t)]V_free(x,t)

式中，V_thr临界速度，ΔV为自由流和阻塞流转换之间的带宽，c_free为自由流速度，c_cong消散波速度。t_i，x_i，v_i分别为已知检测点所在的时刻、位置、速度。t，x，V分别为任意非检测点所在的时刻、位置、速度。σ为空间平滑宽度，τ是时间平滑宽度。

(5)计算检测点的估计速度，具体为：针对某一检测点，通过其余检测点的真实速度，结合步骤(4)所述的非线性插值算法，计算出检测点的估计速度

(6)根据步骤(2)中检测点的真实速度v(x_i,t_i)和步骤(4)中非检测点的估计速度，按照现有方法，构建虚拟轨迹线，对旅行时间进行估计，得到估计旅行时间以上塘-中河高架望江路上匝道，登云路下匝道所匹配的车辆为例，构建其虚拟轨迹线如图2所示。具体如下：

(6.1)给定足够小的时间增量Δt和空间增量Δx；

(6.2)如果一辆车从点A(t_A,x_A)出发，那么该辆车的虚拟轨迹线可以通过重复以下步骤构建；

(6.3)假设车辆当前的位置是点(t,x)，采用所使用的插值算法计算出该点的速度v(t,x)；

(6.4)如果v(t,x)＝0，给时间一个增量即t'＝t+Δt，更新当前位置至x'＝x，进入第(6.6)步；

(6.5)否则以v(t,x)为斜率，更新至下一点的位置即t'＝t+Δt，x'＝x+Δx，然后对新点进行速度的插值计算；

(6.6)检查x'≥x_F，若为否，则转入第(6.3)步；

(6.7)否则，计算车辆在x_F处时的时间，T_F＝t+(t'-t)*(x_F-x)/(x'-x)；

(6.8)车辆的估计旅行时间为T_F-T_A，在此过程中，车辆由A点行驶至F点。

(7)微波数据估计速度的标准均方根误差为

卡口数据估计旅行时间的标准均方根误差为

其中，v(x_i,t_i)为真实速度，为估计速度，T(t_j)为真实旅行时间，为估计旅行时间。

(8)采用交通状态预测滚动优化算法，通过计算其目标函数，对空间平滑宽度σ、时间平滑宽度τ两个参数进行标定：

式中：D_microwave(t₀-T_h→t₀)为从t₀-T_h时刻到t₀时刻的微波数据，D_LPR(t₀-T_h→t₀)为同一时间段的卡口数据，t₀为当前时刻，T_h为历史时间窗口，g为目标函数，即两种数据的标准均方根误差之和。

(9)根据步骤(8)标定后的空间平滑宽度σ、时间平滑宽度τ，采用步骤(4)所述的非线性的插值算法，结合当前时刻采集的真实速度，对该路段上任意位置的下一时刻行驶速度进行预测，进一步根据预测后的速度，采用步骤(6)所述的虚拟轨迹线算法，对下一时刻的旅行时间进行预测。交通状态滚动预测示意如图3所示，采用传统默认值方法和滚动优化算法预测结果如表2所示。

表2默认与滚动优化算法预测结果。

从表2可以看出，滚动优化算法的均方根误差(RMSE)、标准均方根误差(NRMSE)对称绝对百分比误差(SMAPE1、SMAPE2)均低于默认算法，从而验证了滚动优化算法优于一般的默认值方法。

Claims (1)

1.一种基于多源异构融合数据的城市快速路交通状态滚动预测方法，包括步骤如下：

(1)针对交通状态待预测的路段，获取检测点的位置，所述检测点包括布设在该路段的卡口和微波检测器；

(2)获得路段的真实旅行时间T(t_j)和检测点(x_i,t_i)处的真实速度v(x_i,t_i)，其中，t_j为在卡口处获取车辆驶入上匝道的时间，结合该车辆驶出下匝道时间，从而得到车辆的真实旅行时间T(t_j)，根据实测道路长度计算车辆在该路段的空间平均速度，并以该平均速度作为车辆在卡口处的行驶速度，即获得卡口处的真实行驶速度；通过微波探测器获取车辆的时间平均速度，从而得到该路段上微波检测器位置处的行驶速度，即得到微波检测器位置点的真实速度；

(3)将该路段上检测点的行驶速度，按照其空间相对位置进行多源异构数据融合；

(4)在步骤(3)的融合数据基础上，重建时空速度场；根据已知检测点的真实速度v(x_i,t_i)，求该路段上除检测点以外的非检测点的估计速度，具体采用如下非线性插值算法：

V(x,t)＝w(x,t)V_cong(x,t)+[1-w(x,t)]V_free(x,t)

式中，V_thr临界速度，ΔV为自由流和阻塞流转换之间的带宽，c_free为自由流速度，c_cong消散波速度；t_i，x_i，v_i分别为已知检测点所在的时刻、位置、速度；t，x，V分别为任意非检测点所在的时刻、位置、速度；σ为空间平滑宽度，τ是时间平滑宽度；

(5)计算检测点的估计速度，具体为：针对某一检测点，通过其余检测点的真实速度，结合步骤(4)所述的非线性插值算法，计算出检测点的估计速度

(6)根据步骤(2)中检测点的真实速度v(x_i,t_i)和步骤(4)中非检测点的估计速度，构建虚拟轨迹线，对旅行时间进行估计，得到估计旅行时间

(7)微波数据估计速度的标准均方根误差为

卡口数据估计旅行时间的标准均方根误差为

其中，v(x_i,t_i)为真实速度，为估计速度，T(t_j)为真实旅行时间，为估计旅行时间，N为微波检测数据量，M为卡口旅行时间数据量；

(8)采用交通状态预测滚动优化算法，通过计算其目标函数，对空间平滑宽度σ、时间平滑宽度τ两个参数进行标定：

式中：D_microwave(t₀-T_h→t₀)为从t₀-T_h时刻到t₀时刻的微波数据，D_LPR(t₀-T_h→t₀)为同一时间段的卡口数据，t₀为当前时刻，T_h为历史时间窗口，g为目标函数，即两种数据的标准均方根误差之和；

(9)根据步骤(8)标定后的空间平滑宽度σ、时间平滑宽度τ，采用步骤(4)所述的非线性的插值算法，结合当前时刻采集的真实速度，对该路段上任意位置的下一时刻行驶速度进行预测，进一步根据预测后的速度，采用步骤(6)所述的虚拟轨迹线算法，对下一时刻的旅行时间进行预测。