CN111582691B - 基于双层规划的客运枢纽多交通方式的运力匹配方法 - Google Patents
基于双层规划的客运枢纽多交通方式的运力匹配方法 Download PDFInfo
- Publication number
- CN111582691B CN111582691B CN202010350462.1A CN202010350462A CN111582691B CN 111582691 B CN111582691 B CN 111582691B CN 202010350462 A CN202010350462 A CN 202010350462A CN 111582691 B CN111582691 B CN 111582691B
- Authority
- CN
- China
- Prior art keywords
- passenger
- transportation
- model
- layer
- double
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000009826 distribution Methods 0.000 claims abstract description 25
- 230000029305 taxis Effects 0.000 claims abstract description 25
- 230000002068 genetic effect Effects 0.000 claims abstract description 24
- 230000035772 mutation Effects 0.000 claims description 10
- 238000011160 research Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010187 selection method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 210000000349 chromosome Anatomy 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008436 biogenesis Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 230000008303 genetic mechanism Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Bioinformatics & Computational Biology (AREA)
- Genetics & Genomics (AREA)
- Data Mining & Analysis (AREA)
- Educational Administration (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Quality & Reliability (AREA)
- Physiology (AREA)
- Game Theory and Decision Science (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Development Economics (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于双层规划的客运枢纽多交通方式的运力匹配方法。该方法包括:建立客运枢纽的各交通方式的广义费用函数的具体表达式,构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,上层子模型调整大巴车及轨道交通线路的发车间隔及出租车单位时间可服务数量,下层子模型为上层子模型提供客流分配结果;基于遗传算法和MSA算法求解双层规划模型,输出适应度函数值最大的个体作为双层规划模型的求解结果。本发明针对客运枢纽旅客需求精确匹配轨道交通、大巴车、出租车的运力,提高客运枢纽各交通方式运能利用率,提升陆侧交通***的整体集散效率,改善旅客换乘体验。
Description
技术领域
本发明涉及交通运行管理技术领域,尤其涉及一种基于双层规划的客运枢纽多交通方式的运力匹配方法。
背景技术
目前大型客运枢纽仅针对单一交通方式或两种交通方式间的运力匹配开展优化研究,几乎未见三种及以上交通方式之间的运力匹配研究。且大型客运枢纽集散交通运力匹配优化的研究大都集中于铁路客运枢纽,但不同的客运枢纽在旅客构成、出行特性及出行需求等方面存在着差异不可一概而论,对于客运枢纽交通***运力匹配的研究还十分匮乏。
大型客运枢纽是多种交通方式连接的节点,目前已有研究往往较少考虑各交通方式间的竞争与合作关系,忽略了多交通方式间的相互影响作用,因此,急需开发一种将各集散交通子***统一为一个整体进行运力资源优化配置的方法。
发明内容
本发明的实施例提供了一种基于双层规划的客运枢纽多交通方式的运力匹配方法,以克服现有技术的问题。
为了实现上述目的,本发明采取了如下技术方案。
一种基于双层规划的客运枢纽多交通方式的运力匹配方法,包括:
建立客运枢纽的各交通方式的广义费用函数的具体表达式,所述各种交通方式包括轨道交通、大巴车和出租车;
构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,所述上层子模型调整大巴车及轨道交通线路的发车间隔及出租车单位时间可服务数量,所述下层子模型为所述上层子模型提供客流分配结果;
基于遗传算法和MSA算法求解所述双层规划模型,输出适应度函数值最大的个体作为双层规划模型的求解结果。
优选地,所述的建立客运枢纽的各交通方式的广义费用函数的具体表达式,所述各种交通方式包括轨道交通、大巴车和出租车,包括:
客运枢纽的经营成本W1的计算方法为:
式中:ω1,ω2,ω3,ω4——分别表示每辆大巴车、每列轨道交通列车、每辆出租车、每个公共停车位经营成本参数;
[]——取整符号;
T——研究时段总时长;
X——每分钟可服务的出租车数量;
q4——选择私家车出行的旅客总数;
cap——私家车平均载客量;
旅客候车成本W2的计算方法为:
式中:γM——货币成本转化为时间成本的折算系数,这里γM未添加表示交通方式的下标是基于所有旅客对时间与货币费用间感知一致的假设;
综合客运枢纽的经营成本W1和旅客候车成本W2的上层子模型的目标函数为:
式中:[]——取整符号;
α1——客运枢纽的运营部门的经营成本的权重系数;
α2——各交通方式出行成本的权重系数。
优选地,所述的构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,所述上层子模型调整大巴车及轨道交通线路的发车间隔及出租车单位时间可服务数量,所述下层子模型为所述上层子模型提供客流分配结果,包括:
构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,设定所述上层子模型的约束条件包括:
1)公共交通满载率约束,包括轨道交通、大巴车在内的公共交通的每班次车辆满载率均不得超过1,假设旅客到达公共交通站台时间服从均匀分布,则公共交通满载率的约束为:
[]——取整符号;
2)出租车时段供给能力约束,出租车到达客运枢纽站点时的剩余运力不得小于需求:
μX*T≥q3
式中:X——每分钟可服务的出租车数量;
μ——每辆出租车平均服务强度;
3)公共交通发车间隔约束,即发车间隔必须在线路允许的范围内:
Hi,min≤Hi≤Hi,max
基于所述上层子模型的目标条件和约束条件,所述上层子模型表示如下:
所述上层子模型将获得的决策变量赋给所述下层子模型,所述下层子模型按照SUE-logit模型对各交通方式客流量进行分配,达到平衡状态时,起讫点r-s间所有已被选择交通方式的客流量应满足以下平衡条件:
起讫点r-s间选择第i种交通方式的概率:
式中:Pi rs——起讫点r-s间选择第i种交通方式的概率
θ——效用函数与费用函数的对换参数(θ>0)
起讫点r-s间各交通方式分配的客流量:
式中:drs——起讫点r-s间的客流量
每个起讫点r-s间各交通方式客流量之和等于该起讫点的客流量,且各交通方式客流量均为非负:
各起讫点的客流量之和等于研究时段内进港旅客总数:
式中:D——研究时段内进港旅客总数。
优选地,所述的基于遗传算法和MSA算法求解所述双层规划模型,输出适应度函数值最大的个体作为双层规划模型的求解结果,包括:
基于遗传算法和MSA算法通过MATLAB软件求解双层规划模型,包括如下的处理步骤:
Step1变量编码及参数设定,将公共交通线路发车间隔及出租车单位时间可服务数量作为遗传算法中每条染色体的基因进行编码,对初始种群大小pop_size,交叉概率Pc、变异概率Pm和最大迭代次数num_iter进行参数设置;
Step2种群初始化,采用随机初始化的方式产生若干个满足种群规模的个体,每个个体对应由不同的公共交通线路发车间隔及出租车单位时间可服务数量构成的运力配置,每个个体的取值均从其对应的范围内随机产生;
Step3进行客流分配,对每个个体按SUE-logit模型进行客流分配得到平衡态下各交通方式的客流量;
Step4计算适应度函数值,将上层规划目标函数的负数作为适应度函数,计算种群中每个个体的适应度函数值;
Step5选择操作,采用锦标赛选择法经常选择,锦标赛选择法的基本思想为,将n个个体作为一个次种群,每次选取每个次种群中适应度最高的那个个体复制到下一代群体中,一般n取2;
Step6交叉操作,采用单点交叉法,即在满***叉概率Pc的情况下,在随机配对的父代个体中随机选择一个交叉点,将位于交叉点之前或之后的那部分变量进行交换,以形成新的后代个体;
Step7变异操作,采用基本位变异的方法,以变异概率Pm随机指定某一位或某几位基因座上的值做变异运算,形成新的后代个体;
Step8收敛性检查,判断是否满足算法终止条件,以达到最大迭代次数作为终止条件,若满足终止条件,算法结束,输出适应度函数值最大的个体作为双层规划模型的求解结果,所述适应度函数值最大的个体为上层规划中目标函数值最小的求解结果,即旅客候车成本和运营部门经营成本最小的运力匹配方式;若不满足终止条件,跳转至Step3。
由上述本发明的实施例提供的技术方案可以看出,本发明实施例的方法针对客运枢纽旅客需求精确匹配轨道交通、大巴车、出租车的运力,提高客运枢纽各交通方式运能利用率,提升陆侧交通***的整体集散效率,改善旅客换乘体验。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于双层规划的客运枢纽多交通方式的运力匹配方法的处理流程图;
图2为本发明实施例提供的一种基于遗传算法和MSA算法,通过MATLAB软件求解双层规划模型的处理流程图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。
针对空港枢纽客流集散引发的动态供需平衡有效调控问题,为解决轨道交通、大巴车、出租车等多交通方式运力与航空旅客需求难以精准匹配的难题,本发明实施例构建了基于双层规划的多交通方式运力资源匹配模型,并设计求解算法求解模型,所得结果可实现空港枢纽多交通方式运力匹配优化。
本发明实施例提出的一种基于双层规划的客运枢纽多交通方式的运力匹配方法的处理流程图如图1所示,包括如下的处理步骤:
步骤S1:分析了客运枢纽中各种交通方式的广义费用,构建了各交通方式的广义费用函数的具体表达式。上述各种交通方式包括轨道交通、大巴车和出租车,上述客运枢纽可以为空港枢纽、铁路枢纽、汽车客运枢纽和水路客运枢纽等。
步骤S2:构建客运枢纽多交通方式运力匹配的双层规划模型,双层规划模型中上层子模型调整大巴车及轨道交通线路发车间隔及出租车单位时间可服务数量,优化***总成本;下层子模型基于随机用户均衡理论和Logit模型,为上层子模型提供客流分配结果。
步骤S3:基于遗传算法和MSA算法求解上述双层规划模型,输出适应度函数值最大的个体作为双层规划模型的求解结果。
在本发明的一种实施方式中,上述步骤S2中的客运枢纽多交通方式运力匹配的双层规划模型的基本形式如下:
式中,x∈Rn,y∈Rm分别为上下层规划问题的决策变量,F,f:Rn+m→R,G:Rn+m→Rp,g:Rn+m→Rq。
由上式可以看出,双层规划模型包含上、下两个子模型,两个子模型均拥有各自的目标函数和约束条件,并通过决策变量建立相互关系。就独立模型而言,两个子模型都有符合自身目标函数和约束条件的最优解,但由于上层子模型的目标函数和约束条件不仅与上层决策变量有关,还依赖于下层子模型的最优解,而下层子模型的最优解又受上层决策变量的影响,因此需设定相应的算法获得同时符合两个子模型要求的满意解。
基于上述分析,本发明实施例中,上层子模型以***总成本最小为优化目标,优化轨道交通线路发车间隔及出租车单位时间可服务数量,***总成本由各交通方式旅客候车成本和运营部门经营成本两部分构成;下层子模型则基于随机用户均衡(Stochastic UserEquilibrium,SUE)理论和Logit模型,为上层子模型提供各交通方式最优客流分配结果。
上述步骤S1中的各交通方式的广义费用函数的表达方式如下:
(1)客运枢纽的经营成本W1
运营部门经营成本主要包括车辆及设备设施的整备检修、维护更换,以及动力成本、人力劳务等,可通过调研得到,按单位时间单位车辆折算,因此客运枢纽的运营部门的经营成本的计算方法为:
式中:ω1,ω2,ω3,ω4——分别表示每辆大巴车、每列轨道交通列车、每辆出租车、每个公共停车位经营成本参数;
[]——取整符号;
T——研究时段总时长(min);
X——每分钟可服务的出租车数量(辆);
q4——选择私家车出行的旅客总数(人次);
cap——私家车平均载客量(人/辆);
(2)旅客候车成本W2
对客运枢纽各交通方式进行运力匹配优化的最大目的是为了均衡利用运力资源尽快疏散客流,从而降低旅客在客运枢纽内的滞留时间,因此优化目标考虑旅客候车成本:
式中:γM——货币成本转化为时间成本的折算系数,这里γM未添加表示交通方式的下标是基于所有旅客对时间与货币费用间感知一致的假设;
因此,综合客运枢纽的经营成本W1和旅客候车成本W2的上层子模型的目标函数为:
式中:[]——取整符号;
α1——客运枢纽的运营部门经营成本的权重系数;
α2——各交通方式出行成本的权重系数。
(3)上层子模型的约束条件
1)公共交通满载率约束,公共交通(轨道交通、大巴车等)每班次车辆满载率均不得超过1。由于假设旅客到达公共交通站台时间服从均匀分布,则公共交通满载率的约束为:
[]——取整符号;
2)出租车时段供给能力约束,出租车到达客运枢纽站点时的剩余运力不得小于需求:
μX*T≥q3
式中:X——每分钟可服务的出租车数量(辆);
μ——每辆出租车平均服务强度(人次)。
3)公共交通(i=1,2)发车间隔约束,即发车间隔必须在线路允许的范围内:
Hi,min≤Hi≤Hi,max
综上所述,基于所述上层子模型的目标条件和约束条件,上述上层子模型可表示如下:
上述下层子模型的表达方式如下:
在对公共交通线路发车间隔及出租车单位时间可服务数量进行调整时,各交通方式的广义费用将发生变化,从而使分配的客流量发生变化,进而影响***总成本。因此,上层子模型将获得的决策变量赋给下层子模型,下层子模型按照SUE-logit模型对各交通方式客流量进行分配,达到平衡状态时,起讫点r-s间所有已被选择交通方式的客流量应满足以下平衡条件:
起讫点r-s间选择第i种交通方式的概率:
θ——效用函数与费用函数的对换参数(θ>0)
起讫点r-s间各交通方式分配的客流量:
式中:drs——起讫点r-s间的客流量(人次)
每个起讫点r-s间各交通方式客流量之和等于该起讫点的客流量,且各交通方式客流量均为非负:
各起讫点的客流量之和等于研究时段内进港旅客总数:
式中:D——研究时段内进港旅客总数(人次)
在本发明的一种实施方式中,上述步骤S3具体包括:
(1)遗传算法
遗传算法(Genetic Algorithm)也称为进化算法,是一种通过模拟达尔文生物进化论的自然选择、模拟遗传学机理的生物进化过程的搜索最优解的方法。目前它被广泛地应用于自动控制、规划设计、组合优化等多个领域。
遗传算法一般由以下5个部分组成:
①参数的编码与解码:将变量编码为一种编码串,把变量的可行解从其自身的解空间转换到遗传算法所能处理的解空间中去,解码的过程与编码正好相反。
②初始化种群:遗传算法是对群体进行进化操作的,因此需要在开始为其准备一些表示起始搜索点的初始种群数据。
③适应度函数:在遗传算法中,各个个体的优劣程度需通过个体的适应度大小来判断,从而决定该个体是否被遗传。度量个体适应度大小的函数就称为适应度函数,适应度函数是目标函数的映射。
④遗传算子:初始群体组成后,遗传操作的任务就是对群体的个体按照它们对“环境”适应的程度施加一定的操作,从而实现优胜劣汰的进化过程。遗传算法的遗传操作包括选择、交叉和变异等3个基本遗传算子。
⑤算法终止条件:遗传算法的终止条件大致可分为两种,一种是达到给定的最大迭代次数;另一种是在一个种群中,适应度最优的解决方案和适应度最劣的的解决方案之间的差异小于一个给定值。
(2)MSA算法
相继平均算法(Methodof SuccessiveAverages,MSA)是求解随机用户均衡(Stochastic User Equilibrium,SUE)问题中使用最为广泛的算法,它的基本思想是沿下降方向预先确定步长,将上次迭代中各路段的流量与本次迭代中各路段所获得的附加流量进行加权平均,得到本次迭代中各路段的交通流量,当两次迭代的各路段交通流量差值小于规定值时,算法停止。
MSA算法的具体步骤如下[48-49]:
④更新路段流量:沿下降方向按预先确定的步长“相继平均”移动,更新各路段流量:
⑤检查收敛性:若已满足收敛指标要求,则终止算法。否则令n=n+1,转到第②步。
2.2双层规划模型的求解算法设计
在上述双层规划模型中,上层子模型以公共交通(机场大巴、轨道交通)线路发车间隔及出租车单位时间可服务数量为优化变量,下层子模型基于SUE-logit模型进行客流分配。上层子模型的求解结果为下层子模型提供输入参数,下层子模型的客流分配结果为上层子模型提供客流输入数据。上层子模型和下层子模型相互影响,最终获得符合目标函数的最优解。
本发明实施例提供的一种基于遗传算法和MSA算法,通过MATLAB软件求解双层规划模型的处理流程如图2所示,包括如下的处理步骤:
Step1变量编码及参数设定。将公共交通线路发车间隔及出租车单位时间可服务数量作为遗传算法中每条染色体的基因进行编码。参数设置包括初始种群大小pop_size,交叉概率Pc、变异概率Pm以及最大迭代次数num_iter等。
Step2种群初始化。采用随机初始化的方式产生若干个满足种群规模的个体,每个个体对应由不同的公共交通线路发车间隔及出租车单位时间可服务数量构成的运力配置。每个个体的取值均从其对应的范围内随机产生。
Step3进行客流分配。对每个个体按SUE-logit模型进行客流分配得到平衡态下各交通方式的客流量。
Step4计算适应度函数值。将上层规划目标函数的负数作为适应度函数,计算种群中每个个体的适应度函数值。
Step5选择操作。采用锦标赛选择法经常选择,锦标赛选择法的基本思想为,将n个个体作为一个次种群,每次选取每个次种群中适应度最高的那个个体复制到下一代群体中,一般n取2。
Step6交叉操作。采用单点交叉法,即在满***叉概率Pc的情况下,在随机配对的父代个体中随机选择一个交叉点,将位于交叉点之前或之后的那部分变量进行交换,以形成新的后代个体。
Step7变异操作。采用基本位变异的方法,以变异概率Pm随机指定某一位或某几位基因座上的值做变异运算,形成新的后代个体。
Step8收敛性检查。判断是否满足算法终止条件,以达到最大迭代次数作为终止条件。若满足终止条件,算法结束,输出适应度函数值最大的个体作为双层规划模型的求解结果,上述适应度函数值最大的个体为上层规划中目标函数值最小的求解结果,即旅客候车成本和运营部门经营成本最小的运力匹配方式。若不满足终止条件,跳转至Step3。
综上所述,本发明实施例的方法针对客运枢纽旅客需求精确匹配轨道交通、大巴车、出租车的运力,提高客运枢纽各交通方式的运能利用率,提升陆侧交通***的整体集散效率,改善旅客换乘体验。
本发明实施例的技术方案模型构建符合实际条件,求解方法高效准确,能够用于空港枢纽轨道交通、大巴车、出租车多交通方式间的运力资源优化配置。
该方法能本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或***实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置及***实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (1)
1.一种基于双层规划的客运枢纽多交通方式的运力匹配方法,其特征在于,包括:
建立客运枢纽的各交通方式的广义费用函数的具体表达式,所述各种交通方式包括轨道交通、大巴车和出租车;
构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,所述上层子模型调整大巴车及轨道交通线路的发车间隔及出租车单位时间可服务数量,所述下层子模型为所述上层子模型提供客流分配结果;
基于遗传算法和MSA算法求解所述双层规划模型,输出适应度函数值最大的个体作为双层规划模型的求解结果;
所述的建立客运枢纽的各种交通方式的广义费用函数的具体表达式,所述各种交通方式包括轨道交通、大巴车和出租车,包括:
客运枢纽的经营成本W1的计算方法为:
式中:ω1,ω2,ω3,ω4——分别表示每辆大巴车、每列轨道交通列车、每辆出租车、每个公共停车位经营成本参数;
[]——取整符号;
T——研究时段总时长;
X——每分钟可服务的出租车数量;
q4——选择私家车出行的旅客总数;
cap——私家车平均载客量;
旅客候车成本W2的计算方法为:
式中:γM——货币成本转化为时间成本的折算系数,这里γM未添加表示交通方式的下标是基于所有旅客对时间与货币费用间感知一致的假设;
综合客运枢纽的经营成本W1和旅客候车成本W2的上层子模型的目标函数为:
式中:[]——取整符号;
α1——客运枢纽的运营部门的经营成本的权重系数;
α2——各交通方式出行成本的权重系数;
所述的构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,所述上层子模型调整大巴车及轨道交通线路的发车间隔及出租车单位时间可服务数量,所述下层子模型为所述上层子模型提供客流分配结果,包括:
构建客运枢纽多交通方式运力匹配的双层规划模型,该双层规划模型包括上层子模型和下层子模型,设定所述上层子模型的约束条件包括:
1)公共交通满载率约束,包括轨道交通、大巴车在内的公共交通的每班次车辆满载率均不得超过1,旅客到达公共交通站台时间服从均匀分布,则公共交通满载率的约束为:
[]——取整符号;
2)出租车时段供给能力约束,出租车到达客运枢纽站点时的剩余运力不得小于需求:
μX*T≥q3
式中:X——每分钟可服务的出租车数量;
μ——每辆出租车平均服务强度;
3)公共交通发车间隔约束,即发车间隔必须在线路允许的范围内:
Hi,min≤Hi≤Hi,max
基于所述上层子模型的目标条件和约束条件,所述上层子模型表示如下:
所述上层子模型将获得的决策变量赋给所述下层子模型,所述下层子模型按照SUE-logit模型对各交通方式客流量进行分配,达到平衡状态时,起讫点r-s间所有已被选择交通方式的客流量应满足以下平衡条件:
起讫点r-s间选择第i种交通方式的概率:
θ——效用函数与费用函数的对换参数,θ>0
起讫点r-s间各交通方式分配的客流量:
式中:drs——起讫点r-s间的客流量
每个起讫点r-s间各交通方式客流量之和等于该起讫点的客流量,且各交通方式客流量均为非负:
各起讫点的客流量之和等于研究时段内进港旅客总数:
式中:D——研究时段内进港旅客总数;
所述的基于遗传算法和MSA算法求解所述双层规划模型,输出适应度函数值最大的个体作为双层规划模型的求解结果,包括:
基于遗传算法和MSA算法通过MATLAB软件求解双层规划模型,包括如下的处理步骤:
Step1变量编码及参数设定,将公共交通线路发车间隔及出租车单位时间可服务数量作为遗传算法中每条染色体的基因进行编码,对初始种群大小pop_size,交叉概率Pc、变异概率Pm和最大迭代次数num_iter进行参数设置;
Step2种群初始化,采用随机初始化的方式产生若干个满足种群规模的个体,每个个体对应由不同的公共交通线路发车间隔及出租车单位时间可服务数量构成的运力配置,每个个体的取值均从其对应的范围内随机产生;
Step3进行客流分配,对每个个体按SUE-logit模型进行客流分配得到平衡态下各交通方式的客流量;
Step4计算适应度函数值,将上层规划目标函数的负数作为适应度函数,计算种群中每个个体的适应度函数值;
Step5选择操作,采用锦标赛选择法经常选择,锦标赛选择法的基本思想为,将n个个体作为一个次种群,每次选取每个次种群中适应度最高的那个个体复制到下一代群体中,n取2;
Step6交叉操作,采用单点交叉法,即在满***叉概率Pc的情况下,在随机配对的父代个体中随机选择一个交叉点,将位于交叉点之前或之后的那部分变量进行交换,以形成新的后代个体;
Step7变异操作,采用基本位变异的方法,以变异概率Pm随机指定某一位或某几位基因座上的值做变异运算,形成新的后代个体;
Step8收敛性检查,判断是否满足算法终止条件,以达到最大迭代次数作为终止条件,若满足终止条件,算法结束,输出适应度函数值最大的个体作为双层规划模型的求解结果,所述适应度函数值最大的个体为上层规划中目标函数值最小的求解结果,即旅客候车成本和运营部门经营成本最小的运力匹配方式;若不满足终止条件,跳转至Step3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010350462.1A CN111582691B (zh) | 2020-04-28 | 2020-04-28 | 基于双层规划的客运枢纽多交通方式的运力匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010350462.1A CN111582691B (zh) | 2020-04-28 | 2020-04-28 | 基于双层规划的客运枢纽多交通方式的运力匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111582691A CN111582691A (zh) | 2020-08-25 |
CN111582691B true CN111582691B (zh) | 2022-04-22 |
Family
ID=72124575
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010350462.1A Active CN111582691B (zh) | 2020-04-28 | 2020-04-28 | 基于双层规划的客运枢纽多交通方式的运力匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111582691B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111950814A (zh) * | 2020-09-01 | 2020-11-17 | 吉林大学 | 一种基于广义费用函数的客运枢纽换乘流线优化方法 |
CN112419716B (zh) * | 2020-11-13 | 2021-09-28 | 东南大学 | 一种轨道站点换乘影响区内的共享单车设施布局配置方法 |
CN112529288B (zh) * | 2020-12-10 | 2023-10-24 | 五邑大学 | 一种多列车差异化定价策略优化方法及*** |
CN113159499B (zh) * | 2021-03-05 | 2022-05-10 | 北京化工大学 | 一种面向交通枢纽到达客流疏散的多出行模式组合调度方法 |
CN112907144A (zh) * | 2021-03-31 | 2021-06-04 | 中仁车汇科技发展(深圳)有限公司 | 一种密集人群的疏散方法 |
CN113420926A (zh) * | 2021-06-29 | 2021-09-21 | 北京交通大学 | 一种城市枢纽间多交通方式运力协同优化方案编制方法 |
CN113806683B (zh) * | 2021-08-09 | 2024-06-14 | 北京交通大学 | 一种大型体育赛事服务人员需求测算与组织调度方法 |
CN113361984B (zh) * | 2021-08-11 | 2021-12-21 | 北京航空航天大学杭州创新研究院 | 基于互馈模型的空港客流与运力流动态耦合方法和*** |
CN115527369B (zh) * | 2022-09-29 | 2023-06-06 | 北京交通大学 | 一种空港枢纽大面积延误条件下的大客流预警与疏散方法 |
CN115409295B (zh) * | 2022-11-01 | 2023-02-14 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种基于瓶颈分析的公交调度方法、电子设备及存储介质 |
CN115620525B (zh) * | 2022-12-16 | 2023-03-10 | 中国民用航空总局第二研究所 | 基于时变动态贝叶斯网络的短时交通客运需求预测方法 |
CN116432888A (zh) * | 2023-06-15 | 2023-07-14 | 中铁第四勘察设计院集团有限公司 | 基于od对备选路径集合的高铁枢纽流线优化方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107330716A (zh) * | 2017-06-01 | 2017-11-07 | 华南理工大学 | 一种考虑***最优的定制公交定价方法 |
WO2018032808A1 (zh) * | 2016-08-19 | 2018-02-22 | 大连理工大学 | 基于大数据的公交线路时刻表协同优化方法 |
CN109409560A (zh) * | 2018-08-16 | 2019-03-01 | 北京交通大学 | 基于多智能体仿真的城市轨道交通客流诱导方法 |
-
2020
- 2020-04-28 CN CN202010350462.1A patent/CN111582691B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018032808A1 (zh) * | 2016-08-19 | 2018-02-22 | 大连理工大学 | 基于大数据的公交线路时刻表协同优化方法 |
CN107330716A (zh) * | 2017-06-01 | 2017-11-07 | 华南理工大学 | 一种考虑***最优的定制公交定价方法 |
CN109409560A (zh) * | 2018-08-16 | 2019-03-01 | 北京交通大学 | 基于多智能体仿真的城市轨道交通客流诱导方法 |
Non-Patent Citations (3)
Title |
---|
"Simulated Annealing Algorithm for Solving A Bi-level Optimization Model on High-speed Railway Station Location";Xiang Rui, Lu Gongyuan, He Dong;《2010 Third International Conference on Information and Computing》;20101231;全文 * |
"空港枢纽陆侧交通运力匹配与换乘仿真优化研究";宋悦;《中国优秀硕博士学位论文全文数据库(硕士) 工程科技II辑》;20210315;全文 * |
"铁路综合客运枢纽衔接方式运营匹配性评价及优化";南天伟;《中国优秀硕博士学位论文全文数据库(硕士) 工程科技II辑》;20151015;第36-51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111582691A (zh) | 2020-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111582691B (zh) | 基于双层规划的客运枢纽多交通方式的运力匹配方法 | |
CN109409560B (zh) | 基于多智能体仿真的城市轨道交通客流诱导方法 | |
CN107092976B (zh) | 一种多目标模型协同优化多条公交线路发车间隔的方法 | |
CN111862579A (zh) | 一种基于深度强化学习的出租车调度方法及*** | |
CN109118111B (zh) | 尾随间隔限制及起飞时隙分配综合策略管理决策支持*** | |
CN114912736A (zh) | 一种电动公交协调优化调度方法 | |
CN115527369B (zh) | 一种空港枢纽大面积延误条件下的大客流预警与疏散方法 | |
CN109118023B (zh) | 一种公共交通线网优化方法 | |
CN111445048A (zh) | 一种响应型接驳公交分时段协调优化方法 | |
CN111967134A (zh) | 基于地铁共线共车的地下物流***优化控制方法 | |
CN112561249A (zh) | 面向实时需求的城市定制公交调度方法 | |
CN115809729A (zh) | 考虑新增线路的城市轨道交通枢纽换乘设施优化配置方法 | |
CN111126712B (zh) | 一种面向通勤廊道的停车充电换乘***规划方法 | |
CN113420926A (zh) | 一种城市枢纽间多交通方式运力协同优化方案编制方法 | |
Du et al. | Modeling the capacity of multimodal and intermodal urban transportation networks that incorporate emerging travel modes | |
CN115222156A (zh) | 基于分时对偶路网考虑用户需求响应的汽车充电调度方法 | |
CN108830401B (zh) | 基于元胞传输模型的动态拥堵收费最优费率计算方法 | |
CN112598328A (zh) | 卫星厅模式下中转登机口多目标分配的优化方法及*** | |
CN116739213A (zh) | 基于代理模型辅助算法的地铁接驳公交优化方法 | |
CN112633596A (zh) | 地铁列车速度曲线和区间运行时间的一体化优化方法 | |
CN115481777A (zh) | 面向多线路公交动态时刻表协同仿真优化方法、设备及介质 | |
CN114254822B (zh) | 一种基于物联网技术的无人机配送网络优化模型及其求解算法 | |
CN116562538A (zh) | 一种考虑轨道交通接驳的微循环公交调度优化方法及装置 | |
CN115169696A (zh) | 一种人工与自动驾驶混行下的地铁接驳公交智能调度方法 | |
CN115511336A (zh) | 一种基于多点联合分配的共享单车停车点选址方法及*** |
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 |