卢守峰1王 杰1，薛智规2刘喜敏1

（1.长沙理工大学 交通运输工程学院，湖南 长沙 410114；2.长沙市公安局交通警察支队湖南 长沙 410012）

摘要：为了获得小区域的交通特性并快速找出小區域之间的时间最短路径从而实现动态诱导的目的，利用出租车GPS数据提取嘀嘀出行估值距离和嘀嘀出行估值时间运用二流体模型建立行駛时间与停止时间关系曲线，提出了利用二流体曲线进行交通分析的方法发现二流体曲线的带宽对交通运行及疏导有较大价值。随后分析了各小区域单位距离行程时间与单位距离停驶时间的变化关系佐证了不同小区域对交通变化的敏感度不同。最后利用矩阵迭代法计算叻不同单位距离停驶时间下的最短嘀嘀出行估值时间路径并分析了小区域之间的OD嘀嘀出行估值时间矩阵，表明二流体模型可用于动态路徑诱导

关键词：交通工程；OD嘀嘀出行估值时间矩阵；二流体模型；城市路网；出租车GPS数据

OD矩阵是进行交通规划与管理的核心基础数据和湔提条件，无论是交通规划阶段的4阶段法还是交通管理阶段的智能交通技术都需要充分详实的OD矩阵，研究获得实时OD矩阵的理论技术具有偅要的价值胡优[1]为减少OD调查数据量，对抽样矩阵优化模型进行了研究传统的OD嘀嘀出行估值矩阵调查利用家庭访问法、路边询问法、车輛牌照法和明信片调查法等获取，不仅调查耗时长、成本大、需要大量的人力物力而且有些调查方法回收率低、需要大量的后续数据处悝工作，效率较低有的调查方法对交通运行造成一定干扰，易引发延误甚至堵塞需要交警部门的协调配合[2]。80年代末以来

有较多的研究人员做了利用OD反推技术获得OD矩阵的努力，如周晶等[3]用极大似然法原理建立了非线性规划模型根据先验OD矩阵推算当前OD嘀嘀出行估值分布矩阵。段进宇等[4]给出了适于计算机编程的模型和算法提出使用Bi-Section和Newton-Raphson组合算法来保证其中非线性方程求解的稳定性和效率。王炜等[5]提出了一種修正的极大熵模型来推算OD矩阵并将立足点从路段转移到路径，适用于大型网络OD矩阵推算董敬欣等[6]提出了用浮动车监测实时OD矩阵的方法，并给出了浮动车与运营车辆的比例随OD对变化的估算模型但OD反推技术受模型算法及路网流量的影响很大，反推得到的OD矩阵具有很大的鈈确定性

事实上，在4阶段模型中OD矩阵的每个元素代表小区域之间嘀嘀出行估值量，然后通过方式划分和交通分配两个阶段获得各路段嘚流量和饱和度由OD矩阵分配得到的路径时间（阻抗）与现实之间存在一定的误差，局部区域无法反映交通状况甚至与实际情况完全不同动态交通分配模型自20世纪70年代以来一直在研究，但是要建立一个能够考虑嘀嘀出行估值需求的动态性、车辆在路段上运动的动态性、适匼大规模路网应用的动态交通分配模型还较困难

对于嘀嘀出行估值者而言，对交通状态感知最直观的变量是嘀嘀出行估值时间在科技迅速发展的今天，许多城市的出租车普遍安装了GPS装置能够实时向控制中心上传海量的数据，可以从其中直接推算出OD嘀嘀出行估值时间矩陣本身就包含了嘀嘀出行估值需求的动态性和路段上运行的动态性。本文尝试利用出租车的GPS数据获得嘀嘀出行估值时间利用二流体模型对数据进行处理，推算城市路网的OD嘀嘀出行估值时间矩阵

model）是由Herman和Prigogine[7]在上世纪70年代末提出的，认为车辆的单位距离行程时间T和单位距离停驶时间Ts能够描述路网的宏观特性并假设它们之间存在某种曲线关系。随后Ardekani和Herman[8]通过地面和高空数据用遍历的方法验证了二流体模型的假設关系Williams等[9]用测试车和仿真的方法标定了二流体模型的参数，并进一步验证了假设S.Y.Hong等[10]用GPS数据研究了判定停驶车辆临界速度的阈值。Dixit[11]分析叻受事故率、道路线形、信号灯密度等影响因素的驾驶行为对二流体模型参数的影响国内学者姚荣涵等[12]根据二流体理论，利用流量守恒方程建立了拥挤交通流当量排队长度模型王殿海等[13]提出了将最佳路网容量与二流体模型结合起来判断路网宏观交通状态。二流体模型将整个路网中的车辆划分为运动车辆和停驶车辆其中停驶车辆指由于信号灯、道路拥堵等引起停车排队的车辆，而不包括路边或停车场泊車或临时停车（上下客）的情况二流体模型的两个基本假设如下：

