第十二届中关村青联杯全国研究生数学建模大赛面向节能的单多列车优化决策问题9

参赛密码（由组委会填写）第十二届“中关村青联杯”全国研究生数学建模竞赛学校海军装备研究院/海军工程大学参赛队号K0032队员姓名1.蒋靖伟2.曹进一3.齐振参赛密码（由组委会填写）-1-第十二届“中关村青联杯”全国研究生数学建模竞赛题目面向节能的单/多列车优化决策问题摘要：本文围绕单/多列车优化决策问题，在合理假设的基础上，利用多岛遗传优化算法和NSGA-Ⅱ多目标优化算法给出了单列车单站点、单列车多站点、多列车多站点的能耗最低运行线路的优化决策，并分析处理了列车发生延误时的优化控制问题。针对问题一（1），建立了单列车单区间节能优化模型。首先通过将时间分段-离散的方法，建立了能耗积分方程的数值求解方法，并制定了末端制动策略使得末端速度在规定时间、规定距离上减小为0。在此基础上，建立了以能耗最低为优化目标，分段数、各分段时间间隔、各段运行工况为决策变量，满足速度、加速度等约束条件的优化模型。通过多岛遗传算法，对模型进行求解，得到A6-A7段能耗为3.37×107J。针对问题一（2），建立了单列车多区间节能优化模型。首先通过理论推导，将时间-最低能耗曲线转换为以最少时间、最低能耗为双目标优化问题的Pareto前端解集，利用NSGA-Ⅱ多目标优化算法分别得到了A6-A7站，A7-A8站Pareto前端解集。其次，在各自能耗-时间Pareto前端解集中，利用多岛遗传算法，对时间分配进行优化建模，得到A6-A7段运行时间117s，A7-A8段运行时间103s，总能耗为6.8×107J。针对问题二（1），建立了多列车全区间节能优化模型，在总能耗一定的情况下，再生能源越多，则总能量越少。基于此，本文首先求解单个列车在整个区间段上的最少能耗，这是对于问题一（2）的推广，区别仅在于将停站时间计入运行时间，没有本质上的区别，本文采用将停站看作除去牵引、巡航、惰行和制动在外的第5种工况，采用与问题一（2）相同的策略，求得单列车在整个运行区间（A1-A14）上的最低能耗，其它车辆采用相同的运行方式。其次，建立了针对发车时间的可再生能源最大模型，为了简-2-化模型，将车辆分成20组，认为各组之间运行方式相同。对于同一组内，以再生能源最多为目标函数建立优化模型，通过多岛遗传算法进行求解各车发车间隔。最终得到：通过合理分配发车间隔，可再生能源可占总耗能约17.5%。针对问题二（2），在问题二（1）的基础上，建立了多车全区间全时间段节能优化模型。将整个时间段分为2个高峰阶段和3个非高峰阶段，其中早高峰分配辆车，晚高峰分配辆车，非高峰阶段车辆平均分配。对于5个阶段，每个段内再按照题二（1）中原则对车辆进行分组，各组内采用和题二（1）相同的模型求解发车间隔，最终早高峰可节能37%左右，晚高峰节能27%左右。针对问题三（1），建立了延误节能追赶模型，在问题二（1）中求出的各段最低能耗-最短时间Pareto解集基础上，通过最短追赶时间-最低能耗多目标优化模型，求出最节能的追赶方案。通过两个算例分析发现，针对减少追赶时间和减少能耗这一对矛盾的目标，不同站点之间应采用不同的策略。其中对于A4段10s的延误，车辆应尽快恢复正常运行。针对问题三（2），建立了概率延误节能追赶模型，通过后车受前车延误影响将延误概率参数纳入题三（1）建立的模型中，并从理论上进行了分析。此外，本文还以问题一（1）为例对设计参数进行了敏感性分析，得到分段时间间隔和分段时刻对结果有最为明显的正效应影响。本文下一步将在丰富最低能耗-最短时间Pareto解集、引入延误概率参数后的方案控制方面进行研究。关键词：列车节能，多岛遗传，NSGA-Ⅱ多目标优化，Pareto前端解-3-一问题重述轨道交通系统的能耗是指列车牵引、通风空调、电梯、照明、给排水、弱电等设备产生的能耗。在低碳环保、节能减排日益受到关注的情况下，针对减少列车牵引能耗的列车运行优化控制近年来成为轨道交通领域的重要研究方向。本题给出列车运行过程、列车动力学模型、运行时间和运行能耗的关系以及再生能量利用原理。并给定列车参数和线路参数。根据已知内容，需要解决的问题如下：问题一：单列车节能运行优化控制问题1、计算寻找一条列车从6A站出发到达7A站的最节能运行的速度距离曲线，其中两车站间的运行时间为110秒。2、计算寻找一条列车从6A站出发到达8A站的最节能运行的速度距离曲线，其中要求列车在7A车站停站45秒，6A站和8A站间总运行时间规定为220秒（不包括停站时间）。问题二：多列车节能运行优化控制问题1、当100列列车以间隔199,...,Hhh从A1站出发，建立优化模型并寻找使所有列车运行总能耗最低的间隔H。2、重新制定运行图和相应的速度距离曲线，考虑高峰时间（早高峰7200秒至12600秒，晚高峰43200至50400秒）发车间隔不大于2.5分钟且不小于2分钟，其余时间发车间隔不小于5分钟，每天240列。问题三：列车延误后运行优化控制问题1、若列车i在车站jA延误ijDT（10秒）发车，找出在确保安全的前提下，首先使所有后续列车尽快恢复正点运行，其次恢复期间耗能最少的列车运行曲线。2、随机变量延误ijDT发车，在尽快恢复正点运行，恢复期间耗能最少的目标函数下，给出列车运行曲线。-4-二符号说明符号符号说明()cn列车运行状态函数N列车站间分段数M站间使用工况的类型数量i列车进入第N段后运行的时间it列车运行时刻t离散时间步长itQait时刻列车牵引加速度itait时刻列车实际加速度itS计算距离，是列车到刚通过的一站的距离()iet列车在第N段中it时刻的能耗1itit时刻牵引加速度与最大加速度百分比2itit时刻制动加速度与最大加速度百分比T第N段列车总运行时间itvit时刻列车运行速度TS站间总距离-5-三单列车节能运行优化控制问题3.1问题分析问题一（1）要求我们建立速度距离曲线的数学模型，制定列车在6A站到7A站运行110s耗能最少的方案。列车发动机耗能与运行工况密切相关，在四种运行工况（牵引、巡航、惰行和制动）中，牵引阶段发动机耗能，巡航阶段发动机是否耗能取决于列车当时受到的总阻力。总阻力大于0时，列车需要牵引，发动机耗能；总阻力小于或等于0时，列车需要制动，发动机不耗能。单质点模型中，列车运动符合牛顿运动学定律，根据题目所给“列车参数”和“线路参数”，可以得到列车牵引力，列车运行总阻力和列车制动力等参数。分析不同阶段列车的受力情况，建立列车动力学模型，得到列车在不同工况下的发动机能耗。节能运行的关键在于列车在行驶过程中工况的交替使用，单列车节能运行优化控制问题的本质是一个单目标优化问题。列车运行的总能耗最小为目标函数，运行工况、两站之间列车运行工况的阶段数等参数为决策变量，列车的启止速度、不同路段的限速、最大加减速度等为约束条件。采用多岛遗传算法MIGA作为优化策略，对发动机的总能耗结果进行全局寻优，确定能耗最小条件下，列车的运行策略，得到最节能运行的速度距离曲线。由于列车的运行策略受到不同路段限速和坡度等参数的影响，很难直接得到连续的速度-距离曲线公式。在算法的搜索过程中，为了快速的寻找最优解，对问题进行离散化处理，得到的数值解能够有效地解决实际问题。问题一（2）要求建模计算出列车从6A站出发到达8A站的最节能运行的速度距离曲线，其中列车在7A车站停站45秒，6A站和8A站间总运行时间规定为220秒（不包括停站时间）。相比问题一（1），问题一（2）只是增加了一段站间路程，并在途中7A站进行了停留。问题一（2）是在问题一（1）模型的基础上，将目标函数替换为两段路程能耗的和最小，约束条件为总运行时间一定。根据题干中图5列车站间运行时间与能耗的关系曲线，我们确定列车运行时间与最低能耗存在类似的关系。由于在运行时间与最低能耗的关系曲线下方不存在可行解，可以认为其实质上就是两个目标的Pareto前端解集[1]，于是我们将求解列车运行时间与能耗关系曲线的问题转化为两个目标的优化问题。我们采用改进型非劣分类遗传算法（NSGA-Ⅱ）求解多目标优化问题，得到站间各自的Pareto最优前端解集。利用该解集，分别对站间各自的运行时间进行遍历，保证运行总时间一定，使得总能耗最低的时间分配方案及站间运行策略，就是问题一（2）的解。-6-3.2模型假设1、列车在运行过程中不制动减速，只在进站停车时进行制动。2、列车采用制动工况减速进站停车，不采用惰行工况进站停车。3.3单列车单区间节能优化模型3.3.1模型的建立根据题目所给的《附录：路线参数》，可以确定如图3.1所示的列车参考坐标系：计算公里标A1A14AiAi+1起始公里标0公里标22903始发站终点站公里标175公里标第i站第i+1站计算距离运行方向图3.1列车参考坐标系图中，1A站为始发站，14A站为终点站，列车由始发站1A向终点站14A运行，1A位于公里标22903m处，14A位于公里标175m处，起始公里标0位于终点站右侧。其中计算公里标（m）是到起点的距离，计算距离（m）是到刚通过的一站的距离。根据公里标得到A6站到A7站的距离是1354m。在两车站间运行时间一定的条件下，计算寻找列车从6A站出发到达7A站的最节能运行的速度距离曲线，问题的本质是制定一种列车在约束条件下的运行策略，使得发动机的总能耗最低。建立的单列车单区间节能优化模型的如下：目标函数：0minmin(,())TEetctdt约束条件：max()()vtvt(0)()0vvT0()1354TTiSvtdt决策变量：()ct这里决策变量()ct是泛函，表示时刻与列车能耗有关的运行状态，T是列车达到下一站的运行时间。目标函数表示列车从6A站出发到达7A站的总能耗最低，是t时刻能耗在[0,]T时间段内的积分，其中t时刻能耗是时间t和决策变量的函数。约束条件为：1、列车的运行速度小于限制速度；2、列车在每站间的启止速度为0；3、列车从出发至停止进站时行进的距离TS为6A站到7A站的距离。-7-这里由于很难求得耗能积分的解析解，为了得到数值解，我们对时间进行了离散。由于列车在不同时刻的工况不同，首先将列车的运行时间进行分段，保证每段时间内列车处于相同的工况。然后对每一段时间以步长为1秒进行离散，通过对每秒钟的能耗进行求和，得到发动机总能耗的数值解，具体的分段离散过程在模型的求解中会有详细叙述。为了满足在规定时间点、规定距离上使得末端速度减小为0，需要对末端制动方案进行设计。因此，模型求解将从目标函数的离散求解、优化模型求解和末端制动三个方面进行阐述。3.3.2模型的求解3.3.2.1目标函数的离散化处理1、全区间上耗能积分数值求解由于目标函数连续积分求解析解比较困难，而列车在运行过程中通常会采用牵引到接近限制速度后，交替使用惰行、巡航、牵引三种工况，直至接近下一车站采用制动进站停车。我们首先对积分的时间区间进行分段，在时间轴上将区间分为N段，每一段里列车处于相同的运行工况，N就是列车从6A站运行至7A站过程中使用工况的总数量。分段后决策变量()ct就可以看做是n的函数，即()cn。表示列车在第n段的运行状态，这里的1,2,...,nN。规定第1段起始时刻从6A站出发，列车处于牵引工况；第N段截止时刻到达6A站，列车处于制动工况；中间的2N段，列车交替使用惰行、巡航、牵引三种工况。nt是列车在第n段的初始时间，1nnnttt是列车在第n段运行的时间，则nTNt。然后对每一段的时间进行离散，得到能耗离散后的表达式如下：11011(,())(())(())nnnntNNTtitnnitEetctdtetdtett（3-1）其中，11(())nnNttnetdt是将总运行时间分成N段，1()nnttetdt是列车在第N段的能耗；11(())nntNinitett是对第N段的能耗以1秒为间隔进行离散，()iet是列车在第N段中it时-8-刻的能耗。1t是离散后的时间步长，inttit，i是列车进入第N段后运行的时间。这样，我们就将连