CN102419905A

CN102419905A - 基于交通波理论的高速公路事故的交通影响范围确定方法

Info

Publication number: CN102419905A
Application number: CN2011102312926A
Authority: CN
Inventors: 余贵珍; 王云鹏; 刘玉敏; 鲁光泉; 田大新
Original assignee: Beihang University
Current assignee: Beihang University
Priority date: 2011-08-12
Filing date: 2011-08-12
Publication date: 2012-04-18
Anticipated expiration: 2031-08-12
Also published as: CN102419905B

Abstract

本发明公开了一种基于交通波理论的高速公路事故影响范围的确定方法，用于在无出入口匝道的高速公路上发生交通事故后，确定交通事故的时空影响范围。本发明方法首先确定高速公路在不同交通状况下的车流量和车流密度，绘制高速公路流量一密度图，然后运用交通波理论，分析从事故开始到事故影响消除过程中交通波的聚积和消散过程，计算各阶段的交通波速度，绘制交通事故时空演化图，最后确定事故在各阶段的影响时间范围和空间范围。本发明方法简便、稳定性好、可靠性高，可实时、准确预测交通事故影响范围，向高速公路管理部门提供稳定、及时、可靠的交通事故信息，从而提高交通事故快速处置措施的实施效率。

Description

基于交通波理论的高速公路事故的交通影响范围确定方法

技术领域

本发明属于交通信息及控制领域，涉及一种智能交通技术领域的方法，尤其是一种基于交通波理论的高速公路事故的交通影响范围确定方法。

背景技术

高速公路上的交通事故在给人们带来生命财产损失的同时，也会引起大范围的交通拥堵，增加车辆油耗和废气排放，带来能源消耗和环境污染问题。一旦高速公路上发生交通事故，部分道路就会被占用或者封闭，从而影响道路交通的正常运行，事故发生地点因此便成为交通瓶颈，道路的通行能力降低使之无法满足交通需求，进而导致交通拥堵，增加二次事故发生的可能性。如不及时有效地评估突发事件对高速公路道路交通运营安全的影响，并让高速公路上的驾驶员获悉影响的时间，这种影响往往会从点扩展到面，波及临近的高速公路，进而影响到整个路网的运输效率和营运安全，严重导致大面积路网拥堵甚至瘫痪。事故影响分析是高速公路系统分析中的一个重要方面，它是实行有效的交通控制、事故快速处理、交通诱导及公路规划等的基础。对事故时空影响范围进行预算估计，采取有针对性的交通事故处置措施和诱导、控制、疏导等交通管制措施，并通过高速公路信息发布系统告知高速公路上的行驶者，能避免交通事故造成的交通拥堵现象进一步恶化，有效减轻驾驶员的不舒适度，降低油耗，减小行驶时间损失，为高速公路交通设施与控制方案改善提供重要的定量依据。

但是，现有的交通事故影响研究存在一些不足，还无法有效应用于高速公路，主要包括几下几方面：(1)目前对事故地点引发的交通影响进行分析的研究较少，研究重点主要集中在对事故的特点、成因分析和预防对策上；(2)现阶段的事故影响范围研究主要集中在城市道路中交叉口处的交通事故研究上。由于在高速公路上的车辆运行速度明显高于交叉口处的车辆，并且高速公路上的运行成分单一，交通流特性十分明显，所以适合交叉口处的事故影响范围估算方法并不适合高速公路上的事故影响估算；(3)目前对高速公路交通事故的影响分析主要是确定性排队分析方法，这种方法虽然分析简单，使用方便，但是忽略了交通波的扩散速度对交通流各参数的影响，计算结果与实际偏差较大；(4)由于交通检测技术的局限性，在国际上尚无直接检测车辆停车、排队和延误的方法，传统的检测方法多利用视频检测技术对道路上交通流量、车流速度和排队长度进行事后处置，对道路基础设施依赖性高，无法实时准确地利用检测技术直接确定交通事故的时空影响范围。

发明内容

本发明的目的在于针对现有对高速公路事故的交通影响范围研究中存在的不足及实际需要，提出一种基于交通波理论的高速公路事故的交通影响范围确定方法。

本发明的一种基于交通波理论的高速公路事故的交通影响范围确定方法，包括如下步骤：

步骤一、确定高速公路上的交通流量和车流密度，绘制两相邻车道的交通流量-密度图；

步骤二、分析交通事故过程各阶段的交通波，并确定各阶段的交通波的速度；

步骤三、根据各阶段的交通波速度绘制出交通事故的时空演化图；

步骤四、确定两股交通波相遇时的事故影响范围、相遇时间点和事故影响消散的时间点，然后确定各时间阶段交通事故的空间影响范围、最长影响距离和最长影响时间。

本发明的优点与积极效果在于：

(1)与现有技术相比，本发明所需交通参量简单、易于获取，对高速公路基础设施条件依赖性低，避开了一般道路交通检测手段不完善、可靠性低等缺陷，为准确计算事故拥堵排队长度提供了可靠数据；

(2)该方法利用交通波理论对交通事故发生点上游车辆排队的聚积和消散过程进行分析，确定交通事故的实时空间影响范围和最大排队长度的计算方法，具有计算简便，运算速度快，可靠性高等优点，为交通事故的时空影响范围的确定奠定了坚实基础；

(3)通过对高速公路上交通事故发生点上游车辆排队的聚积和消散过程的分析，对高速公路上事故的时空影响范围进行了快速、准确地估算，为事故的安全态势分析和应急处置策略提供有效依据，为高速公路的交通控制提供了可靠信息材料。

附图说明

图1为本发明的事故交通影响范围确定方法的整体流程图；

图2为本发明实施例的交通流量-密度图；

图3为本发明实施例中交通事故发生初期的产生的交通波的示意图；

图4为本发明实施例中lane 1被清理后的初期产生的交通波示意图；

图5为本发明实施例中lane 1被清理后的后期产生的交通波示意图；

图6为本发明实施例中两车道清理完毕后的初期产生的交通波示意图；

图7为本发明实施例中两车道清理完毕后的后期产生的交通波示意图；

图8为本发明实施例在交通流量-密度图中事故发生到影响消除各阶段交通波示意图；

图9为本发明实施例的交通事故时空演化图。

具体实施方式

下面将结合附图和实施例对本发明作进一步的详细说明。

如图1所示，一种基于交通波理论的高速公路交通事故的交通影响范围确定方法，包括如下步骤：

第一步：确定高速公路上的交通流量和车流密度。

根据高速公路设计资料，确定高速公路的通行能力、拥挤密度等固定交通参量。利用高速公路道路上布设的固定交通信息检测器例如地感线圈检测器，检测高速公路上车辆经过时线圈检测器磁场变化，检测器据此计算出交通流量、车辆速度、车辆时间占有率和长度等交通参数，获取道路交通信息，对所获取交通信息数据进行处理，剔除散点后，计算出高速公路上不同交通状况下的交通流量和车流密度，得出正常交通状况下各车道上的交通流量-密度图。具体包括如下步骤：

(1)获取道路交通信息：

根据高速公路设计资料，确定高速公路的最大通行能力q_max、允许的最大行驶速度v_max、拥挤情况下的车流密度k_jam等。

采用地感线圈检测器获取道路交通信息。设从事故发生地点上游L范围内没有出口和进口，并且认为L足够长，车队没有延伸过这一路段，还设停车线下游的交通流顺畅，没有阻塞现象。上游车道的交通流量稳定，恒为一定值q₁。在高速公路上每隔距离S安置检测器，根据检测器所检测的数据，计算获得如下数据：道路交通流量、车辆速度、车辆时间占有率和长度等。对同一路段进行多时段的反复多次数据采集，统计分析计算当前道路上的交通流量和车流密度，并记录发生交通事故、天气不良、道路施工等不良交通环境下的交通事件。

(2)处理数据，绘制图形：