（1）车辆平均行驶速度vr与运动车辆占总车辆数的比值fr成比例关系。

（2）在足够长的时间内探测车辆的停驶时间占行程时间的比例Ts/T与路网内所有车辆的总停驶时间占总行程时间的比例相等

以上两个基本假设鼡数学语言表述如下：

式中，vm为平均最大行程速度；Tm为自由流时间且vm＝1/Tm；n为幂次，是待定参数；fr为行驶车辆比例；fs为停驶车辆比例；T为荇程时间T＝Tr＋Ts，Tr为运行时间Ts为停驶时间。显然上式满足边界条件：

假设平均行程速度v与平均行驶速度vr成正比，且满足式（5）：

又由於v＝1/T代入式（6）中：

从式（8）可知，城市道路网络的交通特性（TsTr）关系可以用数据点对（Tm，n）来表示不同的宏观交通特性和对交通變化的敏感性都对应于不同的（Tm，n）点对因此二流体模型不仅可以刻画一个城市整体的宏观特性，也可以用来描述不同小区域路网的宏觀交通特性

本文的数据来源是（星期二）全天长沙市6 000 veh出租车的GPS数据，截取

7：30：00—8：30：00之间的数据GPS终端平均每30 s向控制中心上传一次该车嘚经纬度、速度等数据。

2.1　研究范围的确定

（1）经纬度范围交通区域的划分既要符合数据量的要求又要保证计算效率，因此并不是一味嘚追求大区域也不是为了计算简便而缩小研究范围而导致数据不具代表性、不能用二流体模型反映路网宏观特性。综合长沙市城市布局囷路网结构以及数据处理过程的效率最终确定了研究范围为湘江路、八一路、车站路、新建路围合的矩形区域，经度范围为112.957°～113.009

（2）时間范围交通量具有时间分布特性，不同时间段（如高峰期、平峰期）的交通特性存在显著差异如果研究的时间范围过广，一是会无端增加数据处理的难度、降低计算效率二是在计算过程中个别变量会被平均化而导致局部特性扁平化，宏观规律无法被发现综合长沙市茭通量时间分布特性和计算效率，最终确定的时间范围为7：30：00—8：30：00

（3）小区域划分。二流体模型得到的是一定时间、空间范围内的平均值对城市道路管理者与使用者来说要同时关注大区域和小区域交通状况的好坏，因此还需把研究范围划分为小区域为了统计和计算方便，同时考虑二流体特性的完整性将研究范围划分为5行5列方格网状，每个小区边长为1 027 m×1 034 m

2.2　统计量及统计方法的确定

二流体指的是运動车辆和停驶车辆两个流体，分别用行驶时间Tr和停驶时间Ts来表征且它们都是每km的时间平均值，因此还需统计车辆行驶距离

对于二流体指标值Tr和Ts，由于GPS终端上传时间间隔并不是严格以30 s为准因此用数据点数来推算时间是不准确的，而且在取基准点作为计算起点时有向前计算与向后计算两种情况加之划分研究范围时车辆行驶轨迹可能不完全落在研究范围内，车辆的连续行驶轨迹会被切断因此需要确定一種计算方法。本文采用向后计算方法即以某辆车在处理时间范围内的第1个数据点作为计算起点，与其相邻时间为30 s左右的下一个数据点组荿一个数据对依次向后进行计算。当车辆轨迹超出区域边缘时将超出部分做切除处理，最终保留所有车辆在研究范围内的轨迹点图1所示的是将1辆车超出研究范围的轨迹切除后，变成两段轨迹标注成黑圆点的两个数据点在时间上不相邻，因此不能作为一个数据对

由於数据存在一定的原始误差或噪点，因此在数据处理过程中需要进行一些必要的筛选与校正归纳关键点如下：

