徐鹏 张婷
摘 要：开展区域公路网突发事件影响范围的分析和研究，对突发事件发生后应急措施的制定和效用评估具有重要的现实意义。基于交通波原理，通过对公路突发事件引起的拥堵传播过程的规律和特点进行分析，将拥堵传播过程划分为不同时区。考虑到下游分流和限速措施对不同时区拥堵传播的影响，结合数学分析的方法，设计了一个公路突发事件的时空影响范围模型。经过实际数据计算验证，该模型能较好地估算公路突發事件时空影响范围。
关键词：公路突发事件;事件影响范围;交通波;拥堵
中图分类号：U125
文献标志码：A
公路网作为我国交通运输系统中的重要分支，其承担着城市间互相协作和发挥城市经济辐射作用的重任。因此，对区域公路网突发事件影响范围进行分析和研究，在理论层面和现实层面都具有重大意义。与城市道路相比，公路网的交通流在时间上的变化不同，没有明显的早晚高峰;公路网交通流大小和交叉口密集度介于城市道路和高速公路之间，常见车辆主要是机动车，少有非机动车混行，可以从交通波理论角度构建适合于公路网的拥堵扩散分析模型。公路突发事件发生后，上游会及时采取拥堵信息发布等限速措施缓解道路拥堵。突发事件影响范围受上游缓堵作用干预，掌握突发事件本身产生的集结波、道路清障产生的消散波以及上游缓堵措施产生的缓堵波之间相互作用下公路受突发事件的影响情况，对后续交通应急组织的有效实施至关重要。
突发事件影响范围一直是国内外学者重要研究的内容。大量学者从不同的理论出发，研究出各种不同适用范围的模型。陈绍宽等[1]将时间维度引入到传统的Morans I指数中，分析了城市道路交通状态的时空变化规律。CHIEN等[2]通过交通波理论，分析了不同交通场景的道路事件的延误情况，并基于此建立了一种计算突发事件总延误时间的模型，但是计算得出的总延误与交通仿真软件仿真结果相比存在很大误差。YANG等[3]将皮尔逊相关系数应用到时间自相关分析与交通时空相关性分析，研究了能对道路状态进行较好预测并能及时应对道路突发事件的模型。王云鹏[4]通过分析高速公路发生事件后上游交通流的排队和消散过程，建立了多车道封锁情况下高速公路的事件影响时空范围模型，但是该模型没有考虑实际道路情况中上游由于拥堵信息发布产生的交通波对拥堵的缓解反应，模型计算的事件总影响时间与实际相比会偏大。王建军[5]基于流体力学车流理论，研究了在干预措施作用下的突发事件对道路交通状态的改变过程，并提出干预措施作用的3种情况的车流波波速模型以及干预措施解除后的车速变化模型。金书鑫[6]将拓扑关系理论应用到区域高速路网交通事故影响区划分中，按道路节点重要度对事故处理先后顺序进行划分从而达到提高事故处理效率的目的。何雅琴[7]借用水波原理提出了交通影响系数的概念，设计了一个基于时间和距离的交通影响系数模型，并计算出突发事件综合影响程度。HOJATI等[8]研究了道路上各种因素作用下交通事故影响的持续时间模型，为缩短事故影响时间的措施提供依据。李扬等[9]对新增车辆的通行拥堵进行预测，使用K-Medoids聚类算法将交通流运行状态划分为顺畅、阻滞、拥堵三类，引入交通流特征参数构建累积Logistic回归模型量化新增车辆对路段运行状态的影响，最后基于支持向量回归机预测新增车辆通行时间。
目前国内外面向公路网模型的研究较少，忽略了公路网的特点在突发事件引起的交通拥堵传播范围中的影响，对公路网的事件拥堵特征研究不具备良好参考性，故本文设计了一种针对公路网突发事件的时空影响范围模型。
1 基于交通波理论的双车道拥堵消散过程
当区域公路网发生突发事件时，事件发生点处的交通通行能力发生骤减而形成一个虚拟瓶颈;当事件发生点占用车道导致交通供给能力下降以至不能满足交通需求，在采取疏导措施之前，拥堵向事件所在路段上游蔓延并扩大到周围路网[10]。由于公路网没有明显的高峰时间段，故可以将交通车流看作匀速的，交通波沿道路匀速移动。从空间上看，道路出现拥堵后，突发事件影响线会随着车辆排队的增加向上游运动。交通拥堵扩散示意如图1所示。由图1可见，当拥堵的集结波从瓶颈处向上游扩散时，更多路段的通行能力和速度也随着事件影响线的运动而下降。
按交通波理论，当公路上车流因突发事件发生而引起交通密度变化时，会在车流中产生密度不同的两部分界面的波的传播[11]。设密度不同的两部分界面沿空间传播的波速为U（K1，K2），则波速的一般表达式为
U（K1，K2）=（Q2-Q1）（K2-K1）。（1）
式中：Q1，Q2分别为上游交通量、下游交通量，veh/h;K1，K2分别为上游交通密度、下游交通密度，veh/km。
拥堵传播一段时间后，由于拥堵信息传播，上游交通分流的作用或者车速控制的作用，从分流点形成一个缓堵波并向上游传播，会对事件影响线的移动产生影响。采取事件点清障措施后，道路通行能力恢复，消散波开始向上游传播，同时事件影响线将向下游移动。因此，本文研究在缓堵波的干涉下，突发事件产生的集结波和道路清障产生的消散波所确定的公路突发事件时空影响范围。
本文以公路为研究对象，以双车道发生突发事件为例，考虑到事件点车道清障需要分次处理，如果在整个双车道路段以及时间组成的时空区域内研究交通波的全程，则事件影响的空间和时间范围可以分为3个时区：（1）事件发生后，两条车道都被封锁，事件发生点上游形成拥堵并逐渐蔓延;（2）其中一条车道清障完成，该车道交通流逐渐恢复;（3）另一条车道清障完成，事件点断面逐渐恢复正常通行能力[12]。
没有缓堵波作用时的双车道拥堵影响过程如图2所示。由图2可以看出：事件发生时刻t=0，发生位置的横断面L=0（L为事故影响线距离），此时双车道都被封锁，产生从事件点出发向上游传播的波速为W12的集结波;在t1时刻，两条道路中任一条完成清障工作，其产生从事件点出发向上游传播的消散波，波速为W23;两者在t2时刻汇集产生波速为W13的新消散波，此后在两条道路上只有波速为W13的交通波;在t3时刻，另一条道路完成清障工作，其产生从事件点出发向上游传播的消散波，波速为W34;在t4时刻，波速为W34的消散波与波速为W13的交通波汇集成新的消散波，波速为W14;直到tmax时刻，道路上的交通状态完全恢复，事件影响得以消除[9]。
当事件发生后的某时刻，在上游某一点由于拥堵信息的发布使得车流量减少而产生一个向上游传播的缓堵波时，若该缓堵波和以上的集结波或消散波相遇，则会对上述的拥堵影响过程产生一种良性的干扰，加速拥堵消散过程[13]。缓堵波与其他交通波相遇的情况可以分为3种时区，分别是与W12、W13、W14相遇，这3种情况对解除事件影响的过程有所差异。
2 控制波干涉作用下时空影响范围模型
2.1 缓堵波在时区1与事件影响线相遇
假设在突发事件发生前，车道路段上的交通流是均匀的，交通状态描述为（V0、K0、Q0），表示为初始状态，其中：V0为初始车流速度;K0为初始车流密度;Q0为初始交通流量[14]。缓堵波和集结波在事发路段的传播过程如图3所示。ND段和OD段分别表示缓堵波和集结波在相遇之前的运动过程，DE段表示两波相遇之后形成新的交通波的传播过程[9]。
在tn时刻，上游距离事发点Ln处产生一缓堵波，此时Ln处下游的交通状态为（Vn、Kn、Qn），则该缓堵波的波速Uh为
Uh=Q0-QnK0-Kn，（2）
设该缓堵波向下游运动与集结波相遇的时刻为t2，则由方程
Uh×（t2-tn）+W12×t2=Ln（3）
得
t2=Ln+Uh×tnUh+W12，（4）
两波相遇点距离事件发生点的距离为Lmax （即为事件最长影响距离），之后两者形成波速为Wh的新消散波，其到达事件点的时刻为tmax（即为事件最长影响时间），则
Lmax=W12×t2，（5）
tmax=t2+LmaxWh，（6）
由以上可得事件的时空影响域为
L（t）=W12×t，0 Lmax-Wh×（t-t2），t2 2.2 缓堵波在时区2与事件影响线相遇 突发事件发生前交通状态为（V0、K0、Q0），表示为初始状态，其中：V0为初始车流速度;K0为初始车流密度;Q0为初始交通流量。缓堵波和交通波在事发路段的传播过程如图4所示。 在t1时刻，双车道其中一條车道完成清障工作;在t2时刻，波速为W12和W23的两波相遇形成波速为W13的新交通波，并在t2时刻达到最大拥堵长度Lmax。 Lmax=W23×W12W23-W12×t1，（8） t2=t1+LmaxW23，（9） 在tn时刻，上游距离事发点Ln处产生一缓堵波，波速为Un。NE段表示该缓堵波与波速为W13的交通波在相遇之前的运动过程，EF段表示该缓堵波与波速为W13的交通波在t3时刻相遇之后形成的波速为Wn的新交通波的传播过程。由方程 Un×（t3-tn）+[Lmax-W13×（t3-t2）]=Ln（10） 得 t3=Ln-Lmax+Un×tn-W13×t2Un-W13，（11） 到tmax时刻，波速为Wn的交通波传播到达事件发生点结束，tmax即为事件最长影响时间。 tmax=t3+Lmax-W13×（t3-t2）Wn，（12） 由以上可得事件的时空影响域为 L（t）=W12×t，0 Lmax-W13×（t-t2），t2 Lmax-W13×（t3-t2）-Wn×（t-t3），t3 。（13） 2.3 缓堵波在时区3与事件影响线相遇 突发事件发生前交通状态为（V0、K0、Q0），表示为初始状态，其中：V0为初始车流速度;K0为初始车流密度;Q0为初始交通流量。缓堵波和交通波在事发路段的传播过程如图5所示。 在t1时刻，双车道其中一条车道完成清障工作;在t2时刻，波速为W12和W23的两波相遇形成波速为W13的新交通波，并在t2时刻达到最大拥堵长度Lmax;在t3时刻，另外一条车道完成清障工作;在t4时刻，波速为W13和W34的两波相遇形成波速为W14新交通波。t2和t3算法与之前相同，本节不再重复计算。由方程 Lmax=W34 ×（t4 -t3 ） + W13 ×（t4 -t2 ）（14） 得 t4=Lmax+W34×t3+W13×t2W34+W13，（15） 在tn时刻，上游距离事发点Ln处产生一波速为Un的缓堵波。EP段表示该缓堵波和波速为W14的交通波相遇之前的运动过程，PF段表示该缓堵波和波速为W14交通波在t5时刻相遇之后形成的波速为Wn的新交通波的传播过程。由方程 Ln-Lmax=Un×（t5-tn）-W13×（tn-t2）-W34×（t5-t4）（16） 得 t5=Ln-Lmax+（Un+W13）×tn-W13×t2-W34×t4Un-W34，（17） 最长影响时间tmax为 tmax=t5+Ln-Un×（t5-tn）Wn，（18） 由以上可得出该情况下事件的时空影响域为 [7] 何雅琴. 基于水波原理的城市突發事件交通影响程度研究[J]. 交通运输系统工程与信息， 2017， 17（5）： [8] HOJATI A T， FERREIRA L， WASHINGTON S， et al. Hazard based models for freeway traffic incident duration[J]. Accident Analysis & Prevention， 2013， 52（12）： [9] 李扬，胡尧，商明菊，等. 新增车辆通行拥堵预测模型[J]. 贵州大学学报（自然科学版）， 2019，36 （5）： [10]AL-DEEK H， GARIB A， RADAWN A E. New method for es-timating freeway incident congestion[J]. Transportation Research Record， 1995， 1494： [11]吴正. 高速公路突发事件和干涉车流波的非线性数学模型[J]. 西安公路交通大学学报， 2001， 21（2）： [12]BRESSAN A， NGUYEN K T. Optima and equilibria for traffic flow on networks with backward propagating queues [J]. Networks & Heterogeneous Media， 2017， 10（4）： [13]SHEU J B， CHOU Y H， CHEN A. Stochastic modeling and real-time estimation of incident effects on surface street traffic congestion [J]. Intelligent Transportation Systems， 2014， 28（5）： [14]胡晓健， 王炜， 陆建. 道路交通意外事件影响范围确定方法[J]. 东南大学学报， 2007， 37（5）： [15]丛浩哲， 王俊骅， 童世鑫. 高速公路网络交通突发事件辐射范围预测模型[J]. 同济大学学报， 2012， 40（3）： [16]IRAWAN M Z，SUMI T，MUNAWAR A. Implementation of the 1997 indonesian highway capacity manual（MKJI） volume delay function[J]. Journal of the Eastern Asia Society for Transportation Studies， 2010， 8： [17]曹志远， 郭忠印， 张起森， 等. 高速公路重大交通事故时空影响范围研究[J]. 公路工程， 2011， 36（6）： （责任编辑：周晓南） The Time Space Influence Range Model of Highway Emergency under Evacuating Measures XU Peng， ZHANG Ting* （College of Civil and Transportation Engineering， Hohai University， Nanjing 210098， China） Abstract： It is of great practical significance to carry out analysis and research on the impact scope of regional road network emergencies， which has important implications for formulating and evaluating emergency measures after emergencies. Based on the principle of traffic wave， this paper divides the congestion propagation process into different time zones by analyzing the law and characteristics of the congestion propagation process caused by highway emergency. Considering the influence of downstream diversion and speed limit measures on congestion propagation in different time zones， using the method of mathematical analysis， a time and space impact range model of highway emergencies was designed. After the actual data calculation and verification， the model can better estimate the time impact and space range of highway emergencies. Key words： highway emergency event; event impact range; traffic wave; congestion
151-156.
171-181.
21-27.
30-39.
77-87.
717-748.
669-674.
934-939.
414-422.
350-360.
55-73.
