考虑交通量随机波动的随机用户均衡配流模型
卫翀1,2, 邵春福1
1.北京交通大学 城市交通复杂系统理论与技术教育部重点实验室,北京 100044
2.首都世界城市顺畅交通协同创新中心,北京 100124

作者简介:卫翀(1981-),男,副教授,博士.研究方向:交通网络建模与分析.E-mail:chwei@bjtu.edu.cn

摘要

研究了随机用户均衡路网中交通量的随机波动现象。基于随机效用理论,运用贝叶斯定理构建了可以反映随机用户均衡路网中路径路段交通量概率分布的交通配流模型,避免了交通量的随机性在交通配流模型中被忽略的问题。侧重分析了路网中驾驶员的随机路径选择行为以及道路通行能力的不确定性对交通量波动的影响。设计了一种马科夫链-蒙特卡洛算法用来根据交通量概率分布估计路径路段交通量的统计特征。在一个中等规模路网上对提案模型进行了验证,结果表明:模型可以有效地评估出路径选择行为以及通行能力不确定性等因素对交通量波动的影响。

关键词: 交通运输系统工程; 随机用户均衡; 交通量概率分布; 贝叶斯定理; 随机效用理论
中图分类号:U491 文献标志码:A 文章编号:1671-5497(2015)05-1408-06
Stochastic user equilibrium traffic assignment model consisting fluctuations in traffic flows
WEI Chong1,2, SHAO Chun-fu1
1.MOE Key Laboratory for Transportation Complex Systems Theory and Technology, Beijing Jiaotong University, Beijing 100044, China
2.Center of Cooperative Innovation for Beijing Metropolitan Transportation, Beijing 100124, China
Abstract

The stochastic nature of traffic flow in stochastic user equilibrium networks was investigated. Based on the random utility theory, a stochastic traffic assignment model that can reflect the probability distribution of traffic flow was constructed using Baye's theorem. Driver's route choice behavior and the uncertainty in road capacity were taken into account in the model as the sources of fluctuation in traffic flow. A Markov chain Monte Carlo algorithm was also developed in order to estimate the characteristics of traffic flows from the probability distribution. The proposed model was tested on a medium-sized network. Results show that the model can be used to evaluate the effects of route choice behavior and uncertainty in road capacity on the fluctuations in traffic flows.

Keyword: engineering of communication and transportation; stochastic user equilibrium; probability distribution of traffic flows; Baye's theorem; random utility theory
0 引 言

路网中交通量的随机波动与路网运营质量紧密相关[1], 是建立交通配流模型时需要考虑到的一个重要现象。在以往的研究中, 研究者们提出了多种可以在交通配流模型中描述交通量随机波动性的建模方法。例如, Cascetta[2]以及Watling[3]提出了基于马科夫链理论的逐日学习型交通配流模型。Asakura等[1, 4], Shao等[5]从OD间交通需求随机波动的角度出发讨论了路网中交通量的随机波动性。Hazelton[6, 7]基于Hammersley-Clifford定理构建了一种遵循随机效用理论并可以输出路径路段交通量统计特征的CSUE(Conditional stochastic user equilibrium)模型, 但是该模型并不能给出交通量概率分布的具体形式, 并且该模型的计算量会随交通需求量的增加而增加。Watling[8]基于渐进均衡条件建立了广义随机用户均衡模型(GSUE)并应用泰勒展开讨论了模型的解法, 但是该模型同样不能明确交通量概率分布的具体形式。Wei等[9]讨论了均衡网络中交通量概率分布的情况, 但是没有考虑路网属性的不确定性与交通量随机波动性之间的关系。国内也有学者从交通系统演进的角度探讨了交通量的随机特性[10, 11], 但是这些研究的目的并不在于探索均衡网络中交通量的概率分布。在O-D需求确定的随机用户均衡路网中路径交通量的随机性主要源于驾驶员的随机路径选择行为和路网本身属性的不确定性, 并且交通量实质上是随机变量, 其概率分布的形式应在交通配流模型中得以明确体现。

本文提出了一种基于随机效用理论, 考虑到了随机用户均衡网路中交通量随机波动性的交通配流模型, 并开发算法根据模型反映出的交通量概率分布估计交通量的统计特征。

1 交通量的随机波动性
1.1 通行能力确定的情况

下面将从随机效用理论出发逐步推导出随机用户均衡路网中交通量的概率分布。基于随机效用理论, 可用下式来描述路径选择问题:

P(subi|ci)=P(Uci> Uhh, hci)=P(Vci+εci> Vh+εhh, hci)(1)

