基于跟驰特性的智能网联车混合交通流轨迹重构

蒋阳升; 刘梦; 王思琛; 高宽; 姚志洪

doi:10.3969/j.issn.0258-2724.20200735

基于跟驰特性的智能网联车混合交通流轨迹重构

doi: 10.3969/j.issn.0258-2724.20200735

蒋阳升^{1, 2,},
刘梦^{1, 2},
王思琛^{1, 2},
高宽^{1, 2},
姚志洪^{1, 2}

1.
西南交通大学交通运输与物流学院，四川成都 611756
2.
西南交通大学综合交通大数据应用技术国家工程实验室，四川成都 611756

基金项目: 国家自然科学基金（52002339）；四川省科技计划（2021YJ0535，2020YFH0026）；广西壮族自治区科技计划（2021AA01007AA）；中央高校基本科研业务费专项资金（2682021CX058）

详细信息

作者简介:
蒋阳升（1976—），男，教授，博士，研究方向为交通系统优化，E-mail：jiangyangsheng@swjtu.cn

中图分类号: U491
计量
- 文章访问数: 934
- HTML全文浏览量: 452
- PDF下载量: 45
- 被引次数: 0
出版历程
- 收稿日期: 2020-11-10
- 修回日期: 2021-01-13
- 网络出版日期: 2021-04-15
- 刊出日期: 2021-04-15

Trajectory Reconstruction for Traffic Flow Mixed withConnected Automated Vehicles Based on Car-Following Characteristics

JIANG Yangsheng^{1, 2
,},
LIU Meng^{1, 2},
WANG Sichen^{1, 2},
GAO Kuan^{1, 2},
YAO Zhihong^{1, 2}

1.
School of Transportation and Logistics，Southwest Jiaotong University，Chengdu 611756，China
2.
National Engineering Laboratory of Integrated Transportation Big Data Application Technology，Southwest Jiaotong University，Chengdu 611756，China

摘要

摘要:
车辆轨迹数据蕴含着丰富的时空交通信息，是交通状态估计的基础数据之一. 为解决现有数据采集环境难以获得全样本车辆轨迹的问题，面向智能网联环境，构建了混合交通流全样本车辆轨迹重构模型. 首先，分析了智能网联环境下混合交通流的车辆构成及其轨迹数据采集环境；然后，提出了基于智能驾驶员跟驰模型的车辆轨迹重构模型，实现了对插入轨迹数量、轨迹位置和速度等参数的估计；最后，设计仿真试验验证了模型在不同交通流密度和智能网联车（connected automated vehicle，CAV）渗透率条件下的适用性. 试验结果表明：CAV和网联人工驾驶车（connected vehicle，CV）的渗透率为8%和20%时，该车辆轨迹重构模型在不同交通流密度下均能重构84%以上的车辆轨迹；重构轨迹准确性随着CAV和CV渗透率的增加而提高；当交通密度为70辆/km，且CAV渗透率仅为4%的情况下，模型也能重构82%的车辆轨迹.
- 智能交通 /
- 跟驰模型 /
- 智能网联车 /
- 混合交通流
Abstract:
Vehicle trajectory data contains massive spatial-temporal traffic information, which is one of the necessary data for traffic state estimation. To solve the problem that it's difficult to obtain the fully sampled vehicular trajectory in the existing data collection environment, oriented to the connected and automated environment, a fully sampled trajectory reconstruction model of mixed traffic flow is proposed . Firstly, vehicle composition and trajectory data collection environment of mixed traffic flow with the connected automated vehicle (CAV) are analyzed. Then, a vehicle trajectory reconstruction model is proposed based on intelligent driver car-following model. Based on this, the number of inserted trajectories, trajectory position and speed are estimated. Finally, numerical simulation is designed to investigate the influence of traffic density and penetration rate of CAVs. Results show that, when the penetration rates of CAV and connected vehicle (CV) at 8% and 20%, respectively, the model can reconstruct more than 84% vehicular trajectories under different traffic densities. The accuracy of the reconstructed trajectories increases with the increase in penetration rates of CAV and CV. Besides, when the traffic density is 70 veh/km and the penetration rate of CAV at a low level of 4%, the proposed model can reconstruct 82% vehicular trajectories.
- intelligent transportation system (ITS) /
- car-following model /
- connected automated vehicle /
- mixed traffic flow

HTML全文

随着无人驾驶技术的发展，未来道路上的智能网联车（connected automated vehicle，CAV）将逐渐普及. 但最新研究表明^[1]，预计2045年道路上L4级别的CAV渗透率仅能达到24.8%. 因此，未来很长一段时间内，道路上将普遍呈现由常规车（regular vehicle，RV）、网联人工驾驶车（connected vehicle，CV）和CAV组成的混合交通流. 网联车（CV和CAV）能够提供海量的车辆轨迹数据，而这些数据蕴含了丰富的时空交通信息，是交通状态估计的基础数据之一. 但由于RV无法提供自身的轨迹数据，在智能网联环境下全样本车辆轨迹数据仍难以获得. 因此，如何利用已有的轨迹数据来重构所有车辆的轨迹具有重要的意义^[2].

