具体实施方式
本发明提出的基于固定检测器的城市道路交通事故影响范围确定方法为:
1)用城市道路组合交通流来描述车辆到达的特征,计算车辆到达率,
2)利用设置在路段上的固定检测器,检测车辆离去率,
3)根据车辆的到达率和离去率,建立基于固定检测器的排队长度测定方法,
4)分析确定交通事故造成拥堵的扩散过程,
5)确定交通事故的影响范围。
根据以上基本思路,提出基于固定检测器的城市道路交通事故影响范围的具体确定方法:
1、车辆到达率的计算方法:
设:城市道路网完全由信号控制交叉口组成,路段上的车辆到达为组合交通流,由稳定流和离散流组成,组合交通流具有周期性特征,且组合交通流周期c1等于上游交叉口的信号周期c;稳定流由上游交叉口直行相位放行的连续行驶车辆构成,稳定流持续时间g1等于上游交叉口直行相位绿灯时间g,稳定流的到达服从泊松分布;离散流由上游交叉口左转和右转相位放行的车辆构成,离散流由上游交叉口左转和右转相位放行的车辆构成,离散流持续时间d1等于上游交叉口 周期时间减去直行绿灯相位时间的剩余时间d,且d=c-g,离散流的到达服从泊松分布;设周期开始时间为上游交叉口直行相位开始时间,则路段上车辆到达率为:
其中,λ1为稳定流的到达率,λ2为离散流的到达率,组合交通流在道路空间上的分布如附图2所示。
2、固定检测器检测车辆的离去率
在路段上设置固定检测器S,检测器检测所得的数据为交叉口车辆的离去率,检测器检测得到的交通流参数如图3所示。
3、基于固定检测器的排队长度计算方法为:
3.1)交通事故造成的拥堵-消散过程如附图4所示。t
0时刻交通事故发生,拥堵扩散持续时间包括三个部分:第一部分为车辆到达时组合交通流所在周期剩余时间Δc
1;第二部分为累计周期时间
为t
0到任意时刻t内包括的完整组合交通流周期的个数;第三部分为拥堵扩散结束时组合交通流所在周期剩余时间Δc
2,且有0<Δc
1,Δc
2<c1,拥堵扩散持续时间内的任意时刻t为:
3.2)根据拥堵扩散持续到t时刻累计到达的车辆数D(t)和累计离去的车辆数T(t)计算路段上车辆的排队长度,任意时刻t路段上车辆的排队长度L(t)为: 其中,N1和N2分别为车辆到达时组合交通流所在周期剩余时间Δc1内到达的车辆数和拥堵扩散结束时组合交通流所在周期剩余时间Δc2内到达的车辆数;设固定检测器检测的时间间隔为Δt,qi为第i个时间间隔内固定检测器检测到的车辆离去率,i为固定检 测器检测过的时间间隔的序号,n为到任意时刻t时固定检测器检测过的次数;kj为路段阻塞密度,一般取为100-150辆车/公里;
①确定N1:当Δc1>c1-g1时,N1=λ1(Δc1-c1+g1)+λ2(c1-g1),其中λ1(Δc1-c1+g1)为车辆到达时组合交通流所在周期剩余时间Δc1内以稳定流形式到达的车辆数,λ2(c1-g1)为车辆到达时组合交通流所在周期剩余时间Δc1内以离散流形式到达的车辆数;当Δc1<c1-g1时,N1=λ2Δc1,此时,车辆到达时组合交通流所在周期剩余时间Δc1内仅包括离散流,且到达的车辆数为λ2Δc1;
②确定N2:当Δc2<g1时,N2=λ1Δc2,此时,拥堵扩散结束时组合交通流所在周期剩余时间Δc2内仅包括稳定流,且到达的车辆数为λ1Δc2;当Δc2>g1时,N2=λ1g1+λ2(Δc2-g1),其中,λ1g1为拥堵扩散结束时组合交通流所在周期剩余时间Δc2内以稳定流形式到达的车辆数,λ2(Δc2-g1)为拥堵扩散结束时组合交通流所在周期剩余时间Δc2内以离散流形式到达的车辆数;
3.3)设固定检测器设置点S到事故点的距离为l
0,采样间隔为Δt,则基于固定检测器的排队长度为:
其中,t
S为拥堵排队扩散到固定检测器设置点S的时刻,排队长度L(t)计算示意图如图5所示。
4、交通事故造成拥堵的扩散过程
4.1)设交通事故发生在东西方向路段上,交通事故造成的拥堵在路段上扩散,当排队长度L(t)大于事故点到最近的上游交叉口的距离LAB时,拥堵排队会影响到距离事故点最近的上游交叉口;
4.2)当拥堵扩散到距离事故点最近的上游交叉口时该交叉口东西方向直行为绿灯相位:
如图6所示,东进口方向的车辆开始滞留在交叉口内部,西进口方向的车辆通行顺畅,南、北进口车辆在停车线处排队等待,
任意时刻t时交叉口北进口从停车线开始计算的排队长度为:
任意时刻t时交叉口东进口从停车线开始计算的排队长度为:
任意时刻t时交叉口南进口从停车线开始计算的排队长度为:
其中,NN1、NE1、NS1分别为北进口、东进口、南进口车辆到达时组合交通流所在周期剩余时间Δc1内到达的车辆数;NN2、NE2、NS2分别为北进口、东进口、南进口拥堵扩散结束时组合交通流所在周期剩余时间Δc2内到达的车辆数;qNi、qEi、qSi分别为第i个时间间隔内设置在北出口、东出口和南出口的固定检测器检测到的车辆离去率;td为拥堵排队扩散到最近的上游交叉口的时刻;lN和ls为td时刻北进口和南进口方向已有的排队长度;d1和d2分别为交叉口东西方向和南北方向的宽度;
东西方向绿灯相位结束,转为红灯相位,如图7所示,此时,交叉口处于死锁状态,西进口方向开始出现排队,任意时刻t时西进口方向排队长度为:
其中,Δg为拥堵排队扩散到上游交叉口时东西方向绿灯相位的剩余时间;NW1为西进口车辆到达时组合交通流所在周期剩余时间Δc1内到达的车辆数;NW2为西进口拥堵扩散结束时组合交通流所在周期剩余时间Δc2内到达的车辆数;qWi分别为第i个时间间隔内设置在西出口的固定检测器检测到的车辆离去率;
4.3)拥堵扩散到距离事故点最近的上游交叉口时该交叉口东西方向直行为 红灯相位:
①交通事故路段剩余通行能力大于车辆到达率时,拥堵路段排队长度逐渐减小,拥堵停止继续扩散;
②交通事故路段剩余通行能力小于车辆到达率时:
如图8所示,在上一东西方向直行绿灯相位时东进口方向车辆已经进入该交叉口内部,南、北进口方向车辆在绿灯相位也无法继续通行,排队等待,东、西进口车辆由于红灯状态而在停车线处排队等待,本周期时间结束,东西方向变为下一信号周期的绿灯相位,交叉口进入死锁状态,如图9所示,则:
任意时刻t时东进口从停车线开始计算的排队长度为:
任意时刻t时南进口从停车线开始计算的排队长度为:
任意时刻t时西进口从停车线开始计算的排队长度为:
任意时刻t时北进口从停车线开始计算的排队长度为:
其中,lW为td时刻西进口方向已有的排队长度。
5、确定交通事故的影响范围
5.1)交通事故发生时,根据交通事故类型确定交通事故清理时间,按照《道路交通事故处理办法》相关规定,将交通事故分为轻微事故、一般事故、重大事故和特大事故四类,第一类轻微交通事故的清理时间一般为20-30分钟;第二类一般事故的清理时间为30-60分钟,第三类重大事故的清理时间为60-90分钟,第四类特大事故的清理时间为90-120分钟,从而可以确定交通事故清理结束的时刻t1;
5.2)假设如下路网:由三条东西方向道路和三条南北方向道路构成的方格网式道路网,距离事故点最近的上游交叉口为O交叉口,距离O交叉口最近的东、南、西、北四个方向的交叉口分别为E1、S1、W1、N1交叉口,O交叉口与E1、S1、W1、N1四个交叉口之间的路段长度分别为lO-E1、lO-S1、lO-W1、lO-N1, 距离E1交叉口最近的北、南两个方向的交叉口分别为E2、E3交叉口,距离W1交叉口最近的北、南两个方向的交叉口分别为W2、W3交叉口,假设路网如图10所示;
5.3)根据3.3)中基于固定检测器的排队长度测定方法,计算任意时刻t时事故点的拥堵排队长度L(t),将事故点的拥堵排队长度L(t)和事故点与O交叉口之间的距离LAB进行比较,如果L(t)<LAB,则交通事故造成的拥堵排队不会扩散到距离事故点最近的上游O交叉口,事故影响范围仅限于此交通事故发生路段,此时的交通事故影响范围即为事故点拥堵排队长度;如果L(t)>LAB,则交通事故造成的拥堵排队就会扩散到O交叉口,转5.4);
5.4)根据拥堵排队扩散到O交叉口时该交叉口东西方向直行为绿灯相位还是红灯相位,采用4.2)或4.3)中拥堵排队长度测定方法,计算O交叉口东、南、西、北四个方向进口的排队长度LEO(t)、LSO(t)、LWO(t)、LNO(t),将其和O交叉口与E1、S1、W1、N1四个交叉口之间的路段长度分别进行比较,如果LEO(t)<lO-E1,则O交叉口东进口方向的拥堵排队不会扩散到E1交叉口;如果LEO(t)>lO-E1,则O交叉口东进口方向的拥堵排队会扩散到E1交叉口;南、西、北三个方向的确定方法与此相同;
5.5)重复以上计算判断过程,直至t1时刻,交通事故现场清理完毕,交通事故发生路段通行能力恢复正常,停止以上计算判断过程,此时,确定交通事故影响范围为拥堵排队扩散到的所有交叉口及其各方向排队共同构成的区域范围。
以图10所示路网,运用VISSIM仿真软件模拟交通事故造成拥堵的扩散过程。
示例:一个完全由信号控制交叉口构成的标准方格路网。交叉口基础数据如表1所示:
表1交叉口基础数据表
模拟交通事故在500s时刻发生在W1交叉口与O交叉口之间的路段上,事故点距上游交叉口144m,由于固定检测器距上游交叉口124m,故l0=20m。交通事故造成W1交叉口与O交叉口之间的路段通行能力下降,可用检测器实时检测。设交通事故拥堵扩散时间为20分钟,模拟数据记录时间为500s-1700s。模拟与计算的排队长度随时间变化如附图11所示。
模拟中排队车辆显示为白色,运动车辆显示为蓝色。不同时刻交通事故的模拟影响范围与计算影响范围的如表2所示。
表2事故影响范围数据表
用排队长度总和S描述交通事故影响范围,S1为模拟得到的交通事故影响范围,S2为计算得到的交通事故影响范围,相对误差σ=|S1-S2|/S1。各时刻两者的误差如表3所示:
表3模拟与计算误差表
时间(s) |
1300 |
1440 |
1560 |
1700 |
相对误差 |
13.5% |
14.9% |
13.0% |
11.1% |
基于固定检测器计算得到的交通事故影响范围与模拟结果相符,说明本发明设计的方法能较好地估计交通事故影响范围,从而为交通控制管理措施的实施提供理论依据。