式中:P(subi|ci)是一个条件概率, ci为驾驶员i的路径选择结果, subi表示驾驶员i服从随机用户行为(Stochastic user behavior)[7]。随机用户行为的定义如下:驾驶员选择了一条他或她所感知的效用最大的路径的这样一种行为被称为随机用户行为。应注意式(1)中条件概率的前提条件仅仅表明驾驶员的路径选择结果是ci; 另一方面subi仅仅表示了驾驶员选择路径的行为规则, 但是并没有包含路径选择结果的信息。如式(1)所示, 驾驶员i的路径选择结果是ci这一条件已给定的前提下subi可以被等价地表达为 Uci> Uhh, hci, 其中Uh为路径h的随机效用, 随机效用的可观测部分Vh由路径行程时间决定, ε h为驾驶员对路径行程时间的感知误差。

接下来在交通网络的背景下进一步扩展式(1)。为了方便起见, 用I来表示网络中所有驾驶员集合; 用向量c表示cjjI, 即网络中所有驾驶员的路径选择结果; 用R表示网络中所有路径的集合; 用向量f表示路径交通量集合。参照式(1), 当所有驾驶员的路径选择结果c为给定条件时, 驾驶员i服从随机用户行为这一事件的概率可以表示为:

P(subi|c)=P(Vci(c)+εci> Vh(c)+εhh, hci)(2)

式中:Vh(c)为根据所有驾驶员的路径选择结果所确定下来的可观测效用。具体来说, Vh(c)包含了下面一系列的关系:由于路径交通量是路径选择结果的叠加, 所以根据c可以计算出路网中的路径交通量f, 并且这个结果是唯一确定的; 根据路径交通量又可以计算出每条路径的行程时间, 路径的行程时间进而又决定了Vh的值。

在随机用户行为的基础上, Hazelton[7]指出随机用户均衡可以表达为:当路网上所有的驾驶员都服从随机用户行为时, 路网处于随机用户均衡状态。如果用sue来表示路网处于随机用户均衡状态, 那么sue可以等价地表示为subiiI。当所有驾驶员的路径选择结果c作为给定条件时, 路网处于用户均衡状态这一事件的概率为:

P(sue|c)=iIP(subi|c)=p(r|c)fr(c)3

式中:fr(c)为根据所有驾驶员的路径选择结果c所确定的路径r的交通量; p(r|c)为在路径选择结果c作为条件给定的情况下, 路径r的效用为最大这一事件发生的概率。

p(r|c)=P(Vr(c)+εr> Vh(c)+εhh, hr)(4)

需要补充说明的是, 式(3)中假设所有驾驶员都是同质的(homogeneous), 也就是说如果ci=cm。式(3)所显示的概率实际上是路径选择结果的似然函数, 应用贝叶斯定理可以从似然函数反过来推导出在路网已处于用户均衡状态这一条件给定的情况下路网上驾驶员路径选择结果是c的概率:

P(c|sue)=P(sue|c)P(c)P(sue)5

式中:驾驶员的路径选择结果是随机变量, 从贝叶斯统计学的观点来看, P(c)可以被理解为是路径选择结果的先验概率, 这个先验概率是一个均匀分布; P(sue)表示网络处于用户随机均衡的概率, 这一概率的值并不受路径选择结果c的影响。基于以上两点原因, 式(5)可以被改写为:

P(c|sue)P(sue|c)(6)

对于一个交通配流问题, 研究的重心在于分析交通量的分布情况, 而不仅仅是掌握所有驾驶员的路径选择行为结果。交通量实质上是驾驶员路径选择结果的叠加, 所以可以通过式(6)推导出网络已处于随机用户均衡状态下路径交通量的概率分布。在开始推导之前, 应该先注意到如果所有驾驶员的路径选择结果c已经给定, 路网中所有路径的交通量可以被唯一确定下来。但是反过来, 如果网络中路径交通量给定, 与之对应的路径选择结果c并不是唯一确定的, 而是有M种可能的情况。M的计算公式为:

M=c:f(c)=f1=nNqn!rRnfr!7

式中:f(c)为根据c算出的路径交通量向量; ∀ c:f c=f表示所有可以使等式f(c)=f成立的路径选择结果; N为网络中OD对的集合; Rn为第n对OD之间的路径集合; qn为第n对OD之间的交通需求。

根据式(6)(7), 网络处于随机用户均衡状态下路径交通量为f的概率可以表示为:

P(f|sue)=c:f(c)=fP(c|sue)c:f(c)=fP(sue|c)=c:f(c)=fp(r|c)fr(c)=nNqn!rRnfr!p(r|f)fr8