近年来，有许多学者对车辆轨迹重构问题进行了探讨，主要为基于经典交通流理论和基于统计概率模型的方法. 这些方法按轨迹重构数据来源可分为三大类：定点检测器、移动检测器和多源数据. Coifman^[3]基于双线圈检测器数据和交通流理论，通过在时空图中估计路段行程时间重构车辆轨迹. Rao等^[4]利用车牌识别数据构建了一种车辆路径重构方法，进一步可基于路径对交通网络的历史OD （original-destination）进行估计. Hao等^[5]、Wan等^[6]、Shan等^[7]基于概率论的思想，分别采用随机模型、期望最大化算法（expectation-maximum，EM）及最大似然估计，从微观角度重构浮动车任意相邻采样点之间的运动轨迹. Li等^[8]提出了基于K最近邻回归算法利用不完整轨迹来扩充轨迹样本量. Xie等^[9]基于修正的交通信号配时方案、车辆通过时间和车辆行程时间构建微观仿真环境，实现了对干道车辆轨迹重构. Mehran等^[10]通过融合秒级浮动车数据、音视频数据（audio video interleaved，AVI）及信号控制参数，提出了一种基于变分理论城市干道车辆运行轨迹估计方法，并利用东京一条包含五个交叉口的单行道实测数据进行了方法验证和误差分析. 唐克双等^[11]通过融合视频和定点检测器流量数据，基于交通流理论和交通仿真思想，构建了一种不依赖高频浮动车数据，适用于多车道的轨迹重构方法. 综上分析可知，现有研究均能够对车辆轨迹数据进行有效重构，但均只考虑了人工驾驶环境，未考虑未来智能网联混合交通流的情况. 随着智能网联技术的发展，CAV除了能提供自身轨迹数据外，还能通过感知设备（摄像头和雷达）提供自身一定范围内的车辆轨迹数据. 因此，CAV采集的移动传感数据将会给轨迹重构提供一种新的思路. 如何充分考虑混合交通流环境下不同车辆类型提供的数据，构建全样本车辆轨迹重构模型，是未来CAV环境下车辆轨迹重构亟待解决的关键问题.

为解决该问题，本文拟研究混有RV、CV和CAV的混合交通流全样本车辆轨迹重构方法. 充分利用CAV和CV提供的轨迹数据信息，提出基于跟驰特性的混合交通流轨迹重构模型，并通过仿真试验验证模型在不同交通流密度和CAV渗透率条件下的合理性与有效性.

1. 跟驰特性分析

在本文研究中，不考虑车辆换道的影响，选取较为简单的场景来阐述轨迹重构思想. 在微观交通流分析中，一般采用跟驰模型描述车辆之间的跟驰行为. 该模型认为当前车辆是否加减速与本身速度、前后车速度差和位置差有关. 因此，在已知部分轨迹的情况下，采用该模型可分析相邻轨迹中存在其他车辆轨迹的可能性.

跟驰模型主要可分为五种类型：刺激反应模型、安全间距模型、社会力模型、优化速度模型和低阶线性模型^[12]，均可作为本文轨迹重构的基础模型. 而作为社会力模型的一种，智能驾驶员模型（intelligent driver model，IDM）能够很好地描述熟练驾驶员的驾驶习惯，应用范围更广. 因此，本文以IDM为例描述车辆之间的跟驰特性，如式（1）所示.

$\begin{split} & {a_p}\left( t \right) = {a_{\rm{M}}}\left[ {1 - {{\left( {\dfrac{{{v_p}\left( t \right)}}{{{v_{\rm{f}}}}}} \right)}^\delta } - } \right. \\&\quad \left. {{{\left( {\dfrac{{{s_0} + T{v_p}\left( t \right) + \dfrac{{{v_p}\left( t \right)\left[ {{v_p}\left( t \right) - {v_{p - 1}}\left( t \right)} \right]}}{{2\sqrt {{a_{\rm{M}}}b} }}}}{{{x_{p - 1}}\left( t \right) - {x_p}\left( t \right)}}} \right)}^2}} \right] \end{split},$

(1)

式中： ${a_p}(t)$ 为时刻t车辆 $p$ 的加速度； ${v_{\rm{f}}}$ 为期望速度； ${a_{\rm{M}}}$ 为最大加速度；T为安全车头时距； $b$ 为舒适减速度的绝对值； ${s_0}$ 为最小停车间距； ${v_p}(t)$ 为时刻t车辆 $p$ 的速度； ${x_p}(t)$ 为时刻t车辆 $p$ 的位置； $\delta$ 为加速度指数（一般取4）.

2. 基于IDM的车辆状态估计

2.1 数据环境分析

本文研究的混合交通流中有三种类型的车辆：RV、CV和CAV. 其中，CAV除提供自身的轨迹数据外，还可通过感知设备获取周围一定范围内的车辆位置和状态信息，因此，能提供周围一定范围内车辆的轨迹数据. 混合交通流中不同车辆组成如图1所示. CAV探测范围内的车辆信息均可采集，不论是RV还是CV. CV由于不具备感知设备，所以只能提供自身的轨迹数据. 而RV由于其不具备或未打开联网功能，所以既不能提供自身的轨迹数据，也不能提供周围车辆的轨迹数据.

图 1 混合交通流中不同车辆组成

Figure 1. Composition of different types of vehicles in mixed traffic flow

下载: 全尺寸图片幻灯片