剔除发生交通事故或者灾害性天气等不良交通环境下路段上交通流量和车辆密度数据。选取两相邻车道作为研究对象，计算出无交通事件发生情况下的平均交通流量和相应的车辆密度，并绘制交通流量-密度图，如图2所示。

图2为本发明实施例中某段高速公路的交通流量-密度图。横坐标表示车辆密度，纵坐标表示交通流量。

分别表示单条车道和两条车道的最大交通流量，

分别表示单条车道和两条车道在最大通行能力时的车流密度，

分别表示单条车道和两条车道的拥挤车流密度。

第二步：分析交通事故过程各阶段的交通波，并确定各阶段的交通波的速度。

本发明运用交通波理论，分析从事故开始到事故影响消除过程中交通波的积聚和消散过程，并计算各阶段的交通波速度。

交通波理论是交通研究中进行分析、设计、仿真和决策的重要手段，它描述了同向运动不同状态的两股交通流相遇时状态的转移过程，当道路交通因为交通状况的改变而发生车辆密度的变化时，不同状态的交通流如q₁，k₁，v与q₂，k₂，v₂，相遇并产生传播速度为u_w＝(q₂-q₁)/(k₂-k₁)的交通波，进而定量分析车流的积聚和消散过程。其中，q₁、q₂分别表示两股交通流的交通流量，k₁、k₂分别表示两股交通流的车辆密度，v₁、v₂分别表示两股交通流的速度。

下面给出一个实施例进行具体说明。

本发明实施例，在t_A时刻时，高速公路某横断面处发生一起交通事故，事故发生后将该方向上的所有交通即刻封锁并进行事故处理；在时间t_B时刻，一条车道完成清理工作，该车道上的车辆恢复行驶；在时间t_C时刻，另外一条车道也清理完毕，道路全面通车，所有车辆恢复正常行驶。

(1)首先用交通波理论分析自事故发生到事故影响消散的过程中各阶段的交通流量及车辆密度变化。

可将交通波传播过程分为以下几个阶段：

1)当t＝t_A时，在高速公路某横断面上发生交通事故，该方向上的所有交通即时被封锁，交通状况如图3所示，发生事故前后的交通状况分别为自由流状况下的交通流q₁，k₁和道路阻塞状况下的交通流q₂，k₂相遇，并产生波速为W₂₁的交通波。

2)当t＝t_B时，第一条车道lane 1被清理，部分车辆恢复行驶，第一条车道被清理后事故地点的交通流q₃，k₃与阻塞状况下的交通流q₂，k₂相遇产生波速为W₃₂的交通波。第一条车道恢复初期，道路上有两股交通波，波速分别为W₂₁和W₃₂，如图4所示，此状态持续到波速为W₂₁和W₃₂的两股交通波相遇时刻为止。

3)第一条车道交通恢复后期，即从波速为W₂₁和W₃₂的两股交通波相遇时刻开始到两车道清理前的时间，波速为W₂₁和W₃₂的两股交通波相遇，产生波速为W₃₁的交通波，如图5所示。

4)当t＝t_C时，全部道路清理完毕，两车道lane 1和lane 2均被清理后的交通流q₄，k₄遇第一条车道被清理后的交通流q₃，k₃，产生波速为W₄₃的交通波。初期阶段，道路上有波速为W₃₁和W₄₃的两股交通波，此状态持续到波速为W₄₃和W₃₁的两股交通波相遇时刻为止，如图6所示。

5)从波速为W₃₁和W₄₃的两股两波相遇时刻开始到交通恢复正常为止。交通波W₃₁和W₄₃相遇产生波速为W₄₁的交通波。如图7所示，道路上仅有波速为W₄₁的交通波，波速为W₄₁的交通波完成后，交通恢复正常。

(2)计算交通波速度：