式中:

p(r|f)=P(Vr(f)+εr> Vh(f)+εhh, hci)(9)

Vr(f)的可观测效用是由与路径交通量对应的路径行程时间所确定的。

需要指出的是, 式(8)所示的概率模型与Sheffi[12]提出的多项式概率模型是完全不同的。Sheffi模型中路径行程时间是恒定不变的, 不受路径交通量大小的影响, 所以其模型并非本文所讨论的均衡模型范畴。

1.2 通行能力不确定的情况

式(8)反映了随机路径选择行为与路径交通量随机波动性之间的关系。本文接下来进一步考察路网属性的不确定性对路径交通量随机波动性的影响。这里仅讨论路段通行能力的不确定性这一因素。但是提出的模型具有普遍性, 其他因素可以通过相同的方法来考虑。

s来表示路段通行能力并且设通行能力服从正态分布。用μ σ 分别表示该正态分布的均值和标准差。在交通配流问题中, sueμ σ 被视为给定的条件, 配流结果可以用条件概率P(f|sue, μ , σ )来表述。推导P(f|sue, μ , σ )和推导P(f|sue)的思路类似, 可以从P(sue|c, s)开始讨论。该概率表示如下:

P(sue|c, s)=p(r|c, s)fr(c)10

式中:p(r|c, s)表示在cs作为给定条件的情况下, 路径r的效用为最大这一事件发生的概率, 其中路径的可观测效用由cs决定的。

通过P(sue|c, s)可以进一步得到P(f, s|sue, μ , σ ), 有:

P(f, s|sue, μ, σ)=c:f(c)=cP(c, s|sue, μ, σ)=c:f(c)=cP(sue|c, s)P(c)P(s|μ, σ)P(μ, σ)P(sue, μ, σ)c:f(c)=cp(r|c, s)fr(c)P(s|μ, σ)=nNqn!rRnfr!p(r|f, s)frP(s|μ, σ)(11)

式中:P(μ , σ )和P(sue, μ , σ )为常数, 因为这些概率的值都不受fs影响; P(sue|c, s)已经由式(10)给出; p(r|f, s)表示在fs作为条件给定的情况下, 路径r的效用为最大这一事件发生的概率; P(s|μ , σ )表示给定均值和标准差的情况下路段通行能力为s的概率。

P(f, s|sue, μ , σ )可以最终推导出当sueμ σ 为给定条件时, 路网中路径交通量为f的概率:

P(f|sue, μ, σ)=sP(f, s|sue, μ, σ)dssnNqn!rRnfr!p(r|f, s)frP(s|μ, σ)ds12

式(12)给出了路径交通量f所对应概率密度函数的一般形式, 概率密度函数的具体表达将随着与p(r|f, s)相对应的路径选择模型的不同而变化。通过概率密度函数的表达式实际上并不能直观地反映出路径路段交通量的均值和方差等特征, 因此本文将通过马科夫链-蒙特卡洛(Markov chain Monte Carlo)算法来根据推导出的概率密度函数估计出交通量的均值和方差等特征。

2 交通量的统计特征

根据已知概率密度函数可以计算随机变量的统计特征。式(12)反映了路径交通量这一随机变量所对应的概率密度函数。虽然该概率密度函数的形式较为复杂, 但可用马科夫链-蒙特卡洛算法[13, 14]来估计交通量的统计特征值。马科夫链-蒙特卡洛算法的数学基础是马科夫过程, Andrieu等[15]详细介绍了该算法的基本原理。应用该算法可以产生一系列来自概率分布P(f|sue, μ , σ )的随机样本。通过汇总这些随机样本便可以得到f的统计特征值。然而, 马科夫链-蒙特卡洛算法在运行过程中需要计算式(12)所示的概率密度函数, 因此会产生频繁的积分计算。本文通过将数据扩张技术(Data augmentation)与马科夫链-蒙特卡洛算法结合使用来避免积分计算。所设计的算法并不直接从概率分布P(f|sue, μ , σ )中抽取f的随机样本, 而是从式(11)所示的联合概率分布P(f, s|sue, μ , σ )中抽取fs的随机样本, 然后用抽出的f的随机样本来估计出路径交通量的均值和方差等统计特征。

以下用ftst表示从P(f, s|sue, μ , σ )中抽取的第t个样本。用马科夫链-蒙特卡洛算法抽出ftst的过程可简要概括如下:

步骤1 抽出一个随机数f* , 再从多维正态分布中抽取一个随机样本s* , 该多维正态分布均值为st-1, 协方差矩阵是一个维数与st-1相对应的单位矩阵。

