CN115409388A - 多车型定制公交运营优化方法 - Google Patents
多车型定制公交运营优化方法 Download PDFInfo
- Publication number
- CN115409388A CN115409388A CN202211069012.0A CN202211069012A CN115409388A CN 115409388 A CN115409388 A CN 115409388A CN 202211069012 A CN202211069012 A CN 202211069012A CN 115409388 A CN115409388 A CN 115409388A
- Authority
- CN
- China
- Prior art keywords
- vehicle
- belongs
- cost
- constraints
- path
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005457 optimization Methods 0.000 title claims abstract description 32
- 238000004422 calculation algorithm Methods 0.000 claims description 39
- 230000009977 dual effect Effects 0.000 claims description 13
- 230000002068 genetic effect Effects 0.000 claims description 12
- 238000007689 inspection Methods 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 108090000623 proteins and genes Proteins 0.000 description 4
- 230000004044 response Effects 0.000 description 3
- 230000035772 mutation Effects 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 101150067055 minC gene Proteins 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/02—Reservations, e.g. for tickets, services or events
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明涉及城市公共交通领域,公开了一种多车型定制公交运营优化方法,包括:获取定制公交预约信息;根据预约信息,确定多车型定制公交车队配置和路径优化的目标函数与约束条件,其中,目标函数是最小化运营成本和超时惩罚成本,约束条件包括:所有的乘客都需要被服务、车辆服务里程和服务时间限制、每种类型车辆的数量限制和乘客时间窗限制;根据目标函数和约束条件,建立混合整数非线性规划模型;对混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案。从而解决了现有技术无法解决的多车型定制公交系统运营优化问题,而且有利于企业降低运营成本。
Description
技术领域
本发明涉及城市公共交通领域,具体涉及一种多车型定制公交运营优化方法。
背景技术
随着互联网、云计算和人工智能等新技术的发展,居民日常出行方式正在发生深刻变化。定制公交车作为一种新兴的交通出行方式,提供需求响应、高效率和高质量的大众服务。自20世纪70年代需求响应型公交服务问世以来[1],定制公交已在全球多个城市开通运营,如里斯本和旧金山[2]。在我国,定制公交服务已经推广到北京、河北、山西等29个省份,运营线路5400余条,年客运量接近1.8亿人次[3]。我们需要充分利用大数据、智能技术,提供需求响应式服务,“人少坐小车,人多坐中巴”,用多样化、定制化的公交服务,进一步满足市民出行的多样化需求。但是,现有技术中关于定制公交的运营方案是针对单一车型的定制公交的车队管理和路径规划,因此,对于多种车型的定制公交运营不适用。
参考文献
[1]Rahimi,M.,Amirgholy,M.,Gonzales,E.J..System modeling of demandresponsive transportation services:evaluating cost efficiency of service andcoordinated taxi usage.Transp.Res.Part E Logist.Transp.Rev.2018,112:66-83.
[2]Martínez,L.M.,Viegas,J.M.,Eiro,T..Formulating a new expressminibus service design problem as aclusteringproblem.Transport.Sci.2015,49(1):85-98..
[3]孙宏阳.大城市公交如何吸引客流?定制公交已在京冀等29个省份开通[N].北京日报,2021-11-12.
发明内容
为了解决上述技术问题,本发明提供一种多车型定制公交运营优化方法,适用于多种车型的定制公交的运营。
本发明公开了一种多车型定制公交运营优化方法,包括:
获取定制公交预约信息;所述预约信息包括上车站点、下车站点、期望上车时间和期望下车时间;
根据所述预约信息,确定多车型定制公交车队配置和路径优化的目标函数与约束条件,其中,所述目标函数是最小化运营成本和超时惩罚成本,所述约束条件包括:所有的乘客都需要被服务、车辆服务里程和服务时间限制、每种类型车辆的数量限制和乘客时间窗限制;
根据所述目标函数和约束条件,建立混合整数非线性规划模型;
对所述混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案。
优选的,所述对所述混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案,具体包括:
将所述混合整数非线性规划模型中的非线性目标进行线性化处理,得到混合整数线性规划模型,所述非线性目标包括运营成本和超时惩罚成本;
使用Dantzig-Wolfe(D-W)分解所述混合整数线性规划模型进行分解重构,引入新的变量,得到限制主问题模型和子问题模型;
采用Cplex求解器求解小规模子问题,采用启发式算法求解大规模子问题,优化异构车队配置、定制公交路线和时刻表以及订单指派方案。
优选的,所述混合整数非线性规划模型包括:目标函数及所述目标函数的约束条件,具体包括:
所述目标函数如下述公式:
minCS=CF+CO+CP, (1)
其中,固定运营成本CF取决于使用了多少辆车,表述为:
可变运营成本CO取决于车辆类型及线路长度,表述为:
惩罚成本CP定义为超出乘客r可接受的最晚到达时间的惩罚成本,表述为:
所述约束条件如下述公式(5)~(27)所示:
xkii=0,k∈K,i∈N (8)
ATkj≤ATki+HTki+tij+Z(1-xkij),k∈K,i,j∈V (9)
ATkj≥ATki+HTki+tij-Z(1-xkij),k∈K,i,j∈V (10)
ATkd(r)≥ATko(r)-Z(1-αrk),k∈K,r∈R (11)
ATko(r)≤(EDTr+tmax)αrk+Z(1-αrk),r∈R,k∈K (16)
(EDTr-tmax)αrk≤ATko(r)+HTko(r),r∈R,k∈K (17)
zka-Z(1-xksi)≤zki,i∈N,k∈K (20)
zka+Z(1-xksi)≥zki,i∈N,k∈K (21)
zks=0,k∈K (24)
xkij,ymk,αrk∈{0,1};ATki,HTki,zki≥0,k∈K,i,j∈V,r∈R,m∈M (27)
其中,公式(1)~公式(27)所用数学符号的含义如下:
N:候选站点集合,i∈N;
s,t:虚拟起点和终点的索引;
V:所有节点的集合,包括候选站和虚拟节点,即:N∪{s,t};
M:定制公交的类型,m∈M;
K:车辆集合,k∈K;
R:订单集合,r∈R;
capm:m型车辆的容量,即座位数,m∈M;
dij:站点i和站点j之间的距离,i,j∈V;
tij:站点i和站点j之间的行驶时间,i,j∈V;
o(r):订单r的上车站点,r∈R;
d(r):订单r的下车站点,r∈R;
EDTr:订单r的期望上车时间,r∈R
EATr:订单r的期望下车时间,r∈R;
tmax:时间偏差阈值;
LATr:订单r最晚可接受到达时间,即LATr=EATr+tmax,r∈R;
Fm:m型车辆的最大可用数量,m∈M;
cf:一条定制公交线路的固定运营成本;
cm:m型车辆的单位距离运营成本,m∈M;
cp:在LATr之前未到达的单位超时时间的惩罚成本,r∈R;
CF:总固定运营成本;
CO:总可变运营成本;
CP:总超时惩罚成本;
CS:系统总成本;
Z:一个极大的正整数;
xkij:决策变量,若车辆k经过弧(i,j)则取1,否则取0,k∈K,i,j∈V;
ymk:决策变量,若车辆k的类型为m则取1,否则取0,k∈K,m∈M;
αrk:决策变量,若订单r分配给车辆k则取1,否则取0,r∈R,k∈K;
ATki;决策变量,车k到达站点i的时间,i∈V,k∈K;
HTki:决策变量,车k在站点i的等待时间,i∈V,k∈K;
zki:决策变量,车k到达站点i时所装载的乘客人数(即订单数),i∈V,k∈K;
其中,约束(5)-(8)是线路约束,用于确保每辆车可以匹配可行的路线;
约束(9)-(12)是车辆的到达时间限制,确保所有车辆能够有序、合理地完成运输任务;
约束(13)-(17)是订单指派约束,确保每个订单被指派到可行的车辆线路与车辆;约束(18)-(24)是车辆的容量限制约束,确保装载乘客的数量不会超过车辆的容量。约束(25)是车辆的最大行驶时间约束;
约束(26)是车队规模限制;约束条件(27)定义了决策变量的取值范围。
优选的,所述将所述混合整数非线性规划模型中的非线性目标进行线性化处理,得到混合整数线性规划模型,具体包括:
分别对所述运营成本和对所述超时惩罚成本进行线性化处理。
优选的,所述对所述运营成本进行线性化处理,具体包括:
所述对所述超时惩罚成本进行线性化处理,具体包括:
定义一个辅助变量actTrk∈R+∪{0},表示订单r在车辆k上的实际下车时间;
actTrk满足以下约束:
actTrk≥ATkd(r)-Z(1-αrk),k∈K,r∈R, (32)
actTrk≤ATkd(r)+Z(1-αrk),k∈K,r∈R, (33)
actTrk≤Zαrk,k∈K,r∈R, (34)
将所述超时惩罚成本CP用actTrk表示为:
优选的,所述使用Dantzig-Wolfe(D-W)分解所述混合整数线性规划模型进行分解重构,引入新的变量,得到限制主问题模型和子问题模型,具体包括:
引入新的集合Ω为定制公交的所有可行路径的集合,ω∈Ω;其中,Ωm表示使用m型车辆的路径集合;新的变量λωk为0-1变量,如果路径ω分配给车辆k则取1,否则取0;
将所述混合整数线性规划模型等价转化为如下的限制主问题模型:
主问题的目标函数为:
主问题约束条件为:
λωk∈{0,1},ω∈Ω,k∈K (40)
其中,目标(36)是最小化总成本,包括固定运营成本和可变运营成本以及超时惩罚成本;约束(37)表示订单必须全部被服务;约束(38)确保每辆车至多分配一条路线;约束(39)限制了m型车辆的可用数量;约束(40)为变量的取值范围;
将所述混合整数线性规划模型分解后得到子问题模型:
子问题目标函数如公式(41)所示:
xii=0,i∈N (45)
ATj≤ATi+HTi+tij+Z(1-xij),i,j∈V (46)
ATj≥ATi+HTi+tij-Z(1-xij),i,j∈V (47)
ATd(r)≥ATo(r)-Z(1-αr),r∈R (48)
ATo(r)≤(EDTr+tmax)αr+Z(1-αr),r∈R (52)
(EDTr-tmax)αr≤ATo(r)+HTo(r),r∈R (53)
za-Z(1-xsi)≤zi,i∈N (56)
za+Z(1-xsi)≥zi,i∈N (57)
zks=0,k∈K (60)
xij,ym,αr∈{0,1},ATi,HTi,zi≥0,i,j∈V,r∈R,m∈M (62)
约束(42)-(63)的含义与所述约束(5)-(12)和(14)-(25)的含义相同;根据上述子问题的最优解xij得到的一条车辆的路径ω,添加到路径集合Ω中;或,根据最优解xij得到可变运营成本根据ATi和αr的值,得到路径ω所服务乘客的超时惩罚成本将和作为目标系数添加到主问题的式(36)中;αr和ym分别作为约束(37)和(39)的系数。
优选的,所述启发式算法包括遗传算法,蚁群算法、禁忌搜索。
优选的,所述采用所述遗传算法求解大规模子问题,具体包括:
生成定制公交初始的可行路径的集合Ω,得到每条路径所使用的车型以及可变运营成本和超时惩罚成本;
调用求解器lpsolve对所述主问题进行求解,并计算所述主问题决策变量的对偶变量,使用对偶变量对子问题目标函数进行更新;其中,所述子问题的目标函数对应的值称之为主问题的检验数;
使用遗传算法对所述子问题进行求解,得到检验数和一条可行的路径以及它的时刻表,该路径为列生成方法中一个潜在的列;同时,得到这条路径所服务的订单以及所使用的车型;
基于子问题的解,判断检验数是否符合列生成的检验标准,若符合标准则向所述主问题中添加列,即将获得的路径作为一个元素添加到集合Ω中,并调用求解器lpsolve对所述主问题进行求解,并计算所述主问题决策变量的对偶变量,使用对偶变量对子问题目标函数进行更新;若不符合标准,则表明已经不存在有价值的列,则对当前所述主问题进行求解,得到的解即为原问题的最优解,输出定制公交的路线以及时刻表,还有每条路线所使用的车型和其服务的订单。
本发明与现有技术相比,具有以下优点:
本发明公开了一种多车型定制公交运营优化方法,通过将多类型定制公交的情况引入到模型中,允许车辆在站点有等待时间,将车队组成、订单-路线分配、路径规划、时刻表优化综合考虑,以最小化系统总成本为目标,综合考虑线路运营成本以及服务订单的超时惩罚成本,将当前的异构车队的使用和路径优化问题视为传统路径规划问题的进一步拓展,能够为定制公交系统提供决策支持与系统规划,适用于目前定制公交系统运营企业中多类型定制公交的实际情况。
附图说明
图1为本发明实施例提供的编码方式;
图2为本发明实施例提供的交叉规则;
图3为本发明实施例提供的变异规则;
图4为本发明实施例提供的路线及时刻表和等待时间结果;
图5为本发明实施例提供的多车型定制公交运营优化方法流程示意图。
具体实施方式
本发明提供了一种基于列生成和启发式算法相结合的多车型定制公交运营优化方法,采用精确算法与启发式算法相结合的方式解决多车型定制公交系统中的运营优化问题,降低定制公交运营企业的运营成本,具有决策可靠性高且计算稳定性强的优点。
为了使本发明的目的、技术方案及优点更加清楚,以下通过具体实施,对本发明进一步详细说明。
本发明实施例提供的一种多车型定制公交运营优化方法,包括:
获取定制公交预约信息;所述预约信息包括上车站点、下车站点、期望上车时间和期望下车时间;根据所述预约信息,确定多车型定制公交车队配置和路径优化的目标函数与约束条件,其中,所述目标函数是最小化运营成本和超时惩罚成本,所述约束条件包括:所有的乘客都需要被服务、车辆服务里程和服务时间限制、每种类型车辆的数量限制和乘客时间窗限制;根据所述目标函数和约束条件,建立混合整数非线性规划模型;对所述混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案。
所述对所述混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案,具体包括:
将所述混合整数非线性规划模型中的非线性目标进行线性化处理,得到混合整数线性规划模型,所述非线性目标包括运营成本和超时惩罚成本;
使用Dantzig-Wolfe(D-W)分解所述混合整数线性规划模型进行分解重构,引入新的变量,得到限制主问题模型和子问题模型;
采用Cplex求解器求解小规模子问题,采用启发式算法求解大规模子问题,优化异构车队配置、定制公交路线和时刻表以及订单指派方案。具体的,如图5所示,该方法可以包括下述步骤:
步骤S1:收集定制公交用户终端上传的定制公交预约信息,包括上下车站点以及期望上下车时间,并进行数据预处理;
步骤S2:根据收集到的预约信息,确定多车型定制公交车队配置和路径优化的目标函数与约束条件,其中,目标函数是最小化运营成本和超时惩罚成本之和,约束条件主要包括:所有的乘客都需要被服务、车辆服务里程和服务时间限制、每种类型车辆的数量限制、乘客时间窗限制等;
步骤S3:根据所述目标函数和约束条件,建立混合整数非线性规划模型;
步骤S4:将混合整数非线性规划模型中的非线性目标进行线性化处理,包括运营成本和超时惩罚成本两部分,得到混合整数线性规划模型;
步骤S5:使用Dantzig-Wolfe(D-W)分解对S4中得到的混合整数线性规划模型进行分解重构,引入新的变量,得到限制主问题模型和子问题模型;
步骤S6:分别采用Cplex求解器和启发式算法求解小规模子问题与大规模子问题,优化异构车队配置、定制公交路线和时刻表以及订单指派方案。
在一个实施例中,当定制公交运营系统为用户提供7个备选站点进行选择时,上述步骤S1:收集定制公交用户终端上传的定制公交预约信息,包括上下车站点以及期望上下车时间,设定时间为60个单位,其中前20个订单的信息如表1所示,具体包括:
表1用户上传的预定信息(前20)
对数据进行预处理,即加入虚拟起始点后,上下车站点序号发生变化,所以前20个订单的信息处理如表2所示,具体包括:
表2用户上传的预定信息处理后
在一个实施例中,上述步骤S2:根据收集到的预约信息,确定多车型定制公交车队配置和路径优化的目标函数与约束条件,其中,目标函数是最小化运营成本和超时惩罚成本之和,约束条件主要包括:所有的乘客都需要被服务、车辆服务里程和服务时间限制、每种类型车辆的数量限制、乘客时间窗限制等;
在一个实施例中,上述步骤S3:根据所述目标和约束条件,建立混合整数非线性规划模型,最小化定制公交运营系统的运营成本和超时惩罚成本;其中,所述混合整数非线性规划模型包括:目标函数及其约束条件,具体包括:
目标函数如下述公式(1):
min CS=CF+CO+CP, (1)
其中,公式(1)使多车型的定制公交的路线固定运营成本CF、可变运营成本CO和未在规定时间内送达的惩罚成本CP之和CS最小化。其中固定运营成本CF取决于使用了多少辆车,可表述为:
可变运营成本CO取决于车辆类型及线路长度,可表述为:
惩罚成本CP定义为超出乘客r可接受的最晚到达时间的惩罚成本,可表述为:
所述约束条件如下述公式(5)~(27)所示:
xkii=0,k∈K,i∈N (8)
ATkj≤ATki+HTki+tij+Z(1-xkij),k∈K,i,j∈V (9)
ATkj≥ATki+HTki+tij-Z(1-xkij),k∈K,i,j∈V (10)
ATkd(r)≥ATko(r)-Z(1-αrk),k∈K,r∈R (11)
ATko(r)≤(EDTr+tmax)αrk+Z(1-αrk),r∈R,k∈K (16)
(EDTr-tmax)αrk≤ATko(r)+HTko(r),r∈R,k∈K (17)
zka-Z(1-xksi)≤zki,i∈N,k∈K (20)
zka+Z(1-xksi)≥zki,i∈N,k∈K (21)
zks=0,k∈K (24)
xkij,ymk,αrk∈{0,1};ATki,HTki,zki≥0,k∈K,i,j∈V,r∈R,m∈M (27)
其中,公式(1)~公式(27)所用数学符号的含义如下:
N:候选站点集合,i∈N;
s,t:虚拟起点和终点的索引;
V:所有节点的集合,包括候选站和虚拟节点,即:N∪{s,t};
M:定制公交的类型,m∈M;
K:车辆集合,k∈K;
R:订单集合,r∈R;
capm:m型车辆的容量,即座位数,m∈M;
dij:站点i和站点j之间的距离,i,j∈V;
tij:站点i和站点j之间的行驶时间,i,j∈V;
o(r):订单r的上车站点,r∈R;
d(r):订单r的下车站点,r∈R;
EDTr:订单r的期望上车时间,r∈R
EATr:订单r的期望下车时间,r∈R;
tmax:时间偏差阈值;
LATr:订单r最晚可接受到达时间,即LATr=EATr+tmax,r∈R;
Fm:m型车辆的最大可用数量,m∈M;
cf:一条定制公交线路的固定运营成本;
cm:m型车辆的单位距离运营成本,m∈M;
cp:在LATr之前未到达的单位超时时间的惩罚成本,r∈R;
CF:总固定运营成本;
CO:总可变运营成本;
CP:总超时惩罚成本;
CS:系统总成本;
Z:一个极大的正整数;
xkij:决策变量,若车辆k经过弧(i,j)则取1,否则取0,k∈K,i,j∈V;
ymk:决策变量,若车辆k的类型为m则取1,否则取0,k∈K,m∈M;
αrk:决策变量,若订单r分配给车辆k则取1,否则取0,r∈R,k∈K;
ATki;决策变量,车k到达站点i的时间,i∈V,k∈K;
HTki:决策变量,车k在站点i的等待时间,i∈V,k∈K;
zki:决策变量,车k到达站点i时所装载的乘客人数(即订单数),i∈V,k∈K。
其中,约束(5)-(8)是线路约束,确保每辆车可以匹配可行的路线。具体来说,约束(5)确保每辆车必须从虚拟起点s出发并到达虚拟终点t;约束条件(6)是流平衡约束,约束(7)保证每辆车一次最多可以访问一个站点;约束(8)显示了当i=j时xkij的取值。
约束(9)-(12)是车辆的到达时间限制,确保所有车辆能够有序、合理地完成运输任务。其中,约束条件(9)和(10)确定每辆车在每个车站的到达时间,约束(11)确保车辆k在订单r目的地的到达时间始终不小于车辆k在其始发站o(r)的到达时间;约束(12)限制了车辆的最大等待时间,即每辆车在每个站点的等待时间不能超过指定时间。
约束(13)-(17)是订单指派约束,确保每个订单被指派到可行的车辆线路与车辆。其中,约束(13)确保每个订单都被服务;约束(14)和(15)保证只有当车辆k经过订单r的上下车站时,才有可能将订单r分配给车辆k;约束(16)和(17)将车辆和订单的服务时间联系起来,即车辆k到达订单r的上车站点o(r)的时间在该订单可接受的时间窗口内,即[EDTr-tmax,EDTr+tmax],订单r才可能分配给车辆k。
约束(18)-(24)是车辆的容量限制约束,确保装载乘客的数量不会超过车辆的容量。其中,约束条件(18)-(19)表明车辆k在i站和j站所装载的订单数量的关系;同时,约束(20)-(21)保证在虚拟节点s处分配给车辆k的订单数等于在其第一个车站i处分配给该车的订单数;约束(22)表示m型车辆k容量约束;约束(23)表明每辆车只能使用一种车型;约束(24)将虚拟节点处每辆车上的订单数量限制为0。
约束(26)是车队规模限制。由于资源有限,给出了每类型车辆数量的上限。
约束条件(27)定义了决策变量的取值范围。
在一个实施例中,上述步骤S4:将混合整数非线性规划模型中的非线性部分线性化,包括运营成本和超时惩罚成本两部分的线性化过程,具体包括:
同样的,我们定义一个辅助变量actTrk∈R+∪{0},表示乘客r在车辆k上的实际下车时间。所以actTrk可以表述为以下约束:
actTrk≥ATkd(r)-Z(1-αrk),k∈K,r∈R, (32)
actTrk≤ATkd(r)+Z(1-αrk),k∈K,r∈R, (33)
actTrk≤Zαrk,k∈K,r∈R, (34)
由上可知,当αrk=1时,actTrk=ATkd(r);当αrk=0时,actTrk=0。同样地,将超时惩罚成本CP用actTrk表示为:
因此,在不改变任何解空间的情况下,将非线性目标函数(3)和(4)转换为等效线性形式。然后,可以使用列生成算法直接求解重新构造的模型,以获得原始非线性模型的最优解。
在一个实施例中,上述步骤S5:使用D-W分解对所述混合整数线性规划模型进行分解重构,引入新的变量得到重构的规划模型,包括限制主问题模型和子问题模型,具体包括:
步骤S51:引入新的集合Ω为定制公交的所有可行路径的集合,ω∈Ω;其中Ωm表示使用m型定制公交的路径集合。新的变量λωk为0-1变量,如果路径ω分配给车辆k则取1,否则取0;
步骤S52:将所述混合整数线性规划模型等价转化为如下的限制主问题模型:
主问题的目标函数为:
主问题约束条件为:
λωk∈{0,1},ω∈Ω,k∈K (40)
目标(36)是最小化总成本,包括固定运营成本和可变运营成本以及超时惩罚成本;约束(37)表示订单必须全部被服务;约束(38)确保每辆车至多分配一条路线;约束(39)限制了m型车辆的可用数量;约束(40)为变量的取值范围。
步骤S53:将所述混合整数线性规划模型分解后得到子问题模型:
子问题目标函数如公式(41)所示:
xii=0,i∈N (45)
ATj≤ATi+HTi+tij+Z(1-xij),i,j∈V (46)
ATj≥ATi+HTi+tij-Z(1-xij),i,j∈V (47)
ATd(r)≥ATo(r)-Z(1-αr),r∈R (48)
ATo(r)≤(EDTr+tmax)αr+Z(1-αr),r∈R (52)
(EDTr-tmax)αr≤ATo(r)+HTo(r),r∈R (53)
za-Z(1-xsi)≤zi,i∈N (56)
za+Z(1-xsi)≥zi,i∈N (57)
zks=0,k∈K (60)
xij,ym,αr∈{0,1},ATi,HTi,zi≥0,i,j∈V,r∈R,m∈M (62)
约束(42)-(63)的含义与S3中约束(5)-(12)和(14)-(25)的含义相同。根据上述子问题的最优解xij可以得到的一条车辆的路径ω,并添加到路径集合Ω中;还可以根据最优解xij得到可变运营成本同时根据ATi和αr的最优值,得到路径ω所服务乘客的超时惩罚成本最后将和作为目标系数添加到主问题的式(36)中;αr和ym分别作为约束(37)和(39)的系数。
在一个实施例中,上述步骤S6:分别采用Cplex求解器和启发式算法对所述重构的规划模型进行求解,得到异构车队配置、定制公交路线和时刻表以及订单指派方案,其中,启发式算法包括遗传算法,蚁群算法、禁忌搜索等。其中小规模算例的子问题采用Cplex商业求解器直接求解,大规模算例的子问题用启发式算法求解,以遗传算法为例,具体包括:
具体步骤如下:
步骤S61:生成定制公交初始的可行路径的集合Ω,并得到每条路径所使用的车型以及可变运营成本和超时惩罚成本;
本步骤中,采用每条路径仅服务一名乘客,以每个乘客的上下车站点作为每条路线的站点,乘客的期望上车时间作为车辆在起点的出发时间,所产生的路径添加到集合Ω中,得到限制主问题。
步骤S62:调用求解器lpsolve对所述主问题进行求解,并计算其决策变量的对偶变量,使用对偶变量对子问题目标函数进行更新;其中,所述子问题的目标函数对应的值称之为主问题的检验数;
步骤S63:使用启发式算法对所述子问题进行求解,得到目标值(即检验数)和一条可行的路径以及它的时刻表,该路径为列生成方法中一个潜在的列;同时还可以得到这条路径所服务的乘客以及所使用的车型;
步骤S631:数据初始化,将所有预定订单的数据,备选站点和系统的数据,定制公交的相关参数输入到算法中并设置遗传算法的参数,其中:
所有预定订单的数据:上车站点、期望上车时间、下车站点、期望下车时间,如表1所示;
备选站点和系统的相关数据:车辆可到达的7个备选站点和两个虚拟站点(其中序号1和9为虚拟站点)以及站点之间的行驶时间(如表3所示)、时间偏差阈值设为2个时间单位、在每个站点的最大停留时间设为3个时间单位、超时单位时间的惩罚成本设为10、每辆车的固定成本设为60
表3备选站点之间的行驶时间
车辆的相关参数:车辆的类型、每种类型的可用车辆数量以及座位数、每种类型车辆的单位距离运营成本、最长行驶时间,如表4所示;
表4定制公交的相关信息
遗传算法设计的参数:种群数量sizepop为10、最大迭代次数maxG为20、交叉概率Pc为0.9、变异概率Pm为0.2。
步骤S632:订单路线分配,将订单尽可能的分配到可行的路径上,则订单-路线分配算法如下:
a)已知初始路径ω、路径上的站点数l、到达ω上各站点i的时间ATi、在ω上各站点的停留时间HTi(初始全为0)、在ω上各站点的装载乘客数量initZi(初始全为0)、ω所需的行驶时间Time;
c)若在(b)中的范围内,则判断订单r的下车站是否在路径ω上,若在,则直接将订单r分配给路径ω,更新initZi,并根据以下公式(63)更新站点i的车辆等待时间:
HTi=max{HTi,max(EDTr-(ATi+tmax),0)} (63)
其中,max(EDTr-(ATi+tmax),0)是为了判断订单r的期望上车时间是否超过了ATi+tmax,若超过,则存在等待时间,则需要更新HTi,若未超过,则等待时间为0;若订单r的下车站点d(r)不在路径ω上,则根据公式(64):
Dn=djd(r)+dd(r)k-djk,j,k∈[i+1,l],n∈[i,l-1] (64)
选择其中min(Dn),并判断min(Dn)+Time是否超过最大行驶时间若超过,则该订单不能分配给路径ω,若未超过,则该订单可以分配给路径ω,并将下车站点d(r)插入到增加距离最少即min(Dn)的路径ω点n的后面,并更新AT、HT、initZ和Time;
d)重复(b)-(c),当路径ω上所有的点都被遍历,或达到车辆的最大容量或达到车辆最大行驶时间时,结束迭代,转到(e);
e)找到下车点在路径ω且上车点不在的订单集合DR,对DR中的订单进行循环,同样根据(b)-(c)的原理,在容量和行驶时间允许的情况下,考虑是否分配订单并插入订单起点o(r);
f)更新路径ω的可变运营成本、所服务订单的超时惩罚成本以及所需要采用的车型;
步骤S633:遗传算法个体编码:采用实数编码的方式,如下图1所示,每个个体由三行数据构成,其中第一行表示路径ω所经过的站点,基因序列1表示虚拟起点,最后一个基因序列表示虚拟终点,第二行表示到达各个站点的时间AT,第三行表示在各个站点的停留时间HT;
步骤S634:生成初始种群,每个个体生成的规则:随机选择一个订单r,根据Dijkstra算法得到订单r上下车站点之间的最短路作为初始路径,再根据S632,分配订单并更新路径、AT、HT和车辆类型,得到和
步骤S635:适应度函数,设置为子问题的目标函数;
步骤S636:选择,采用轮盘赌的方法选择两个个体;
步骤S637:交叉,采用单点交叉的方式,找到两条路径中实际长度较小的值A,生成小于A且大于1的随机正整数值a,将个体中第一行和第二行序列a之前的基因序列进行交换,并修复AT,且将HT的值设为0,对于重复的站点,删除重复的后一个站点,若两条路径的长度分别为6、5,则随机生成小于5大于1的正整数,若为4,则交换方式如图2所示,然后再根据S632,分配订单并更新路径、AT、HT和车辆类型,得到和计算适应度:
步骤S638:变异,同样采用单体变异规则,生成小于路径长度且大于1的正整数a,将序列a的基因随机变为不在路径中的站点序号,如图3所示,并修复AT,且将HT的值设为0,然后再根据S632,分配订单并更新路径、AT、HT和车辆类型,得到和计算适应度:
步骤S639:判断是否满足终止条件,即是否达到最大迭代次数,否则重复S636-S639。
步骤S64:基于子问题的解,判断检验数是否符合列生成的检验标准,即对于本系统是否大于等于0:若符合标准则向所述主问题中添加列,即将获得的路径作为一个元素添加到集合Ω中,并返回步骤S62;若不符合标准,则表明已经不存在有价值的列,则对当前所述主问题进行求解,得到的解即为原问题的最优解,输出定制公交的路线以及时刻表,还有每条路线所使用的车型和其服务的乘客。
按照上述步骤分别对前10个订单和20个对异构车队确定和路径规划模型进行求解,得到多车型定制公交的车队规模以及对订单的响应情况和路径规划方案如表5所示。
表5定制公交运营系统结果示例
注释:CG:列生成;GA:遗传算法
由表5可以看出Cplex求解器可以在短时间内直接求解小规模算例,“列生成+启发式”相结合的算法可以在短时间内求解大规模算例,并得到较优解,对于运营企业是适用的。其中10个订单的路线相关结果如图4所示:
其中;路径1服务序号为1、2、3、4的乘客,使用1类型的定制公交;路径2服务序号为5、6、7、8、9、10的乘客,使用2类型的定制公交。
本发明公开了一种基于列生成和启发式算法相结合的多车型定制公交运营优化方法,基于目前定制公交系统运营企业中多类型定制公交的实际情况,以往针对单一车型的定制公交的车队管理和路径规划的研究不再适用。本发明将多类型定制公交的情况引入到模型中,允许车辆在站点有等待时间,将车队组成、订单-路线分配、路径规划、时刻表优化综合考虑,具有实际使用价值。
本发明公开的基于列生成和启发式算法相结合的多车型定制公交运营优化方法,以最小化系统总成本为目标,综合考虑线路的固定运营成本、可变运营成本以及服务乘客的超时惩罚成本,将当前的异构车队的使用和路径优化问题视为传统路径规划问题的进一步拓展,能够为定制公交系统提供决策支持与系统规划,降低定制公交运营商的系统总成本。
本发明公开的基于列生成和启发式算法相结合的多车型定制公交运营优化方法,在对模型进行求解时将精确算法中的D-W分解技术、列生成技术与启发式算法中的启发式算法进行结合,能够在保证解决方案精确度的基础上提升计算速度。
Claims (8)
1.一种多车型定制公交运营优化方法,其特征在于,包括:
获取定制公交预约信息;所述预约信息包括上车站点、下车站点、期望上车时间和期望下车时间;
根据所述预约信息,确定多车型定制公交车队配置和路径优化的目标函数与约束条件,其中,所述目标函数是最小化运营成本和超时惩罚成本,所述约束条件包括:所有的乘客都需要被服务、车辆服务里程和服务时间限制、每种类型车辆的数量限制和乘客时间窗限制;
根据所述目标函数和约束条件,建立混合整数非线性规划模型;
对所述混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案。
2.根据权利要求1所述的多车型定制公交运营优化方法,所述对所述混合整数非线性规划模型进行处理,得到异构车队配置、定制公交路线和时刻表以及订单指派方案,具体包括:
将所述混合整数非线性规划模型中的非线性目标进行线性化处理,得到混合整数线性规划模型,所述非线性目标包括运营成本和超时惩罚成本;
使用Dantzig-Wolfe(D-W)分解所述混合整数线性规划模型进行分解重构,引入新的变量,得到限制主问题模型和子问题模型;
采用Cplex求解器求解小规模子问题,采用启发式算法求解大规模子问题,优化异构车队配置、定制公交路线和时刻表以及订单指派方案。
3.根据权利要求1所述的多车型定制公交运营优化方法,所述混合整数非线性规划模型包括:目标函数及所述目标函数的约束条件,具体包括:
所述目标函数如下述公式:
min CS=CF+CO+CP, (1)
其中,固定运营成本CF取决于使用了多少辆车,表述为:
可变运营成本CO取决于车辆类型及线路长度,表述为:
惩罚成本CP定义为超出乘客r可接受的最晚到达时间的惩罚成本,表述为:
所述约束条件如下述公式(5)~(27)所示:
xkii=0,k∈K,i∈N (8)
ATkj≤ATki+HTki+tij+Z(1-xkij),k∈K,i,j∈V (9)
ATkj≥ATki+HTki+tij-Z(1-xkij),k∈K,i,j∈V (10)
ATkd(r)≥ATko(r)-Z(1-αrk),k∈K,r∈R (11)
ATko(r)≤(EDTr+tmax)αrk+Z(1-αrk),r∈R,k∈K (16)
(EDTr-tmax)αrk≤ATko(r)+HTko(r),r∈R,k∈K (17)
zka-Z(1-xksi)≤zki,i∈N,k∈K (20)
zka+Z(1-xksi)≥zki,i∈N,k∈K (21)
zks=0,k∈K (24)
xkij,ymk,αrk∈{0,1};ATki,HTki,zki≥0,k∈K,i,j∈V,r∈R,m∈M (27)
其中,公式(1)~公式(27)所用数学符号的含义如下:
N:候选站点集合,i∈N;
s,t:虚拟起点和终点的索引;
V:所有节点的集合,包括候选站和虚拟节点,即:N∪{s,t};
M:定制公交的类型,m∈M;
K:车辆集合,k∈K;
R:订单集合,r∈R;
capm:m型车辆的容量,即座位数,m∈M;
dij:站点i和站点j之间的距离,i,j∈V;
tij:站点i和站点j之间的行驶时间,i,j∈V;
o(r):订单r的上车站点,r∈R;
d(r):订单r的下车站点,r∈R;
EDTr:订单r的期望上车时间,r∈R
EATr:订单r的期望下车时间,r∈R;
tmax:时间偏差阈值;
LATr:订单r最晚可接受到达时间,即LATr=EATr+tmax,r∈R;
Fm:m型车辆的最大可用数量,m∈M;
cf:一条定制公交线路的固定运营成本;
cm:m型车辆的单位距离运营成本,m∈M;
cp:在LATr之前未到达的单位超时时间的惩罚成本,r∈R;
CF:总固定运营成本;
CO:总可变运营成本;
CP:总超时惩罚成本;
CS:系统总成本;
Z:一个极大的正整数;
xkij:决策变量,若车辆k经过弧(i,j)则取1,否则取0,k∈K,i,j∈V;
ymk:决策变量,若车辆k的类型为m则取1,否则取0,k∈K,m∈M;
αrk:决策变量,若订单r分配给车辆k则取1,否则取0,r∈R,k∈K;
ATki;决策变量,车k到达站点i的时间,i∈V,k∈K;
HTki:决策变量,车k在站点i的等待时间,i∈V,k∈K;
zki:决策变量,车k到达站点i时所装载的乘客人数(即订单数),i∈V,k∈K;
其中,约束(5)-(8)是线路约束,用于确保每辆车可以匹配可行的路线;
约束(9)-(12)是车辆的到达时间限制,确保所有车辆能够有序、合理地完成运输任务;约束(13)-(17)是订单指派约束,确保每个订单被指派到可行的车辆线路与车辆;约束(18)-(24)是车辆的容量限制约束,确保装载乘客的数量不会超过车辆的容量;约束(25)是车辆的最大行驶时间约束;约束(26)是车队规模限制;约束条件(27)定义了决策变量的取值范围。
4.根据权利要求1多车型定制公交运营优化方法,所述将所述混合整数非线性规划模型中的非线性目标进行线性化处理,得到混合整数线性规划模型,具体包括:
分别对所述运营成本和对所述超时惩罚成本进行线性化处理。
5.根据权利要求1多车型定制公交运营优化方法,所述对所述运营成本进行线性化处理,具体包括:
所述对所述超时惩罚成本进行线性化处理,具体包括:
定义一个辅助变量actTrk∈R+∪{0},表示订单r在车辆k上的实际下车时间;
actTrk满足以下约束:
actTrk≥ATkd(r)-Z(1-αrk),k∈K,r∈R, (32)
actTrk≤ATkd(r)+Z(1-αrk),k∈K,r∈R, (33)
actTrk≤Zαrk,k∈K,r∈R, (34)
将所述超时惩罚成本CP用actTrk表示为:
6.根据权利要求3所述的多车型定制公交运营优化方法,所述使用Dantzig-Wolfe(D-W)分解所述混合整数线性规划模型进行分解重构,引入新的变量,得到限制主问题模型和子问题模型,具体包括:
引入新的集合Ω为定制公交的所有可行路径的集合,ω∈Ω;其中,Ωm表示使用m型车辆的路径集合;新的变量λωk为0-1变量,如果路径ω分配给车辆k则取1,否则取0;
将所述混合整数线性规划模型等价转化为如下的限制主问题模型:
主问题的目标函数为:
主问题约束条件为:
λωk∈{0,1},ω∈Ω,k∈K (40)
其中,目标(36)是最小化总成本,包括固定运营成本和可变运营成本以及超时惩罚成本;约束(37)表示订单必须全部被服务;约束(38)确保每辆车至多分配一条路线;约束(39)限制了m型车辆的可用数量;约束(40)为变量的取值范围;
将所述混合整数线性规划模型分解后得到子问题模型:
子问题目标函数如公式(41)所示:
xii=0,i∈N (45)
ATj≤ATi+HTi+tij+Z(1-xij),i,j∈V (46)
ATj≥ATi+HTi+tij-Z(1-xij),i,j∈V (47)
ATd(r)≥ATo(r)-Z(1-αr),r∈R (48)
ATo(r)≤(EDTr+tmax)αr+Z(1-αr),r∈R (52)
(EDTr-tmax)αr≤ATo(r)+HTo(r),r∈R (53)
za-Z(1-xsi)≤zi,i∈N (56)
za+Z(1-xsi)≥zi,i∈N (57)
zks=0,k∈K (60)
xij,ym,αr∈{0,1},ATi,HTi,zi≥0,i,j∈V,r∈R,m∈M (62)
7.根据权利要求2所述的多车型定制公交运营优化方法,所述启发式算法包括遗传算法,蚁群算法、禁忌搜索。
8.根据权利要求7所述的多车型定制公交运营优化方法,所述采用所述遗传算法求解大规模子问题,具体包括:
生成定制公交初始的可行路径的集合Ω,得到每条路径所使用的车型以及可变运营成本和超时惩罚成本;
调用求解器lpsolve对所述主问题进行求解,并计算所述主问题决策变量的对偶变量,使用对偶变量对子问题目标函数进行更新;其中,所述子问题的目标函数对应的值称之为主问题的检验数;
使用遗传算法对所述子问题进行求解,得到检验数和一条可行的路径以及它的时刻表,该路径为列生成方法中一个潜在的列;同时,得到这条路径所服务的订单以及所使用的车型;
基于子问题的解,判断检验数是否符合列生成的检验标准,若符合标准则向所述主问题中添加列,即将获得的路径作为一个元素添加到集合Ω中,并调用求解器lpsolve对所述主问题进行求解,并计算所述主问题决策变量的对偶变量,使用对偶变量对子问题目标函数进行更新;若不符合标准,则表明已经不存在有价值的列,则对当前所述主问题进行求解,得到的解即为原问题的最优解,输出定制公交的路线以及时刻表,还有每条路线所使用的车型和其服务的订单。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211069012.0A CN115409388B (zh) | 2022-09-02 | 2022-09-02 | 多车型定制公交运营优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211069012.0A CN115409388B (zh) | 2022-09-02 | 2022-09-02 | 多车型定制公交运营优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115409388A true CN115409388A (zh) | 2022-11-29 |
CN115409388B CN115409388B (zh) | 2023-04-07 |
Family
ID=84163551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211069012.0A Active CN115409388B (zh) | 2022-09-02 | 2022-09-02 | 多车型定制公交运营优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115409388B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115641704A (zh) * | 2022-12-26 | 2023-01-24 | 东风悦享科技有限公司 | 一种智能公交调度方法及系统 |
CN116757339A (zh) * | 2023-08-21 | 2023-09-15 | 倍施特科技(集团)股份有限公司 | 一种基于定制客运的运营优化方法及定制运输装置 |
CN117391564A (zh) * | 2023-10-19 | 2024-01-12 | 北京氢远质投新能源汽车有限公司 | 一种新能源物流车辆补能调度数据模型及调度优化方法 |
CN117455212A (zh) * | 2023-12-26 | 2024-01-26 | 武汉元光科技有限公司 | 一种响应公交的派单方法以及相关设备 |
CN117575292A (zh) * | 2024-01-15 | 2024-02-20 | 西南交通大学 | 基于人机协同决策的需求响应公交车辆柔性调度优化方法 |
CN117829713A (zh) * | 2024-03-05 | 2024-04-05 | 深圳市久通物联科技股份有限公司 | 一种基于混合整数线性规划的成品油智慧运输方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630440A (zh) * | 2009-06-01 | 2010-01-20 | 北京交通大学 | 常规公共交通衔接城市轨道交通的运营协调优化方法与系统 |
CN106886603A (zh) * | 2017-03-03 | 2017-06-23 | 东南大学 | 一种需求响应资源组合优化的分层分布式体系架构和方法 |
CN107330559A (zh) * | 2017-07-03 | 2017-11-07 | 华南理工大学 | 一种多起终点多车型混载的定制公交线路规划方法 |
US20180032964A1 (en) * | 2016-08-01 | 2018-02-01 | Nec Europe Ltd. | Transportation system and method for allocating frequencies of transit services therein |
CN112700029A (zh) * | 2020-12-03 | 2021-04-23 | 北京交通大学 | 一种基于仿真优化框架的定制公交规划方法 |
CN113159499A (zh) * | 2021-03-05 | 2021-07-23 | 北京化工大学 | 一种面向交通枢纽到达客流疏散的多出行模式组合调度方法 |
-
2022
- 2022-09-02 CN CN202211069012.0A patent/CN115409388B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630440A (zh) * | 2009-06-01 | 2010-01-20 | 北京交通大学 | 常规公共交通衔接城市轨道交通的运营协调优化方法与系统 |
US20180032964A1 (en) * | 2016-08-01 | 2018-02-01 | Nec Europe Ltd. | Transportation system and method for allocating frequencies of transit services therein |
CN106886603A (zh) * | 2017-03-03 | 2017-06-23 | 东南大学 | 一种需求响应资源组合优化的分层分布式体系架构和方法 |
CN107330559A (zh) * | 2017-07-03 | 2017-11-07 | 华南理工大学 | 一种多起终点多车型混载的定制公交线路规划方法 |
CN112700029A (zh) * | 2020-12-03 | 2021-04-23 | 北京交通大学 | 一种基于仿真优化框架的定制公交规划方法 |
CN113159499A (zh) * | 2021-03-05 | 2021-07-23 | 北京化工大学 | 一种面向交通枢纽到达客流疏散的多出行模式组合调度方法 |
Non-Patent Citations (1)
Title |
---|
郭晓俊: "基于需求响应的实时定制公交系统研究" * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115641704A (zh) * | 2022-12-26 | 2023-01-24 | 东风悦享科技有限公司 | 一种智能公交调度方法及系统 |
CN116757339A (zh) * | 2023-08-21 | 2023-09-15 | 倍施特科技(集团)股份有限公司 | 一种基于定制客运的运营优化方法及定制运输装置 |
CN116757339B (zh) * | 2023-08-21 | 2023-11-14 | 倍施特科技(集团)股份有限公司 | 一种基于定制客运的运营优化方法及定制运输装置 |
CN117391564A (zh) * | 2023-10-19 | 2024-01-12 | 北京氢远质投新能源汽车有限公司 | 一种新能源物流车辆补能调度数据模型及调度优化方法 |
CN117391564B (zh) * | 2023-10-19 | 2024-04-19 | 北京氢远质投新能源汽车有限公司 | 一种新能源物流车辆补能调度数据模型及调度优化方法 |
CN117455212A (zh) * | 2023-12-26 | 2024-01-26 | 武汉元光科技有限公司 | 一种响应公交的派单方法以及相关设备 |
CN117455212B (zh) * | 2023-12-26 | 2024-03-26 | 武汉元光科技有限公司 | 一种响应公交的派单方法以及相关设备 |
CN117575292A (zh) * | 2024-01-15 | 2024-02-20 | 西南交通大学 | 基于人机协同决策的需求响应公交车辆柔性调度优化方法 |
CN117575292B (zh) * | 2024-01-15 | 2024-03-26 | 西南交通大学 | 基于人机协同决策的需求响应公交车辆柔性调度优化方法 |
CN117829713A (zh) * | 2024-03-05 | 2024-04-05 | 深圳市久通物联科技股份有限公司 | 一种基于混合整数线性规划的成品油智慧运输方法 |
CN117829713B (zh) * | 2024-03-05 | 2024-05-10 | 深圳市久通物联科技股份有限公司 | 一种基于混合整数线性规划的成品油智慧运输方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115409388B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115409388B (zh) | 多车型定制公交运营优化方法 | |
CN109753694B (zh) | 基于全过程出行感知时间的中小城市公交线网设计方法 | |
CN102867408B (zh) | 一种公交出行路线的选择方法 | |
CN107944611B (zh) | 面向个性化出行需求的跨模式组合出行规划方法 | |
Chien et al. | Optimal feeder bus routes on irregular street networks | |
Chow et al. | Activity-based market equilibrium for capacitated multimodal transport systems | |
Chen et al. | Customized bus route design with pickup and delivery and time windows: Model, case study and comparative analysis | |
CN110245377B (zh) | 一种出行方案推荐方法及推荐系统 | |
CN108038578B (zh) | 基于需求预测和中心辐射网络的公共自行车静态调度方法 | |
Durán-Micco et al. | A survey on the transit network design and frequency setting problem | |
CN112561249B (zh) | 面向实时需求的城市定制公交调度方法 | |
CN110909434A (zh) | 一种低碳导向下的城市公共交通干线网络设计方法 | |
CN112309119B (zh) | 一种城市交通系统容量分析优化方法 | |
CN111144618A (zh) | 一种基于两阶段优化模型的需求响应型定制公交网络规划方法 | |
Sayarshad et al. | Optimizing dynamic switching between fixed and flexible transit services with an idle-vehicle relocation strategy and reductions in emissions | |
Oluwadare et al. | Solving school bus routing problem using genetic algorithm-based model | |
CN115936240B (zh) | 共享单车需求预测与投放调度方法 | |
Carlier et al. | A supernetwork approach towards multimodal travel modeling | |
Wei et al. | Optimal integrated model for feeder transit route design and frequency-setting problem with stop selection | |
Ma et al. | Dynamic vehicle routing problem for flexible buses considering stochastic requests | |
CN111882915A (zh) | 采用复合算法及交互模型的应需定制公交线路规划方法 | |
Zgraggen et al. | A model predictive control scheme for intermodal autonomous mobility-on-demand | |
Musolino et al. | A modelling framework to simulate paths and routes choices of freight vehicles in sub-urban areas | |
Liu et al. | Exploring the impact of spatiotemporal granularity on the demand prediction of dynamic ride-hailing | |
Jelen et al. | A multi-agent system for context-aware electric vehicle fleet routing: A step towards more sustainable urban operations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |