交通分布--四阶段法预测第二阶段ijqijqji主要内容概述增长系数法重力模型法机会模型最大熵法重点了解(自学)难点交通分布预测是把交通的发生与吸引量预测获得的各小区的出行量转换成小区之间的出行交换量(即出行分布量),即PA矩阵。一、出行分布量指分区i与分区j之间平均单位时间内的出行量。即qij和qji。二、出行分布矩阵二维表(矩阵)。图矩阵中各个组成元素及其关系。第一节概述PA12jn合计12in合计…...…...1P2PiPQ1A2AiA…...…...…...…...…...…...nPnA发生交通量吸引交通量生成交通量ijqijiijjijjjiiiijjjijiqqAPQqAqP三、分布量预测1、已知现状PA表和预测的Pi和Aj;两个分区之间的交通阻抗矩阵.2、方法增长率/系数法:假定现在和将来的交通分布模式变化不大,在此基础上预测对象区域目标年的PA交通量。(需要有完整的现状OD表)综合法:从分布交通量的实际分析中,剖析PA分布规律,并用数学模型表现此规律,然后用实测数据标定模型参数,最后用标定的模型预测分布交通量。(方法包括重力模型法、介入机会模型法、最大熵模型法等,应用范围更广)第一节概述第二节增长系数法包括两大类:增长函数法和Furness约束条件法一、增长系数法的思路原理:假设在现状分布交通量给定的情况下,预测将来的分布交通量。1、已知现状PA表(qij0,Pi0、Aj0),规划年产生量Pi和吸引量Aj。令k=02、计算各分区第k次发生增长率、吸引增长率3、设f(Fpi,Faj)为增长函数,计算第(k+1)次预测值kiikpiPPFkjjkajAAF),(1kajkpikijkijFFfqq第二节增长系数法4、检验预测结果计算新的产生量Pik+1和吸引量Ajk+1,计算新的增长率:看是否满足误差范围ε。若是,qijk+1即为所求;若不是,令k=k+1,转至第3步,进行迭代。二、各种增长系数法根据增长函数的种类不同,提出各种具体的增长系数法。1、常增长系数法qij的增长仅与i区的产生量增长率或与生成量增长率有关,增长函数为:常数,常)(ajpiFFf11kiikpiPPF11kjjkajAAF1DO123现状值将来值1422820235412203233825现状值910928将来值25182265求将来OD分布交通量(单位:万次)。设定收敛标准ε=3%DO123增长率预测将来值11055120258.3336.66712036.259.3759.375125增长率1.1760.7931.046预测将来值21.2522.7121.0465第二节增长系数法1、常增长系数法特点:由于发生量与吸引量不对称,故预测精度不高,是一种最粗糙的方法,甚至不能保证迭代过程一定能收敛。2、平均增长率法i,j小区之间的分布交通量的增长系数是i小区出行发生量增长和j小区出行吸引量增长的平均值,即:常数,常)(ajpiFFf)+(,平ajpi21)(FFFFfajpi特点:公式简明,易于计算;但迭代步数较多,收敛慢,计算精度低,但仍被广泛使用;随着计算机的发展,正逐渐被D法和F法所取代。DO123增长率将来值111.33.850.99520.126.26.67.212037.47.79.81.00424.9增长率1.0040.9941将来值24.918.12265平均增长率法,迭代第6次的最后结果第二节增长系数法3、Detroit法(底特律法,简称D法)jjjjjjiiajpiajpiDAAAAPPQQFFFFf0000///)(,该方法是Detroit(底特律)市1956年规划首次被开发利用,J.D.Carol提出的。特点:考虑的因素较平均增长率法全面;但同样是收敛速度慢,需要多次迭代才能求得将来年的分布交通量。DO123现状值将来值1422820235412203233825现状值910928将来值25182265求将来OD分布交通量(单位:万次)。设定收敛标准ε=3%第二节增长系数法4、Fratar法(福莱特法)1954年T.J.Fratar提出的。分布交通量的增长率是用出行发生量误差修正量和出行吸引量误差修正量的组合平均值。发生/吸引位置系数特点:收敛速度快、计算精度高、误差小,满足相同的精度条件下迭代次数少,广泛应用;但其计算过程较复杂(说明每次收敛受到更多的约束),一般通过计算机编程实现,或通过专门的交通规划软件计算。2)(kajkpikajkpikajkpiFLLFFFFf,ikpikijkjkajjkajkijkikpiFqALFqPL,DO123现状值将来值1422820235412203233825现状值910928将来值25182265求将来OD分布交通量(单位:万次)。设定收敛标准ε=3%DO123增长率将来值111.33.850.99520.126.16.87.112037.57.59.91.00424.9增长率1.0040.9941将来值24.918.12265F法,迭代第2次的最后结果第二节增长系数法二、各种增长系数法5、Furness约束条件法(佛尼斯法)假设i,j小区间分布交通量qij的增长系数与i小区的发生增长系数ui和j小区的吸引增长系数vj都有关系。但这两个系数不是简单等于产生量或吸引量的增长率。采用循环迭代方法,求解uik与vjk,直到该系数满足设定的收敛标准为止。特点:方法计算相对简单,收敛速度相对较快,也适合编程获得预测结果。jiiijjiijijjijijijAuqvqPvquq00DO123现状值将来值1422820235412203233825现状值910928将来值25182265求将来OD分布交通量(单位:万次)。设定收敛标准ε=3%DO123增长率将来值111.343.744.9612026.116.757.13119.9937.567.519.92124.99增长率111将来值25182265Furness法,迭代第4次的最后结果第二节增长系数法6、增长系数法的特点优点:1)结构简单、实用,不需要交通小区之间的距离和时间2)适用于小时交通量或日交通量等的预测,也可以获得各种交通目的的PA交通量3)适合短期的、局部的预测(费用不高,不用作时间、距离的调查)缺点:1)必须有所有小区的PA交通量2)难以考虑网络或网络结构的变化对出行的影响3)对稀疏矩阵、部分矩阵、交通小区之间交通量较小时(即空元素或零元素),适用性比较差。4)将来交通量仅用一个增长系数表示,缺乏合理性第三节重力模型法交通阻抗???最早出现的重力模型是Casey在1955年提出。一、简单模型(无约束重力模型)1、模型两个交通小区的交通分布量与两个交通小区的出行发生量与吸引量成正比,与交通小区之间的交通阻抗成反比。改进:α,通常称为潜能系数,一般在0.5~1.0间取值;多数情况下,可取经验值α==1或0.5。R广义费用(又叫交通阻抗),可以是出行时间、距离、费用、票价、油耗或指标的综合。(模型构造上不能保证预测值满足交通守恒约束条件)2ijjiijRAPKqijjiijRAPKq费用可以用距离、时间或货币来表示。一般地,可通过一个统一的度量来表示,通常将这个度量称为出行的广义费用。它是典型的出行属性的线性函数,该函数的系数代表着属性的重要程度(权重)。如交通方式k的广义费用的一种可能表达形式为:jijnijvijtijwijijaFatatatataR654321式中:twij-到达或离开车站的步行时间;ttij-在车站的等待时间;tvij-i到j在车旅行时间;tnij-i、j间发生的换乘时间;Fij-i、j间出行需要的费用;φj-i到j的节点费用(特别是停车);δ-随机变量,代表未被包含的属性,如安全、舒适和方便性等;a1,a2,…a6-各费用元素的权重;它们将各种属性转换成合适的度量单位,如货币或时间。第三节重力模型法一、简单模型(无约束重力模型)2、参数的标定采用线性回归方法标定模型中的参数。两边取自然对数:此处Pi,Aj,Rij,qij可从现状调查数据中取若干个分区作为样本,待标定参数有lnK、α、β、-γ,用多元线性回归方法。若预先取α==1,则问题简化为:只须用一元线性回归方法就能标定参数K和γ。ijjiijRAPKqlnlnlnlnlnijjiijRAPKqijjiijRAPKqDO123合计117.07.04.028.027.038.06.051.034.05.017.026.0合计28.050.027.0105.0166.539.3预测值90.336.9预测值38.691.936.012318.017.022.0217.015.023.0322.023.07.00ijRijjiijRAPKqijjiijRAPKq)(1步:建立模型;2步:标定模型参数;88.972.418.9180.275.5237.946.2359.618.843.976.0138.8183.2354.3141.1678.6可以认为求得的qij是第1次近似值,然后通过增长系数法的迭代计算使得两者一致。DO123合计123合计166.539.3预测值90.336.9预测值38.691.936.012314.09.011.029.08.012.0311.012.04.0ijR3步:运用标定的模型预测交通分布。455.1173.1)(124.0ijjiijRAPq第三节重力模型法一、简单模型(无约束重力模型)3、误差说明2)阻抗问题:分区内的出行,阻抗不易确定。若Rij接近于0,则出行分布量将产生偏大的估计;其解决方法是:1)首要原因是这类模型不遵循交通守恒的约束;解决的办法是增加约束条件。一是qii不用引力模型,用其他方法,作为第1次近似值。如以分区规模和交通服务条件作自变量建立回归分析模型;二是修改分布阻抗,采用更一般的函数关系f(Rij);第三节重力模型法常见的交通阻抗函数有:ijijbRaRf1)(Rijf幂型指数型复合型半钟型指数型较常用;复合型用于特定的交通方式的出行分布中;离散型:m为费用等级r阻抗平均值δmij取0,1mijmmijrRf)(阻抗函数的选择??ijijbRijijbRijijijeRRfeRfRRf)()()(第三节重力模型法二、单约束引力模型(A.M.Voorhees乌尔希斯重力模型)用上述阻抗函数代替阻抗项,假定α=β=1,则:考虑分布量的和与发生交通量一致,模型为:满足K的引力模型都叫单约束引力模型。同理,可以定义一个列约束系数:)(ijjiijRfAKPqjijjRfAK)(1iijiRfPK)(1)()(ijjjijjiijRfARfAPq第三节重力模型法美国公路局模型(BPR)在上述修正单约束重力模型的基础上导入反映小区i和小区j之间固有关系的调整系数kij(也叫地域间结合度),模型为:预测时,Kij在将来不做变化,直接使用。ijijijjjijjiijKKRfARfAPq)()(第三节重力模型法三、双约束引力模型满足上述两个约束条件的模型为双约束引力模型。jijjjiRfAKK)(1iijiijRfPKK)(1)(ijjijiijRfAPKKq四、重力模型参数的标定一般情况下,由于研究区域内广义费用Rij不完全,阻抗参数γ等必须单独标定;若采用组合阻抗函数,还需考虑参数的灵活性。第三节重力模型法参数的标定的过程(多采用试算法):1.假设参数初值;2.利用现状数据和模型,得出理论交通分布;利用交通分布,计算某一控制指标的理论值;3.利用现状交通分布,计算控制指标的现状值,如平均交通阻抗;4.计算理论值与实际值两者之间的相对误差,看是否满足一定的误差范围(如3%);5.满足误差范围,初值可取;否则,分析理论值与实际值关系,调整初