根据检测器记录的数据点及事故发生时间t_A、第一条车道清理时间t_B、第二条车道清理时间t_C等，在交通流量-密度图中找出各交通状态下的交通流量q_i和车辆密度k_i：事故前自由流状态下事故发生地的交通流量和交通密度分别为q₁，k₁；事故发生后事故地点的交通流量和交通密度分别为q₂，k₂；第一条车道lane 1被清理后事故地点的交通流量和交通密度分别为q₃，k₃；两条车道均被清理后事故地点的交通流量和交通密度分别为q₄，k₄。进而计算出上述分析中各阶段的交通波速度W_ij。

计算公式为：

W_{ij} = \frac{q_{j} - q_{i}}{k_{j} - k_{i}}

(q_i，k_i)和(q_j，k_j)分别为t_i和t_j时刻的交通状态。本实施例在交通流量---密度图中事故发生到影响消除各阶段交通波示意图如图8所示。

图8中横坐标表示车辆密度，纵坐标表示交通流量。点1表示交通事故发生时刻的交通状态，点2表示交通事故发生时刻后道路交通全部中断的交通状态，点3表示第一条道路清理后的交通状态，点4表示全部道路清理后的交通状态，两点连线的斜率表示交通波的传播速度。

第三步：绘制交通事故时空演化图，具体方法是：作时间-距离图，以事故发生点A为起始点，各时间段中的交通波速度W_ij为线段斜率，结合事故发生时间、第一车道清理完成时间、第二车道清理完成时间，绘制事故发生到事故影响消散过程中车队队尾的时间-距离图，即为交通事故时空演化图，如图9所示。

根据交通波的传播特性，可以通过图9所示的交通事故时空演化图描述事故发生后的车辆排队过程。事故发生在t_A时刻，距离为L的横断面上，两条车道全部封堵，产生波速为W₂₁的交通波；在t_B时刻，其中一条道路清理，产生波速为W₃₂的交通波，并在t_D时刻与波速为W₂₁的交通波相遇，形成新的波速为W₃₁的交通波；在t_C时刻，两条车道全部清理，产生波速为W₄₃的交通波；在t_E时刻波速为W₄₃的交通波与波速为W₃₁的交通波相遇，形成新的波速为W₄₁的交通波；在t_F时刻，事故的影响全部消除，交通恢复正常。

第四步：确定交通事故各阶段时空影响范围。

根据绘制的交通事故的时空演化图，建立两车道的事故交通影响范围确定模型，通过计算两股交通波相遇时的事故影响范围和相遇时间点，利用解析几何方法，进而计算各时段事故的时间影响范围和空间影响范围。具体包括如下步骤：

(1)确定两股交通波相遇时的事故影响范围L_DG、L_EH：

L_{DG} = \frac{W_{32} \times W_{21}}{W_{32} - W_{21}} \times (t_{B} - t_{A})

L_{EH} = \frac{W_{43}}{(W_{43} + W_{31})} [\frac{W_{21} \times (W_{32} + W_{31})}{(W_{32} - W_{21})} \times (t_{B} - t_{A}) - W_{31} \times (t_{C} - t_{B})]

如图9所示，L_DG表示波速分别为W₂₁和W₃₂的两股交通波相遇时的事故影响范围；L_EH表示波速为W₃₁和W₄₃的两股交通波相遇时的事故影响范围。

(2)确定两波相遇时的时间节点t_D、t_E，事故影响消散的时间点t_F：

t_{D} = t_{B} + \frac{L_{DG}}{W_{32}}

t_{E} = t_{C} + \frac{L_{EH}}{W_{43}}

t_{F} = t_{C} + (\frac{1}{W_{43}} + \frac{1}{W_{41}}) \times L_{EH}

(3)确定各时段事故的空间影响范围L(t)：

