CN113177295B - 一种地铁网络列车时刻表快速编制方法 - Google Patents

一种地铁网络列车时刻表快速编制方法 Download PDF

Info

Publication number
CN113177295B
CN113177295B CN202110370870.8A CN202110370870A CN113177295B CN 113177295 B CN113177295 B CN 113177295B CN 202110370870 A CN202110370870 A CN 202110370870A CN 113177295 B CN113177295 B CN 113177295B
Authority
CN
China
Prior art keywords
train
line
station
sub
time
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
Application number
CN202110370870.8A
Other languages
English (en)
Other versions
CN113177295A (zh
Inventor
李树凯
胡雨婷
杨立兴
许永成
王悉
戚建国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jiulian Technology Co ltd
Beijing Jiaotong University
Original Assignee
Beijing Jiulian Technology Co ltd
Beijing Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Jiulian Technology Co ltd, Beijing Jiaotong University filed Critical Beijing Jiulian Technology Co ltd
Priority to CN202110370870.8A priority Critical patent/CN113177295B/zh
Publication of CN113177295A publication Critical patent/CN113177295A/zh
Application granted granted Critical
Publication of CN113177295B publication Critical patent/CN113177295B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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"
    • G06Q50/40—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/02—CAD in a network environment, e.g. collaborative CAD or distributed simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD

Abstract

本发明属于城市轨道交通控制技术领域，涉及一种地铁网络列车时刻表快速编制方法，包括：S1：设置目标地铁网络的基础参数；S2：设置各线路区间的运行时间、停站时间、发车间隔和旅行时间的范围；采集现实数据，设置各换乘车站站台间的换乘走行时间；S3：采集地铁网络的客流信息；S4：采用广义Benders分解算法进行编制；S5：获得大型地铁网络的协同优化列车时刻表和线路间列车接续情况。本发明采用广义Benders分解算法将原问题分解为更易求解的子问题和主问题，具有较好的求解效率，且能够最大程度减少乘客等待时间，提高地铁系统的换乘效率和服务水平。

Description

一种地铁网络列车时刻表快速编制方法
技术领域
本发明属于城市轨道交通控制技术领域，涉及一种地铁网络列车时刻表快速编制方法，尤其涉及一种地铁网络下基于广义Benders分解算法的地铁网络列车时刻表快速编制计算方法。
背景技术
随着全世界城市化和机动化进程的加剧，城市中交通拥堵、意外事故和尾气排放所引起的交通污染等问题越来越严重。为缓解城市交通拥堵，倡导低碳环保出行，优先发展城市公共交通已成为国内外广泛采用的重要举措。作为公共交通的重要组成部分，城市轨道交通具有运量大、速度快、节能环保和准时可靠等特点，而受到各大城市的青睐。随着居民出行的急剧增加和出行需求的多样化，在城市轨道交通网络中的乘客为完成出行，往往从起始点到目的地，有时需要完成一次(或多次)不同线路之间的接驳转运，即换乘。对城市轨道交通换乘站来说，换乘客流比例较高，而换乘效率直接影响乘客出行时间和乘客对地铁系统服务水平的满意度。高质量的换乘可以减少出行时间，增强地铁的吸引力，而不合理的换乘则对出行的便捷性、舒适性和出行需求产生不利影响。随着北京和上海等城市轨道交通网络的快速扩大，网络运营提高了地铁的可达性和覆盖范围，同时也存在大量的换乘客流需求。如果线路间列车运行计划不相互协调，会给乘客换乘带来不便，从而产生较长的等待时间，降低换乘的衔接性和流畅性。因此，研究地铁网络中列车衔接和时刻表优化问题，对地铁运营管理具有重要的指导作用和现实意义。
列车时刻表优化是轨道交通网络运营管理的重要组成部分，通过调整线路上列车的到发时间等，使不同线路上的列车同步到达换乘站，或换乘乘客下车后走行至换乘站台赶上接续列车。通过优化后的时刻表，乘客可以体验无缝换乘衔接和良好的地铁系统服务水平，从而实现换乘乘客和进站乘客等待时间的最小化。
面对具有大量线路、车站和列车的大规模地铁网络，往往需要较长的求解时间来编制时刻表。而在城市轨道交通运营中，由于工作日、节假日、大型活动和疫情爆发期间的客流特性不同，需根据历史或预测的客流特征来编制新时刻表，这给地铁运营公司带来了巨大的计算负担。因此，针对大规模地铁网络的列车时刻表编制问题提出一种计算效率高的优化算法是很有必要的。
本发明提供一种地铁网络列车时刻表快速编制计算方法，以期最大程度减少换乘乘客和进站乘客的总等待时间，并提高线路间换乘效率和客运服务水平。
发明内容
面对大规模地铁网络问题时计算负荷大的情况，本发明公开了一种新的基于广义Benders分解算法的时刻表快速编制计算方法，通过在合理范围内调整两列接续列车的到发时间，从而达到快速计算得到协同优化的时刻表，最大程度减少进站乘客和换乘乘客等待时间的目标。
为了解决上述问题，本发明提出一种基于广义Benders分解算法的地铁网络列车时刻表快速编制方法，具体技术方案如下：
一种地铁网络列车时刻表快速编制方法，包括以下步骤：
S1：设置目标地铁网络的基础参数；
所述基础参数包括：线路数量、每条线路的车站数和目标地铁网络中不同线路间换乘车站情况等；
S2：设置各线路区间的运行时间、停站时间、发车间隔和旅行时间的范围；
采集现实数据，设置各换乘车站站台间的换乘走行时间；
S3：采集地铁网络的客流信息；
所述地铁网络的客流信息包括：换乘车站的换乘乘客人数和进站的乘客人数；
S4：为了快速编制大规模地铁网络列车时刻表，采用广义Benders分解算法进行编制，具体包括以下步骤：
S41：给定计划列车数量和时域；
S42：根据目标大型地铁网络(即大规模地铁网络)相关参数和客流信息等，建立时刻表协同优化模型；
S43：利用广义Benders分解算法，将步骤S42构建的时刻表协同优化模型，分解为容易求解的子问题和主问题，并迭代求解子问题和主问题，直至收敛至最优解；
S5：由步骤S4得到最优解，获得大型地铁网络的协同优化列车时刻表和线路间列车接续情况。
在上述技术方案的基础上，步骤S42的具体步骤如下：
S421、构建时刻表线性约束：
利用公式(1)得到各线路上每列车到达各车站的时间，利用公式(2)得到各线路上每列车离开各车站的时间，
其中：i是列车标号；s是车站标号；m是线路标号；M表示：目标大型地铁网络中所有线路标号的集合；Sm表示：目标大型地铁网络中线路m中所有车站标号的集合；Nm表示：目标大型地铁网络中线路m中所有列车标号的集合；表示：在线路m上，列车i到达车站s的时间；/>表示：在线路m上，列车i离开车站s-1的时间；/>表示：在线路m上，列车i从车站s-1到车站s的运行时间；/>表示：在线路m上，列车i离开车站s的时间；/>表示：在线路m上，列车i在车站s的停站时间；
为满足运营需求，同一车站内两列相邻列车间的发车间隔时间需要满足如式(3)所述约束，
其中，和/>分别是线路m上车站s所规定的两列相邻列车间最小发车间隔时间和最大发车间隔时间；
为保证列车i从离开起始站到到达终点站的时间在合理的总旅行时间内，列车i的总旅行时间需要满足如式(4)所述约束，
其中，是线路m上列车i到达终点站的时间；/>是线路m上列车i离开起始站的时间；/>和/>分别是线路m上规定的最小旅行时间和最大旅行时间；
考虑到运营安全和服务水平，线路m上列车i在车站s的停站时间需要满足如式(5)所述约束，从车站s到车站s+1的运行时间/>需要满足如式(6)所述约束，
其中：和/>分别是线路m上车站s规定的最小和最大停站时间；/>和/>分别是线路m上车站s规定的最小和最大区间运行时间；
为保证所有列车在计划周期H内完成行程，是线路m上最后一列车到达终点站的时间，需满足如式(7)所述约束：
S422、计算换乘等待时间线性约束：
线路m和线路m′交汇于s站，若在换乘站s实现线路m至线路m′的换乘，则需要满足乘客乘坐线路m′上列车i′到达后，换乘走行至线路m的换乘站台时，列车还未驶离换乘站s；
为准确表示列车在换乘站的接续，引入0-1变量满足如式(8)所述约束，
其中，M是足够大的正数，是在换乘站内两条线路站台间的走行时间；T是目标大型地铁网络中所有换乘站的集合，例如：s是线路m′到m的换乘站；注意，如果/>那么/>意味着线路m′上列车i′到达换乘站s足够早，而线路m上接续列车i离开足够晚，则列车i从时间上来看能够接续成功；如果/>那么/>这时情况则相反，/>是乘客从线路m′上的列车i′到达换乘站s后最早能够换乘成功的时间，而/>是线路m上接续列车i的发车时间；
基于对的定义，引出/>满足如式(9)所述约束，
其中，基于所有乘客都将登上第一列接续列车的假设，表示两条相交线路在换乘站实际的列车衔接情况；如果线路m上的列车i，而不是列车i-1，是线路m′上列车i′的第一列接续列车，则得到/>因此，由式(9)分别推出
根据上述约束，如果线路m上的列车i是线路m′上列车i′的第一列接续列车，则和/>满足式(10)，
换乘站s内在两列车间换乘乘客的等待时间计算如约束式(11)所示，如果线路m′上的列车i′成功接续线路m上列车i，则换乘等待时间为列车i′到其他列车的换乘等待时间则为0，
S423、最小化乘客等待时间的目标函数：
为了最小化乘客的等待时间，所述等待时间包括：换乘乘客和从地铁网络外进入车站乘客的等待时间，构建目标函数如式(12)所示，
其中，是相交线路m′和m上列车i′和i在换乘站s的换乘需求，/>是进入车站s等待线路m上列车i的平均乘客到达率；目标函数中的第一项表示所有换乘乘客的等待时间，而第二项表示从地铁网络外进入车站乘客的等待时间；因最大程度减少换乘等待时间会增加列车的衔接数量，延长停站时间和发车间隔时间，从而影响不需要换乘乘客的出行体验；预先给定的权重ρ1和ρ2可以平衡换乘乘客和进入车站乘客的等待时间；
S424、建立时刻表协同优化模型：
结合上述目标函数和约束，构建时刻表协同优化模型如式(13)所示，
上述时刻表协同优化模型是一个混合整数非线性规划问题，称为：时刻表编制问题，其约束是线性的，而目标函数是非线性的。
在上述技术方案的基础上，包含大量0-1变量的混合整数非线性规划问题较复杂，难以在较短时间内求解包含大量线路、车站的大规模地铁网络问题。为减轻大规模问题的计算负担，提出：以广义Benders分解算法来代替传统的集中式求解方法。
基于广义Benders分解算法将步骤S42构建的时刻表编制问题分为两部分：子问题(SP)和主问题(MP)；将一些复杂变量特别是整数变量暂时固定为给定值，剩余部分称为子问题较容易解决；如果子问题不可行，则需要引入松弛变量，使得子问题可行；然后将子问题求解出的其他决策变量值和对偶变量值构建Benders割，而迭代增加的Benders割和时刻表编制问题中仅涉及复杂变量的约束构成主问题，优化求解复杂变量；最后，将主问题求出的复杂变量值代入子问题再次求解，重复以上过程直至收敛至最优解；
所述Benders割包括：可行割和最优割；
Benders割函数包括：可行割函数和最优割函数。
在上述技术方案的基础上，所述步骤S43具体包括以下步骤：
S431、推导与求解子问题：
根据广义Benders分解算法，在时刻表协同优化模型，即式(13)中，定义两列车是否可以衔接的0-1决策变量作为被赋值的复杂变量，求解子问题得到其他决策变量和/>
首先，将全部设为0或者根据原始时刻表的列车衔接设置初始值/>但固定值可能会导致子问题不可行，故将人工变量/>和/>引入式(8)所示的耦合约束以及目标罚函数来使其可行，而不改变最优解。构建子问题如式(14)所示，
其中，如果对子问题可行，则/>固定/>后，我们得到的混合整数非线性规划子问题远比原始的混合整数非线性规划问题容易处理，通过运用Karush-Kuhn-Tucker(KKT)条件得到对偶变量值，式(14)中加粗的/>和分别是对应约束的对偶变量，将对偶变量作为对偶乘子建立拉格朗日对偶问题；对偶变量的定义域取决于其相关约束的性质；比如/>R代表实数。
然后，结合子问题求出的最优值和/>等，求解下列方程组(15)得到上述对偶变量值，式(14)运用KKT条件后得到方程组(15)，
因整数变量固定后，子问题是原问题的限制问题，故其目标函数的最优值是原问题的上限。最后，将第k次迭代的子问题的目标函数值写作Zk，其中，k＝1，2，...，K；K是最大迭代次数，构建上限UBk如式(16)所示，
UBk＝Zk，k＝1，...，K (16)；
S432、构建Benders割函数和主问题：
一些整数变量的给定值可能对耦合约束(8)不可行，从而使得子问题不可行，因此需基于实际情况，向主问题加入防止问题不可行且不改变最优解的可行割。给定以及子问题中求得的对偶乘子，构建Benders割函数如下，其中可行割函数如式(17)和最优割函数如式(18)，
将所有的Benders割一起构成主问题，主问题的目的在于：通过求解最小化μ(即使得μ最小)；在分解算法的k次迭代中，如果子问题不可行(即式(14)中的人工变量和/>不等于0)，将新的可行割/>加进主问题中；如果子问题能求得最优解(即式(14)中的人工变量/>和/>等于0)，将新的最优割/>加进主问题中；可行割和最优割的数量随着每次迭代而增加，因此，构建包括Kf个可行割和Kp个最优割的主问题如式(19)所示，
其中，Kf+Kp＝K，且μ是标量，其界限考虑实际情况给出，这里没有设置其下限。
应注意，主问题是原问题的松弛问题，它从下接近原问题的目标函数。因此，对第k次迭代，式(19)的目标函数值是原问题目标函数的下界，如式(20)所示，
LBk＝μk，k＝1，...，K (20)
S433、基于广义Benders分解算法求解，具体步骤如下：
第一步：初始化，设置迭代计数k＝1和容差值ε，将整数变量赋初值/>
第二步：求解子问题：
如果固定后，子问题可行，那么求解子问题得到目标函数值Zk、决策变量和对偶变量/>否则，如果固定/>后，子问题不可行，则通过求解松弛子问题得到子问题的解，更新UBk；
第三步：生成Benders割：利用子问题的解，如果子问题不可行，按公式(17)生成可行割将新的可行割/>加进主问题中；如果子问题能求得最优解，按公式(18)生成最优割/>将新的最优割/>加进主问题中；
第四步：求解主问题；
求解式(19)，得到μk和整数变量的新解更新LBk；
第五步：检查收敛：
计算上界UBk和下界LBk，如果UBk-LBk≤ε，则算法终止；否则迭代计数k←k+1，继续第二-五步。
在上述技术方案的基础上，在所述广义Benders分解算法中，当一些复杂决策变量固定时，有下列三种情况：
①子问题可行且有界，只有一个最优解，将最优割加至主问题；
②子问题可行但无界，因原问题也无界，故算法停止(Geoffrion，1972)；
③子问题不可行，将可行割加至主问题。
在上述技术方案的基础上，如果子问题不可行是指：式(14)中的人工变量和不等于0；如果子问题能求得最优解是指：式(14)中的人工变量/>和/>等于0。
本发明的有益技术效果如下：
本发明通过历史或预测的客流信息，针对具有较多线路和车站的大规模地铁网络，提出了一种地铁网络列车时刻表快速编制方法，通过调整列车的到发时间来优化列车间衔接情况，采用广义Benders分解算法将原问题分解为更易求解的子问题和主问题，具有较好的求解效率，且能够最大程度减少乘客等待时间，提高地铁系统的换乘效率和服务水平。
附图说明
本发明有如下附图：
图1示出本发明提供的城市地铁线路示意图。
具体实施方式
为了更清楚地说明本发明，下面结合优选实例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解，下面所具体描述的内容是说明性的，而非限制性的，不应以此限制本发明的保护范围。
为证明所提出的模型和求解方法的实用性，我们对实际大型地铁网络进行试验。图1中的网络是北京地铁网络的一部分，包括15条运营线和51个换乘站，其运行方向在图中用箭头标记。
试验考虑的规划周期是早上8点至10点，共150列车从线路首站发出，每条线10列车。每条线的发车间隔时间的上下限列在了表1中，每个车站的列车停站时间是[10，60]s。相邻车站间的列车运行时间是[120，300]s。另外，根据北京地铁网络实际运营数据设置相交线路m和m′在换乘站s的换乘走行时间根据北京地铁网络的工作日实际客流记录，设置乘客需求。而目标函数中的权重系数设为ρ1＝1，ρ2＝0.2。
表1试验的北京地铁网络中每条线的发车间隔的上下限列表
首先，我们尝试用CPLEX求解器直接求解，但没有得到一个可接受的解。而运用Benders分解算法的情况下，目标函数值为4.82×106s，换乘等待时间为2.13×106s，进站乘客等待时间为2.69×106s。CPU的运行时间是395s，远比CPLEX求解器的计算时间短。计算得到对比表2，利用所提出模型优化时刻表前后的结果。目标函数值，换乘等待时间和进站乘客等待时间均得到显著改善，最明显的是与优化前相比，乘客换乘等待时间降低了20.52％。通过轻微调整列车的到发时间和衔接情况，乘客总等待时间明显减少。
表2列车时刻表优化前后的结果对比表
目标函数(s) 换乘等待时间(s) 进站乘客等待时间(s)
优化前 5.67×106 2.68×106 2.99×106
优化后 4.82×106 2.13×106 2.69×106
提升 14.99％ 20.52％ 10.03％
特别地，表2展示了地铁网络中所有换乘车站的换乘等待时间和相应的改进之处。计算优化前后每个换乘车站上10列车的总换乘等待时间的偏差作为每个换乘节点的改进之处。譬如，从1号线至2号线的建国门站换乘节点在优化前后的换乘等待时间分别为23650s和10240s，故其提升了13410s。而少量换乘节点上换乘等待时间的增加，可以减少整个地铁网络上总的乘客换乘等待时间。因为每条线上每列车在首站发车时间的限制，共有8个换乘节点在规划周期内不能接续成功，一些线路上在前几列车里的乘客可能会等待较长时间。
综上，本发明公开的基于广义Benders分解算法的城市轨道交通列车时刻表快速编制方法，构建了一种考虑列车衔接的时刻表协同优化模型，并运用Benders分解算法将原模型分解为容易求解的子问题和主问题两部分，从而能够快速求解得到协调优化后的列车时刻表，最大程度减少换乘乘客和进站乘客的等待时间，大大提高运营服务水平和乘客出行满意度。
显然，本发明的上述实施例仅仅是为清楚地说明本发明所做的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引申出的显而易见的变化或变动仍处于本发明的保护范围之列。
本说明书中未做详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (5)

1.一种地铁网络列车时刻表快速编制方法，其特征在于，包括以下步骤：
S1：设置目标地铁网络的基础参数；
所述基础参数包括：线路数量、每条线路的车站数和目标地铁网络中不同线路间换乘车站情况；
S2：设置各线路区间的运行时间、停站时间、发车间隔和旅行时间的范围；
采集现实数据，设置各换乘车站站台间的换乘走行时间；
S3：采集地铁网络的客流信息；
所述地铁网络的客流信息包括：换乘车站的换乘乘客人数和进站的乘客人数；
S4：采用广义Benders分解算法进行编制，具体包括以下步骤：
S41：给定计划列车数量和时域；
S42：根据目标大型地铁网络相关参数和客流信息，建立时刻表协同优化模型；
S43：利用广义Benders分解算法，将步骤S42构建的时刻表协同优化模型，分解为子问题和主问题，并迭代求解子问题和主问题，直至收敛至最优解；
S5：由步骤S4得到最优解，获得大型地铁网络的协同优化列车时刻表和线路间列车接续情况；
步骤S42的具体步骤如下：
S421、构建时刻表线性约束：
利用公式(1)得到各线路上每列车到达各车站的时间，利用公式(2)得到各线路上每列车离开各车站的时间，
其中：i是列车标号；s是车站标号；m是线路标号；M表示：目标大型地铁网络中所有线路标号的集合；Sm表示：目标大型地铁网络中线路m中所有车站标号的集合；Nm表示：目标大型地铁网络中线路m中所有列车标号的集合；表示：在线路m上，列车i到达车站s的时间；表示：在线路m上，列车i离开车站s-1的时间；/>表示：在线路m上，列车i从车站s-1到车站s的运行时间；/>表示：在线路m上，列车i离开车站s的时间；/>表示：在线路m上，列车i在车站s的停站时间；
为满足运营需求，同一车站内两列相邻列车间的发车间隔时间需要满足如式(3)约束，
其中，和/>分别是线路m上车站s所规定的两列相邻列车间最小发车间隔时间和最大发车间隔时间；
为保证列车i从离开起始站到到达终点站的时间在合理的总旅行时间内，列车i的总旅行时间需要满足如式(4)约束，
其中，是线路m上列车i到达终点站的时间；/>是线路m上列车i离开起始站的时间；/>和/>分别是线路m上规定的最小旅行时间和最大旅行时间；
考虑到运营安全和服务水平，线路m上列车i在车站s的停站时间需要满足如式(5)约束，从车站s到车站s+1的运行时间/>需要满足如式(6)约束，
其中：和/>分别是线路m上车站s规定的最小和最大停站时间；/>和/>分别是线路m上车站s规定的最小和最大区间运行时间；
为保证所有列车在计划周期H内完成行程，是线路m上最后一列车到达终点站的时间，需满足如式(7)约束：
S422、计算换乘等待时间线性约束：
线路m和线路m'交汇于s站，若在换乘站s实现线路m至线路m'的换乘，则需要满足乘客乘坐线路m'上列车i'到达后，换乘走行至线路m的换乘站台时，列车还未驶离换乘站s；
为准确表示列车在换乘站的接续，引入0-1变量满足如式(8)约束，
其中，M是足够大的正数，是在换乘站内两条线路站台间的走行时间；T是目标大型地铁网络中所有换乘站的集合；如果/>那么/>意味着线路m'上列车i'到达换乘站s早，而线路m上接续列车i离开晚，则列车i从时间上来看能够接续成功；如果/>那么/>这时情况则相反，/>是乘客从线路m'上的列车i'到达换乘站s后最早能够换乘成功的时间，而/>是线路m上接续列车i的发车时间；
基于对的定义，引出/>满足如式(9)约束，
其中，基于所有乘客都将登上第一列接续列车的假设，表示两条相交线路在换乘站实际的列车衔接情况；如果线路m上的列车i，而不是列车i-1，是线路m'上列车i'的第一列接续列车，则得到/>因此，由式(9)分别推出
根据上述式(9)约束，如果线路m上的列车i是线路m'上列车i'的第一列接续列车，则和/>满足式(10)，
换乘站s内在两列车间换乘乘客的等待时间计算如约束式(11)所示，如果线路m'上的列车i'成功接续线路m上列车i，则换乘等待时间为列车i'到其他列车的换乘等待时间则为0，
S423、最小化乘客等待时间的目标函数：
为了最小化乘客的等待时间，所述等待时间包括：换乘乘客和从地铁网络外进入车站乘客的等待时间，构建目标函数如式(12)所示，
其中，是相交线路m'和m上列车i'和i在换乘站s的换乘需求，/>是进入车站s等待线路m上列车i的平均乘客到达率；目标函数中的第一项表示所有换乘乘客的等待时间，而第二项表示从地铁网络外进入车站乘客的等待时间；因最大程度减少换乘等待时间会增加列车的衔接数量，延长停站时间和发车间隔时间，从而影响不需要换乘乘客的出行体验；预先给定的权重ρ1和ρ2平衡换乘乘客和进入车站乘客的等待时间；
S424、建立时刻表协同优化模型：
结合上述目标函数和约束，构建时刻表协同优化模型如式(13)所示，
上述时刻表协同优化模型是一个混合整数非线性规划问题，称为：时刻表编制问题，其约束是线性的，而目标函数是非线性的。
2.如权利要求1所述的地铁网络列车时刻表快速编制方法，其特征在于：基于广义Benders分解算法将步骤S42构建的时刻表编制问题分为两部分：子问题和主问题；将整数变量暂时固定为给定值，剩余部分称为子问题；如果子问题不可行，则引入松弛变量，使得子问题可行；然后将子问题求解出的其他决策变量值和对偶变量值构建Benders割，而迭代增加的Benders割和时刻表编制问题中仅涉及复杂变量的约束构成主问题，优化求解复杂变量；最后，将主问题求出的复杂变量值代入子问题再次求解，重复以上过程直至收敛至最优解；
所述Benders割包括：可行割和最优割；
Benders割函数包括：可行割函数和最优割函数。
3.如权利要求2所述的地铁网络列车时刻表快速编制方法，其特征在于：所述步骤S43具体包括以下步骤：
S431、推导与求解子问题：
根据广义Benders分解算法，在时刻表协同优化模型中，定义两列车是否衔接的0-1决策变量作为被赋值的复杂变量，求解子问题得到其他决策变量/>和
首先，将全部设为0或者根据原始时刻表的列车衔接设置初始值/>将人工变量/>和/>引入式(8)所示的耦合约束以及目标罚函数，构建子问题如式(14)所示，
其中，如果对子问题可行，则/>通过运用Karush-Kuhn-Tucker条件得到对偶变量值，式(14)中加粗的/> 和/>分别是对应约束的对偶变量，将对偶变量作为对偶乘子建立拉格朗日对偶问题；对偶变量的定义域取决于其相关约束的性质；
然后，结合子问题求出的最优值和/>，求解下列方程组(15)得到上述对偶变量值，
最后，将第k次迭代的子问题的目标函数值写作Zk，其中，k＝1,2,…,K；K是最大迭代次数，构建上限UBk如式(16)所示，
UBk＝Zk,k＝1,...,K (16)；
S432、构建Benders割函数和主问题：
给定以及子问题中求得的对偶乘子，构建Benders割函数如下，其中可行割函数如式(17)和最优割函数如式(18)，
将所有的Benders割一起构成主问题，主问题的目的在于：通过求解最小化μ；在k次迭代中，如果子问题不可行，将新的可行割/>加进主问题中；如果子问题能求得最优解，将新的最优割/>加进主问题中；可行割和最优割的数量随着每次迭代而增加，因此，构建包括Kf个可行割和Kp个最优割的主问题如式(19)所示，
其中，Kf+Kp＝K，且μ是标量；
对第k次迭代，式(19)的目标函数值是原问题目标函数的下界，如式(20)所示，
LBk＝μk,k＝1,...,K (20)
S433、基于广义Benders分解算法求解，具体步骤如下：
第一步：初始化，设置迭代计数k＝1和容差值ε，将整数变量赋初值/>
第二步：求解子问题：
如果固定后，子问题可行，那么求解子问题得到目标函数值Zk、决策变量和对偶变量/>否则，如果固定/>后，子问题不可行，则通过求解松弛子问题得到子问题的解，更新UBk；
第三步：生成Benders割：利用子问题的解，如果子问题不可行，按公式(17)生成可行割将新的可行割/>加进主问题中；如果子问题能求得最优解，按公式(18)生成最优割/>将新的最优割/>加进主问题中；
第四步：求解主问题；
求解式(19)，得到μk和整数变量的新解更新LBk；
第五步：检查收敛：
计算上界UBk和下界LBk，如果UBk-LBk≤ε，则算法终止；否则迭代计数k←k+1，继续第二-五步。
4.如权利要求3所述的地铁网络列车时刻表快速编制方法，其特征在于：在所述广义Benders分解算法中，当复杂决策变量固定时，有下列三种情况：
①子问题可行且有界，只有一个最优解，将最优割加至主问题；
②子问题可行但无界，因原问题也无界，故算法停止；
③子问题不可行，将可行割加至主问题。
5.如权利要求4所述的地铁网络列车时刻表快速编制方法，其特征在于：如果子问题不可行是指：式(14)中的人工变量和/>不等于0；如果子问题能求得最优解是指：式(14)中的人工变量/>和/>等于0。
CN202110370870.8A 2021-04-07 2021-04-07 一种地铁网络列车时刻表快速编制方法 Active CN113177295B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110370870.8A CN113177295B (zh) 2021-04-07 2021-04-07 一种地铁网络列车时刻表快速编制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110370870.8A CN113177295B (zh) 2021-04-07 2021-04-07 一种地铁网络列车时刻表快速编制方法

Publications (2)

Publication Number Publication Date
CN113177295A CN113177295A (zh) 2021-07-27
CN113177295B true CN113177295B (zh) 2024-02-27

Family

ID=76923158

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110370870.8A Active CN113177295B (zh) 2021-04-07 2021-04-07 一种地铁网络列车时刻表快速编制方法

Country Status (1)

Country Link
CN (1) CN113177295B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115619082B (zh) * 2022-12-02 2023-03-28 天津大学 一种人机协同混流装配线平衡问题的解决方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101853440A (zh) * 2010-05-11 2010-10-06 北京交通大学 列车运行计划自动编制的实现方法及其系统
CN102169524A (zh) * 2010-02-26 2011-08-31 同济大学 城市轨道交通网络客流分布的分阶段多路径模型算法
CN102737129A (zh) * 2012-06-20 2012-10-17 东南大学 一种公交车辆区域行车时刻表编制的方法
WO2017045294A1 (zh) * 2015-09-17 2017-03-23 华南理工大学 一种城市常规公交线网设计方法
CN111967134A (zh) * 2020-07-10 2020-11-20 北京交通大学 基于地铁共线共车的地下物流系统优化控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102169524A (zh) * 2010-02-26 2011-08-31 同济大学 城市轨道交通网络客流分布的分阶段多路径模型算法
CN101853440A (zh) * 2010-05-11 2010-10-06 北京交通大学 列车运行计划自动编制的实现方法及其系统
CN102737129A (zh) * 2012-06-20 2012-10-17 东南大学 一种公交车辆区域行车时刻表编制的方法
WO2017045294A1 (zh) * 2015-09-17 2017-03-23 华南理工大学 一种城市常规公交线网设计方法
CN111967134A (zh) * 2020-07-10 2020-11-20 北京交通大学 基于地铁共线共车的地下物流系统优化控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《基于机器学习的重载列车智能驾驶方法研究》;王悉;《中国博士学位论文全文数据库（工程科技Ⅱ辑）》;全文 *
《铁路列车时刻表与停站方案协同优化研究》;戚建国;《中国博士学位论文全文数据库（工程科技Ⅱ辑）》;全文 *
杨艳.《 地铁线网末班列车时刻表协同编制优化》.《山东交通学院学报》.2020,全文. *

Also Published As

Publication number Publication date
CN113177295A (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
Kang et al. Last train timetabling optimization and bus bridging service management in urban railway transit networks
Kang et al. A case study on the coordination of last trains for the Beijing subway network
Kuppusamy et al. Deep learning based energy efficient optimal timetable rescheduling model for intelligent metro transportation systems
CN110135755B (zh) 一种综合优化片区城乡公交时刻表编制与车辆调度的方法
Yang et al. Analysis of energy consumption reduction in metro systems using rolling stop-skipping patterns
CN109508751B (zh) 高速铁路列车晚点时间预测的深度神经网络模型建模方法
CN105678411A (zh) 一种旅客列车运行方案图编制方法
CN110533219B (zh) 城市轨道交通末班列车时刻表优化方法
CN114312926B (zh) 一种城市轨道交通列车运行调整方案优化方法和系统
CN113177295B (zh) 一种地铁网络列车时刻表快速编制方法
Wang et al. Optimization models for high-speed train unit routing problems
CN111815166A (zh) 基于客流分配的城市公共交通系统的协同优化方法和装置
Zhang et al. Coupling analysis of passenger and train flows for a large-scale urban rail transit system
Zhou et al. Collaborative optimization of energy-efficient train schedule and train circulation plan for urban rail
Li et al. Joint optimization of delay-recovery and energy-saving in a metro system: A case study from China
CN111985814A (zh) 间歇式供电的城际列车开行方案的优化方法及系统
Zhang et al. Real-time energy saving optimization method for urban rail transit train timetable under delay condition
Li et al. Optimal track utilization in electric multiple unit maintenance depots
Zhao et al. Integrated optimization of demand-driven timetable, train formation plan and rolling stock circulation with variable running times and dwell times
Zhang et al. A Multistage Decision Optimization Approach for Train Timetable Rescheduling Under Uncertain Disruptions in a High-Speed Railway Network
Liu et al. Spatiotemporal synchronous coupling algorithm for urban rail transit timetables design under dynamic passenger demand
CN114537479A (zh) 考虑列车开行方案周内波动的动车组运用计划编制方法
CN114021291A (zh) 一种城市轨道交通线网限流方案仿真评价建模方法
CN113715875A (zh) 一种面向初始延误与区间限速的列车运行调控方法和系统
Lai et al. Collaborative optimization model for the design and operation of feeder bus routes based on urban metro

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