发明名称 一种城市轨道交通随机客流加载方法
摘要 本发明公开了一种城市轨道交通随机客流加载方法,包括以下步骤:1.1、初始化,将一天平均划分为多个时间窗,并加载基础数据,同时设定时间窗总数为n,设定标记时间窗的变量t的初始值为1,设定路段拥挤费用的初始值CN<sub>1</sub>=0;1.2、读取第t时间窗的OD分布量,基于基本有效路径集合,计算各路径的首末班车时间,判断该时间与第t时间窗的关系,生成动态有效路径集合;1.3、基于动态有效路径集合和城市轨道交通乘客出行路径选择模型,确定OD间各有效路径选择比例,获得路径流量;1.4、统计各路段的客流量,更新路段拥挤费用CN<sub>t+1</sub>;1.5、终止判断,如果t&lt;n,则t=t+1,从步骤1.2继续迭代;否则,该方法终止。
申请公布号 CN103218670B 申请公布日期 2017.02.08
申请号 CN201310093784.2 申请日期 2013.03.22
申请人 北京交通大学 发明人 姚恩建;张永生;潘龙;杨扬
分类号 G06Q10/04(2012.01)I;G06Q50/30(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京市商泰律师事务所 11255 代理人 毛燕生
主权项 一种城市轨道交通随机客流加载方法,其特征在于,包括以下步骤:1.1、初始化,将一天平均划分为多个时间窗,并加载路网拓扑数据、列车时刻表、各时间窗下OD分布量,同时设定时间窗总数为n,设定标记时间窗的变量t的初始值为1,设定路段拥挤费用的初始值CN<sub>1</sub>=0;1.2、读取第t时间窗的OD分布量,基于基本有效路径集合,计算各路径的首末班车时间,判断该时间与第t时间窗的关系,生成动态有效路径集合RSET<sub>t</sub>;其中,基本有效路径集合的构建过程见步骤2.1‑2.4,而第t时间窗下动态有效路径集合RSET<sub>t</sub>的生成过程如下:路径的首班车时间大于第t时间窗的终止时刻或路径的末班车时间小于第t时间窗的起始时刻的路径不考虑;路径的首班车时间处于第t时间窗之间的路径按正常通行路径考虑;路径的末班车时间处于第t时间窗之间的路径以末班车时间为节点细化为小时间窗,并根据小时间窗所占时间窗的比重确定小时间窗下OD分布量,该路径在该小时间窗下算作正常通行路径;所有正常通行路径均为有效路径,且在第t时间窗下的有效路径组成的集合即为动态有效路径集合RSET<sub>t</sub>;1.3、基于动态有效路径集合RSET<sub>t</sub>和城市轨道交通乘客出行路径选择模型,计算第t时间窗下OD分布量在各有效路径上的分配比例,将第t时间窗下的OD分布量乘以该比例,即获得第t时间窗下各路径的路径流量;所述路径流量的计算公式如下:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>ql</mi><mi>k</mi><mrow><mi>r</mi><mi>s</mi></mrow></msubsup><mo>=</mo><mfrac><mrow><mi>exp</mi><mrow><mo>(</mo><msubsup><mi>V</mi><mi>k</mi><mrow><mi>r</mi><mi>s</mi></mrow></msubsup><mo>)</mo></mrow></mrow><mrow><mo>&Sigma;</mo><mi>exp</mi><mrow><mo>(</mo><msubsup><mi>V</mi><mi>k</mi><mrow><mi>r</mi><mi>s</mi></mrow></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>&times;</mo><msub><mi>q</mi><mrow><mi>r</mi><mi>s</mi></mrow></msub></mrow>]]></math><img file="FDA0001057232830000011.GIF" wi="494" he="143" /></maths>式中,<img file="FDA0001057232830000012.GIF" wi="77" he="63" />为从起点r到终点s的OD之间第k条有效路径的路径流量,q<sub>rs</sub>为从r到s的分时OD分布量,<img file="FDA0001057232830000013.GIF" wi="69" he="55" />为从r到s的OD之间第k条有效路径的路径选择效用的确定项,该确定项同步骤1.2中基本有效路径集合构建时所建立的城市轨道交通乘客出行路径选择效用函数的确定项;1.4、基于路径流量、路径的路段构成和列车时刻表,更新路段拥挤费用CN<sub>t+1</sub>,路段拥挤费用为路段满载率与路段乘车时间的乘积,即CN<sub>t+1</sub>=qplink<sub>t</sub>*TRlink,该值将作为下一次迭代时在步骤1.3中计算路径流量的输入数据;其中,CN<sub>t+1</sub>为在第t+1个时间窗时更新的路段拥挤费用,单位为小时;qplink<sub>t</sub>为第t时间窗的路段满载率;TRlink为路段的乘车时间,单位为小时;1.5、终止判断,如果t&lt;n,则t=t+1,从步骤1.2继续迭代;否则,该方法终止;所述步骤1.2中基本有效路径集合的构建包括以下步骤:2.1、构建城市轨道交通乘客出行路径选择效用函数,该函数的确定项的计算公式如下:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>k</mi><mrow><mi>r</mi><mi>s</mi></mrow></msubsup><mo>=</mo><mo>-</mo><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></msubsup><msub><mi>&beta;</mi><mi>i</mi></msub><msubsup><mi>X</mi><mrow><mi>i</mi><mi>k</mi></mrow><mrow><mi>r</mi><mi>s</mi></mrow></msubsup></mrow>]]></math><img file="FDA0001057232830000021.GIF" wi="454" he="95" /></maths>式中,<img file="FDA0001057232830000022.GIF" wi="68" he="62" />为从r到s的OD之间第k条路径的路径选择效用的确定项;<img file="FDA0001057232830000023.GIF" wi="76" he="63" />为从r到s的OD之间第k条路径的第i个特征属性值;β<sub>i</sub>为相应特征属性的参数;m为特征属性i的个数;上述的特征属性值不仅包含描述城市轨道交通服务水平的变量,即路径的乘车时间、票价、进出站走行时间、站台等车时间、换乘走行时间、换乘等车时间、换乘次数和拥挤费用,还包含描述路网拓扑和出行路径的非直线程度对乘客路径选择影响的变量,即角度费用变量,其中,路径拥挤费用的取值为该路径各组成路段的路段拥挤费用CN<sub>t</sub>之和除以路段条数;角度费用是对乘客出行路径每一路段走向偏离从起点到终点的直达方向的惩罚,所述角度费用的计算公式如下:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>AC</mi><mi>k</mi><mrow><mi>r</mi><mi>s</mi></mrow></msubsup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>L</mi><mi>p</mi></msub><mo>*</mo><mi>t</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>p</mi></msub><mo>/</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001057232830000024.GIF" wi="502" he="111" /></maths>式中:<img file="FDA0001057232830000025.GIF" wi="118" he="70" />为从r到s的OD之间第k条路径的角度费用;L<sub>p</sub>为该条路径的第p条路段的长度;θ<sub>p</sub>为该条路径的第p条路段偏离从起点到终点的直达方向的夹角,N为该条路径的路段条数;2.2、基于等效乘车时间系数,将其他因素均转化为等效乘车时间;2.3、确定乘客对时间的容忍阈值;2.4、利用Double‑sweep搜索算法确定基本有效路径集合。
地址 100044 北京市海淀区西直门外上园村3号