步骤2 计算

τ=P(f* , s* |sue, μ, σ)P(ft-1|f* )P(ft-1, st-1|sue, μ, σ)P(f* |ft-1)13

步骤3 以min(1, τ )的概率接受f* , s* 为抽取的第t个样本, 即以min(1, τ )的概率设ft=f* , st=s* , 以1-min(1, τ )的概率设ft=ft-1, st=st-1

在抽出第一个样本ftst之前需要设定一个初始值f0s0, 这个初始值并不影响抽样过程, 可以任意设定。在上述算法中f* 是来自一个多项式概率分布的随机样本。该分布的参数由ft-1决定, 其概率密度函数为:

P(f* |ft-1)=nNqn!rRnfr* !frt-1qnfr* 14

算法中P(ft-1|f* )的值可通过交换式(14)中ft-1f* 的位置后进行计算。

重复上述步骤t次后产生的结果f1, s1, f2, s2, …, ft, st是从式(11)所对应的概率分布中抽取出的一组随机样本。

3 数值计算示例

本文用文献[9]中的含有13个节点的路网测试了提出的交通配流模型。路网中有3对O-D对, 分别是点1到点2, 点1到点3, 和点4到点2。三对O-D对之间的交通需求分别是50, 30和30。本文中用下面的Logit模型来计算p(r|f, s)的值并用β 表示Logit模型的参数。在实际应用中β 也可被设置为随机变量[16]

β =0.35, μ =40, σ =5。应用前文描述的马科夫链-蒙特卡洛算法, 从概率分布P(f, s|sue, μ , σ )中抽取一系列路径交通量f的随机样本。图1图2用抽出的随机样本的直方图展示了通过节点1, 12, 8, 2的路径和通过节点4, 5, 9, 10, 11, 2的路径所承担交通量的概率分布情况。从这些输出结果可以看出:提案交通配流模型可以完整地描述路网中交通量的随机波动性, 而不是像传统的配流模型那样仅输出交通量的均值。基于马科夫链-蒙特卡洛算法抽取的随机样本还可以进一步估计出交通量的统计特征, 其中交通量的方差具体反映了交通量的随机波动性。

模型参数β σ 实际上分别反映了路径选择行为的随机性和路网属性的不确定性。其中β 为路径选择模型的参数, 该参数的取值越大, 表明驾驶员对路径行程时间的感知误差越小, 即路径选择的随机性越小。又由于路径选择结果与交通量之间存在叠加关系, 所以β 的值对交通量的方差是有影响的。为了考察这一点, 设μ =40, σ =5, 分别计算了β 等于0.2, 0.4, 0.6, 0.8, 1和1.2时路网中路段交通量方差的平均值。各个路段交通量的方差可以用马科夫链-蒙特卡洛算法所抽取的随机样本来计算。图3显示了计算结果。从图中可以看到:随着β 取值的增大, 或者说随着驾驶员对路径行程时间的感知误差的减小, 路段交通量方差的平均值将逐渐减小。

在交通配流模型中σ 表示路段的通行能力。σ 变大会使路径行程时间的不确定性变大, 与之相关联的驾驶员路径选择结果的随机性也会随之加大。

图1 经过点1-12-8-2的路径所承担交通量的随机样本直方图Fig.1 Histogram of simulated samples for the route that visits node 1-12-8-2

图2 经过点4-5-9-10-11-2的路径所承担交通量的随机样本的直方图Fig.2 Histogram of simulated samples for the route that visits node 4-5-9-10-11-2

图3 β 的取值对路段交通量方差的影响Fig.3 Effect of β on the variances of link traffic flows

为了考察σ 的取值对路段交通量随机波动性的影响, 设β =0.35, μ =40, 分别计算了σ 等于2, 4, 6, 8和10时路网中路段交通量方差的平均值。图4为计算结果。从图中可以看出:随着通行能力不确定性的增大, 路段交通量方差的平均值也随之增大。从σ 对路径选择结果的影响以及路径选择结果与交通量的关系可以看出, 这种交通量方差逐渐增大现象本质上是由于路径选择结果的不确定性相应变大引起的。

图4 σ 的取值对路段交通量方差的影响Fig.4 Effect of σ on the variances of link traffic flows

4 结束语

