发明名称 一种综合成本优化的城市公共交通多模式站点停靠方法
摘要 本发明公开了一种综合成本优化的城市公共交通多模式站点停靠方法,该停靠方法包括以下步骤:第1步:采集乘客出行信息、车辆运营信息和沿线站点信息;第2步:首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数,然后,建立对应目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:第3步:将第1步中采集的乘客出行信息、车辆运营信息和沿线站点信息,代入第2步建立的目标函数和约束中,得到发车频率<i>h</i>和停站模式<i>O</i>,获得综合成本最低的多模式站点停靠方式。该方法考虑公交通行能力、客流、时间点、发车间隔、满载率约束以及停站模式,使得乘客出行成本和公交运营成本到达最低。
申请公布号 CN102831767B 申请公布日期 2014.07.09
申请号 CN201210286822.1 申请日期 2012.08.13
申请人 东南大学 发明人 陈峻;舒蕾;王昊
分类号 G08G1/00(2006.01)I 主分类号 G08G1/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,该停靠方法包括以下步骤:第1步:采集乘客出行信息、车辆运营信息和沿线站点信息;第2步:首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数,该目标函数如式(1)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>min</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><msub><mi>C</mi><mn>1</mn></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfenced open='(' close=')'><mtable><mtr><mtd><mfrac><mrow><msub><mi>r</mi><mi>j</mi></msub><msub><mi>h</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>h</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><mn>2</mn></mfrac><mo>+</mo></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>o</mi></msubsup><msub><mi>S</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>h</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>C</mi><mn>2</mn></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfenced open='{' close='}'><mtable><mtr><mtd><msub><mi>L</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>&theta;</mi><mi>i</mi></msub><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>o</mi></msubsup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>L</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>c</mi><mo>)</mo></mrow><mo>+</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><msub><mi>L</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>c</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mrow><mo>(</mo><msub><mi>L</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mi>u</mi></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><msub><mi>f</mi><mn>3</mn></msub><mo>=</mo><msub><mi>C</mi><mn>3</mn></msub><mo>&CenterDot;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>[</mo><msub><mi>&theta;</mi><mi>i</mi></msub><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><mo>&CenterDot;</mo><mi>c</mi><mo>+</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>+</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mi>u</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000450610010000011.GIF" wi="1529" he="695" /></maths>式(1)在式(1)中:z表示综合成本,min(z)表示最小综合成本,f<sub>1</sub>表示乘客出行成本,f<sub>2</sub>表示车辆运营成本,C<sub>1</sub>表示乘客单位等待时间成本,单位:元/min;C<sub>2</sub>表示乘客单位在车时间成本,单位:元/min;C<sub>3</sub>表示乘客单位运营成本,单位:元/min;i表示公交车辆的标记,i=1、2、…、M,M为大于等于3的整数,j表示公交车站的标记,j=1、2、…、N,N为大于等于3的整数,r<sub>j</sub>表示第j个公交车站的平均乘客到达率,人/min;h<sub>i,j</sub>表示经过公交车站j的公交车辆i‑1和公交车辆i之间的车头间距;<img file="FDA0000450610010000012.GIF" wi="74" he="75" />表示“0‑1”变量,对于第O种运营模式,公交车辆i在公交车站j停车时,<img file="FDA0000450610010000013.GIF" wi="77" he="72" />取1;公交车辆i在公交车站j不停时,<img file="FDA0000450610010000014.GIF" wi="75" he="75" />取0;其中,O表示公交车辆的停站模式,当O=1时,表示公交车辆站站停,当O=2时,表示公交车辆为区间车,当O=3时,表示公交车辆为大站快车模式;S<sub>i‑1,j</sub>表示公交车辆i‑1离开公交车站j时,公交车站剩余的乘客总数;L<sub>i,j</sub>表示公交车辆i离开公交车站j时的车上人数;θ<sub>i</sub>表示公交车辆i在线路上受到干扰修正系数;t<sub>i,j</sub>表示公交车辆i在公交车站j‑1和公交车站j之间的运行时间;t表示公交车辆进行开门动作和关门动作的总时间,单位:秒;c表示平均每辆公交车辆的加减速时间,单位:秒;A<sub>i,j</sub>表示公交车辆i到达公交车站j时的下车人数;B<sub>i,j</sub>表示公交车辆i到达公交车站j时的上车人数;u表示单个乘客平均上下车时间,单位:秒;然后,建立对应式(1)目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:公交通行能力的约束为:在研究时间段T内,每小时每个公交车站的公交车辆运送的乘客数必须不大于公交车站供该条线路使用的停车位的客运通行能力,如式(2)和式(3)所示,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA0000450610010000021.GIF" wi="301" he="143" /></maths>     式(2)<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>W</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA0000450610010000022.GIF" wi="320" he="135" /></maths>     式(3)式(2)中,P<sub>j</sub>为公交车站j的客运通行能力;式(3)中,W<sub>i,jk</sub>为公交车辆i运送从公交车站j到公交车站k的乘客数量;客流约束如式(4)至式(8)所示:L<sub>i,j</sub>=L<sub>i,j‑1</sub>+B<sub>i,j</sub>‑A<sub>i,j</sub>     式(4)<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>W</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow><mi>O</mi></msubsup></mrow>]]></math><img file="FDA0000450610010000023.GIF" wi="362" he="116" /></maths>     式(5)<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>W</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup></mrow>]]></math><img file="FDA0000450610010000024.GIF" wi="365" he="126" /></maths>     式(6)W<sub>i,jk</sub>=S<sub>i‑1,jk</sub>+r<sub>j</sub>h<sub>i</sub>     式(7)<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow></msub><mo>=</mo><msub><mi>S</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>jk</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow><mi>O</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>h</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>jk</mi></mrow><mi>O</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000450610010000025.GIF" wi="653" he="87" /></maths>     式(8)其中:L<sub>i,j‑1</sub>表示公交车辆i离开公交车站j‑1时的车上人数;k表示公交车站的标记,k=1、2、…、N,N为大于等于3的整数;W<sub>i,jk</sub>为公交车辆i运送从公交车站j到公交车站k的乘客数量;<img file="FDA0000450610010000031.GIF" wi="94" he="77" />表示“0‑1”变量,对于第O种运营模式,公交车辆i在公交车站j和公交车站k均停车时,<img file="FDA0000450610010000032.GIF" wi="88" he="77" />取1,否则,<img file="FDA0000450610010000033.GIF" wi="93" he="77" />取0;S<sub>i‑1,jk</sub>表示公交车辆i‑1离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;h<sub>i</sub>表示发车时间间隔,即第i‑1辆车和第i辆车在首站的发车时间间隔;r<sub>j,k</sub>表示需要乘坐到公交车站k的乘客在公交车站j的平均乘客到达率,1≤j≤k≤N,S<sub>i,jk</sub>为公交车辆i离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;S<sub>i,j</sub>为公交车辆i离开公交车站j时,公交车站j剩余的乘客数量;时间约束如式(9)至式(13)所示,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>a</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>d</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>&theta;</mi><mi>i</mi></msub><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mi>O</mi></msubsup><mo>+</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><mo>)</mo></mrow><mi>c</mi></mrow>]]></math><img file="FDA0000450610010000034.GIF" wi="667" he="94" /></maths>     式(9)<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>a</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><mo>&CenterDot;</mo><msub><mi>T</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000450610010000035.GIF" wi="353" he="85" /></maths>     式(10)h<sub>i,j</sub>=d<sub>i,j</sub>‑d<sub>i‑1,j</sub>     式(11)<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>h</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>h</mi><mi>i</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>j</mi></munderover><mo>[</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>+</mo><mi>c</mi><msub><mi>T</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>j</mi></munderover><mo>[</mo><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mi>O</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>+</mo><mi>c</mi><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000450610010000036.GIF" wi="739" he="112" /></maths>     式(12)O·h<sub>i</sub>=T     式(13)其中,a<sub>i,j</sub>为公交车辆i到达公交车站j的时刻;d<sub>i,j‑1</sub>为公交车辆i离开公交车站j‑1的时刻;d<sub>i‑1,j</sub>为公交车辆i‑1离开公交车站j的时刻;d<sub>i,j</sub>为公交车辆i离开公交车站j的时刻;<img file="FDA0000450610010000037.GIF" wi="105" he="78" />表示“0‑1”变量,对于第O种运营模式,公交车辆i在公交车站j‑1停车时,<img file="FDA0000450610010000038.GIF" wi="105" he="76" />取1;公交车辆i在公交车站j‑1不停时,<img file="FDA0000450610010000039.GIF" wi="104" he="76" />取0;<img file="FDA00004506100100000310.GIF" wi="103" he="76" />表示“0‑1”变量,对于第O种运营模式,公交车辆i‑1在公交车站j停车时,<img file="FDA00004506100100000311.GIF" wi="105" he="78" />取1;公交车辆i‑1在公交车站j不停时,<img file="FDA00004506100100000312.GIF" wi="102" he="78" />取0;T<sub>0</sub>为公交车辆停站时间,T<sub>0</sub>=2t+max(A<sub>i,j</sub>,B<sub>i,j</sub>)u,T为研究时间段,单位:min;发车间隔约束如式(14)所示:h<sub>min</sub>≤h<sub>i</sub>≤h<sub>max</sub>   式(14)其中:h<sub>min</sub>为最小的发车间隔,单位:min;h<sub>max</sub>为最大的发车间隔,单位:min,满载率约束如式(15)所示:L<sub>i,j</sub>/P<sub>max</sub>≤1   式(15)式中:P<sub>max</sub>为每辆公交车辆的额定最大载客量,人/车;第3步:将第1步中采集的乘客出行信息、车辆运营信息和沿线站点信息,代入第2步建立的目标函数和约束中,得到发车频率h和停站模式O,获得综合成本最低的多模式站点停靠方式。
地址 210096 江苏省南京市四牌楼2号