一种区域路网交通状态宏观基本图的绘制方法与流程

本发明涉及城市路网管理领域，尤其是涉及一种区域路网交通状态宏观基本图的绘制方法。

背景技术：

城市路网作为现实存在的复杂网络，优化其单个“节点”或“边”的管理方法通常不能达到优化系统的目的，有时甚至会适得其反，采用优化更大范围思路的区域交通管理政策更容易取得较好的效果。

区域交通管理政策的实施需要相关依据作为参照，为此需要一种新的评价区域交通状态的方法。宏观基本图(macroscopicfundamentaldiagram，下简称mfd)的出现为评价区域交通状态提供了很好的技术选择。

mfd是使用区域路网的平均流量、密度绘制的交通流基本图，其思想起源于1969年godfrey对现实网络中mfd性质的探究，并由geroliminis和daganzo在2008年对稳定mfd的存在性给出了证明。mfd的实际意义是在实际交通网络可以容纳的车辆数是有限的，并且系统达到最大服务能力时的车辆数不是系统能容纳的最大车辆数，控制系统中的车辆数(平均密度)可以显著改善系统的性能，而控制标准可以通过mfd得到。

scats系统(sydneycoordinatedadaptivetrafficsystem，悉尼自适应交通控制系统)是我国使用较为广泛的城市交通信号控制系统，在城市交通管理中具有重要作用。理论上绘制区域mfd需要路段的长度及各时段的密度、流量，但装备scats系统的城市一般不具备获得大规模的地面道路密度的条件，scats系统只有在停止线才有单检测线圈，而单线圈是得不到密度的，目前，国外虽有使用scats数据做mfd的实证先例，但他们要么用了额外的数据估算得到了密度，要么脱离mfd的定义，转而从mfd的稳定存在，在局部拥挤时存在“滞后回圈”现象等mfd性质出发，做出满足这类性质的图线。

技术实现要素：

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于scats线圈数据、实时准确的区域路网交通状态宏观基本图的绘制方法。

本发明的目的可以通过以下技术方案来实现：

一种区域路网交通状态宏观基本图的绘制方法，包括以下步骤：

1)通过scats线圈获取交叉口车道的流量和饱和度数据，并进行数据清洗，用以去除异常线圈数据和异常时段数据；

2)使用scats数据计算占用率；

3)对占用率进行修正；

4)将一定时间段内的占用率和流量进行加权平均；

5)用加权平均后的占用率及流量绘制路网交通状态宏观基本图。

所述的步骤1)中，判断是否属于异常线圈数据的方法具体为：

当n/naspect＜0.9且n0/naspect＜0.01时，则判定为异常线圈数据并删除，

其中，naspect为期望数据量，t为时段总长，ti为单个线圈第i条记录的周期时长，n为单个线圈的数据量，n0为线圈流量数据中含有不大于0的数据条数。

所述的步骤1)中，判断是否属于异常时段数据的方法具体为：

当交叉口处于空闲状态，即交叉口的5min流量小于15，并且scats线圈数据出现流量为0但饱和度不为0的数据时，则被判定为异常。

所述的步骤2)中，占用率occupancy的计算式为：

其中，gds为相位时长，包括绿灯时间和绿间隔时间，ts为线圈空闲时间，nvg为车辆数，tsmf为最大流量时平均每辆车的空闲时间，ds为车道饱和度。

所述的步骤3)中，对以下三种情况的占用率需要进行修正，包括：

(1)当交叉口5min流量小于15，且h>0.3时，则有：

其中，occupancy_new为修正后的占用率，occupancy_opt为最大流量时线圈的占用率，flow为5min的折算流量，mf为最大流量，h为下游排队的时间占相位时间的比例；

(2)当交叉口5min流量大于15，且h>0.3时，则有：

其中，dsnew为修正后的车道饱和度，m为三角形基本图中的常数，取值为3～5，kp为最大流量下单车占用的时间；

(3)当ds>0.95且k>1.3时，则有

其中，k为拥挤造成的scats数据换算占用率偏离理论值的程度值。

所述的步骤5)具体包括以下步骤：

以5min的平均占用率为横坐标，平均流量为纵坐标，绘制散点图，即得到区域路网的宏观基本图。

所述的宏观基本图中出现的平均流量的峰值反映当前的道路系统能够达到的最优状态，对应的平均占用率为采取区域管制的阈值。

与现有技术相比，本发明具有以下优点：

本发明仅采用scats线圈采集到的数据绘制mfd，不需要额外的交通检测数据，数据获取更为简单、方便,mfd的形状可以反映区域路网的运行状况，其中的占用率可用于实时反映区域路网的交通拥挤情况，并且有明确的阈值用于判断是否需要采取管控措施，对于城市交通管理，具有很强的参考意义。

附图说明

图1为本发明的方法流程图。

图2为2016年6月上海市部分线圈修正前后的scats线圈占用率-流量图，其中，图(2a)为修正前的scats线圈占用率-流量图，图(2b)为修正后的scats线圈占用率-流量图。