从以上分析可知，在该混合交通流环境下，主要需要对CAV探测范围之外的RV进行轨迹重构. 混合交通流轨迹时空图如所示. CAV两辆（CAV1、CAV2），CV三辆（CV1、CV2、CV3）；将阴影区域A中前一辆CAV探测到的最后一辆车记为 ${{\rm{V}}_{n - 1}}$ ，后面一辆CV记为 ${{\rm{V}}_n}$ ；红色虚线表示CAV轨迹；蓝色点线表示CV轨迹；黑色实线表示RV轨迹；阴影区域A和B即为需要重构轨迹的区域. 在跟驰模型中，本车的跟驰行为主要与前车有关，即前车的轨迹直接影响本车的轨迹. 以阴影区域A为例，若区域A中无其他车辆， ${{\rm{V}}_n}$ 会跟随 ${{\rm{V}}_{n - 1}}$ ，本车 ${{\rm{V}}_n}$ 与前车 ${{\rm{V}}_{n - 1}}$ 构成跟驰行为，可采用跟驰模型进行描述. 但若区域A中还存在未被观测的RV，则本车 ${{\rm{V}}_n}$ 用跟驰模型计算的加速度会与实际加速度存在较大差异. 因此，构建轨迹重构模型试图寻找区域A内的最佳轨迹条数，来实现这种差异的最小化.

图 2 混合交通流轨迹时空图

Figure 2. Spatio-temporal diagram of mixed traffic flow trajectories

下载: 全尺寸图片幻灯片

CV提供的轨迹数据可分为是否在CAV的检测范围内. 当CV在CAV的检测范围内，如图2中的CV1和CV3，这两辆车的轨迹可直接由CAV观测，故该车提供的数据有无对整个模型的求解没有影响. 当CV不在CAV检测范围内时，如图2中的CV2，此时需要利用该轨迹和CAV提供的轨迹重构阴影区域A和B中的车辆轨迹.

2.2 估计插入车辆的速度

车辆时空轨迹示意如所示，图中： ${A_1}{\text{,}} {A_2}{\text{,}}\cdots{\text{,}} {A_k}$ 为时空区域； $k$ 为时间总步数； $\Delta t$ 为时间步长； ${x_n}\left( {\Delta t} \right)$ 为 $\Delta t$ 时 ${{\rm{V}}_n}$ 的位置. 考虑到交通流稳定时，前后车的速度差别较小. 因此，第k个时空区域 ${A_k}$ 中车辆 ${{\rm{V}}_n}$ 的速度 $v_n$ 的估计如式（2）所示^[13].

图 3 车辆轨迹示意

Figure 3. Schematic diagram of vehicle trajectories

下载: 全尺寸图片幻灯片

$\begin{split} & v_n = \frac{{d\left( {{A_k}} \right)}}{{T\left( {{A_k}} \right)}} =\\ &\quad \frac{{ {{x_{n - 1}}\left( {t + \Delta t} \right) - {x_{{{n}} - 1}}(t)} + {{x_{{n}}}\left( {t + \Delta t} \right) - {x_{{n}}}(t)} }}{{2\Delta t}} \end{split},$

(2)

式中： $d\left( {{A_k}} \right)$ 和 $T\left( {{A_k}} \right)$ 分别为区域 ${A_k}$ 的位置差和时间差.

由于插入的RV在区域 ${A_k}$ 中，因此，该车在时刻 $k\Delta t$ 的速度可以用区域 ${A_k}$ 的速度 $v\left( {{A_k}} \right)$ 进行估计.

2.3 估计插入车辆的位置

在任意时刻，插入车辆的位置都需满足与前后车之间的安全间距. 因此，将V_n前插入的第i辆车记为 ${{\rm{V}}_{n,i}}$ ， ${{\rm{V}}_{n,i}}$ 的位置x_n，i只能出现在一定的区间范围内，如式（3）所示.

${x_n}\left( t \right) + {s_0} + l \leqslant {x_{n{\text{,}}i}}\left( t \right) \leqslant {x_{n - 1}}\left( t \right) - {s_0} - l,$

(3)

式中：l为车身长度.

当获得插入车辆的速度和可能位置区间后，可将 ${{\rm{V}}_{n,i}}$ 和 ${{\rm{V}}_n}$ 分别作为前车和后车，而 ${{\rm{V}}_{n,i}}$ 插入位置 ${x_{n,i}}\left( t \right)$ 会对后车 ${{\rm{V}}_n}$ 产生影响，这种影响会反映到后车的加速度上. 因此，可以根据IDM代入式（1）计算 ${{\rm{V}}_n}$ 的加速度 ${{\hat a}_n}\left( {{x_{n,i}}\left( t \right)} \right)$ .

当 ${x_{n,i}}\left( t \right)$ 在可行域内取不同的值时， ${{\rm{V}}_n}$ 会得到不同的加速度. 而当估计加速度 ${\hat a_n}\left( {{x_{n,i}}(t)} \right)$ 与实际加速度 ${a_n}\left( t \right)$ 最接近时，对应的 ${x_{n,i}}\left( t \right)$ 即为要插入车辆 ${{\rm{V}}_{n,i}}$ 的位置.

因此，可构建目标函数求得最优的插入车辆位置，同时考虑到车辆不同时刻的位置即为车辆的轨迹，可选用 ${{\rm{V}}_n}$ 的实际加速度和估计加速度之间均方根误差最小为优化目标，如式（4）所示.

$\min \sqrt {\frac{1}{k}{{\sum\limits_{t = 0}^{k\Delta t} {\left[ {{{\hat a}_n}\left( {{x_{n,i}}(t)} \right) - {a_n}\left( t \right)} \right]^2} }}}.$

(4)

