卢晓珊 黄海军 尚华艳
(北京航空航天大学 经济管理学院,北京100191)(首都经济贸易大学信息学院,北京100070)
可变信息板(VMS,Variable Message Signs)是城市智能交通系统(ITS,Intelligent Transportation System)[1]的重要组成部分,通常安装在高速公路上,用来发布道路车流量、突发事故或临时管制等信息,引导驾驶员进行合理的路径选择.VMS有助于优化车流的时空分布,提高交通运输网络的运行效率[2].
仿真模拟方法被广泛应用于VMS的效果评价.利用元胞自动机,文献[3-4]在1个OD对、2条路径的简单网络中,分别研究了基于实时和预测拥挤系数的VMS反馈策略对流量分布的影响.文献[5]利用元胞传输模型(CTM,Cell Transmission Model)研究了VMS时间反馈信息对重复拥挤和非重复拥挤2类车流的影响.除了信息发布策略的影响,VMS的诱导效果还受到VMS数量、选址、驾驶员反应、多种信息装置等影响.例如,文献[6-7]研究了驾驶员对信息的服从率和路径转移问题.
但到目前为止,研究VMS选址的文献并不多.文献[8]应用CTM研究了VMS位置的优化问题,但是没有考虑多种信息装置同时存在的情况.实际路网中一定比例的车辆还配备先进的出行者信息系统(ATIS,Advanced Traveler Information System).ATIS与VMS的信息内容不尽相同,拥有ATIS的驾驶员通常能够获得更多实时道路信息.当ATIS与VMS并存时,出行者的路径选择行为受到不同信息诱导装置的影响[9-10],相关的研究非常少.最近,文献[11]通过建立优化模型,研究了随机事故和ATIS环境下的VMS选址问题.但以上研究假设ATIS仅在出发时刻提供当前最短路径信息,与智能交通系统的实时性要求有一定差距,因此有必要在ATIS提供实时交通信息的环境下,进一步研究VMS的最优选址问题.
路段传输模型(LTM,Link Transmission Model)由 Yperman于2007 年提出[12].该模型假设路段内的车辆运动满足运动波理论,能够较好地模拟激波、排队形成以及排队消散.与传统的元胞自动机和元胞传输模型(CTM)相比,LTM模型可以将整个路段看成一个元胞,计算速度更快,能够应用到较大规模网络上,且便于实现所有节点的动态流量加载.
本文将出行者分为有ATIS装置的驾驶员和无ATIS装置的驾驶员2类.ATIS提供路段的实时路况(包括事故信息)和实时预测行驶时间,VMS提供有效路段上的事故信息.以系统总拥挤延迟最小为目标,应用LTM模型,确定不同ATIS渗透率下的VMS最优位置.本研究对于多种交通信息共同作用下的驾驶员路径选择行为和VMS选址具有启发价值,所建立的LTM模型也可以扩展到较大规模的路网上.
考虑1个OD对、2条路径构成的简单路网,如图1所示,路段 L1和L2构成路径R1,路段 L1和L3构成路径R2.从起点1出发的车辆可以选择路径R1或R2到达终点3.在某个时刻路段L2上发生交通事故.将研究时域[0,T]划分为K个时段,每个时段的长度为 Δt,即 T=KΔt,并假设 T足够大,使得所有车辆均能驶离网络.假设如下:
1)根据有无ATIS装置将驾驶员分为2类.假设ATIS渗透率为α,即安装ATIS的车辆比例为α,未安装ATIS的车辆比例为1-α.路段L2和L3上的实时路况(包括事故信息)和预测行驶时间通过ATIS提供给驾驶员参考.
2)路段L1上某处设置VMS信息板,能够发布下游路段L2和L3的交通事故等路况信息.未安装ATIS装置的驾驶员都能看到并完全相信VMS板上的事故信息.
3)在未获得事故信息前,所有驾驶员按照logit模型进行路径选择,假设有ATIS装置的驾驶员对路径行驶时间的感知系数为θ1,无ATIS装置的驾驶员对路径行驶时间的感知系数为θ2,且θ1≥θ2.
4)驾驶员通过VMS板或车载ATIS装置获得事故信息后,均避开事故路段出行.安装ATIS的车辆能够实时地获得事故信息,而未安装ATIS的车辆只有到达VMS信息板后方能获得.
图1 简单交通网络示意图
因此,在事故处理时段内,经过VMS板时事故信息尚未发布的无ATIS驾驶员按照logit模型选择路径,有ATIS驾驶员和经过VMS板时获得事故信息的无ATIS驾驶员均避开事故路段出行;事故处理完毕后,经过VMS板时获知事故正在发生的无ATIS驾驶员仍然按照有事故情况来选择路径,有ATIS驾驶员和经过VMS板时获得事故解除信息的无ATIS驾驶员将按照logit模型进行路径选择.
通过感应线圈和监控设备,控制中心获得路段流入率、流出率、交通事故等路况信息.计算中心根据这些数据,可以预测每个时段进入路段的车辆行驶时间,并反馈给控制中心;控制中心将该信息发送到每一辆载有ATIS装置的车内信息板上,并动态刷新所提供的信息.有ATIS的驾驶员根据ATIS实时预测信息对路径选择做出决策.
文献[5]提出的时间反馈策略是该时刻离开路段的车辆在该路段上的行驶时间.实际上,对于此刻进入路段的车辆而言,此信息已滞后不能完全反映当前的路况.为真正预测该路段的行驶时间,根据文献[13]的方法计算每个时段每条路段上的实际行驶时间,即时段[t,t+Δt]内进入路段a的车辆在该路段上的行驶时间为
其中,ra(t)是时段[t,t+Δt]内进入路段 a 的车辆在该路段上的排队长度,ra(t)=max{ra(t-Δt)+Δt(ua(t)-(t)),0}是自由流时间;(t)是时段[t,t+Δt]内路段 a出口的流出能力;ua(t)是时段[t,t+Δt]内路段 a上的车辆流入率.
图1为网络构建LTM示意图,如图2所示.其中路段L1被拆分成L'1和L″1两条路段,路段L2被拆分成 L'2和 L″2两条路段.节点 0,1',2'和 4 为虚拟节点,1'和2'分别表示可变信息板安装位置和事故发生点.路段L0和L4是虚拟路段,L0的发送能力为单位时间进入路网的交通需求,L4的接收能力无穷大.假设时刻ts节点2'处发生了严重交通事故,事故持续时间为td,事故处理时间内路段L'2出口发送能力保持不变,路段L″2入口接收能力立即退化为零.
图2 构建的路段传输模型
类似于 Newell简化运动波理论[14],本文的LTM模型使用三角形基本图(图3).该基本图由自由流速度(v),最大通行能力(qmax)以及堵塞密度(ρjam)确定.激波向后传播速度(w)和临界密度(ρcrit)也可通过这3个参数计算得到,即 ρcrit=qmax/v和 w=qmaxv/(qmax-ρjamv).
图3 三角形基本图
1.2.1 路段模型
LTM模型包括路段发送能力和接收能力2部分[12].假设 Sa(t)表示路段 a 在时段[t,t+ Δt]内的发送能力,它由路段上游进入的车辆数和路段出口流出能力共同决定;Ra(t)表示路段a在时段[t,t+Δt]内的接收能力,它由路段下游的车辆数和路段入口流入能力共同决定.可以得到公式:
其中,N(x,t)表示到t时刻为止累计通过位置x的车辆数,La,va分别表示路段 a 的入口位置、出口位置、长度、自由流速度;(t)分别表示时段[t,t+Δt]内入口流入能力及出口流出能力.
1.2.2 节点模型
LTM包含5类节点:普通节点,起点,终点,分流节点,会聚节点.ub(t)为节点i下游路段b∈Bi在时段[t,t+Δt]内的流入率,它由上游路段的发送能力以及路段b与其它路段的连接情况共同决定.主要考虑以下3种情形,其它情形可相应转化.
1)普通节点(图4a).路段 b在时段[t,t+Δt]内的流入率为
2)发散节点(图4b).未获得事故信息时,所有驾驶员根据logit模型选择路径,路段b在时段[t,t+Δt]内的流入率(uATISb(t)和 uVMSb(t)分别表示安装和未安装ATIS装置的车辆)为
其中,τb(t)和 τc(t)分别表示时段[t,t+ Δt]内进入路段b和c的车辆在这2条路段上的预测行驶时间.
3)会聚节点(图4c).路段 b在时段[t,t+Δt]内的流入率为
图4 节点的3种类型
1.2.3 累计流量更新
考虑路网各节点的流量加载.累计流量的更新可参考文献[12].本文LTM模型(图2)只有起点1发生需求加载,故各路段的累计流入量和累计流出量可按照下列公式更新:
设路段b是路段a在路径p的下游路段,则对路径累计流量更新如下:
如果车辆不能在自由流时间内流出路段,则表明路段存在拥挤,且车辆发生了拥挤延迟[15].于是,路段a的拥挤延迟计算如下:
研究时段内的系统总拥挤延迟为
假定路段L1,L2和L3上流量-密度关系均满足图3所示的三角形基本图,模拟参数设置如下:路段L1,L2,L3分别有3条、2条、2条车道,长度分别为4km,3km,3.6km;3条道路的最大通行能力均为qmax=2340 veh/(h·ln),堵塞密度ρjam=65 veh/(km·ln),自由流速度 v=72 km/h.安装ATIS的驾驶员对行驶时间的感知系数θ1=1.0,无ATIS的驾驶员的感知系数 θ2=0.1.假设前400 s车辆的交通需求是1.95 veh/s,之后变成零.事故地点在距离路段L2起点1.4 km处,发生时间ts=300 s,持续时间td=45 s.假设ATIS和VMS信息发布均不发生时间延迟.时段长度Δt=1 s.
图5给出了不同ATIS渗透率情形下VMS的最优选址使得系统总延迟最小和相应的系统总拥挤延迟的变化.可以发现,VMS的位置变化存在一个临界值:当ATIS渗透率小于50%时,随着ATIS渗透率的增加,VMS的最优选址越来越远离路口;当ATIS渗透率大于50%时,VMS的最优选址越来越靠近路口.相应的系统总拥挤延迟先变小后变大,当ATIS渗透率为50%时系统总延迟取得最小值.下面分析VMS最优选址和相应的系统总拥挤延迟的变化原因.
图5 不同ATIS渗透率下的VMS最优选址和相应的系统总拥挤延迟
首先分析不同ATIS渗透率下VMS位置变化的原因.当ATIS渗透率为30%(小于50%)时,VMS的最优选址在距离路段L1起点2.42 km处(图5),对应的路段L2和L3上的累计车流量变化如图6a所示.可以发现,事故发生时由于安装ATIS的车辆能够获得实时事故信息,故避开事故路段L2,选择路段L3;没安装ATIS并且在事故发生前经过VMS板的车辆继续选择路段L2.事故解除后不久,有ATIS的驾驶员转而选择路段L2,而此时从VMS板上获得事故信息的无ATIS的驾驶员恰好到达路口,改选路段L3.可见,在事故发生后的一个连续时间段内,利用L2和L3这2条路段同时疏导车流,能够最小化系统总拥挤延迟.
图6 ATIS渗透率为30%时,VMS在路段L1上不同位置对路段L2和L3上累计车流量的影响
当VMS安装在ATIS渗透率为零时的最佳位置,即距离路段L1起点3.98 km处时(如图6b),事故发生后从VMS板上获得事故信息的无ATIS驾驶员立即到达路口,此时3类驾驶员都选择路段L3.
类似地,事故解除后不久所有驾驶员都选择路段L2.由于L2和L3的通行能力都小于L1,故必然会在路口处形成排队增加系统总拥挤延迟.因此,当ATIS渗透率小于50%时,随着渗透率的增加,为使从VMS板上获得事故信息的无ATIS的驾驶员延迟到达路口,VMS的最优选址应越来越远离路口.
当ATIS渗透率为80%(大于50%),VMS的最优选址距离路段L1起点3.2 km(图5),对应的路段L2和L3上累计车流量变化见图7a.
类似ATIS渗透率为30%时的情况,在事故发生后的一个连续时间段内,利用L2和L3这2条路段同时疏导车流,能够最小化系统总拥挤延迟.
当VMS安装在ATIS渗透率为50%时的最佳位置,即距离路段L1起点2.32 km处时(见图7b),虽然事故发生时情况与图7a类似,但是由于VMS距离路口较远,获得VMS板上事故信息的无ATIS驾驶员没有及时到达路口,使得当事故解除后不久,2类驾驶员都选择路段L2,在路口处形成排队.获得VMS板上事故信息的无ATIS驾驶员需等待排队消散后才能进入L3路段,增加了系统总拥挤延迟.
图7 ATIS渗透率为80%时,VMS在路段L1上不同位置对路段L2和L3上累计车流量的影响
因此,当ATIS渗透率大于50%时,随着渗透率的增加,为使从VMS板上获得事故信息的无ATIS的驾驶员及时到达路口,VMS的最优选址应越来越靠近路口.
图8显示了不同ATIS渗透率下的车辆时空演化图.
图8 不同ATIS渗透率下的车辆演化时空图
当ATIS渗透率为零时(图8a),在事故发生时段选择路段L3的车辆较多,造成大量的车辆在路段L1排队,增加了系统总拥挤延迟.
当ATIS渗透率为90%时(图8b),安装有ATIS装置的车辆过多,引发“信息过度”现象[16],大部分驾驶员在事故发生时都选择路段L3出行,也导致系统总拥挤延迟增大.
当ATIS渗透率为50%时(图8c),事故发生时段选择路段L2和L3的车辆数相当,事故结束后也能充分利用这2条路段疏导车流,故系统总拥挤延迟能够达到最小.
本文应用路段传输模型,研究了ATIS环境下的VMS选址问题.研究发现,ATIS渗透率存在一个临界值,低于该值时,随着ATIS渗透率的增加,VMS的最优选址应该远离路口;高于该值时,VMS的最优选址应该靠近路口.并且,随着ATIS渗透率的增加,系统总拥挤延迟先变小后变大,在某个ATIS渗透率下系统总拥挤延迟取得最小值.本文揭示了ATIS市场渗透率在VMS位置优化中的重要作用,具有实践指导意义.但实际工作中,还需要根据ATIS的渗透率、反馈策略[17]、路网结构[18]、事故发生的时空特征以及驾驶员的风险态度等因素,综合确定VMS的最佳位置.
References)
[1]杨东凯,吴今培,张其善.智能交通系统(ITS)的发展及其模型化研究[J].北京航空航天大学学报,2000,26(1):22 -25 Yang Dongkai,Wu Jinpei,Zhang Qishan.Evolution of intelligent transport system(ITS)and its model research[J].Journal of Beijing University of Aeronautics and Astronautics,2000,26(1):22-25(in Chinese)
[2]Chorus C G,Molin E J E,Wee B V.Use and effects of advanced traveler information services(ATIS):a review of the literature[J].Transport Reviews,2006,26(2):127 -149
[3]Wang W X,Wang B H,Zheng W C,et al.Advanced information feedback in intelligent traffic systems[J].Phys Rev E,2005,72(6):6702-6707
[4]Dong C F,Ma X,Wang G W,et al.Prediction feedback in intelligent traffic systems[J].Physica A,2009,388(21):4651 -4657
[5]尚华艳,黄海军,高自友.基于元胞传输模型的实时交通信息设计[J].北京航空航天大学学报,2008,34(2):234-238 Shang Huayan,Huang Haijun,Gao Ziyou.Design real-time traffic information by cell transmission model[J].Journal of Beijing U-niversity of Aeronautics and Astronautics,2008,34(2):234 -238(in Chinese)
[6]Chatterjee K,Hounsell N B,Firmin P E,et al.Driver response to variable message sign information in London[J].Transportation Research Part C,2002,10(2):149 -169
[7]Dia H,Panwai S.Modelling drivers’ compliance and route choice behaviour in response to travel information[J].Nonlinear Dynamics,2007,49(4):493 -509
[8]尚华艳,黄海军,高自友.基于元胞自动传输模型的可变信息标志选址问题研究[J].物理学报,2007,56(8):4342 -4347 Shang Huayan,Huang Haijun,Gao Ziyou.Locating the variable message signs by cell transmission model[J].Acta Physica Sinica,2007,56(8):4342 -4347(in Chinese)
[9]Bogers E A I,Viti F,Hoogendoorn S P.Joint modeling of advanced travel information service,habit,and learning impacts on route choice by laboratory simulator experiments[J].Transportation Research Record,2005,1926:189 -197
[10]李志纯,黄海军.先进的旅行者信息系统对出行者选择行为的影响研究[J].公路交通科技,2005,22(2):95 -99 Li Zhichun,Huang Haijun.Modeling the impacts of advanced traveler information systems on travelers’travel choice behaviors[J].Journal of Highway and Transportation Research and Development,2005,22(2):95 -99(in Chinese)
[11]Chiu Y C,Huynh N.Location configuration design for dynamic message signs under stochastic incident and ATIS scenarios[J].Transportation Research Part C,2007,15(1):33 -50
[12]Yperman I.The link transmission model for dynamic network loading[D].Leuven,Belgium:Department of Civil Engineering,Faculty of Engineering,Katholieke Universiteit,2007
[13]Huang H J,William H K L.Modeling and solving the dynamic user equilibrium route and departure time choice problem in network with queues[J].Transportation Research Part B,2002,36(3):253 -273
[14]Newell G F.A simplified theory of kinematic waves in highway traffic,Part I:General theory,Part II:Queuing at freeway bottlenecks,Part III:Multi-destination flows[J].Transportation Research Part B,1993,27(4):281 -313
[15]龙建成,高自友,赵小梅.基于路段传输模型的道路出口渠化[J].吉林大学学报:工学版,2009,39(增刊2):41 -46 Long Jiancheng,Gao Ziyou,Zhao Xiaomei.Channelization at road exit based on link transmission model[J].Journal of Jilin University:Engineering and Technology Edition,2009,39(Spplement 2):41-46(in Chinese)
[16]FernandezL J E,de Cea Ch J,Valverde G G.Effect of advanced traveler information systems and road pricing in a network with non-recurrent congestion[J].Transportation Research Part A,2009,43(5):481 -499
[17]Chen B,Xie Y B,Tong W,et al.A comprehensive study of advanced information feedbacks in real-time intelligent traffic systems[J].Physica A,2012,391(8):2730 -2739
[18]Gao S,Huang H.Real-time traveler information for optimal adaptive routing in stochastic time-dependent networks[J].Transportation Research Part C,2012,21(1):196 -213