L (t) = \{\begin{matrix} W_{21} \times (t - t_{A}) & t_{A} \leq t \leq t_{D} \\ \frac{(W_{21} - W_{31}) \times W_{32}}{W_{32} - W_{21}} \times (t_{B} - t_{A}) + W_{31} \times (t - t_{A}) & t_{D} < t \leq t_{E} \\ \frac{W_{43}}{W_{43} - W_{31}} [(t_{C} - t_{B}) \times W_{31} - \frac{W_{32} - W_{31}}{W_{21} - W_{32}} \times (t_{B} - t_{A}) \times W_{21}] + W_{41} \times t & t_{E} < t \leq t_{F} \end{matrix}

(4)确定事故最长影响距离L_max和事故最长影响时间T_max：

L_{\max} = \frac{W_{32} \times W_{21}}{W_{32} - W_{21}} \times (t_{B} - t_{A})

T_max＝t_F-t_A+α

其中，α为事故影响时间修正系数，该参数值的大小可根据交通事故影响时间的实际观测值和模型理论值的比较结果进行标定，一般α取值范围为[0.04599，0.092]。

实现上述方法仅需要在高速公路断面上安置检测器。检测器利用感应式信号控制系统的现有设施，系统性采集高精度的实时交通数据，并把实时采集的车流信息发送到后台服务器数据库中，进行交通道路交通路段交通流量的事实状态和拥堵实时提取、排队长度计算、交通诱导发布和交通事件的自动检测。

通过本方法及其实施装置，在不增加执行机构的基础上，只需安装较少检测器即可实现对高速公路上交通事故的时空影响范围进行分析计算，帮助驾驶员在拥堵状况下及时采取必要措施，合理选择行车路线，避开拥堵，减少了交通延误、堵塞和事故，从而可大幅度地提高路网的通行能力，使路网畅通，降低能源消耗，减少排放量，改善环境污染，提高汽车运输生产率和经济效益。

Claims

1.一种基于交通波理论的高速公路事故影响范围的确定方法，其特征在于，所述方法包括如下步骤：

2.根据权利要求1所述的一种基于交通波理论的高速公路事故影响范围的确定方法，其特征在于，所述的步骤一具体包括：

步骤1.1：获取道路交通信息：根据高速公路设计资料，确定高速公路的固定交通参量，主要包括：最大通行能力、允许的最大行驶速度和拥挤情况下的车流密度；在高速公路上每隔距离S安置检测器，利用检测器对同一路段进行多时段的数据采集，根据每次采集的数据获取该路段的车辆流量和车辆速度，并记录发生在不良交通环境下的交通事件；

步骤1.2：统计不良交通环境下该路段的车辆流量和车辆密度数据，在步骤1.1中获取的数据中剔除不良交通环境下的数据，然后选取两相邻车道作为研究对象，计算出无交通事件发生情况下的平均车辆流量和相应的车辆密度，并绘制出两相邻车道的交通流量-密度图。

3.根据权利要求1所述的一种基于交通波理论的高速公路事故影响范围的确定方法，其特征在于，所述的步骤二具体包括：

步骤2.1：用交通波理论分析自事故发生到事故影响消散的过程中各阶段的车辆流量及车辆密度变化，具体各阶段为：

1)t_A时刻，在高速公路某横断面上发生交通事故，该方向上的所有交通即时被封锁，发生事故前的自由流状况下的交通流和发生事故后的道路交通全部阻塞状况下的交通流相遇，产生波速为W₂₁的交通波；

2)t_B时刻，第一条车道被清理，部分车辆恢复行驶，恢复初期，第一条车道被清理后事故地点的交通流与道路交通全部阻塞状况下的交通流相遇产生波速为W₃₂的交通波，此时道路上存在波速分别为W₂₁和W₃₂的两股交通波；

3)在波速分别为W₂₁和W₃₂的两股交通波相遇时刻开始，道路上仅有波速为W₃₁的交通波；

4)t_C时刻，两条车道都清理完毕，两车道均被清理后的交通流与第一条车道被清理后的交通流相遇，产生波速为W₄₃的交通波，此时道路上有波速为W₃₁和W₄₃的两股交通波；

5)在波速为W₃₁和W₄₃的两股交通波相遇时刻开始，道路上仅有波速为W₄₁的交通波，在波速为W₄₁的交通波完成后，交通恢复正常；

步骤2.2：在交通流量-密度图中找出各交通状态下的交通流量q_i和车辆密度k_i，然后根据下式确定各阶段的交通波的速度W_ij：