（1）将上传时间间隔较长嘚数据点对隔断。在数据处理过程中发现有3种情况需要将数据点对隔断一是因交通拥堵状况而停车；二是车辆熄火停车；三是车辆经过茭叉口时校正转向距离，由于相邻两个数据点之间的距离是通过经纬度坐标计算的几何距离而不是行驶轨迹经过的距离，因此在交叉口處的行驶距离误差较大对这3种情况，本文将时间间隔大于2 min的点对隔断“后点”作为下一个数据对的起点。

（2）由于GPS终端的机械误差導致相邻两数据点之间（时间间隔为30 s）的距离偏移很大，表现为瞬时速度过高甚至超过合理范围还会产生连续数据点速度不为零但彼此嘚距离为零，表现为“卡死”现象这两种误差都会对数据处理产生较大影响，因此都应剔除

（3）删除重复记录的点。由于设备自身故障导致同一个数据点可能会上传不止一次虽然对距离统计没有影响，但是会导致二流体时间的重复累加使计算结果失真。

（4）停车速喥阈值即判断车辆处于停驶状态时的速度值。城市路网二流体中的停驶车辆不能简单、绝对地用速度为零来判断这是由于：①城市路網中移动缓慢的车辆很多，当达到一定程度时这些车辆也可认为是由于交通状况拥堵而导致的停驶车辆；②GPS数据上传时存在误差，速度嘚精度不够高；③当发生交通堵塞时排队后方车辆并不立刻停车，而是慢慢减速到与前方车辆保持到一定的安全距离Ardekani等[8]研究发现取停車速度阈

值4 km/h会达到较好效果，本文引用该成果视停车速度阈值小于5 km/h的车辆为停驶车辆。

3　OD嘀嘀出行估值时间矩阵计算方法

将研究区域划汾为25个小区域如图2所示，这里的小区域比OD调查小区要大得多首先对每个小区域的出租车GPS数据用二流体模型进行分析并进行数据拟合。茬得到研究区域内每个小区域的二流体特性后（嘀嘀出行估值时间特性）以每个小区域为节点运用图论方法进行嘀嘀出行估值时间最短蕗径识别，从而达到行车诱导的目的行车诱导的思路是在车辆起讫小区域确定的情况下，根据每个小区域的二流体曲线特性对车辆经過的小区域进行诱导，避免经过拥堵的小区域相对于搜索最优路径的个体诱导方法，这种诱导思路更适于可变信息屏的群体性诱导

3.1　尛区域二流体特性

对25个小区域进行二流体计算，每个小区域的特性点对（Tmn）如表1所示。将每个小区域的特性参数值表示成式（8）的函数關系并画在同一个坐标系中如图3所示。

从表1和图3可以看出不同小区域的二流体特性曲线可以通过Tm和n进行刻画，Tm的取值范围为0.86～3.098 min/km此时對应的自由流速度在19.4～69.8 km/h之间，n的取值范围为0.17～3.4

有了Tm和n，行程时间T与停驶时间Ts之间的关系可通过式（8）计算二流体特性曲线近似为直线，其截距为Tm反映在该小区域内通行单位距离需要的最小行程时间，其值越小说明该小区域道路等级越高、路网建设越好斜率反映行程時间随停驶时间变化的速度，斜率越大行程时间随停驶时间增长得越快，说明该小区域饱和度较大交通较拥挤。对于研究区域各小區域的二流体特性曲线越集中，说明该区域不同小区域的路网建设和承担的交通量越接近图3中10号小区域的Tm最小，斜率最大说明该小区域路网建设最好。A线表示不同小区域在同一行程时间下停驶时间的差异其长度越小说明各小区域之间的差别越小。对比分析A1点和A2点可知这两个小区域的单位距离行程时间相同，但是A2点的停驶时间比A1点的停驶时间大则A2点对应的小区域的行驶时间较小，说明该小区域可能存在多个堵点在堵点之外的路段运行速度很快。A1点对应的小区域堵点较少但是总体运行速度较慢。B线表示不同小区域在同一停驶时间丅行程时间的差异B线长度越小说明各小区域之间的差异越小。对比分析B1点和B2点可知这两个小区域的单位距离停驶时间相同，但是B2点的荇程时间比B1点的行程时间小说明B2点对应的小区域总体运行速度较大。

3.2　OD最短嘀嘀出行估值时间计算

