CN102831767A - 一种综合成本优化的城市公共交通多模式站点停靠方法 - Google Patents
一种综合成本优化的城市公共交通多模式站点停靠方法 Download PDFInfo
- Publication number
- CN102831767A CN102831767A CN2012102868221A CN201210286822A CN102831767A CN 102831767 A CN102831767 A CN 102831767A CN 2012102868221 A CN2012102868221 A CN 2012102868221A CN 201210286822 A CN201210286822 A CN 201210286822A CN 102831767 A CN102831767 A CN 102831767A
- Authority
- CN
- China
- Prior art keywords
- bus stop
- transit vehicle
- formula
- passenger
- public transit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Traffic Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种综合成本优化的城市公共交通多模式站点停靠方法,该停靠方法包括以下步骤:第1步:采集乘客出行信息、车辆运营信息和沿线站点信息;第2步:首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数,然后,建立对应目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:第3步:将第1步中采集的乘客出行信息、车辆运营信息和沿线站点信息,代入第2步建立的目标函数和约束中,得到发车频率h和停站模式O,获得综合成本最低的多模式站点停靠方式。该方法考虑公交通行能力、客流、时间点、发车间隔、满载率约束以及停站模式,使得乘客出行成本和公交运营成本到达最低。
Description
技术领域
本发明属于城市公共交通领域,具体来说,涉及一种综合成本优化的城市公共交通多模式站点停靠方法。
背景技术
公共交通因其自身大容量的特点成为缓解城市拥堵的必要措施,受到各个城市的青睐。公交停站模式的合理选择能够在一定程度上提升公交运营质量,对减少乘客出行成本、降低公交公司运营成本、缓解运力与运量之间的矛盾、改善城市交通环境和缓解客流拥挤等方面具有显著的作用和相当重要的意义。
公交停站模式主要包括站站停、大站快车和区间车三种,如附图1所示,其中:站站停模式提供正常的公交上下客服务,大站快车和区间车为长距离出行和重点区域内的乘客出行服务,公交停站次数减少,运营速度提高,车辆周转时间缩短。目前常规公交多采用站站停模式,在少数线路辅助区间车模式,多模式停站未成为常态。该种模式带来公交车站延误尤为明显,乘客出行等待时间长,影响居民选择公交出行。若将三种模式组合,对于乘客,不同出行目的的乘客得到分类服务,降低互相影响,减少出行中的停留时间,缩短总出行时间;对于运营企业,停站次数减少,加减速频率降低,有效减少尾气排放,达到较高的燃油使用效率,降低***运营成本。
发明内容
技术问题:本发明所要解决的技术问题是:提供一种综合成本优化的城市公共交通多模式站点停靠方法,该方法考虑公交通行能力、客流、时间点、发车间隔、满载率约束以及停站模式,使得乘客出行成本和公交运营成本到达最低。
技术方案:为解决上述技术问题,本发明采用的综合成本优化的城市公共交通多模式站点停靠方法,该停靠方法包括以下步骤:
第1步:采集乘客出行信息、车辆运营信息和沿线站点信息;
第2步:首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数,该目标函数如式(1)所示:
在式(1)中:z表示综合成本,min(z)表示最小综合成本,f1表示乘客出行成本,f2表示车辆运营成本,C1表示乘客单位等待时间成本,单位:元/min;C2表示乘客单位在车时间成本,单位:元/min;C3表示乘客单位运营成本,单位:元/min;
i表示公交车辆的标记,i=1、2、…、M,M为大于等于3的整数,
j表示公交车站的标记,j=1、2、…、N,N为大于等于3的整数,
rj表示第j个公交车站的平均乘客到达率,人/min;
hi,j表示经过公交车站j的公交车辆i-1和公交车辆i之间的车头间距;
表示“0-1”变量,对于第O种运营模式,公交车辆i在公交车站j停车时,取1;公交车辆i在公交车站j不停时,取0;其中,O表示公交车辆的停站模式,当O=1时,表示公交车辆站站停,当O=2时,表示公交车辆为区间车,当O=3时,表示公交车辆为大站快车模式;
Si-1,j表示公交车辆i-1离开公交车站j时,公交车站剩余的乘客总数;
Li,j表示公交车辆i离开公交车站j时的车上人数;
θi表示公交车辆i在线路上受到干扰修正系数;
ti,j表示公交车辆i在公交车站j-1和公交车站j之间的运行时间;
t表示公交车辆进行开门动作和关门动作的总时间,单位:秒;
c表示平均每辆公交车辆的加减速时间,单位:秒;
Ai,j表示公交车辆i到达公交车站j时的下车人数;
Bi,j表示公交车辆i到达公交车站j时的上车人数;
u表示单个乘客平均上下车时间,单位:秒;
然后,建立对应式(1)目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:
公交通行能力的约束为:在研究时间段T内,每小时每个公交车站的公交车辆运送的乘客数必须不大于公交车站供该条线路使用的停车位的客运通行能力,如式(2)和式(3)所示,
式(2)中,Pj为公交车站j的客运通行能力;式(3)中,Wi,jk为公交车辆i运送从公交车站j到公交车站k的乘客数量;
客流约束如式(4)至式(8)所示:
Li,j=Li,j-1+Bi,j-Ai,j 式(4)
Wi,jk=Si-1,jk+rjhi 式(7)
其中:Li,j-1表示公交车辆i离开公交车站j-1时的车上人数;
k表示公交车站的标记,k=1、2、…、N,N为大于等于3的整数;
Wi,jk为公交车辆i运送从公交车站j到公交车站k的乘客数量;
Si-1,jk表示公交车辆i-1离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;
hi表示发车时间间隔,即第i-1辆车和第i辆车在首站的发车时间间隔;
rj,k表示需要乘坐到公交车站k的乘客在公交车站j的平均乘客到达率,1≤j≤k≤N,
Si,jk为公交车辆i离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;
Si,j为公交车辆i离开公交车站j时,公交车站j剩余的乘客数量;
时间约束如式(9)至式(13)所示,
hi,j=di,j-di-1,j 式(11)
O·hi=T 式(13)
其中,ai,j为公交车辆i到达公交车站j的时刻;
di,j-1为公交车辆i离开公交车站j-1的时刻;
di-1,j为公交车辆i-1离开公交车站j的时刻;
di,j为公交车辆i离开公交车站j的时刻;
T0为公交车辆停站时间,T0=2t+max(Ai,j,Bi,j)u,
T为研究时间段,单位:min;
发车间隔约束如式(14)所示:
hmin≤hi≤hmax 式(14)
其中:hmin为最小的发车间隔,单位:min;hmax为最大的发车间隔,单位:min,
满载率约束如式(15)所示:
Li,j/Pmax≤1 式(15)
式中:Pmax为每辆公交车辆的额定最大载客量,人/车;
第3步:将第1步中采集的乘客出行信息、车辆运营信息和沿线站点信息,代入第2步建立的目标函数和约束中,得到发车频率h和停站模式O,获得综合成本最低的多模式站点停靠方式。
有益效果:本发明充分考虑停站模式的选择对乘客出行以及车辆运营两个方面的综合成本的影响,并且将公交站点通行能力作为其中的一个约束和检验条件考虑,弥补可能因多模式停站导致通行能力不足进而引发公交延误的缺陷。本发明提出的综合成本优化的城市公共交通多模式中途站点停靠方法,以综合成本最低为目标,进行某时段内城市公共交通多模式中途站点停靠控制,进行站站停、大站快车和区间车的组合优化。本发明可以在调查客流,乘客及站点信息的基础上,通过建立综合成本最低的优化模型,并应用遗传算法进行求解,给出最优的发车频率以及停站模式的组合选择。实例表明,在客流特征明显时采用多种模式组合的停站模式,优化后成本较优化前有明显降低、停站时间显著减少、延误减少,具有一定的可行性和合理性。
附图说明
图1为本发明的流程图。
图2为本发明中的三种公交车辆的停站模式示意图。
图3为本发明实施例中北京公交运通105线路车站乘客上、下车人数图。
具体实施方式
下面结合附图,对本发明的技术方案进行详细的说明。
如图1所示,本发明的一种综合成本优化的城市公共交通多模式站点停靠方法,包括以下步骤:
第1步:采集乘客出行信息、车辆运营信息和沿线站点信息.
在第1步中,乘客出行信息包含在各个公交车站,上车人数Bi,j、下车人数Ai,j、乘客单位等待时间成本C1、乘客单位在车时间成本C2、乘客单位运营成本C3和单个乘客平均上下车时间u。其中,上车人数Bi,j、下车人数Ai,j和单个乘客平均上下车时间u,采集方法均可采用现有技术,包含摄像法、跟车调查法等,详细内容可参阅参考文献(例如,王炜、过秀成著《交通工程学》)。乘客单位等待时间成本C1和乘客单位在车时间成本C2可以利用“生产法”进行计算,即时间价值C为年度人均GDP除以每年工作时间,人均GDP可以通过查阅当地年鉴获得。乘客单位运营成本C3可以根据公交公司调查公布的数据获得。
车辆运营信息包含有站间运行时间ti,j和车辆加减速时间c。站间运行时间ti,j和车辆加减速时间c采用平均值。这两个参数采集方法可采用现有技术,包含基于GPS的信息采集方法、基于移动通信的信息采集方法,基于问卷调查的信息采集技术,摄像法、跟车调查法等,详细内容可参阅参考文献(例如,王炜、过秀成著《交通工程学》)。
沿线站点信息包含站点总数N、绿信比g/c、清空时间tc、平均停靠时间td和每车额定最大载客量Pmax。其中,站点总数N可以查阅公交公司公布的数据获得。每车额定最大载客量Pmax根据不同车型查阅设备说明获得。绿信比g/c、清空时间tc和平均停靠时间td可以通过现有技术中的摄像法、基于GPS的信息采集方法等获得,具体可参阅参考文献(例如:吴兵和李晔著《交通管理和控制》,王炜和过秀成著《交通工程学》)。
第2步:首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数,该目标函数如式(1)所示:
在式(1)中:z表示综合成本,min(z)表示最小综合成本,f1表示乘客出行成本,f2表示车辆运营成本,C1表示乘客单位等待时间成本,单位:元/min;C2表示乘客单位在车时间成本,单位:元/min;C3表示乘客单位运营成本,单位:元/min;
i表示公交车辆的标记,i=1、2、…、M,M为大于等于3的整数,
j表示公交车站的标记,j=1、2、…、N,N为大于等于3的整数,
rj表示第j个公交车站的平均乘客到达率,人/min;
hi,j表示经过公交车站j的公交车辆i-1和公交车辆i之间的车头间距;
表示“0-1”变量,对于第O种运营模式,公交车辆i在公交车站j停车时,取1;公交车辆i在公交车站j不停时,取0;其中,O表示公交车辆的停站模式,当O=1时,表示公交车辆站站停,当O=2时,表示公交车辆为区间车,当O=3时,表示公交车辆为大站快车模式;
Si-1,j表示公交车辆i-1离开公交车站j时,公交车站剩余的乘客总数;
Li,j表示公交车辆i离开公交车站j时的车上人数;
θi表示公交车辆i在线路上受到干扰修正系数;
ti,j表示公交车辆i在公交车站j-1和公交车站j之间的运行时间;
t表示公交车辆进行开门动作和关门动作的总时间,单位:秒;
c表示平均每辆公交车辆的加减速时间,单位:秒;
Ai,j表示公交车辆i到达公交车站j时的下车人数;
Bi,j表示公交车辆i到达公交车站j时的上车人数;
u表示单个乘客平均上下车时间,单位:秒;
然后,建立对应式(1)目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:
公交通行能力的约束为:在研究时间段T内,每小时每个公交车站的公交车辆运送的乘客数必须不大于公交车站供该条线路使用的停车位的客运通行能力,如式(2)和式(3)所示,
式(2)中,Pj为公交车站j的客运通行能力;式(3)中,Wi,jk为公交车辆i运送从公交车站j到公交车站k的乘客数量;
客流约束如式(4)至式(8)所示:
Li,j=Li,j-1+Bi,j-Ai,j 式(4)
Wi,jk=Si-1,jk+rjhi 式(7)
其中:Li,j-1表示公交车辆i离开公交车站j-1时的车上人数;
k表示公交车站的标记,k=1、2、…、N,N为大于等于3的整数;
Wi,jk为公交车辆i运送从公交车站j到公交车站k的乘客数量;
Si-1,jk表示公交车辆i-1离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;
hi表示发车时间间隔,即第i-1辆车和第i辆车在首站的发车时间间隔;
rj,k表示需要乘坐到公交车站k的乘客在公交车站j的平均乘客到达率,1≤j≤k≤N;
Si,jk为公交车辆i离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;
Si,j为公交车辆i离开公交车站j时,公交车站j剩余的乘客数量;
时间约束如式(9)至式(13)所示,
hi,j=di,j-di-1,j 式(11)
O·hi=T 式(13)
其中,ai,j为公交车辆i到达公交车站j的时刻;
di,j-1为公交车辆i离开公交车站j-1的时刻;
di-1,j为公交车辆i-1离开公交车站j的时刻;
di,j为公交车辆i离开公交车站j的时刻;
T0为公交车辆停站时间,T0=2t+max(Ai,j,Bi,j)u,
T为研究时间段,单位:min;
发车间隔约束如式(14)所示:
hmin≤hi≤hmax 式(14)
其中:hmin为最小的发车间隔,单位:min;hmax为最大的发车间隔,单位:min,
满载率约束如式(15)所示:
Li,j/Pmax≤1 式(15)
式中:Pmax为每辆公交车辆的额定最大载客量,人/车。
在步骤2中,t、c和u均为常数。O表示公交车辆的停站模式,如图2所示。当O=1时,表示公交车辆站站停,对应图2中的上图。当O=2时,表示公交车辆为区间车,对应图2中的下图。当O=3时,表示公交车辆为大站快车模式,对应图2中间的图。θi表示公交车辆i在线路上受到干扰修正系数,包括车辆之间的互相影响和拥堵,对运行时间产生影响。公交通行能力采用公交通行能力和服务质量手册(美国交通运输研究委员会编著,杨晓光,腾靖等译.公共交通通行能力和服务质量手册(原著第二版)[M].北京:中国建筑工业出版社,2010.15-16.)中的方法测算,具体如下:
每小时公交停靠站的通行能力(BS)为:
其中:Bs表示停靠站的公交通行能力,单位:辆公交车/小时;
Bl表示单个车位的公交车通行能力,单位:辆公交车/小时;
Ncl表示有效的公交车位数;
3600表示代表单位时间,即1小时;
g/c表示绿信比(有效绿灯时间占总信号周期的比例);
tc表示清空时间,单位:秒;
td表示平均停靠时间,单位:秒;
Z表示满足期望进站失败率的标准正态变量;
cv表示停靠时间波动系数。
单条公交线路一般占用一个固定的停车位,取有效车位数Ncl为1进行计算,每小时公交车站的客运通行能力则可以通过以下公式进行计算:
P=min{Pmaxf(PHF),PmaxBs(PHF)}
式中:P表示客运通行能力,单位:人/小时;
Pmax表示每车额定最大载客量,单位:人/车;
f表示额定发车频率,单位:车/小时;
PHF表示高峰小时系数。
第3步:将第1步中采集的乘客出行信息、车辆运营信息和沿线站点信息,代入第2步建立的目标函数和约束中,得到发车频率h和停站模式O,获得综合成本最低的多模式站点停靠方式。
在第3步中,对目标函数和约束的求解方法采用遗传算法、搜索算法,或者下降法。这些算法为现有技术。例如可以参阅由邢文训,谢金星著的《现代优化计算方法》。对目标函数和约束的求解采用遗传算法,可以参阅参考文献(王小平,曹立明《遗传算法——理论、应用与软件实现》西安:西安交通大学出版社,2002.)遗传算法具体包括以下各步骤:
31)设置参数:分别设置种群大小、选择概率、交叉概率、变异概率、循环次数和发车次数的最小值;
32)进行初始化:采用完全随机的方法产生初始种群,即在优化过程中随机选取n个点作为初始解,对各公交车站选用停站模式的乘客进行二进制编码;
33)计算目标函数值:根据初始化产生的发车次数,由时间段除以发车次数得到发车频率hi,根据已知参数T0和c,由时间约束中公式(12)测得hi,j;已知车站乘客到达率情况rj,k和初始车站剩余乘客数,由客流约束条件中公式(7)测得运送乘客的数量;根据式(5)和式(6)测得上、下车乘客数量,从而确定目标函数值;
34)计算个体适应值:根据步骤33)测算的目标函数值,采用适应转换方法在适应值函数上加上惩罚函数,计算出个体适应值;
35)按照步骤34)计算的个体适应值大小对种群进行适应度由大到小的排序,然后采用轮盘赌法选择进入下一子代种群的个体;
36)按照交叉概率和变异概率对染色体进行交叉、变异操作,交叉采用双亲遗传法,变异通过赋予每一个基因相对较小的变异概率Pm完成;
37)重复步骤33)-步骤36),直到达到设置的循环次数,从而完成一个发车次数的循环;
38)将发车次数加1,重复步骤33)-步骤37),产生模型最优解。
作为优选,在步骤31)中,交叉概率Pm=0.8;变异概率Pc=0.005;种群大小取20;循环次数取100,发车次数的最小值取3。
在步骤35)中,轮盘赌法为现有技术,具体可以参阅玄光男,陈润伟著的《遗传算法与优化工程》。轮盘赌法按照以下步骤实施:首先计算个体适应值并将种群中个体适应值叠加,得到总适应值其中,fx表示个体适应值,F表示总适应值,X为种群中的个体个数;其次,用每个个体的适应值除以总适应度值得到个体被选择的概率,即其中,px表示概率;再次,计算个体的累计概率构造一个轮盘;最后,进行轮盘选择即产生一个[0,1]区间内的随机数,若小于或等于个体x的累积概率且大于个体x-1的累积概率,选择个体x进入子代种群。
为了验证本发明方法的有效性,采用本发明的方法对北京运通105线路进行综合成本优化的多模式中途站点停靠设置,具体步骤如下:
第1步:采集乘客出行信息、车辆运营信息和沿线站点信息,具体如下:
运通105线路沟通南北,穿越中关村科技园区等大型电子卖场,途径北京交通大学、北京理工大学、人民大学、清华大学、北京大学和北京体育大学等高校,路过圆明园景点,线路里程15.495km,途径25个站点。现行发车间隔是12min,不区分平峰和高峰。采用车辆拥有座位数39个,额定最大载客量为60人,高峰小时系数取0.85。
研究时间段T为晚高峰时段1h;线路上车站总数N=25;乘客单位等待时间成本C1=0.4元/min;乘客单位在车时间成本C2=0.2元/min;乘客单位运营成本C3=1元/min;车辆加减速时间c=40s;单个乘客平均上下车时间u=5s;t忽略不计。
关键站点通行能力相关参数如表2所示。
表2关键站点通行能力相关参数
站点 | g(s) | C(s) | tc(s) | td(s) | Cv | Za | g/c |
中关村南 | 75 | 140 | 10 | T0 | 0.6 | 1.44 | 0.54 |
海淀黄庄北 | 83.8 | 160 | 8.5 | T0 | 0.6 | 1.44 | 0.52 |
人民大学 | 83.8 | 160 | 5 | T0 | 0.6 | 1.44 | 0.52 |
各个站点乘客上下车人数如图3所示。根据实际测得的每个车站上下车人数,按比例分配,得到站点乘客到达率。
线路考虑干扰系数,根据实测情况的包括站点距离和运行时间如表3所示。
表3线路基本情况
不同停站模式的车辆根据各站点客流情况确定的车站停靠方式如表4所示。
表4不同车辆停站方式
第2步:建立综合成本最优的规划模型:
首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数:
然后,建立对应式(1)目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:
按照遗传算法对模型求解,具体步骤如下:
31)设置参数:种群大小取20、交叉概率Pm=0.8、变异概率Pc=0.005、循环次数取100、发车次数的最小值为3,选择概率在轮盘赌法中计算获得。
32)进行初始化,采用完全随机的方法产生初始种群,即在优化过程中随机选取n个点作为初始解,对各公交车站选用停站模式的乘客进行二进制编码。采用可变长度的遗传算法编码,停站模式组合总的长度为发车次数,时间段确定根据发车次数即可得出发车间隔。
33)计算目标函数值:根据初始化产生的发车次数,由时间段除以发车次数得到发车频率hi,根据已知参数T0和c,由时间约束中公式(12)测得hi,j;已知车站乘客到达率情况rj,k和初始车站剩余乘客数,由客流约束条件中公式(7)测得运送乘客的数量;根据式(5)和式(6)测得上、下车乘客数量,从而确定目标函数值;
34)计算个体适应值。根据目标函数值,采用适应转换方法在适应值函数上加上惩罚函数,计算出个体适应值;
35)按照步骤34)计算的个体适应值大小对种群进行适应度由大到小的排序,然后采用轮盘赌法选择进入下一子代种群的个体;
36)按照交叉概率和变异概率对染色体进行交叉、变异操作,交叉采用双亲遗传法,变异通过赋予每一个基因相对较小的变异概率Pc完成,
37)重复步骤33)-步骤36),直到达到设置的循环次数,从而完成一个发车次数的循环;
38)将发车次数加1,重复步骤33)-步骤37),产生模型最优解。
根据计算获得的最优解对结果进行分析:
计算得最终结果为发10次车,发车间隔为6min,5辆站站停模式,1辆区间车和4辆大站快车。得到结果如表6所示,1代表站站停模式,2代表区间车模式,3代表大站快车模式。
表5计算结果
优化之前每小时成本为23352,优化后目标值为11419,较优化前综合成本降低51.1%,优化效果明显。将优化前在中关村南站的停靠时间和优化后的时间比较,可以发现:优化前运通105发车间隔为12min,共有5辆车,平均停站时间为38s;优化后发车间隔为6min,共有10辆车,平均停站时间为21.1s,在车站停靠时间减少。优化前停站时间不平稳,时间波动较大;优化后停站时间较为平稳,多在20s-30s之间。
Claims (7)
1.一种综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,该停靠方法包括以下步骤:
第1步:采集乘客出行信息、车辆运营信息和沿线站点信息;
第2步:首先,建立基于乘客出行成本和车辆运营成本的综合成本最低的目标函数,该目标函数如式(1)所示:
在式(1)中:z表示综合成本,min(z)表示最小综合成本,f1表示乘客出行成本,f2表示车辆运营成本,C1表示乘客单位等待时间成本,单位:元/min;C2表示乘客单位在车时间成本,单位:元/min;C3表示乘客单位运营成本,单位:元/min;
i表示公交车辆的标记,i=1、2、…、M,M为大于等于3的整数,
j表示公交车站的标记,j=1、2、…、N,N为大于等于3的整数,
rj表示第j个公交车站的平均乘客到达率,人/min;
hi,j表示经过公交车站j的公交车辆i-1和公交车辆i之间的车头间距;
表示“0-1”变量,对于第O种运营模式,公交车辆i在公交车站j停车时,取1;公交车辆i在公交车站j不停时,取0;其中,O表示公交车辆的停站模式,当O=1时,表示公交车辆站站停,当O=2时,表示公交车辆为区间车,当O=3时,表示公交车辆为大站快车模式;
Si-1,j表示公交车辆i-1离开公交车站j时,公交车站剩余的乘客总数;
Li,j表示公交车辆i离开公交车站j时的车上人数;
θi表示公交车辆i在线路上受到干扰修正系数;
ti,j表示公交车辆i在公交车站j-1和公交车站j之间的运行时间;
t表示公交车辆进行开门动作和关门动作的总时间,单位:秒;
c表示平均每辆公交车辆的加减速时间,单位:秒;
Ai,j表示公交车辆i到达公交车站j时的下车人数;
Bi,j表示公交车辆i到达公交车站j时的上车人数;
u表示单个乘客平均上下车时间,单位:秒;
然后,建立对应式(1)目标函数的约束条件式,包含公交通行能力、客流、时间点、发车间隔及满载率五个方面的约束:
公交通行能力的约束为:在研究时间段T内,每小时每个公交车站的公交车辆运送的乘客数必须不大于公交车站供该条线路使用的停车位的客运通行能力,如式(2)和式(3)所示,
式(2)中,Pj为公交车站j的客运通行能力;式(3)中,Wi,jk为公交车辆i运送从公交车站j到公交车站k的乘客数量;
客流约束如式(4)至式(8)所示:
Li,j=Li,j-1+Bi,j-Ai,j 式(4)
Wi,jk=Si-1,jk+rjhi 式(7)
其中:Li,j-1表示公交车辆i离开公交车站j-1时的车上人数;
k表示公交车站的标记,k=1、2、…、N,N为大于等于3的整数;
Wi,jk为公交车辆i运送从公交车站j到公交车站k的乘客数量;
Si-1,jk表示公交车辆i-1离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;
hi表示发车时间间隔,即第i-1辆车和第i辆车在首站的发车时间间隔;
rj,k表示需要乘坐到公交车站k的乘客在公交车站j的平均乘客到达率,1≤j≤k≤N,
Si,jk为公交车辆i离开公交车站j时,公交车站j未上车的需要乘坐到k站的乘客数量;
Si,j为公交车辆i离开公交车站j时,公交车站j剩余的乘客数量;
时间约束如式(9)至式(13)所示,
hi,j=di,j-di-1,j 式(11)
O·hi=T 式(13)
其中,ai,j为公交车辆i到达公交车站j的时刻;
di,j-1为公交车辆i离开公交车站j-1的时刻;
di-1,j为公交车辆i-1离开公交车站j的时刻;
di,j为公交车辆i离开公交车站j的时刻;
T0为公交车辆停站时间,T0=2t+max(Ai,j,Bi,j)u,
T为研究时间段,单位:min;
发车间隔约束如式(14)所示:
hmin≤hi≤hmax 式(14)
其中:hmin为最小的发车间隔,单位:min;hmax为最大的发车间隔,单位:min,
满载率约束如式(15)所示:
Li,j/Pmax≤1 式(15)
式中:Pmax为每辆公交车辆的额定最大载客量,人/车;
第3步:将第1步中采集的乘客出行信息、车辆运营信息和沿线站点信息,代入第2步建立的目标函数和约束中,得到发车频率h和停站模式O,获得综合成本最低的多模式站点停靠方式。
2.根据权利要求1所述的综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,所述的第1步中,乘客出行信息包含在各个公交车站,上车人数Bi,j、下车人数Ai,j、乘客单位等待时间成本C1、乘客单位在车时间成本C2、乘客单位运营成本C3和单个乘客平均上下车时间u。
3.根据权利要求1所述的综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,所述的第1步中,车辆运营信息包含有站间运行时间ti,j和车辆加减速时间c。
4.根据权利要求1所述的综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,所述的第1步中,沿线站点信息包含站点总数N、绿信比g/c、清空时间tc、平均停靠时间td和每车额定最大载客量Pmax。
5.根据权利要求1所述的综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,所述的第3步中,对目标函数和约束的求解方法采用遗传算法、搜索算法,或者下降法。
6.根据权利要求5所述的综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,所述的第3步中,对目标函数和约束的求解采用遗传算法,具体包括以下各步骤:
31)设置参数:分别设置种群大小、选择概率、交叉概率、变异概率、循环次数和发车次数的最小值;
32)进行初始化:采用完全随机的方法产生初始种群,即在优化过程中随机选取n个点作为初始解,对各公交车站选用停站模式的乘客进行二进制编码;
33)计算目标函数值:根据初始化产生的发车次数,由时间段除以发车次数得到发车频率hi,根据已知参数T0和c,由时间约束中公式(12)测得hi,j;已知车站乘客到达率情况rj,k和初始车站剩余乘客数,由客流约束条件中公式(7)测得运送乘客的数量;根据式(5)和式(6)测得上、下车乘客数量,从而确定目标函数值;
34)计算个体适应值:根据步骤33)测算的目标函数值,采用适应转换方法在适应值函数上加上惩罚函数,计算出个体适应值;
35)按照步骤34)计算的个体适应值大小对种群进行适应度由大到小的排序,然后采用轮盘赌法选择进入下一子代种群的个体;
36)按照交叉概率和变异概率对染色体进行交叉、变异操作,交叉采用双亲遗传法,变异通过赋予每一个基因相对较小的变异概率Pm完成;
37)重复步骤33)-步骤36),直到达到设置的循环次数,从而完成一个发车次数的循环;
38)将发车次数加1,重复步骤33)-步骤37),产生模型最优解。
7.根据权利要求6所述的综合成本优化的城市公共交通多模式站点停靠方法,其特征在于,所述的步骤31)中,交叉概率Pm=0.8;变异概率Pc=0.005;种群大小取20;循环次数取100,发车次数的最小值取3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210286822.1A CN102831767B (zh) | 2012-08-13 | 2012-08-13 | 一种综合成本优化的城市公共交通多模式站点停靠方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210286822.1A CN102831767B (zh) | 2012-08-13 | 2012-08-13 | 一种综合成本优化的城市公共交通多模式站点停靠方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102831767A true CN102831767A (zh) | 2012-12-19 |
CN102831767B CN102831767B (zh) | 2014-07-09 |
Family
ID=47334872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210286822.1A Expired - Fee Related CN102831767B (zh) | 2012-08-13 | 2012-08-13 | 一种综合成本优化的城市公共交通多模式站点停靠方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102831767B (zh) |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104464293A (zh) * | 2014-12-15 | 2015-03-25 | 大连理工大学 | 一种公交专用道上公交停靠站泊位数确定方法 |
CN104504638A (zh) * | 2014-12-18 | 2015-04-08 | 南京信息工程大学 | 一种公交车发车时刻表制定方法 |
CN105185105A (zh) * | 2015-06-26 | 2015-12-23 | 东南大学 | 基于车辆gps和公交ic卡数据的公交换乘识别方法 |
CN105448082A (zh) * | 2015-12-30 | 2016-03-30 | 清华大学 | 一种可变间隔发车的快速公交车组合调度方法 |
CN106056909A (zh) * | 2016-06-29 | 2016-10-26 | 东南大学 | 一种中小城市快速公交阈值确定方法 |
CN106297273A (zh) * | 2016-09-29 | 2017-01-04 | 百度在线网络技术(北京)有限公司 | 班车路线的处理方法及装置 |
CN106781432A (zh) * | 2016-11-16 | 2017-05-31 | 青岛海信网络科技股份有限公司 | 一种公交车辆调度的方法及装置 |
CN107092976A (zh) * | 2017-03-28 | 2017-08-25 | 东南大学 | 一种多目标模型协同优化多条公交线路发车间隔的方法 |
CN107516147A (zh) * | 2017-08-10 | 2017-12-26 | 中南大学 | 一种高速铁路本线列车开行方案优化方法及其*** |
CN107766969A (zh) * | 2017-09-28 | 2018-03-06 | 东南大学 | 一种基于地铁服务能力瓶颈区段识别的大站快线布设方法 |
CN108009685A (zh) * | 2017-12-07 | 2018-05-08 | 厦门理工学院 | 一种用于城市公共交通的智能移动终端 |
CN108648441A (zh) * | 2018-06-05 | 2018-10-12 | 武汉理工大学 | 基于站点群的动态公交运行方法 |
CN109544969A (zh) * | 2019-01-09 | 2019-03-29 | 张鸿青 | 公交车智能跳站运行方法及装置 |
CN110046845A (zh) * | 2019-03-06 | 2019-07-23 | 青岛日日顺物流有限公司 | 一种场区车辆装载方法和装载*** |
CN110288829A (zh) * | 2019-06-19 | 2019-09-27 | 河海大学 | 一种港湾式公交停靠站停车泊位数量设置判别方法 |
CN110533227A (zh) * | 2019-08-08 | 2019-12-03 | 东南大学 | 一种确定可变线路式公交固定站点和动态站点的方法 |
CN110782694A (zh) * | 2019-10-31 | 2020-02-11 | 同济大学 | 一种网联公交车到站停靠管理优化方法 |
CN110796877A (zh) * | 2019-05-14 | 2020-02-14 | 广州学塾加软件科技有限公司 | 一种面向单向公交线路的交通信号控制与公交调度协同控制方法 |
CN110941691A (zh) * | 2020-02-19 | 2020-03-31 | 武汉元光科技有限公司 | 一种顺风车搭乘撮合方法及装置 |
CN110969851A (zh) * | 2019-12-09 | 2020-04-07 | 浙江顿悟科技有限公司 | 一种智慧城市公交*** |
CN111063209A (zh) * | 2019-12-23 | 2020-04-24 | 北京航空航天大学 | 配合区间车的自动驾驶公交组合调度优化方法 |
CN111667087A (zh) * | 2019-03-08 | 2020-09-15 | 南京农业大学 | 一种考虑污染排放的公交跳站运行方法 |
CN112101751A (zh) * | 2020-08-31 | 2020-12-18 | 华南理工大学 | 一种公共交通优化方法、***、装置及存储介质 |
CN112509357A (zh) * | 2020-11-19 | 2021-03-16 | 北京清研宏达信息科技有限公司 | 一种基于遗传算法的公交发车时间间隔优化方法 |
CN113781787A (zh) * | 2021-11-15 | 2021-12-10 | 深圳市都市交通规划设计研究院有限公司 | 公交发车时刻表生成方法及*** |
CN114758522A (zh) * | 2022-04-25 | 2022-07-15 | 驭势科技(北京)有限公司 | 一种车辆的月台停靠方法、装置、设备及介质 |
CN114897445A (zh) * | 2022-07-12 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | 公交线网的停靠点调整优化方法、装置及可读存储介质 |
CN115424437A (zh) * | 2022-08-22 | 2022-12-02 | 厦门筑成信创城市规划设计有限公司 | 站点响应式公交调度方法及装置 |
CN115457793A (zh) * | 2022-08-31 | 2022-12-09 | 东风悦享科技有限公司 | 一种应用于公交调度中始发站非均匀发车的方法及*** |
CN115860298A (zh) * | 2022-12-05 | 2023-03-28 | 四川警察学院 | 公交中间站停靠线路优化方法、装置、设备及存储介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110246356B (zh) * | 2019-06-12 | 2021-05-04 | 东南大学 | 一种确定可变线路公交配置车辆数和类型的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090254405A1 (en) * | 2008-04-08 | 2009-10-08 | Benjamin Leslie Hollis | Simultaneous vehicle routing, vehicle scheduling, and crew scheduling |
CN101572011A (zh) * | 2009-06-10 | 2009-11-04 | 上海理工大学 | 城市公交智能化调度管理***及方法 |
CN101630440A (zh) * | 2009-06-01 | 2010-01-20 | 北京交通大学 | 常规公共交通衔接城市轨道交通的运营协调优化方法与*** |
CN101799981A (zh) * | 2010-02-09 | 2010-08-11 | 华南理工大学 | 多模式公共交通区域调度控制方法 |
CN101807222A (zh) * | 2010-02-09 | 2010-08-18 | 武汉大学 | 基于站点的城市公交线网优化配置方法 |
CN102044149A (zh) * | 2011-01-12 | 2011-05-04 | 北京交通大学 | 一种基于时变客流的城市公交运营协调方法与装置 |
-
2012
- 2012-08-13 CN CN201210286822.1A patent/CN102831767B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090254405A1 (en) * | 2008-04-08 | 2009-10-08 | Benjamin Leslie Hollis | Simultaneous vehicle routing, vehicle scheduling, and crew scheduling |
CN101630440A (zh) * | 2009-06-01 | 2010-01-20 | 北京交通大学 | 常规公共交通衔接城市轨道交通的运营协调优化方法与*** |
CN101572011A (zh) * | 2009-06-10 | 2009-11-04 | 上海理工大学 | 城市公交智能化调度管理***及方法 |
CN101799981A (zh) * | 2010-02-09 | 2010-08-11 | 华南理工大学 | 多模式公共交通区域调度控制方法 |
CN101807222A (zh) * | 2010-02-09 | 2010-08-18 | 武汉大学 | 基于站点的城市公交线网优化配置方法 |
CN102044149A (zh) * | 2011-01-12 | 2011-05-04 | 北京交通大学 | 一种基于时变客流的城市公交运营协调方法与装置 |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104464293A (zh) * | 2014-12-15 | 2015-03-25 | 大连理工大学 | 一种公交专用道上公交停靠站泊位数确定方法 |
CN104504638A (zh) * | 2014-12-18 | 2015-04-08 | 南京信息工程大学 | 一种公交车发车时刻表制定方法 |
CN105185105A (zh) * | 2015-06-26 | 2015-12-23 | 东南大学 | 基于车辆gps和公交ic卡数据的公交换乘识别方法 |
CN105448082A (zh) * | 2015-12-30 | 2016-03-30 | 清华大学 | 一种可变间隔发车的快速公交车组合调度方法 |
CN106056909B (zh) * | 2016-06-29 | 2018-09-07 | 东南大学 | 一种中小城市快速公交阈值确定方法 |
CN106056909A (zh) * | 2016-06-29 | 2016-10-26 | 东南大学 | 一种中小城市快速公交阈值确定方法 |
CN106297273A (zh) * | 2016-09-29 | 2017-01-04 | 百度在线网络技术(北京)有限公司 | 班车路线的处理方法及装置 |
CN106781432A (zh) * | 2016-11-16 | 2017-05-31 | 青岛海信网络科技股份有限公司 | 一种公交车辆调度的方法及装置 |
CN107092976A (zh) * | 2017-03-28 | 2017-08-25 | 东南大学 | 一种多目标模型协同优化多条公交线路发车间隔的方法 |
CN107092976B (zh) * | 2017-03-28 | 2020-06-02 | 东南大学 | 一种多目标模型协同优化多条公交线路发车间隔的方法 |
CN107516147A (zh) * | 2017-08-10 | 2017-12-26 | 中南大学 | 一种高速铁路本线列车开行方案优化方法及其*** |
CN107766969A (zh) * | 2017-09-28 | 2018-03-06 | 东南大学 | 一种基于地铁服务能力瓶颈区段识别的大站快线布设方法 |
CN107766969B (zh) * | 2017-09-28 | 2021-02-19 | 东南大学 | 一种基于地铁服务能力瓶颈区段识别的大站快线布设方法 |
CN108009685A (zh) * | 2017-12-07 | 2018-05-08 | 厦门理工学院 | 一种用于城市公共交通的智能移动终端 |
CN108648441B (zh) * | 2018-06-05 | 2021-08-10 | 武汉理工大学 | 基于站点群的动态公交运行方法 |
CN108648441A (zh) * | 2018-06-05 | 2018-10-12 | 武汉理工大学 | 基于站点群的动态公交运行方法 |
CN109544969A (zh) * | 2019-01-09 | 2019-03-29 | 张鸿青 | 公交车智能跳站运行方法及装置 |
CN109544969B (zh) * | 2019-01-09 | 2020-12-15 | 张鸿青 | 公交车智能跳站运行方法及装置 |
CN110046845A (zh) * | 2019-03-06 | 2019-07-23 | 青岛日日顺物流有限公司 | 一种场区车辆装载方法和装载*** |
CN110046845B (zh) * | 2019-03-06 | 2022-05-20 | 日日顺供应链科技股份有限公司 | 一种场区车辆装载方法和装载*** |
CN111667087A (zh) * | 2019-03-08 | 2020-09-15 | 南京农业大学 | 一种考虑污染排放的公交跳站运行方法 |
CN110796877A (zh) * | 2019-05-14 | 2020-02-14 | 广州学塾加软件科技有限公司 | 一种面向单向公交线路的交通信号控制与公交调度协同控制方法 |
CN110288829A (zh) * | 2019-06-19 | 2019-09-27 | 河海大学 | 一种港湾式公交停靠站停车泊位数量设置判别方法 |
CN110288829B (zh) * | 2019-06-19 | 2020-07-17 | 河海大学 | 一种港湾式公交停靠站停车泊位数量设置判别方法 |
CN110533227A (zh) * | 2019-08-08 | 2019-12-03 | 东南大学 | 一种确定可变线路式公交固定站点和动态站点的方法 |
WO2021083232A1 (zh) * | 2019-10-31 | 2021-05-06 | 同济大学 | 一种网联公交车到站停靠管理优化方法 |
CN110782694A (zh) * | 2019-10-31 | 2020-02-11 | 同济大学 | 一种网联公交车到站停靠管理优化方法 |
CN110969851A (zh) * | 2019-12-09 | 2020-04-07 | 浙江顿悟科技有限公司 | 一种智慧城市公交*** |
CN111063209A (zh) * | 2019-12-23 | 2020-04-24 | 北京航空航天大学 | 配合区间车的自动驾驶公交组合调度优化方法 |
CN110941691A (zh) * | 2020-02-19 | 2020-03-31 | 武汉元光科技有限公司 | 一种顺风车搭乘撮合方法及装置 |
CN112101751A (zh) * | 2020-08-31 | 2020-12-18 | 华南理工大学 | 一种公共交通优化方法、***、装置及存储介质 |
CN112509357A (zh) * | 2020-11-19 | 2021-03-16 | 北京清研宏达信息科技有限公司 | 一种基于遗传算法的公交发车时间间隔优化方法 |
CN113781787A (zh) * | 2021-11-15 | 2021-12-10 | 深圳市都市交通规划设计研究院有限公司 | 公交发车时刻表生成方法及*** |
CN114758522A (zh) * | 2022-04-25 | 2022-07-15 | 驭势科技(北京)有限公司 | 一种车辆的月台停靠方法、装置、设备及介质 |
CN114897445A (zh) * | 2022-07-12 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | 公交线网的停靠点调整优化方法、装置及可读存储介质 |
CN114897445B (zh) * | 2022-07-12 | 2022-11-25 | 深圳市城市交通规划设计研究中心股份有限公司 | 公交线网的停靠点调整优化方法、装置及可读存储介质 |
CN115424437A (zh) * | 2022-08-22 | 2022-12-02 | 厦门筑成信创城市规划设计有限公司 | 站点响应式公交调度方法及装置 |
CN115457793A (zh) * | 2022-08-31 | 2022-12-09 | 东风悦享科技有限公司 | 一种应用于公交调度中始发站非均匀发车的方法及*** |
CN115457793B (zh) * | 2022-08-31 | 2023-06-23 | 东风悦享科技有限公司 | 一种应用于公交调度中始发站非均匀发车的方法及*** |
CN115860298A (zh) * | 2022-12-05 | 2023-03-28 | 四川警察学院 | 公交中间站停靠线路优化方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102831767B (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102831767B (zh) | 一种综合成本优化的城市公共交通多模式站点停靠方法 | |
CN103678917B (zh) | 一种基于模拟退火算法的公交实时到站时间预测方法 | |
CN109543934A (zh) | 城市公交线网的综合指标的评价方法 | |
CN103440422B (zh) | 一种基于带时间窗的到站时间预测的公交晚点恢复方法 | |
CN111754039B (zh) | 纯电动公交线网综合集成优化设计的方法 | |
CN105702035B (zh) | 一种利用历史公交数据评估乘车难易程度的方法 | |
CN112819316B (zh) | 一种综合客运枢纽轨道交通***的枢纽运能识别方法 | |
CN103745089A (zh) | 一种多维公共交通运行指数评价方法 | |
CN102945616B (zh) | 一种公交站点信息的采集方法 | |
CN110239600B (zh) | 面向再生制动能量利用的列车运行控制方法及*** | |
CN103226801A (zh) | 基于多用户配流模型的空港集疏运交通量确定方法 | |
CN112580866B (zh) | 基于全程车与区间车组合调度的公交线路串车优化方法 | |
CN110084397A (zh) | 一种地铁直达线路规划方法 | |
CN110599760A (zh) | 一种多模式交通网络下的出行行为模拟方法 | |
CN111859718B (zh) | 一种区域多制式轨道交通车站拥挤系数计算方法及*** | |
CN112446648B (zh) | 基于非高峰时段的城市轨道交通混合运输方法和*** | |
CN112309119A (zh) | 一种城市交通***容量分析优化方法 | |
CN106274994A (zh) | 一种列车全线到站时间预测方法和*** | |
CN111881939A (zh) | 一种基于聚类算法的共享单车停车区布设方法 | |
CN117401001A (zh) | 复杂工况下的城轨多列车行车调度综合节能控制方法及装置 | |
CN116542404B (zh) | 一种客运枢纽场站抵站旅客接续换乘时间的预测方法 | |
CN103226666B (zh) | 复杂功能互通立体交叉***功效评价方法 | |
CN116307580A (zh) | 运力调度方法、装置、电子设备及存储介质 | |
CN111859717B (zh) | 最小化区域多制式轨道交通旅客拥挤系数的方法及*** | |
CN115049272A (zh) | 一种基于电池交换进行中间站充电的电动公交调度方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140709 Termination date: 20170813 |