综上，该优化模型的目标函数为式（4），约束条件为 ${{x_{n,i}}}$ 和 ${{\hat a}_n}( {{x{_{n,i}(t)}}})$ ，为典型的非线性优化模型，可用MATLAB自带的函数直接进行求解.

3. 基于IDM的车辆轨迹重构

由第2节的分析可知，得出插入车辆后，可由模型计算其插入的位置和速度. 因此，车辆轨迹重构模型即为需要寻找已知轨迹中间未知轨迹的数量，插入车辆轨迹示意如图4所示. 当插入车辆数已知时，可采用第2节中的方法依次求得每辆插入车辆轨迹的具体位置和速度. 考虑到车辆之间的最小安全间距，插入车辆数m的最大可能取值为N，如式（5）所示.

图 4 插入车辆轨迹示意

Figure 4. Schematic diagram of inserted vehicle trajectories

下载: 全尺寸图片幻灯片

$\begin{split} & N = {\rm{ent}}\left( {\rm{min}} \left\{ {\dfrac{{{x_{n - 1}}\left( 0 \right) - {x_n}\left( 0 \right)}}{{{s_0} + l}},\dfrac{{x_{n - 1}}\left( {\Delta t} \right) - {x_n}\left( {\Delta t} \right)}{{{s_0} + l}},}\cdots, \right. \right. \\ &\quad \left. \left. { \dfrac{{{x_{n - 1}}\left( {k\Delta t} \right) - {x_n}\left( {k\Delta t} \right)}}{{{s_0} + l}}}\right\}\right) . \\[-20pt] \end{split}$

(5)