W_{ij} = \frac{q_{j} - q_{i}}{k_{j} - k_{i}} - - - (1)

其中，发生交通事故前的自由流状态下的事故发生地的交通流量和交通密度分别为q₁，k₁；交通事故发生后的道路交通全部阻塞状况下的事故地点的交通流量和交通密度分别为q₂，k₂；第一条车道被清理后的事故地点的交通流量和交通密度分别为q₃，k₃；两车道均被清理后的事故地点的交通流量和交通密度分别为q₄，k₄。

4.根据权利要求1或3所述的一种基于交通波理论的高速公路事故影响范围的确定方法，其特征在于，步骤三中所述的交通事故的时空演化图具体通过下面方法得到：以事故发生点A为起始点，各时间段中的交通波速度W_ij为线段斜率，结合事故发生时间、第一车道清理完成时间、第二车道清理完成时间，绘制事故发生到事故影响消散过程中车队队尾的时间-距离图，所得到的时间-距离图就是该交通事故的时空演化图。

5.根据权利要求1所述的一种基于交通波理论的高速公路事故影响范围的确定方法，其特征在于，所述的步骤四具体包括：

(1)确定两股交通波相遇时的事故影响范围L_DG、L_EH：

L_{DG} = \frac{W_{32} \times W_{21}}{W_{32} - W_{21}} \times (t_{B} - t_{A}) - - - (2)

L_{EH} = \frac{W_{43}}{(W_{43} + W_{31})} [\frac{W_{21} \times (W_{32} + W_{31})}{(W_{32} - W_{21})} \times (t_{B} - t_{A}) - W_{31} \times (t_{C} - t_{B})] - - - (3)

其中，W₂₁表示t_A时刻发生事故前后的交通流所产生的交通波的波速；W₃₂表示t_B时刻第一条车道被清理后事故地点的交通流与t_A时刻发生事故后道路交通全部阻塞状况下的事故地点的交通流相遇所产生的交通波的波速；W₃₁表示波速分别为W₂₁和W₃₂的两股交通波相遇所产生的交通流的波速；W₄₃表示t_C时刻两条车道都清理完毕后的交通流与第一条车道被清理后的交通流相遇所产生的交通波的波速；L_DG表示波速分别为W₂₁和W₃₂的两股交通波相遇时的事故影响范围；L_EH表示波速为W₃₁和W₄₃的两股交通波相遇时的事故影响范围；

(2)确定两股交通波相遇时的时间节点t_D、t_F和事故影响消散的时间点t_F：

t_{D} = t_{B} + \frac{L_{DG}}{W_{32}} - - - (4)

t_{E} = t_{C} + \frac{L_{EH}}{W_{43}} - - - (5)

t_{F} = t_{C} + (\frac{1}{W_{43}} + \frac{1}{W_{41}}) \times L_{EH} - - - (6)

其中，W₄₁表示波速为W₃₁和W₄₃的两股交通波相遇所产生的交通波的波速；

(3)确定各阶段交通事故的空间影响范围L(t)：

L (t) = \{\begin{matrix} W_{21} \times (t - t_{A}) & t_{A} \leq t \leq t_{D} \\ \frac{(W_{21} - W_{31}) \times W_{32}}{W_{32} - W_{21}} \times (t_{B} - t_{A}) + W_{31} \times (t - t_{A}) & t_{D} < t \leq t_{E} \\ \frac{W_{43}}{W_{43} - W_{31}} [(t_{C} - t_{B}) \times W_{31} - \frac{W_{32} - W_{31}}{W_{21} - W_{32}} \times (t_{B} - t_{A}) \times W_{21}] + W_{41} \times t & t_{E} < t \leq t_{F} \end{matrix} - - - (7)

(4)确定交通事故最长影响距离L_max和交通事故最长影响时间T_max：

L_{\max} = \frac{W_{32} \times W_{21}}{W_{32} - W_{21}} \times (t_{B} - t_{A}) - - - (8)

T_max＝t_F-t_A+α (9)

其中，α为事故影响时间修正系数。