图3为2016年5月16日至2016年5月22日周一到周日的三个区域的mfd，其中，图(31a)为静安寺区域周一的mfd，图(32a)为静安寺区域周二的mfd，图(33a)为静安寺区域周三的mfd，图(34a)为静安寺区域周四的mfd，图(35a)为静安寺区域周五的mfd，图(36a)为静安寺区域周六的mfd，图(37a)为静安寺区域周日的mfd，图(31b)为人民广场区域周一的mfd，图(32b)为人民广场区域周二的mfd，图(33b)为人民广场区域周三的mfd，图(34b)为人民广场区域周四的mfd，图(35b)为人民广场区域周五的mfd，图(36b)为人民广场区域周六的mfd，图(37b)为同济大学区域周日的mfd，图(31c)为同济大学区域周一的mfd，图(32c)为同济大学区域周二的mfd，图(33c)为同济大学区域周三的mfd，图(34c)为同济大学区域周四的mfd，图(35c)为同济大学区域周五的mfd，图(36c)为同济大学区域周六的mfd，图(37c)为同济大学区域周日的mfd。

具体实施方式

下面结合附图和具体实施例对本发明进行详细说明。

实施例

本方法利用scats数据(形式如表1所示)，提出了一种绘制区域路网交通状态宏观基本图(mfd)的方法，本发明使用scats数据得到的平均占有率和流量绘制mfd，对scats数据计算得到的占有率进行修正，可以实时反映区域路网的交通拥挤情况，并有明确的阈值用于判断是否需要采取管控措施。

表1scats数据格式

如图1所示，本发明解决其技术问题所采用的技术方案包括以下步骤：

首先，对scats数据进行数据清洗，包括异常线圈数据清洗及异常时段数据清洗，具体为：

(1)异常线圈数据清洗

删除数据缺失较多或长期没有流量的线圈。数据异常严重程度的判断基于期望数据量naspect，计算公式为：

t——研究时段总长，单位s

ti——单个线圈第i条记录的周期时长，单位s

n——单个线圈的数据量

若线圈在研究时段存在数据缺失，那么该线圈的n/naspect将小于1；若线圈长时间流量为0，设线圈流量不大于0有n0条记录，那么n0/naspect将接近于0。取n/naspect＜0.9，n0/naspect＜0.01作为异常线圈的判断标准，删除这些线圈的数据。

(2)异常时段数据清洗

对于流量为0、ds不为0的数据。考虑两种情况：第一种，该现象有可能是交叉口过饱和造成的；第二种，该现象是单纯的数据异常。为排除第一种情况，统计不同时段交叉口5min流量分布作为中间表,选择交叉口5min流量是否小于15作为阈值来分辨交叉口是否空闲，空闲时发生的流量为0、大于0的现象将被识别为异常。当出现流量为0、大于0且对应交叉口5min流量小于15的记录超过70条，

绘制mfd时跳过该时段数据。

接着，使用scats数据计算占用率，计算占用率的公式如下：

ds——scats系统记录的车道饱和度(degreeofsaturation)

gds——相位时长，单位s，包括绿灯时间和绿间隔时间(黄灯时间和全红时间之和)

ts——检测出的线圈空闲时间，单位s

nvg——检测出的车辆数

tsmf——最大流量时平均每辆车的空闲时间，可以用最大流量mf和kp值计算得到，单位s；

再者，对部分情况下的占用率进行修正，如图2所示，具体包括：

设定三类指标：

occupancy_opt——最大流量时线圈的占用率

h——下游排队的时间占相位时间的比例

k——拥挤造成的scats数据换算占用率偏离理论值的程度

tsmf——最大流量时平均每辆车的空闲时间，可以用最大流量mf和kp值计算得到，单位s

mf——最大流量，单位辆/h

flow——5min的折算流量，单位辆/h

对三种情况下的占有率进行修正，得到修正后的占用率(occupancy_new)：

1、对于当前交叉口5min流量小于15，同时h>0.3：

2、交叉口5min流量大于15，同时h>0.3：

其中m为三角形基本图中的常数，一般取3～5

将dsnew重新计算占用率，即得到修正后的占用率。

3、ds>0.95且k>1.3：

然后，将一定时间段的占用率和流量加权平均，为提高mfd的稳定性，将区域路网5min内计算得到的占用率及流量按相位时长加权平均；

最后，用加权平均后的占用率及流量绘制路网交通状态宏观基本图(mfd)。

以5min的平均占用率为横坐标，平均流量为纵坐标，绘制散点图，即得到区域路网的宏观基本图(mfd)，其中出现的流量峰值反映了当前的道路系统可以达到的最优状态，对应的平均占用率为采取区域管制的阈值。

若以时间为横坐标，平均占用率为纵坐标，绘制折线图，可以得到一天中区域路网交通状况的变化情况。

采用上述方法，分析上海三个区域，2016年5月16日至2016年5月22日(一周)的mfd如图3所示。