一种针对可变限速控制的元胞传输仿真模型改进方法
【专利摘要】一种针对可变限速控制的元胞传输仿真模型改进方法。首先检测某段包含进出口匝道的快速道路的交通流密度、流量和速度。采用最优线性拟合法分别拟合数据得到包含通行能力下降和可变限速控制下的交通流基本图。基于实时交通流数据确定各时段主线起点及出入口匝道的交通需求。在各仿真步，采用两类基本图分别计算各元胞内流入流出的最大车辆数，更新各元胞内交通流量、密度及速度。本发明弥补了基本元胞传输模型假设过于简单且缺乏考虑可变限速控制影响的不足，模拟了瓶颈区域通行能力下降、拥堵交通流中时走时停现象及不同驾驶员遵从度下可变限速控制对交通流运行的影响，为基于仿真的计算可变限速控制对交通安全及效率的影响奠定了基础。
【专利说明】
一种针对可变限速控制的元胞传输仿真模型改进方法
技术领域
[0001] 本发明属于交通仿真技术领域，具体涉及一种针对可变限速控制的元胞传输仿真 模型的改进方法。
【背景技术】
[0002] 交通仿真是研究复杂交通问题的重要工具，通过一定数学模型对交通流运行进行 描述，可以动态、准确地模拟交通流运行规律，对交通流时间和空间变化特征予以再现。元 胞传输模型作为一种运动学模型的差分离散形式，被广泛应用于模拟快速道路交通流运行 之中，但其基本假设过于简单，对其实行二次开发可以准确地模拟关键交通流特征与现象。 同时，随着越来越多的快速道路采用可变限速控制这一动态交通控制策略，在交通仿真模 型中能够准确模拟可变限速对交通运行的影响变得愈发重要。因此，提出一种针对可变限 速控制的元胞传输仿真模型的改进方法，能够更有针对性的刻画可变限速控制下的交通流 特性以及一些关键交通流特性。
[0003] 目前宏观仿真模型中，元胞传输模型是可变限速控制中采用较多的仿真模型之 一。但基本元胞传输模型对交通流状态的假设过于简单，无法模拟快速道路上通行能力下 降及时走时停交通流特征，也无法直接模拟不同驾驶员遵从度下可变限速控制对交通流影 响。本发明提出基于基本元胞传输模型的针对可变限速控制的对元胞传输仿真模型的改 进，相比基本元胞传输仿真模型更具有针对性且对交通流状态的假设更加合理，改进后的 模型可以仿真得到一些以往模型无法仿真得到的典型的交通流特征。
【发明内容】
[0004] 本发明要解决的问题是:基本元胞传输模型作为一种被广泛应用于快速道路交通 流运行模拟的模型，能够模拟快速道路上交通拥堵生成与传播过程，但基本模型无法模拟 快速道路上通行能力下降及时走时停交通流特征，也无法直接模拟不同驾驶员遵从度下可 变限速控制对交通流影响。本发明提出一种针对可变限速控制的元胞传输仿真模型改进方 法，分别提出了针对可变限速控制下的元胞传输仿真模型中各步骤中规则修正的方法，建 立了能够模拟瓶颈区域通行能力下降现象、拥堵交通流中时走时停现象以及不同驾驶员遵 从度下可变限速控制对交通流运行影响的元胞传输模型。二次开发后的元胞传输模型弥补 了之前基本元胞传输模型对交通流状态的假设过于简单的缺陷与不足。
[0005] 本发明技术方案为：
[0006] 本发明提出一种科学合理的针对可变限速控制的元胞传输仿真模型的改进方法， 基于基本元胞传输模型各元胞交通流状态更新步骤，依据不同的交通流状态的假设修改模 型中发送能力和接受能力的计算方法，实现了针对可变限速控制下对元胞传输模型的二次 开发，弥补了以往基本元胞传输仿真模型无法模拟一些特有交通现象的缺陷。实例显示，本 发明提出的元胞传输仿真模型改进方法很好的实现了预期效果，比基本元胞传输模型更好 的模拟了交通流特性。
【附图说明】
[0007] 图1为本发明元胞传输仿真模型改进方法流程。
[0008] 图2为典型路段的几何特征与仿真模型示意图。
[0009] 图3为瓶颈元胞中包含通行能力下降的基本图。
[0010] 图4为可变限速控制下交通流基本图。
[0011] 图5为三种条件下车辆传递示意图。
【具体实施方式】
[0012] 本发明是基于快速道路上通行能力下降及时走时停交通流特征，以及不同驾驶员 遵从度下可变限速控制对交通流的影响特性，修改基本元胞传输模型中每个元胞内允许流 入和流出的最大车辆数的计算公式，基于此结果计算各元胞内的交通流量、交通密度和交 通流速度，其中可变限速控制的效果在交通流速度的计算公式修正中得以体现，从而得到 针对可变限速控制的元胞传输仿真模型的三种情况下的改进方法。针对可变限速控制的元 胞传输仿真模型改进方法流程可参考图1。
[0013] 第一步是确定各时间段主线起点、入口匝道以及出口匝道上交通需求。具体方法 为，选定一段长度的快速路段，路段内需包含一个入口匝道和一个出口匝道以及三个交通 流检测器，记录其主线路段车道数。选定数据检测的时间间隔T，仿真时间步长ΔΤ。典型路 段的几何特征与元胞传输仿真模型中元胞位置如图2所示。用设置的交通流检测器以时间 间隔T周期性采集主线起点、入口匝道和出口匝道上的平均交通量，并以时间间隔T周期性 采集各位置交通流平均速度。对检测所得的交通流数据进行处理集计求得各位置交通量作 为交通需求，依据检测到的交通流速度计算自由流车流的平均速度，元胞长度由每时间间 隔内车辆以自由流速度行驶的距离确定。
[0014] 第二步是在某时刻k，依据各元胞内交通流状态及交通流基本图的右枝，计算每个 元胞内允许流入的最大车辆数δ (k)。
[0015] 改进后的元胞传输模型在计算每个元胞内允许流入的最大车辆数时考虑了对瓶 颈区通行能力下降和时走时停的交通流特性仿真。具体步骤为首先确定瓶颈元胞内时走时 停运动波的概率阈值Φ，然后在k时刻产生介于0到1之间(包含0和1)的随机数ζ(1〇,依据下 式判断不同条件相应的瓶颈元胞最大接受车辆数的计算公式：
[0016]
(I)
[0017] 其中，
[0018] δ (k)为瓶颈元胞最大接受车辆数；
[0019] ζ〇〇为介于[0，1]之间随机数；
[0020] Φ为产生时走时停波动的概率阈值；
[0021] VTH为判断交通流是否拥堵的速度阈值(通常取值为50mph或80km/h);
[0022] v(k)为k时刻每个元胞中交通流速度；
[0023] e(k)为k时刻随机部分。
[0024]在公式（1)中当k时刻产生的随机数小于时走时停波动的概率阈值且k时刻元胞中 的交通流速度小于拥堵速度阈值时，表明k时刻元胞中交通流处于拥堵状态且拥堵状态引 发了时走时停现象。因此此时瓶颈元胞最大接受车辆数的计算公式为S(k) = S(k)+e(k);反 之，当k时刻不满足产生的随机数小于时走时停波动的概率阈值且k时刻元胞中的交通流速 度小于拥堵速度阈值这一条件时，瓶颈元胞最大接受车辆数的计算公式为S(k) =δ(1〇。
[0025] 其中δ(1〇和e(k)的计算过程如下，依据交通流数据绘制包含通行能力下降的交通 流基本图如图3所示，其中Qd为通行能力下降后流量，dd为通行能力下降发生时刻的密度阈 值。基于如图3所示的交通流基本图得到δ(1〇的计算公式如下：
[0026] (2)
[0027] 其中，
[0028] δ(1〇为元胞单位时间内能够接受从上游元胞流入的最大车辆数(也称为接受能 力）；
[0029] d(k)为k时刻元胞中交通流密度；
[0030] w为运动波传播速度；
[0031] djam为阻塞密度；
[0032] η为车道数；
[0033] Qvsl为某限速值下最大流量；
[0034] (T jam为瓶颈元胞的拥堵密度。
[0035] 在公式(2)的基础上增加随机部分ε(1〇用于在元胞传输模型中产生交通流波动， 得到e(k)的计算公式如下：
[0036]
(5)
[0037] 其中，
[0038] e(k)为k时刻产生的随机噪声；
[0039] P为决定时走时停交通流波动幅度的参数。
[0040] 第三步是在某时刻k，依据各元胞内交通流状态及如图4所示的可变限速控制下交 通流基本图的左枝，计算每个元胞内允许流出的最大车辆数 〇(k)，图4中V%为限速值，该限 速值下最大流量为Qvsl，对应的密度为dvSL。。
[0041] 改进后的元胞传输模型在计算每个元胞内允许流出的最大车辆数时考虑了驾驶 员对可变限速遵从情况的仿真。具体步骤为，基于如图3所示的交通流基本图的左枝计算单 个元胞单位时间内为下游元胞提供的最大车辆数(也称为发送能力)〇(k)，计算公式如下：
[0042] o(k) =min{V/sL(k) · d(k) · n,Qvsl} (4)
[0043] 其中，
[0044] V^Uk)为限速值为VSL(k)时车辆平均速度，VSUkhVsUlO+Vc^VsUk)为k时刻限 速值，Vo为超速幅度；
[0045] d(k)为k时刻元胞中交通流密度；
[0046] η为车道数；
[0047] Qvsl为某限速值下最大流量。
[0048]第四步为根据相邻元胞内流入流出车辆数及车辆流向，计算各元胞内单位时间的 交通流量值。依据元胞所处的位置不同交通流量的计算方法不同。具体步骤为首先需要判 断当前元胞是主线元胞、入口匝道元胞还是出口匝道元胞，三种情况下的车辆传递示意图 如图5所示。
[0049 ] 1)主线区域，元胞i+1内k时刻的流量qi (k)为上游i元胞能够提供的车辆数〇i (k)和 元胞i+1能够接受的车辆数3i+i(k)的较小值，具体计算公式为
[0050] qi(k) =min{〇i(k)，3i+i(k)} (5)
[0051] 2)入口匝道与主线交汇区域，汇入区之后主线元胞i + 1中的车辆来源于主线上游 元胞i以及入口阻道元胞j。当主线元胞i+Ι在k时刻接受能力大于上游主线元胞i以及阻道 元胞j发送能力之和，即3^1(1〇多〇 1(1〇+〇」(1〇时，主线及阻道元胞流量的计算公式分别为91 (k) =〇i(k)和qj(k) =min{〇j-：L(k)，δ」(1〇};反之，SP5i+1(k)〈〇i(k)+〇j(k)时，主线及阻道元胞 流量的计算公式分别为qi(k)=max{5 i+1(k)-〇j(k) A+1(k) · (1-ρ〇η)}和qj(k)=min{〇j(k)， δ1+1(1〇 ·ρ〇η}，其中？。"为入口匝道汇入主线的最大率，计算方法为入口匝道通行能力与主 线通行能力之比。
[0052] 3)出口匝道与主线交汇区域，主线元胞i中车辆流入主线下游元胞i+Ι和出口匝道 元胞j + Ι中。当出口匝道元胞j + Ι在k时刻接受能力大于离开主线车辆数，即h+1(k)>〇1 (k) · p〇ff时（其中，p〇ff为主线元胞i中离开主线的比例），主线及阻道元胞流量的计算公式 分别为qi(k) = 〇i-：L(k) · (Ι-poff)和qj(k) = 〇i(k) · p〇ff;反之，即5j+1(k)〈〇i(k) · p〇ff时，主 线及阻道元胞流量的计算公式分别为9办）=111；[11{51+1(10,5」+1(10*(11。￡|0/^過}和9」(10 = 5j+i(k) 0
[0053] 第五步为从起点元胞向终点元胞依次传递元胞内交通流量，更新各元胞内交通密 度。具体步骤为确定了 k时刻主线与匝道元胞内流量后，对k+Ι时刻各元胞内车辆数进行更 新，按如下公式计算各元胞内交通流密度：
[0054] di(k+l) = di(k)+A T/(Li · m) · (qi-i(k)_qi(k)) (6)
[0055] 其中，
[0056] △ T为仿真时间步长，等于车辆以自由流速度行驶过单位元胞的时间；
[0057] U为单位元胞的长度。
[0058] 第六步为根据交通流基本图中流量-密度-速度关系同时考虑限速值，计算各元胞 内交通流速度，计算公式如下：
[0059] (7)
[0060] 其中，
[0061] VF为可变限速控制下自由流速度；
[0062] dm为当前限速值下最大流量对应的密度值。
[0063] 下面结合附图对发明的仿真模型改进方法进行了实例演示：
[0064] 假设某2英里长如图2所示几何特征的路段，路段内车道数为4车道，路段内包含一 个入口匝道和一个出口匝道，及3个交通流检测器。在元胞传输模型中对上述参数进行设 置，每个元胞长度设置为0.1英里，出入口匝道元胞位置和真实情况一致。数据检测的时间 间隔T为30秒，仿真时间步长设为6秒。
[0065] 假定仿真运行到时刻k，假设瓶颈元胞内时走时停运动波的概率阈值Φ=0.5Λ时 亥|J产生的随机数ζ (k) = 0.4，k时刻每个元胞中交通流速度v(k) = 9 2km/h，判断交通流是否 拥堵的速度阈值VTH = 80km/h。因此参考公式（1)可知不符合时走时停运动波产生条件，故瓶 颈元胞最大接受车辆数计算公式为S(k) = S(k)。依据实测历史交通流数据绘制该路段交通 流量密度图可知Qd = 5680veh/h，对应的dd为56veh/km/ln，w = 12km/h，瓶颈元胞的拥堵密度 Y jam=110veh/km/ln，k时刻的交通流密度d(k)为62veh/km/ln。由于d(k) >dd，依据公式 (2)可知3(k)=w· (cT jam-d(k)) · n = 12 · (110-62) · 4 = 2304veh〇
[0066] 依据如图4所示的可变限速控制下交通流基本图的左枝，假设限速值VSL = 80km/h， 该限速值下最大流量以51^ = 798〇￥611/11，对应的密度(^[ = 46￥611/1011/111，山￡1111=120￥611/1011/ ln，超速幅度Vci = 5km/h，将上述参数值带入公式(4)得0(k) sL(k) · d(k) · n，QvsL} =min{(80+5) · 62 · 4,7980} =min{21080,7980} =7980veh0
[0067] 假设当前元胞为主线区域第3个元胞，则由公式(5)可知当前元胞k时刻的流量q2 (k) =min{〇2(k)，δ3(1〇 } =min{7980,2304} =2304veh/h，同理计算第4个元胞，则q3(k)= min{〇3(k)，δ4(1〇}，假设q3(k) = 2103veh/h，依据更新规则从起点元胞向终点元胞对每个元 胞依次更新其k时刻的流量，则k+Ι时刻主线第4个元胞内交通流密度d 3(k+1) = d3(k) + Δ T/ (L3 · n3) · (q2(k)-q3(k)) = 62+6/3600/(0.16 · 4) · (2304-2103) =62.52veh/km/ln〇d3(k+ 1) = 62.52veh/km/ln>46veh/km/ln = dvsL，因此由公式(7)可知主线第4个元胞的交通流速 度计算结果如下：
[0068] v(k+l) = (djam(k+l)-d(k+l)) · w/d(k+l)
[0069] =(120-62.52) · 12/62.52 = llkm/h
[0070] 本实例中仅给出了一个元胞在一个仿真时刻的交通流量、密度和速度的仿真计算 过程，其他元胞在其他时刻的仿真计算流程与此类似。
【主权项】
1. 一种针对可变限速控制的元胞传输仿真模型改进方法，其特征是包括W下步骤： 1) 确定各时间段主线起点、入口应道W及出口应道上交通需求，具体步骤包括： 101) 选定一段长度的快速路段并记录其主线路段车道数，在快速道路路段合理位置设 置交通流检测器，选定数据检测的时间间隔T及仿真时间步长ΔΤ; 102) 采用步骤101)中设置的交通流检测器W时间间隔T周期性采集主线起点、入口应 道和出口应道上的平均交通量、密度及速度； 103) 依据步骤102)中检测到的历史交通流数据求得自由流车流的平均速度，元胞长度 由每时间间隔内车辆W自由流速度行驶的距离确定； 2) 在当前仿真时间步中，依据各元胞内交通流状态及考虑通行能力下降的交通流基本 图右枝，计算每个元胞内允许流入的最大车辆数S化），具体步骤包括： 201) 确定瓶颈元胞内时走时停运动波的概率阔值Φ，在k时刻产生介于0到1之间包含0 和1的随机数ζ化），依据不同交通流状态选择相应的瓶颈元胞最大接受车辆数δ化)的计算 方法;当k时刻产生的随机数小于时走时停波动的概率阔值且k时刻元胞中的交通流速度小 于拥堵速度阔值时，表明k时刻元胞中交通流处于拥堵状态且拥堵状态引发了时走时停现 象，此时瓶颈元胞最大接受车辆数的计算公式为δ化)=δ化)+e化）；反之，当k时刻不满足产 生的随机数小于时走时停波动的概率阔值且k时刻元胞中的交通流速度小于拥堵速度阔值 运一条件时，瓶颈元胞最大接受车辆数的计算公式为S化)=δ化）； 202) 依据交通流数据绘制包含通行能力下降的交通流基本图，找到当前快速路段上通 行能力下降后流量Qd和通行能力下降发生时刻的密度阔值dd，则δ化)的计算公式如下：其中， S化)为元胞单位时间内能够接受从上游元胞流入的最大车辆数(也称为接受能力）； cKk)为k时刻元胞中交通流密度； W为运动波传播速度； dja^为阻塞密度； η为车道数； Qvsl为某限速值下最大流量； d ' jam为瓶颈元胞的拥堵密度； 203) 当元胞中出现拥堵且拥堵状态引发了时走时停现象时，随机部分ε化)的计算方法 如下式： e{k)=(p-0{k)-c{k), where cV(7v)G[-l, 1] 其中， ε化)为k时刻产生的随机噪声； ?^为决定时走时停交通流波动幅度的参数； 3) 绘制可变限速控制下交通流基本图，得到限速值VsL，该限速值下最大流量Qvsl和对应 的密度dvsL的取值，依据每个元胞交通流状态将上述参数带入公式Wk)=min{r SL化）· d 化）· n，QvsL}计算每个元胞内允许流出的最大车辆数； 4) 依据步骤3)中计算得到的各元胞内允许流入和流出的最大车辆数，同时基于每个元 胞所处在快速道路上的位置确定元胞中的车辆流向，计算各元胞内单位时间的交通流量 值。每个元胞依据其所处位置分为主线元胞、入口应道元胞和出口应道元胞Ξ种情况，具体 流量计算公式如下： 401) 主线区域，元胞i+1内k时刻的流量qi化)为上游i元胞能够提供的车辆数〇i化)和元 胞i+1能够接受的车辆数Si+i化)的较小值，具体计算公式为qi化)=min{〇i化），δ?+1化)}; 402) 入口应道与主线交汇区域，汇入区之后主线元胞i + 1中的车辆来源于主线上游元 胞iW及入口应道元胞j ;当主线元胞i + 1在k时刻接受能力大于上游主线元胞iW及应道元 胞j发送能力之和，即Sw(k)>〇i(k)+〇^k)时，主线及应道元胞流量的计算公式分别为qi 化）=〇i化)和qj化）=min{〇j-i化），Sj化）};反之，即δ?+1化)<〇i化)+〇j化)时，主线及应道元胞 流量的计算公式分别为qi化）=max{Si+i化)-Oj化），δ?+1化）·（l-p〇n)巧口qj化）=min{〇j化）， 5i+i化）· P。。}，其中P。。为入口应道汇入主线的最大率，计算方法为入口应道通行能力与主 线通行能力之比； 403) 出口应道与主线交汇区域，主线元胞i中车辆流入主线下游元胞i+1和出口应道元 胞j+1中；当出口应道元胞j+1在k时刻接受能力大于离开主线车辆数，即化）化）· Pnff时，其中，Pnff为主线元胞i中离开主线的比例，主线及应道元胞流量的计算公式分别为 qi 化）=〇i-i 化）·（l-p〇ff)和 qj 化）=〇i 化）· Poff;反之，即 Sj+i 化)<〇i 化）· poff 时，主线及应 道元胞流量的计算公式分别为qi化）=min{Si+i化），Sj+i化）·（l-p〇ff)/p〇ffWPq^k) = Sj+i 化)； 5) 在快速道路路段上从起点元胞向终点元胞依次传递元胞内交通流量，更新各元胞内 交通密度。将所得交通流密度与当前限速值下最大流量对应的密度值进行比较，选择合适 公式计算各元胞内交通流速度； 6) 若未到达仿真结束时刻，则进入下一个仿真时间步，回到步骤2);若到达仿真结束时 亥IJ，则结束仿真，输出所需交通流数据。2. 根据权利要求1所述的一种针对可变限速控制的元胞传输仿真模型改进方法，其特 征是所述步骤2)中，依据包含通行能力下降的交通流基本图右枝，找到当前快速路段上通 行能力下降后流量Qd和通行能力下降发生时刻的密度阔值dd;具体而言是绘制无可变限速 控制下交通流密度流量关系散点图，可W发现在某一密度范围内，交通流量值会发生突降， W该密度值为分界点将该散点图划分为左枝和右枝并分别进行线性最优拟合，依据线性拟 合的结果可W分别得到通行能力下降前后的流量及通行能力下降发生时刻对应的密度阔 值，运两个参数值可用于计算各元胞最大接受车辆数。3. 根据权利要求1所述的一种针对可变限速控制的元胞传输仿真模型改进方法，其特 征是所述步骤3)中，V^sL化)为限速值为VsL化)时车辆平均速度，V^sL化）的计算公式为rsL 化)=VSL化)+Vo,其中，VSL化)为k时刻限速值，Vo为超速幅度;驾驶员对可变限速控制的不同 遵从度就是通过调整超速幅度大小来体现的。4. 根据权利要求1所述的一种针对可变限速控制的元胞传输仿真模型改进方法，其特 征是所述步骤5)中，确定了k时刻主线与应道元胞内流量后，对k+1时刻各元胞内车辆数进 行更新时，各元胞内交通流密度的计算公式为di化+l) = di化)+ΔΤ/化i ·η〇 · (qi-i化)-qi 化)），其中，A Τ为仿真时间步长，Li为单位元胞的长度。5.根据权利要求1所述的一种针对可变限速控制的元胞传输仿真模型改进方法，其特 征是所述步骤5)中，根据交通流基本图中流量-密度-速度关系及当前限速值，计算各元胞 内交通流速度的公式如下：其中， Vf为可变限速控制下自由流速度； dvsL为当前限速值下最大流量对应的密度值。
【文档编号】G06F17/50GK105975705SQ201610318492
【公开日】2016年9月28日
【申请日】2016年5月13日
【发明人】李志斌, 刘攀, 王炜, 徐铖铖
【申请人】东南大学