在处于随机用户均衡的道路网中, 路径交通量是随机波动的, 其所对应的概率分布可以运用贝叶斯定理根据随机效用理论推导出来。基于马科夫链-蒙特卡洛算法可以进一步计算出随机用户均衡网络中路径路段交通量的均值和方差等统计特征值。研究结果表明:驾驶员的路径选择行为以及道路通行能力的不确定性均会影响到交通量的随机波动特性, 这种影响在路径交通量的概率密度函数中有明确的体现。数值计算结果显示, 当驾驶员对行程时间的感知误差较大时或道路通行能力的不确定性较高时都会使得路径路段交通量的方差变大。基于随机效用理论所建立的随机用户均衡配流模型可以描述出交通量的随机波动性, 为交通网络可靠性分析提供必要的信息。

The authors have declared that no competing interests exist.

参考文献
[1] Asakura Y, Kashiwadani M. Road network reliability caused by daily fluctuation of traffic flow[C]∥Proceedings of the 19th PTRC Summer Annual Meeting, Brighton UK, 1991: 73-84. [本文引用:2]
[2] Cascetta E. A stochastic process approach to the analysis of temporal dynamics in transportation networks[J]. Transportation Research Part B, 1989, 23(1): 1-17. [本文引用:1]
[3] Watling D. Asymmetric problems and stochastic process models of traffic assignment[J]. Transportation Research Part B, 1996, 30(5): 339-357. [本文引用:1]
[4] Asakura Y. Reliability measures of an origin and destination pair in a deteriorated road network with variable flows[C]∥Proceedings of the 4th Meeting of the EURO Working Group in Transportation, Newcastle UK, 1996: 273-288. [本文引用:1]
[5] Shao H, Lam W H K, Tam M L. A reliability-based stochastic traffic assignment model for network with multiple user classes under uncertainty in demand [J]. Networks and Spatial Economics, 2006, 6(3/4): 173-204. [本文引用:1]
[6] Hazelton M L, Lee S, Polak J W. Stationary states in stochastic process models of traffic assignment: a Markov chain Monte Carlo approach[C]∥Proceedings of the 13th International Symposium on Transportation and Traffic Theory, Lyon, France, 1996: 341-357. [本文引用:1]
[7] Hazelton M L. Some remarks on stochastic user equilibrium[J]. Transportation Research Part B, 1998, 32(2): 101-108. [本文引用:3]
[8] Watling D P. A second order stochastic network equilibrium model. I: theoretical foundation: a second order stochastic network equilibrium model. II: solution method and numerical experiments[J]. Transportation Science, 2002, 36(2): 149-183. [本文引用:1]
[9] Wei C, Asakura Y, Iryo Y. The posterior probability distribution of traffic flow: a new scheme for the assignment of stochastic traffic flow[J]. Tranportmetrica A: Transport Science, 2013, 9(8): 753-771. [本文引用:1]
[10] 刘天亮, 黄海军. 日常择路行为的多智能体模拟[J]. 物理学报, 2007, 56(11): 6321-6325.
Liu Tian-liang, Huang Hai-jun. Multi-agent simulation on day-to-day route choice behavior[J]. Acta Physica Sinica, 2007, 56(11): 6321-6325. [本文引用:1]
[11] 刘诗序, 关宏志, 严海. 预测信息下的驾驶员逐日路径选择行为与系统演化[J]. 北京工业大学学报, 2012, 38(2): 269-274.
Liu Shi-xu, Guan Hong-zhi, Yan Hai. Drivers'day-to-day route choice behavior and system evolution under forecast information[J]. Journal of Beijing University of Technology, 2012, 38(2): 269-274. [本文引用:1]
[12] Sheffi Y. Urban Transportation Networks: Equilibrium Analysis with Mathematical Programming Methods[M]. Englewood Cliffs, NJ: Prentice-Hall, 1985. [本文引用:1]
[13] Metropolis N, Rosenbluth A W, Rosenbluth M N. Equations of state calculations by fast computing machines[J]. Journal of Chemical Physics, 1953, 21(6): 1087-1091. [本文引用:1]
[14] Hastings W K. Monte Carlo sampling methods using Markov chains and their applications[J]. Biometrika, 1970, 57(1): 97-109. [本文引用:1]
[15] Andrieu C, de Freitas N, Doucet A, et al. An introduction to MCMC for machine learning[J]. Machine Learning, 2003, 50(1): 5-43. [本文引用:1]
[16] 罗清玉, 孙宝凤, 吴文静, . 基于Mixed logit模型的拥挤收费下交通方式分担率预测[J]. 吉林大学学报: 工学版, 2010, 40(5): 1230-1234.
Luo Qing-yu, Sun Bao-feng, Wu Wen-jing, et al. Predicting traffic mode split under congestion pricing based on mixed logit model[J]. Journal of Jilin University (Engineering and Technology Edition), 2010, 40(5): 1230-1234. [本文引用:1]