发明名称 一种基于时变客流的城市公交运营协调方法与装置
摘要 本发明公开了城市公交管理技术领域中的一种基于时变客流的城市公交运营协调方法与装置。该装置利用公交动态客流量和公交实时GPS数据,基于公交运营受实时路况随机扰动的实际,以运营商成本、拥挤里程比例和换乘乘客平均候车时间为目标,构建基于公交时变数据的多变量运营协调总表达式,通过基于快速分类的非支配遗传算法寻求最优的动态运营协调方案。本发明有利于改善城市公交系统的资源利用和有效衔接。
申请公布号 CN102044149B 申请公布日期 2012.08.08
申请号 CN201110005702.5 申请日期 2011.01.12
申请人 北京交通大学 发明人 陈旭梅;于雷;郭淑霞;李欣;荆长林
分类号 G08G1/00(2006.01)I 主分类号 G08G1/00(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 黄家俊
主权项 1.一种基于时变客流的城市公交运营协调方法,其特征是该方法包括以下步骤:步骤1:利用公交GPS数据对车辆行程时间进行预测,推算出车辆在换乘站的到站时间间隔;(1)基于指数平滑法的方法预测:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>&alpha;</mi><msub><mi>y</mi><mi>t</mi></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>t</mi></msub></mrow></math>]]></maths>式中:<img file="FDA0000151564210000012.GIF" wi="68" he="52" />为第t+1个行程时间的预测值;y<sub>t</sub>为第t个行程时间的实测值;<img file="FDA0000151564210000013.GIF" wi="40" he="52" />为第t个行程时间的预测值;α为修正系数,且0<α<1;(2)基于时间点指数法的方法预测:<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>t</mi></msub><mo>=</mo><msubsup><mi>y</mi><mi>t</mi><mo>*</mo></msubsup><mo>&times;</mo><msub><mi>I</mi><mi>t</mi></msub></mrow></math>]]></maths>式中:<img file="FDA0000151564210000015.GIF" wi="47" he="68" />为某线路第t个行程时间;<img file="FDA0000151564210000016.GIF" wi="55" he="71" />为第t个行程时间的趋势值;I<sub>t</sub>为t时段的修正系数;步骤2:根据动态客流量,确定公交运营的拥挤里程比例表达式;<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mn>2</mn><mi>t</mi></msubsup><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>Q</mi><mi>vk</mi><mi>t</mi></msubsup><mo>&times;</mo><msub><mi>l</mi><mi>vk</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>w</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>Q</mi><mi>vw</mi><mi>t</mi></msubsup><mo>&times;</mo><msub><mi>l</mi><mi>vw</mi></msub></mrow></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow></math>]]></maths>式中:<img file="FDA0000151564210000021.GIF" wi="47" he="54" />为t时段内公交线路的拥挤里程比例,以百分比表示;k,w为断面编号;n,N分别为拥挤断面和线路总的站点个数;m为线路数;<img file="FDA0000151564210000022.GIF" wi="63" he="54" />为t时段第v条公交线路第k个满载率超过阈值上限的断面客流量,单位人次;<img file="FDA0000151564210000023.GIF" wi="68" he="53" />为t时段第v条公交线路第w个断面的客流量,单位人次;l<sub>vk</sub>为第v条公交线路第k个断面所属的站间距,单位公里;l<sub>vw</sub>为第v条公交线路第w个断面所属的站间距,单位公里;步骤3:根据步骤1得到的车辆在换乘站的到站时间间隔,建立换乘乘客平均候车时间表达式;<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mn>3</mn><mi>t</mi></msubsup><mo>=</mo><msubsup><mi>T</mi><mrow><mi>r</mi><mo>-</mo><mi>b</mi></mrow><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>T</mi><mrow><mi>b</mi><mo>-</mo><mi>r</mi></mrow><mi>t</mi></msubsup></mrow></math>]]></maths>其中:<img file="FDA0000151564210000025.GIF" wi="52" he="62" />为换乘乘客平均候车时间;<img file="FDA0000151564210000026.GIF" wi="78" he="62" />为t时段轨道交通换乘常规公交平均候车时间;<img file="FDA0000151564210000027.GIF" wi="77" he="61" />为t时段常规公交换乘轨道交通候车时间;步骤4:对运营商成本进行分析,进而完成公交调度协调总表达式的建立;模型目标一:<maths num="0005"><![CDATA[<math><mrow><mi>min</mi><msubsup><mi>Z</mi><mn>1</mn><mi>t</mi></msubsup><mo>:</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mn>1</mn><mi>t</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&mu;</mi><mi>vr</mi></msub><mo>&times;</mo><msub><mi>l</mi><mi>vr</mi></msub><mo>&times;</mo><msubsup><mi>n</mi><mi>vr</mi><mi>t</mi></msubsup></mrow></math>]]></maths>模型目标二:<maths num="0007"><![CDATA[<math><mrow><mi>min</mi><msubsup><mi>Z</mi><mn>2</mn><mi>t</mi></msubsup><mo>:</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mn>2</mn><mi>t</mi></msubsup><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>Q</mi><mi>vk</mi><mi>t</mi></msubsup><mo>&times;</mo><msub><mi>l</mi><mi>vk</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>w</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>Q</mi><mi>vw</mi><mi>t</mi></msubsup><mo>&times;</mo><msub><mi>l</mi><mi>vw</mi></msub></mrow></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow></math>]]></maths>模型目标三:<maths num="0009"><![CDATA[<math><mrow><mi>min</mi><msubsup><mi>Z</mi><mn>3</mn><mn>1</mn></msubsup><mo>:</mo></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mn>3</mn><mi>t</mi></msubsup><mo>=</mo><msubsup><mi>T</mi><mrow><mi>r</mi><mo>-</mo><mi>b</mi></mrow><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>T</mi><mrow><mi>b</mi><mo>-</mo><mi>r</mi></mrow><mi>t</mi></msubsup></mrow></math>]]></maths>约束条件:<img file="BDA0000043582210000162.GIF" wi="945" he="622" />因变量:<img file="FDA0000151564210000036.GIF" wi="46" he="54" />为t时段运营商成本,单位元;<img file="FDA0000151564210000037.GIF" wi="47" he="54" />为t时段拥挤里程比例,以百分比表示;<img file="FDA0000151564210000038.GIF" wi="46" he="53" />为换乘乘客平均候车时间,单位秒;自变量:<img file="FDA0000151564210000039.GIF" wi="53" he="53" />为常规公交第v条线路r方向t时段内的发车数,单位辆;参数:r为公交线路走行方向,r=1为上行方向,r=2为下行方向;μ<sub>ir</sub>为公交第i条线路单位距离的成本,单位元/公里;l<sub>vr</sub>为公交第v条线路r方向的线路长度,单位公里;<img file="FDA00001515642100000310.GIF" wi="59" he="55" />为公交第i条线路上行方向t时段内场站内车辆数,单位辆;<img file="FDA0000151564210000041.GIF" wi="64" he="53" />为公交第i条线路下行方向t时段到达上行场站内的车辆数,单位辆;k,w为公交线路断面编号;n为拥挤断面个数;m为协调优化的公交线路数;N为公交线路总的站点数;i,j为站点编号;<img file="FDA0000151564210000042.GIF" wi="42" he="60" />为t时段从站点i上车到站点j下车的客流量,单位人;<img file="FDA0000151564210000043.GIF" wi="62" he="53" />为t时段第v条线路第k个拥挤断面客流量,单位人;<img file="FDA0000151564210000044.GIF" wi="69" he="53" />为t时段第v条线路第w个断面客流量,单位人;l<sub>vk</sub>为第v条线路第k个断面所属的站间距,单位公里;v为公交线路编号;<img file="FDA0000151564210000045.GIF" wi="134" he="69" />为t时段从站点i上车到站点j下车的客流量所能接受的拥挤数;<img file="FDA0000151564210000046.GIF" wi="48" he="54" />为t时段某条公交线路第k个断面的状态,为0-1变量;步骤5:采用基于快速分类的非支配遗传算法对步骤4的公交调度协调总表达式求解,输出所选择线路的运营协调方案;步骤6:通过网络将运营协调方案发送至公共交通网络调度指挥中心。
地址 100044 北京市海淀区西直门外上园村3号北京交通大学科技处