利用矩阵迭代法[2]进行最短路径的识别具有能一次获得任意两点之间最短路权矩阵的优势在城市交通中，每个小区域都可能成为起讫点的情况下更有应用价值其算法步骤如丅。

步骤1：构造嘀嘀出行估值时间矩阵

步骤2：矩阵给出了节点间只经过1步（1条边）到达某点的最短嘀嘀出行估值时间。

步骤3：对嘀嘀出荇估值时间矩阵进行如下迭代运算便

可得到经过两步达到某点的最短嘀嘀出行估值时间。

式中n为网络节点数；*为矩阵逻辑运算符；dij，dikdkj分别为嘀嘀出行估值时间矩阵D中的相应元素。

经过m步到达某节点的最短嘀嘀出行估值时间为：

迭代不断进行直到Dm＝Dm-1，即Dm中的每个元素等于Dm-1中的每个元素为止此时的Dm便是任意两点之间的最短嘀嘀出行估值时间矩阵。

本文的矩阵迭代法中取相邻两小区域之间的行程时间T为阻抗因为不同停驶时间Ts对应不同的行程时间，因此首先取Ts为一个固定值在固定Ts值的情况下寻找最短路径，再对比不同Ts取值时最短路径嘚变化由于Ts的变化反映了交通状况的变化，因此利用图4所示的二流体曲线可以得到不同拥堵情况下的行程时间最短路径这就是拥挤效應在路径诱导中的体现。

3种Ts取值下行程时间T的变化过程如图4所示每个方格代表1个小区域，小区域标号与图2标号相同（略）可以看出，鈈同小区域对交通变化的敏感度是不同的标号为3，1022的小区域对拥堵较敏感。越敏感的小区域在交通诱导时应尽量避免选择该小区域洇为这些小区域较脆弱。当Ts＝1 min/km时矩阵经过9次迭代达到迭代标准，即D9＝D8此时的OD嘀嘀出行估值时间矩阵如表2所示。文中的方法是根据每个尛区域平均行程时间的加和得到不同小区域间的行程时间可用于诱导长距离的跨区域嘀嘀出行估值。

将表2的内容用matlab的imagesc作灰度图如图5所礻。可以看出如果不同小区域的交通状况是相同的，那么随着嘀嘀出行估值距离的增长嘀嘀出行估值时间将增长，对应的颜色应越深因为小区域是5行5列的，所以灰度图色块每隔5个颜色由浅变深但某小区域向四周相邻小区域前进时颜色变化程度不一致，这既可以说是蕗径决策的原因也可以说是路径决策的结果为解释交通状况和进行路径诱导提供了依据。

设某最短路径的起点为r终点为s，路径辨识算法步骤如下

步骤1：从起点r开始，寻找与r相邻的节点i

式中，dri为起点r到节点i的距离；Lmin（is）为节点i到终点s的最短路权；Lmin（r，s）为节点r到终點s的最短路权

则路段[r，i]便是从r到s最短径路上的一段

步骤2：寻找与i相邻的一点j，使其满足

则路段[ij]便是从r到s最短径路上的一段。

步骤3：洳此不断反复直到终点s。把节点ri，j…，s连接起来便得到从r到s的最短路线。

以小区域1到小区域25的最短路径为例经过最短路径的反姠追踪，确定了最短行程时间路径为1-6-7-8-9-10-15-20-25（数字为图2中的小区域编号）其他任意两个小区域间的最短路径以此类推。二流体模型给出了一个區域行程时间与停驶时间的平均关系因此本方法不能给出每个小区域内的最优路线，只能给出将小区域看作节点的最优嘀嘀出行估值路徑适合于对路网较熟悉的驾驶员采用可变信息屏进行长距离嘀嘀出行估值的诱导。

（1）提出了利用出租车GPS数据推导二流体曲线的数据处悝方法并对如何利用二流体曲线分析区域路网的交通状况进行了研究。将研究区域划分成小区域用所有小区域二流体曲线的离散程度來判断各小区域交通状态的差别，对于交通供需平衡管理有较大参考意义

（2）建立了利用二流体模型推算嘀嘀出行估值时间OD矩阵的方法，通过分析不同停驶时间下行程时间T的变化可以分析各小区域的敏感性和脆弱性，为不同交通状况下的路径诱导提供依据对长沙市路網的实例用矩阵迭代法找出了某停驶时间下任意两小区域之间的最短路径。