在车辆 ${{\rm{V}}_n}$ 和车辆 ${{\rm{V}}_{n - 1}}$ 之间依次插入的m辆车( ${{\rm{V}}_{n,1}},{{\rm{V}}_{n,2}},\cdots,{{\rm{V}}}_{n,m}$ ）， ${{\rm{V}}_{n,m}}$ 为距离 ${{\rm{V}}_n}$ 最近的车辆，即为 ${{\rm{V}}_n}$ 的前车. 插入m辆车后，类比式（2），可粗略估计 ${{\rm{V}}_{n,m}}$ 的速度和位置分别如式（6）、（7）所示.

${v_{n,m}}\left( t \right) = {v_n}\left( t \right) - \frac{{ {{v_n}\left( t \right) - {v_{n - 1}}\left( t \right)} }}{{ {m + 1} }},$

(6)

${x_{n,m}}\left( t \right) = {x_n}\left( t \right) - \frac{{{{x_n}\left( t \right) - {x_{n - 1}}\left( t \right)} }}{{{m + 1} }}.$

(7)

在估计插入的速度和位置后，可重构 ${{\rm{V}}_{n,m}}$ 的轨迹. 根据IDM由 ${{\rm{V}}_{n,m}}$ 的轨迹可求出 ${{\rm{V}}_n}$ 的加速度. 插入m辆RV后代入式（1）可计算 ${{\rm{V}}_n}$ 的理论加速度.

此处，选用 ${{\rm{V}}_n}$ 加速度的均方根误差作为衡量标准，则均方根误差最小时对应的 $m$ 即为最佳插入的车辆数，如式（8）所示.

$\mathop {{\rm{arg min}}}\limits_m \sqrt {\frac{1}{k}{{\sum\limits_{t = 0}^{k\Delta t} {\left[ {{{\hat a}_n}\left( {{x_{n,i}}(t)} \right) - {a_n}\left( t \right)} \right]^2} }}}.$

(8)

因此，该问题也为一个非线性优化问题，目标函数为式（8），约束条件为式（1）、（6）、（7）. 综上，当确定最佳插入车辆数量以后，可采用第2节中的方法对每一插入车辆轨迹进行精确地重构.

4. 数值仿真分析

4.1 仿真参数设置

为验证模型的有效性，基于IDM设计元胞自动机进化规则，模拟单车道高速公路交通流运行规律. 所有车辆之间跟驰行为均用IDM描述^[14-15]. 仿真路段长度1000 m，元胞长度0.1 m，仿真时间1000 s，仿真步长1 s. 为使仿真数据更趋向于真实数据，在仿真过程给车辆设置随机减速模拟真实道路上的车流运行状况.

本文共设计两种基准试验方案，第一种方案是将CAV的渗透率恒定为8%，CV渗透率恒定为20%，改变试验中的交通流密度，研究不同交通流密度对模型的影响. 第二种方案将交通流密度设定为70 辆/km，改变CAV和CV的渗透率，研究不同渗透率对模型的影响. 本文IDM的参数^[16]取值为： ${s_0}$ =2 m， ${a_{\rm{M}}}$ =1 m/s²， $T$ =1.5 s， $b$ =2 m/s²， $l$ =4.5 m.

4.2 交通流密度的影响

为研究道路上交通流密度对本文模型结果的影响，将CAV和CV的渗透率分别设定为8%和20%. 首先，可基于仿真试验获取路段的流量-密度-速度关系如图5所示.

图 5 流量-密度-速度关系

Figure 5. Relationship between volume，density and speed of traffic flow

下载: 全尺寸图片幻灯片

由图5可知：仿真路段中的最佳密度为40 辆/km，即当交通流密度大于40 辆/km时，路段上的车辆开始出现拥挤；当交通流密度大于70 辆/km时，车辆的平均速度低于20 km/h，而高速公路中很少出现车辆速度小于20 km/h的情况. 因此，为研究本文模型在不同交通状态下的重构效果，此处将交通流密度设为20、30、40、50、60、70 辆/km. 为避免随机因素的影响，同一交通流密度下采用不同的随机种子进行仿真和轨迹重构，取其平均值作为最终结果. 轨迹重构图像示例如图6所示. 不同交通流密度条件下求得的插入车辆数如表1所示.

图 6 交通流密度为60辆/km时轨迹时空图

Figure 6. Spatio-temporal diagram of the trajectories when traffic flow density is 60 veh/km

下载: 全尺寸图片幻灯片

由表1可知：当交通流密度小于等于40 辆/km时，会出现插入的车辆数小于实际车辆数；而当交通流密度大于40 辆/km时，则出现插入的车辆数大于实际车辆数. 从百分比误差来看，当出现拥堵前，随着交通流密度的增加，模型求得的插入车辆数的误差逐渐减小，插入辆数的百分比误差也逐渐减小；当出现拥堵后，随着交通流密度的增加，模型求得的插入车辆数的误差逐渐增加，插入辆数的百分比误差也逐渐增大. 分析可知，开始拥堵之前，大多数车辆自由行驶，且车辆之间相互干扰较少，所以误差较大. 开始拥堵后，当两辆车之间插入的车辆数较大时，模型求得的插入车辆数会产生一定的误差. 随着交通流密度的增大，车辆之间的间距逐渐减小，出现两辆车之间插入的车辆数较大的情况增多，故插入辆数的误差逐渐增大.

表 1 不同交通流密度条件下插入车辆数估计

Table 1. Estimation of the numbers of inserted vehicles under different traffic densities

交通流密度/ （辆•km⁻¹）	CAV/ 辆	插入车辆/辆	实际车辆/辆	插入误差/辆	插入误差百分比/%
20	2	5.80	6.40	0.60	9.38
30	2	14.05	13.80	0.35	2.54
40	3	15.80	16.05	0.25	1.56
50	4	21.75	21.38	0.31	1.50
60	5	27.85	25.75	2.10	8.16
70	6	31.40	27.25	4.30	15.78
注：表中插入车辆数为多次仿真试验的平均值，因此会出现不是整数的情况；插入误差百分比为插入辆数的误差与实际插入车辆数的百分比.

下载: 导出CSV

| 显示表格

当交通流密度小于60 辆/km时，本文模型求得的插入车辆数的误差小于1辆，误差小于10%. 综上分析可知，本文模型求得的插入车辆数较为准确，适用于高速公路的各种交通状态.

为进一步评估本文模型的轨迹插入效果，选用每一时刻插入车辆的实际和估计位置之间的平均绝对误差（MAE，e_MAE）、平均绝对百分比误差（MAPE，e_MAE）和均方根误差（RMSE，e_MAE）表征重构轨迹的准确性，分别如式（9）～（11）所示.

${e_{\rm{MAE}}} = \frac{1}{k}\sum\limits_{t = 0}^{k\vartriangle t} {\left| {\hat x(t) - x(t)} \right|},$

(9)

${e_{\rm{MAPE}}} = \frac{1}{k}\sum\limits_{t = 0}^{k\Delta t} {\left| {\frac{{\hat x(t) - x(t)}}{{x(t)}}} \right|}\; \times \; {100{\text{%}} },$

(10)

${e_{\rm{RMSE}}} = \sqrt {\frac{1}{k}\sum\limits_{t = 0}^{k\Delta t} {{{\left[ {\hat x(t) - x(t)} \right]}^2}} },$

(11)

式中： $\hat x(t)$ 为插入车辆的估计位置； $x(t)$ 为插入车辆的实际位置.

根据式（9）～（11），可得不同交通流密度下的误差值，如表2所示. 由表2可知：当交通流密度为20 辆/km时，由于需要插入的车辆数过少，所以重构轨迹车辆位置的平均绝对误差相较其他密度较高. 当交通流密度大于20 辆/km时，重构轨迹车辆位置的误差随着交通流密度的增大而增加，但变动幅度不大. 在不同的交通流密度条件下，重构轨迹车辆位置的平均绝对误差的平均值为12.97 m，平均绝对百分比误差的平均值为0.37%，均方根误差的平均值为19.55 m. 从整体上来说，本文模型在不同交通流密度下重构轨迹车辆位置的误差较为稳定.

表 2 不同交通流密度条件下轨迹重构误差

Table 2. Trajectory reconstruction error under different traffic density conditions

交通流密度/ （辆•km⁻¹）	CAV/ 辆	MAE/ m	MAPE/ %	RMSE/ m
20	2	15.02	0.36	18.28
30	2	11.90	0.29	14.30
40	3	12.14	0.31	18.68
50	4	12.30	0.35	20.39
60	5	13.09	0.42	22.19
70	6	13.39	0.52	23.47
平均值		12.97	0.37	19.55

下载: 导出CSV

| 显示表格

4.3 CAV和CV渗透率的影响

CAV和CV提供的数据作为模型的已知数据，研究其渗透率对模型结果的影响尤为重要.

4.3.1 CAV渗透率的影响

为研究不同CAV渗透率对模型结果的影响，将交通流密度设为70 辆/km，CV渗透率设为20%. CAV的检测范围设为100 m. 若车辆在道路上均匀分布，则CAV检测范围包含7辆车. 但考虑到车辆的随机分布和感知设备的误差，此处假设CAV能记录检测范围内车辆数为6. 将CAV渗透率分别设置为4%、6%、8%、10%、12%和14%，研究不同渗透率下的车辆轨迹重构效果，其中，插入车辆数估计误差如表3所示，轨迹重构结果示例如图7所示.

表 3 不同渗透率条件下插入车辆数估计

Table 3. Estimation of the numbers of inserted vehicles under different CAV penetration rates

CAV 渗透率/%	CAV/辆	插入车辆/辆	实际车辆/辆	插入误差/辆	插入误差百分比/%
4	3	45.50	38.70	6.80	17.57
6	4	41.50	35.20	5.80	16.48
8	6	31.40	27.25	4.30	15.78
10	7	29.44	25.44	4.00	15.72
12	8	26.80	23.20	3.60	15.52
14	10	19.60	17.00	2.60	15.29

下载: 导出CSV

| 显示表格

图 7 CAV渗透率为14%时轨迹时空图

Figure 7. Spatio-temporal diagram of the trajectories when CAV penetration rate is 14%

下载: 全尺寸图片幻灯片

由表3可知：当道路上的交通流密度为70 辆/km时，随着CAV渗透率的增加，需要插入的车辆数减少，轨迹重构的难度减小，求出的插入车辆数的误差和百分比误差均逐渐减小；当渗透率只有4%的情况下，模型也能够较好地估计插入轨迹的数量. 因此，本文模型在未来不同渗透率环境下均能很好应用. 进一步，可计算重构轨迹车辆位置误差如表4所示.

表 4 CAV不同渗透率条件下轨迹重构误差

Table 4. Trajectory reconstruction error under different CAV penetration rates

CAV渗透率/%	CAV/辆	MAE/m	MAPE/%	RMSE/m
4	3	16.76	0.68	33.77
6	4	15.69	0.63	28.95
8	6	13.39	0.52	23.47
10	7	12.98	0.51	21.86
12	8	12.23	0.49	19.90
14	10	10.45	0.42	16.45

下载: 导出CSV

| 显示表格

由表4可知：当道路上的交通流密度为70 辆/km时，随着CAV渗透率的增加，重构轨迹车辆位置的误差逐渐减小，即本文模型的准确性逐渐提高. 以MAE为例，当CAV的渗透率在4%～8%范围内时，平均每增加一辆CAV，平均绝对误差减小1.12 m. 当CAV的渗透率在8%～14%范围内时，平均每增加一辆CAV，平均绝对误差减小0.73 m. 因此，当CAV渗透率大于8%时，平均绝对误差趋于稳定.

4.3.2 CV渗透率的影响

将交通流密度设为70 辆/km，CAV渗透率设为8%，研究CV渗透率为16%、18%、20%、22%、24%时轨迹重构结果，如表5所示.

表 5 CV不同渗透率轨迹重构结果

Table 5. Reconstruction results of different CV penetration rates

CV 渗透率/%	CV/辆	插入误差/辆	插入误差百分比/%	MAE/m	MAPE/ %	RMSE/ m
16	11	5.40	0.18	17.46	0.69	32.43
18	13	4.33	0.16	13.98	0.55	24.27
20	14	4.30	0.16	13.39	0.52	23.47
22	15	4.20	0.15	12.93	0.52	21.50
24	17	3.80	0.14	11.50	0.48	19.11

下载: 导出CSV

| 显示表格

由表5可知：随着CV渗透率的提高，插入辆数的误差逐渐减小，插入数量的百分比误差也逐渐减小，即轨迹重构效果不断提升. 此外，车辆位置的误差也越小，即轨迹重构的准确性有所提高.

为研究CAV和CV渗透率对试验结果的影响，以CAV渗透率为横轴，CV渗透率为纵轴，绘制均方根误差热力图，如图8所示.

图 8 均方根误差热力图

Figure 8. Heat map of RMSE

下载: 全尺寸图片幻灯片

由图8可知：均方根误差趋势为随着CAV或CV渗透率的提高，轨迹重构效果逐渐变好. 但与CV渗透率相比，CAV渗透率对结果的影响更大. 分析可知，增加一辆CAV能够获得道路上多辆车的轨迹数据，而增加一辆CV，只能获取该车的轨迹数据. 因此，当CAV渗透率较高时，已知轨迹的车辆较多，此时CV对试验结果的影响较小.

5. 结　论

1）在出现拥堵前，随着交通流密度的增加，本文模型求得的插入车辆数的误差逐渐减小；而在出现拥堵后，随着交通流密度的增加，本文模型求得的插入车辆数的误差逐渐增大.

2）本文模型重构轨迹车辆位置的平均绝对误差在不同的交通流密度条件下均在可接受范围内，即本文模型在不同的交通流密度下均能较为准确地重构所有车辆的轨迹.

3）重构轨迹的准确性随着CAV和CV渗透率的增加而提高. 与CV渗透率相比，CAV渗透率对结果的影响更大. 且当CAV渗透率较低时，本文模型也能较好地重构所有车辆的轨迹.

本文证明了基于跟驰特性的车辆轨迹重构模型的有效性，但选取的场景简单，仅研究了单车道高速公路车辆轨迹重构模型，未来将结合换道模型，研究智能网联环境下多车道车辆轨迹重构模型.

致谢：西南交通大学综合交通大数据应用技术国家工程实验室交大数科创新中心项目（JDSKCXZX202003）和重庆交通大学重庆市交通运输工程重点实验室开放课题（2018TE01）的资助.

图 1 混合交通流中不同车辆组成

Figure 1. Composition of different types of vehicles in mixed traffic flow

下载: 全尺寸图片幻灯片

图 2 混合交通流轨迹时空图

Figure 2. Spatio-temporal diagram of mixed traffic flow trajectories

下载: 全尺寸图片幻灯片

图 3 车辆轨迹示意

Figure 3. Schematic diagram of vehicle trajectories

下载: 全尺寸图片幻灯片

图 4 插入车辆轨迹示意

Figure 4. Schematic diagram of inserted vehicle trajectories

下载: 全尺寸图片幻灯片

图 5 流量-密度-速度关系

Figure 5. Relationship between volume，density and speed of traffic flow

下载: 全尺寸图片幻灯片

图 6 交通流密度为60辆/km时轨迹时空图

Figure 6. Spatio-temporal diagram of the trajectories when traffic flow density is 60 veh/km

下载: 全尺寸图片幻灯片

图 7 CAV渗透率为14%时轨迹时空图

Figure 7. Spatio-temporal diagram of the trajectories when CAV penetration rate is 14%

下载: 全尺寸图片幻灯片

图 8 均方根误差热力图

Figure 8. Heat map of RMSE

下载: 全尺寸图片幻灯片

表 1 不同交通流密度条件下插入车辆数估计

Table 1. Estimation of the numbers of inserted vehicles under different traffic densities

交通流密度/ （辆•km⁻¹）	CAV/ 辆	插入车辆/辆	实际车辆/辆	插入误差/辆	插入误差百分比/%
20	2	5.80	6.40	0.60	9.38
30	2	14.05	13.80	0.35	2.54
40	3	15.80	16.05	0.25	1.56
50	4	21.75	21.38	0.31	1.50
60	5	27.85	25.75	2.10	8.16
70	6	31.40	27.25	4.30	15.78
注：表中插入车辆数为多次仿真试验的平均值，因此会出现不是整数的情况；插入误差百分比为插入辆数的误差与实际插入车辆数的百分比.

下载: 导出CSV

表 2 不同交通流密度条件下轨迹重构误差

Table 2. Trajectory reconstruction error under different traffic density conditions

交通流密度/ （辆•km⁻¹）	CAV/ 辆	MAE/ m	MAPE/ %	RMSE/ m
20	2	15.02	0.36	18.28
30	2	11.90	0.29	14.30
40	3	12.14	0.31	18.68
50	4	12.30	0.35	20.39
60	5	13.09	0.42	22.19
70	6	13.39	0.52	23.47
平均值		12.97	0.37	19.55

下载: 导出CSV

表 3 不同渗透率条件下插入车辆数估计

Table 3. Estimation of the numbers of inserted vehicles under different CAV penetration rates

CAV 渗透率/%	CAV/辆	插入车辆/辆	实际车辆/辆	插入误差/辆	插入误差百分比/%
4	3	45.50	38.70	6.80	17.57
6	4	41.50	35.20	5.80	16.48
8	6	31.40	27.25	4.30	15.78
10	7	29.44	25.44	4.00	15.72
12	8	26.80	23.20	3.60	15.52
14	10	19.60	17.00	2.60	15.29

下载: 导出CSV

表 4 CAV不同渗透率条件下轨迹重构误差

Table 4. Trajectory reconstruction error under different CAV penetration rates

CAV渗透率/%	CAV/辆	MAE/m	MAPE/%	RMSE/m
4	3	16.76	0.68	33.77
6	4	15.69	0.63	28.95
8	6	13.39	0.52	23.47
10	7	12.98	0.51	21.86
12	8	12.23	0.49	19.90
14	10	10.45	0.42	16.45

下载: 导出CSV

表 5 CV不同渗透率轨迹重构结果

Table 5. Reconstruction results of different CV penetration rates

CV 渗透率/%	CV/辆	插入误差/辆	插入误差百分比/%	MAE/m	MAPE/ %	RMSE/ m
16	11	5.40	0.18	17.46	0.69	32.43
18	13	4.33	0.16	13.98	0.55	24.27
20	14	4.30	0.16	13.39	0.52	23.47
22	15	4.20	0.15	12.93	0.52	21.50
24	17	3.80	0.14	11.50	0.48	19.11

下载: 导出CSV

参考文献(16)

[1]	BANSAL P, KOCKELMAN K M. Forecasting Americans’ long-term adoption of connected and autonomous vehicle technologies[J]. Transportation Research Part A:Policy and Practice, 2017, 95: 49-63. doi: 10.1016/j.tra.2016.10.013
[2]	PUNZO V, BORZACCHIELLO M T, CIUFFO B. On the assessment of vehicle trajectory data accuracy and application to the Next Generation SIMulation (NGSIM) program data[J]. Transportation Research Part C:Emerging Technologies, 2011, 19(6): 1243-1262. doi: 10.1016/j.trc.2010.12.007
[3]	COIFMAN B. Estimating travel times and vehicle trajectories on freeways using dual loop detectors[J]. Transportation Research Part A:Policy and Practice, 2002, 36(4): 351-364. doi: 10.1016/S0965-8564(01)00007-6
[4]	RAO W M, WU Y J, XIA J X, et al. Origin-destination pattern estimation based on trajectory reconstruction using automatic license plate recognition data[J]. Transportation Research Part C:Emerging Technologies, 2018, 95: 29-46. doi: 10.1016/j.trc.2018.07.002
[5]	HAO P, BORIBOONSOMSIN K, WU G Y, et al. Modal activity-based stochastic model for estimating vehicle trajectories from sparse mobile sensor data[J]. IEEE Transactions on Intelligent Transportation Systems, 2017, 18(3): 701-711. doi: 10.1109/TITS.2016.2584388
[6]	WAN N F, VAHIDI A, LUCKOW A. Reconstructing maximum likelihood trajectory of probe vehicles between sparse updates[J]. Transportation Research Part C:Emerging Technologies, 2016, 65: 16-30. doi: 10.1016/j.trc.2016.01.010
[7]	SHAN X N, HAO P, CHEN X H, et al. Vehicle energy/emissions estimation based on vehicle trajectory reconstruction using sparse mobile sensor data[J]. IEEE Transactions on Intelligent Transportation Systems, 2019, 20(2): 716-726. doi: 10.1109/TITS.2018.2826571
[8]	LI Z, KLUGER R, HU X B, et al. Reconstructing vehicle trajectories to support travel time estimation[J]. Transportation Research Record:Journal of the Transportation Research Board, 2018, 2672(42): 148-158. doi: 10.1177/0361198118772956
[9]	XIE X, VAN LINT H, VERBRAECK A. A generic data assimilation framework for vehicle trajectory reconstruction on signalized urban arterials using particle filters[J]. Transportation Research Part C:Emerging Technologies, 2018, 92: 364-391. doi: 10.1016/j.trc.2018.05.009
[10]	MEHRAN B, KUWAHARA M, NAZNIN F. Implementing kinematic wave theory to reconstruct vehicle trajectories from fixed and probe sensor data[J]. Transportation Research Part C:Emerging Technologies, 2012, 20(1): 144-163. doi: 10.1016/j.trc.2011.05.006
[11]	唐克双,徐天祥,潘昂,等. 基于定点检测数据的城市干道车辆轨迹重构[J]. 同济大学学报(自然科学版),2016,44(10): 1545-1552. TANG Keshuang, XU Tianxiang, PAN Ang, et al. Signal timing and detector data-based reconstruction of vehicle trajectories on urban arterials[J]. Journal of Tongji University (Natural Science), 2016, 44(10): 1545-1552.
[12]	王昊, 金诚杰. 交通流理论及应用[M]. 北京: 人民交通出版社, 2020: 50-79.
[13]	WANG Y P, WEI L, CHEN P. Trajectory reconstruction for freeway traffic mixed with human-driven vehicles and connected and automated vehicles[J]. Transportation Research Part C: Emerging Technologies, 2020, 111: 135-155. doi: 10.1016/j.trc.2019.12.002
[14]	TALEBPOUR A, MAHMASSANI H S. Influence of connected and autonomous vehicles on traffic flow stability and throughput[J]. Transportation Research Part C: Emerging Technologies, 2016, 71: 143-163. doi: 10.1016/j.trc.2016.07.007
[15]	YAO Z H, XU T R, JIANG Y S, et al. Linear stability analysis of heterogeneous traffic flow considering degradations of connected automated vehicles and reaction time[J]. Physica A:Statistical Mechanics and Its Applications, 2021, 561: 125218.1-125218.15. doi: 10.1016/j.physa.2020.125218
[16]	MILANÉS V, SHLADOVER S E. Modeling cooperative and autonomous adaptive cruise control dynamic responses using experimental data[J]. Transportation Research Part C:Emerging Technologies, 2014, 48: 285-300. doi: 10.1016/j.trc.2014.09.001

施引文献

附加材料(0)

访问统计

点击查看大图

图(8) / 表(5)

计量

文章访问数: 934
HTML全文浏览量: 452
PDF下载量: 45
被引次数: 0

1. 跟驰特性分析
2. 基于IDM的车辆状态估计
2.1 数据环境分析
2.2 估计插入车辆的速度
2.3 估计插入车辆的位置
3. 基于IDM的车辆轨迹重构
4. 数值仿真分析
4.1 仿真参数设置
4.2 交通流密度的影响
4.3 CAV和CV渗透率的影响
5. 结　论

基于跟驰特性的智能网联车混合交通流轨迹重构

doi: 10.3969/j.issn.0258-2724.20200735

作者简介: 蒋阳升（1976—），男，教授，博士，研究方向为交通系统优化，E-mail：jiangyangsheng@swjtu.cn

计量

出版历程

Trajectory Reconstruction for Traffic Flow Mixed withConnected Automated Vehicles Based on Car-Following Characteristics

1. 跟驰特性分析

2. 基于IDM的车辆状态估计

2.1 数据环境分析

2.2 估计插入车辆的速度

2.3 估计插入车辆的位置

3. 基于IDM的车辆轨迹重构

4. 数值仿真分析

4.1 仿真参数设置

4.2 交通流密度的影响

4.3 CAV和CV渗透率的影响

4.3.1 CAV渗透率的影响

4.3.2 CV渗透率的影响

5. 结 论

计量

出版历程

目录

1. 跟驰特性分析

2. 基于IDM的车辆状态估计

2.1 数据环境分析

2.2 估计插入车辆的速度

2.3 估计插入车辆的位置

3. 基于IDM的车辆轨迹重构

4. 数值仿真分析

4.1 仿真参数设置

4.2 交通流密度的影响

4.3 CAV和CV渗透率的影响

5. 结 论

作者简介:
蒋阳升（1976—），男，教授，博士，研究方向为交通系统优化，E-mail：jiangyangsheng@swjtu.cn

5. 结　论

5. 结　论