CN109858154A - 一种基于多目标优化的列车节能运行方法 - Google Patents
一种基于多目标优化的列车节能运行方法 Download PDFInfo
- Publication number
- CN109858154A CN109858154A CN201910099114.9A CN201910099114A CN109858154A CN 109858154 A CN109858154 A CN 109858154A CN 201910099114 A CN201910099114 A CN 201910099114A CN 109858154 A CN109858154 A CN 109858154A
- Authority
- CN
- China
- Prior art keywords
- train
- speed curve
- time
- energy
- individuals
- 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.)
- Pending
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005265 energy consumption Methods 0.000 claims abstract description 23
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 230000002068 genetic effect Effects 0.000 claims abstract description 17
- 239000000126 substance Substances 0.000 claims abstract description 7
- 230000001174 ascending effect Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000005634 sigma model Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Train Traffic Observation, Control, And Security (AREA)
Abstract
本发明公开了一种基于多目标优化的列车节能运行方法,先获取列车运行线路中各个运行区间的限速信息和坡度信息,然后根据列车的各种基本信息构建列车的单质点模型;再利用单质点模型分别计算列车在各个运行区间中的牵引能耗Et、制动能耗Er和惰行能耗Eb,以构建列车的多目标优化模型;再采用多目标多目标遗传算法对目标函数F(x)进行求解,得到各个运行区间的Pareto最优解集;最后根据列车到达当前运行区间的实际运行时间和时刻表规定运行时间的时间差值,从下一运行区间对应的Pareto最优解集中获取最优速度曲线。本发明提供的列车节能运行方法,不仅能使列车运行更为节能,而且还列车晚点或早到问题进行弥补,使列车运行更为智能方便和舒适。
Description
方法领域
本发明属于列车运行优化领域,主要涉及一种基于多目标优化的列车节能运行方法。
背景技术
城市轨道交通是我国重要的交通运输方式,由于其作业频繁、运输量大的特点,列车耗电量巨大,因此对城市轨道交通列车运行节能优化具有重要的意义。现有技术中对列车的节能优化问题已经有了一定的研究。目前列车节能运行方法一般在分析列车运行线路特征和列车限速的基础上,预先计算出列车在其运行线路上各个运行区间的节能控制序列或相应的列车速度曲线,在列车进行区间运行时,控制列车按照对应的列车速度曲线运行便可达到节能的效果。
然而,目前列车节能运行方法大部分都是将列车运行的优化和对列车的控制分开处理,当列车按照预先计算得到的速度曲线运行时,若列车实际运行时间与列车运行时刻表规定的运行时间之间产生误差,那么预先计算得到的各个行车区间之间的速度曲线便不再适用。例如,当列车运行出现延迟时,若继续按照预先计算出的速度曲线运行,会将延迟时间不断累加,导致后续运行区间中列车的运行时间与时刻表上的运行时间之间的时间差不断累积,致使列车运行不仅影响后续列车运行且耗费更多能量。
发明内容
为了解决上述问题,本发明提供了一种能够及时调整列车运行的推荐速度曲线的列车节能运行优化方法。
本发明是按照以下方法实现的:
一种基于多目标优化的列车节能运行优化方法,包括以下步骤:
S1、获取列车运行线路中各个运行区间的限速信息、坡度信息、线路曲线半径信息和线路隧道长度信息;
S2、根据列车的质量、所述限速信息和所述坡度信息、所述线路曲线半径信息和所述线路隧道长度信息;计算所述列车在所述各个运行区间的牵引力、制动力、基本阻力和附加阻力,以构建所述列车的单质点模型;
S3、利用所述单质点模型分别计算所述列车在各个运行区间中的牵引能耗Et、制动能耗Er和惰行能耗Eb,以构建所述列车的多目标优化模型,其中所述多目标优化模型如下:
目标函数为minY=F(x)=f(Ke,Kt,Kc)T,
约束条件为
其中,所述目标函数F(x)用于指示各个运行区间推荐速度曲线x的函数,Ke为列车的运行能耗,Kt为列车的运行时间,Kc为列车的运行舒适度,Ta为时刻表规定运行时间,T为实际运行时间,ta的值表示可以接受的时间误差,σ模型的宽度系数,ak为目标舒适度;
S4、采用多目标遗传算法对所述目标函数F(x)进行求解,得到各个运行区间的Pareto最优解集;其中,所述Pareto最优解集中的每个解用于指示一条推荐速度曲线;
S5、根据所述列车到达当前运行区间的实际运行时间和时刻表规定运行时间的时间差值,从所述下一运行区间对应的Pareto最优解集中获取最优速度曲线,以使所述列车在下一运行区间按照所述最优速度曲线运行;所述最优速度曲线用于表示时刻表规定运行时间为各个区间预设时刻表时间、预设停车时间和所述时间差值的和的推荐速度曲线。
进一步地,所述列车的单质点模型为:
其中,Fp表示控制每单位质量所受到的牵引力,fb表示每单位质量所受的制动力,f0(v)表示每单位质量所受基本阻力,fadd(x)表示列车受到的线路附加阻力,且fadd(x)=0,x1为列车在运行区间中的位置,v表示列车的运行速度,列车所经历的时间为t=t(x1)。
进一步地,所述运行优化方法还包括如下步骤:
利用PID控制结构跟踪所述当前运行区间的最优速度曲线,得到模拟速度曲线;所述模拟速度曲线用于模拟列车运行的所述最优速度曲线;
在列车完成所述当前运行区间的运行时,将所述模拟速度曲线与实际运行速度曲线进行比较,以计算所述最优速度曲线的节能率。
进一步地,在步骤S4中通过如下步骤对所述目标函数F(x)进行求解:
S41、利用多目标遗传算法的生成所述多目标遗传算法的初始种群;所述初始种群为用于指示满足所述多目标函数的约束条件的解的集合,所述初始种群中的每一个个体对应所述多目标优化函数的一个解,所述初始种群的大小为N;
S42、对所述初始种群中的所述个体进行进化选择,并使用快速非支配排序对所述个体进行分层,以使所述Pareto最优解集逼近全局Pareto最优前沿;
S43、计算同层间个体的拥挤距离,使得Pareto最优解集均匀分布在所述全局Pareto最优前沿上;
S44、采用精英保留策略在预设的第一迭代次数内产生子代种群,得到所述多目标函数的Pareto最优前沿的非支配解集,以将所述Pareto最优前沿的非支配解集设定为所述目标函数F(x)的解。
进一步地,在所述步骤S42中通过如下步骤对所述个体进行分层包括如下步骤:
S421、从所述初始种群群体中获取至少两个非支配解集;
S422、针对所述两个非支配解集中的第i个非支配解集,将所述第i个非支配解集设置为第i非支配层Fi,将所述第i非支配层中的所有个体的非支配序值iranki都赋值为i,并将所述第i非支配层中所有个体从所述初始种群中除去;其中,i为整数,且1≤i≤N。
进一步地,在所述步骤S43中所述拥挤距离计算方法包括如下步骤:
S431、将处于同一层的个体的拥挤距离L[i]d初始化为0;
S432、对处于同一层的个体按照第m个目标函数值进行升序排列,m为整数,且,m<N;
S433、将位于排序两端的个体的拥挤距离设定为无穷大;
S434、对于排序中间的个体,采用如下公式计算其拥挤距离:
其中,L[i+1]m为第i+1个个体在第m个目标函数处的取值,和分别表示第m个目标函数的最大值和最小值;
S435、选择拥挤距离在预设数值范围内的个体组进行保留,以使所述多目标遗传算法迭代后得到的初代种群中的个体均匀分布于所述全局Pareto最优前沿上;
进一步地,在所述步骤S44中通过如下步骤得到述多目标函数的Pareto最优前沿的非支配解集:
S441、将进化选择后的初始种群和所述子代种群合并,形成大小为2N的新种群Ri;
S442、对所述新种群Ri采用快速非支配排序分层算法,获取每一个个体的非支配序值irank,并按照个体所在的非支配层Fj计算其对应的拥挤距离;
S443、在预设的第二迭代次数内根据对种群Ri计算得到的每一个个体的非支配序值irank进行择优选择,从中选取N个作为新的父代种群Pi+1,得到所述多目标函数的Pareto最优前沿的非支配解集。
进一步地,通过如下步骤模拟所述最优速度曲线:
对列车的状态空间模型进行线性化,得到线性化后的列车模型状态空间;
所述列车的状态空间模型为:
其中,为列车的加速度,M为列车质量,a、b、c为戴维斯参数,v为列车速度,fg为运行区间的坡道附加阻力;
所述列车模型状态空间为:
其中,为在vo处的偏导,为v在vo处的偏导,为u在vo处的偏导,vo=0;
利用PID控制结构和所述列车模型状态空间跟踪所述最优速度曲线,得到所述模拟速度曲线;
所述PID控制结构的传递函数为Gc(s)=KPKi(s)+Ki(s)/s+Kd(s)s,其中,Kp为比例控制系数,Ki为积分控制系数,Kd为微分控制系数。
与现有技术相比,本发明采用以列车能耗、运行时间、和运行舒适度作为衡量因子,并建立了多目标优化模型,通过对多目标优化模型求解得到列车运行的推荐速度曲线,使得列车按照该推荐速度曲线运行就能实现运行过程中既节能又舒适。本发明还利用利用PID控制结构追踪最优速度曲线的实际运行情况,采取每个运行区间重新优化的方式,对运行区间选取最优速度曲线方式对列车晚点或早到问题进行调整,使列车的运行更为智能、方便和舒适。
附图说明
为了更清楚地说明本发明实施例技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中所述的一种基于多目标优化的列车节能运行优化方法的步骤示意图;
图2为多目标遗传算法原理图;
图3为PID控制***结构图;
图4为列车运行示意图;
图5为Pareto最优解集示意图;
图6原计划速度曲线;
图7列车晚点重挑的速度曲线;
图8列车早到重挑的速度曲线;
图9全线列车上行运行优化结果图;
图10全线列车下行运行优化结果图。
具体实施方式
为了充分地了解本发明的目的、特征和效果,以下将结合附图与具体实施方式对本发明的构思、具体步骤及产生的方法效果作进一步说明。
如图1所示,本发明公开了一种基于多目标优化的列车节能运行优化方法,包括以下步骤:
S1、获取列车运行线路中各个运行区间的限速信息、坡度信息、线路曲线半径信息和线路隧道长度信息;
S2、根据列车的质量、限速信息、坡度信息、线路曲线半径信息和线路隧道长度信息,计算列车在所述各个运行区间的牵引力、制动力、基本阻力和附加阻力,按照《列车牵引计算规程》里的公式参数以及计算方法对列车的运行过程建立了以下的所述列车的单质点模型:
其中,Fp表示控制每单位质量所受到的牵引力,fb表示每单位质量所受的制动力,f0(v)表示每单位质量所受基本阻力,fadd(x)表示列车受到的线路附加阻力,且fadd(x)=0,x1为列车在运行区间中的位置,v表示列车的运行速度,列车所经历的时间为t=t(x1);
S3、利用单质点模型分别计算列车在各个运行区间中的牵引能耗Et、制动能耗Er和惰行能耗Eb,以构建列车的多目标优化模型,其中多目标优化模型如下:
目标函数为minY=F(x)=f(Ke,Kt,Kc)T,
约束条件为
其中,目标函数F(x)用于指示各个运行区间推荐速度曲线x的函数,Ke为列车的运行能耗,Kt为列车的运行时间,Kc为列车的运行舒适度,Ta为时刻表规定运行时间,T为实际运行时间,ta的值表示可以接受的时间误差,σ模型的宽度系数,ak为目标舒适度。S4、如图2所示,采用多目标遗传算法对所述目标函数F(x)进行求解,得到各个运行区间的Pareto最优解集;其中,Pareto最优解集中的每个解用于指示一推荐速度曲线。
其中,对目标函数F(x)求解的步骤如下:
S41、利用多目标遗传算法的生成Pareto最优解集的初始种群;初始种群为用于指示满足多目标函数的约束条件的解的集合,初始种群中的每一个解表示一个个体,初始种群的大小为N;
S42、对初始种群中的个体进行进化选择,并使用快速非支配排序对个体进行分层,以使Pareto最优解集逼近全局Pareto最优前沿;
S43、计算同层间个体的拥挤距离,使得Pareto最优解集均匀分布在全局Pareto最优前沿上,如图5所示;
S44、采用精英保留策略在预设的第一迭代次数内产生子代种群,得到多目标函数的Pareto最优前沿的非支配解集,以将Pareto最优前沿的非支配解集设定为目标函数F(x)的解。
其中:
步骤S42中对个体进行分层的步骤如下:
S421、从初始种群群体中获取至少两个非支配解集;
S422、针对两个非支配解集中的第i个非支配解集,将第i个非支配解集设置为第i非支配层Fi,将第i非支配层中的所有个体的非支配序值iranki都赋值为i,并将第i非支配层中所有个体从初始种群中除去;其中,i为整数,且1≤i≤N。
步骤S42中拥挤距离计算方法包括如下步骤:
S431、将处于同一层的个体的拥挤距离初始化为0;
S432、对处于同一层的个体按照第m个目标函数值进行升序排列,m为整数,且,m<N;
S433、将位于排序两端的个体的拥挤距离设定为无穷大;
S434、对于排序中间的个体,采用如下公式计算其拥挤距离:
其中,L[i+1]m为第i+1个个体在第m个目标函数处的取值,和分别表示第m个目标函数的最大值和最小值;
S435、选择拥挤距离在预设数值范围内的个体组进行保留,以使初代种群中的个体均匀分布于全局Pareto最优前沿上。
步骤S44中得到多目标函数的Pareto最优前沿的非支配解集的步骤如下:
S441、将进化选择后的初始种群和子代种群合并,形成大小为2N的新种群Ri;
S442、对新种群Ri采用快速非支配排序分层算法,获取每一个个体的非支配序值irank,并按照个体所在的非支配层Fj计算其对应的拥挤距离;
S443、在预设的第二迭代次数内根据对种群Ri计算得到的每一个个体的非支配序值irank进行择优选择,从中选取N个作为新的父代种群Pi+1,得到所述多目标函数的Pareto最优前沿的非支配解集。
S5、根据列车到达当前运行区间的实际运行时间和时刻表规定运行时间的时间差值,从下一运行区间对应的Pareto最优解集中获取最优速度曲线,以使列车在下一运行区间按照所述最优速度曲线运行;最优速度曲线用于表示时刻表规定运行时间为各个区间预设时刻表时间、预设停车时间和时间差值的和推荐速度曲线。
其中,如图3所示,在步骤S5中利用PID控制结构跟踪所述当前运行区间的最优速度曲线,得到模拟速度曲线;所述模拟速度曲线用于模拟最优速度曲线。
其中,在步骤S5中,列车到达当前运行区间的实际运行时间和时刻表规定运行时间的时间差值包含早点或晚点。当列车到达当前运行区间的时间运行时间比时刻表规定运行时间晚点Δt秒时,则从下一运行区间的Pareto解集中获取时刻表规定运行时间比预设时刻表时间少Δt秒的推荐速度曲线作为最优速度曲线,以使列车在下一运行区间以该最优速度曲线运行,减小列车在该线路的中的运行误差,提高列车的准时性。
同理,当列车到达当前运行区间的时间运行时间比时刻表规定运行时间早点Δt秒时,则从下一运行区间的Pareto解集中获取时刻表规定运行时间比预设时刻表时间多Δt秒的推荐速度曲线作为最优速度曲线,以使列车在下一运行区间以该最优速度曲线运行。
通过如下步骤模拟最优速度曲线:
对列车的状态空间模型进行线性化,得到线性化后的列车模型状态空间;
列车的状态空间模型为:
其中,为列车的加速度,M为列车质量,a、b、c为戴维斯参数,v为列车速度,fg为运行区间的坡道附加阻力;
列车模型状态空间为:
其中,为在vo处的偏导,为v在vo处的偏导,为u在vo处的偏导,vo=0;
利用PID控制结构和所述列车模型状态空间跟踪最优速度曲线,得到模拟速度曲线;
PID控制结构的传递函数为Gc(s)=KPKi(s)+Ki(s)/s+Kd(s)s,其中,Kp为比例控制系数,Ki为积分控制系数,Kd为微分控制系数。
下面,以具体的例子对本发明的多目标优化的列车节能运行方法的效果进行说明。
设定选取亦庄线亦庄火车站至次渠站为运行实例,两站间的距离为1334米,时刻表上的列车运行时间为110秒,各站的停站时间为45秒。读取到列车运行在该运行区间的牵引能耗为15.67kw/h。下表1是该运行区间的线路坡度信息,表2是列车参数,其中,亦庄火车站的出站1214米内的限制速度为80km/h,列车在其他线路区间的最大限制速度为55km/h。
表1亦庄火车站至次渠的线路坡度
车辆数据:列车的牵引力以及制动力的大小是列车速度的函数,在不同速度下对应不同的值。
其中,列车的牵引力为:
Fp=fp(v)(N/KN)
列车的制动力为:
fb=fb(v)(N/KN)
其中,Fp为列车受到的牵引力,fb表示列车的制动力,v表示列车的运行速度。
在列车空载的条件下,采用北京地铁亦庄线已有的牵引/制动特性曲线,然后经多项式拟合,分别得到其表达式。其中最大牵引力如下式所示:
最大制动力如下式所示:
在上式中,v表示列车运行速度,单位是km/h,FMax表示最大牵引力,单位是KN,BMax表示最大制动力,单位是KN。
表2北京地铁亦庄线列车参数
在使用多目标遗传算法求解时,采用北京地铁亦庄线的亦庄火车站至次渠站的线路数据,程序中种群的最大迭代代数设为100代,选择概率Ps的取值为0.9,交叉概率Pc的取值为0.85,变异概率Pm的取值为0.004。图5中的点为通过多目标遗传算法优化得到的Pareto最优解集,图6为在图5中找的运行时间为110秒所对应的点所对应的推荐速度曲线,即为原推荐速度曲线。
当列车根据预设的推荐速度曲线从亦庄火车站运行到次渠站,在这期间,若列车实际运行时间与列车运行时刻表规定的运行时间有误差,或者因为列车在次渠站多滞留了一些时间(假设比原先时刻表晚了15秒),则预设的从次渠站到下一站的推荐速度曲线便不再适用,若按照预设的推荐速度曲线运行会将延时效果一直累加,导致从次渠站运行到下一站的运行时间与时刻表相差更大。
此时如果列车在生成优化速度曲线时采用的是多目标Pareto寻优,当列车从亦庄火车站运行到次渠站后,在停站期间便可以在预存的Pareto解集中,对应原来的最优解,在时间特性上选一个对应运行时间少15秒的解,并调出该解所对应的列车运行速度曲线,以此保证列车到达下一站时对应的运行时间,能够满足时刻表的要求。
具体地:
如图4所示,若采用图6的原推荐速度曲线运行,但是列车到达亦庄火车站的时间比原先时刻表规定的时间晚了8秒,若不及时调整列车的运行曲线,列车很可能在次渠站比原计划晚到,给列车运行带来很大不便。所以如图6所示的原推荐速度曲线便不再适用,查看该速度曲线对应的列车运行时间为110.05秒,运行能耗为14.12kwh,于是在对应的Pareto解集图5中寻找运行时间段[100,104]秒内的对应解,由于段内不止一个解,所以需要选出运行时间最接近102秒的解,经过计算找出的解其对应运行时间为102.12秒,对应的运行能耗为15.52kwh,其对应的运行推荐速度曲线如图7所示,所以图7为列车晚点重新挑选的推荐速度曲线。
相反地,若采用图6的推荐速度曲线运行,但是列车到达亦庄火车站的时间比原先时刻表规定时间早了8秒,图6的推荐速度曲线便不再适用,查看该推荐速度曲线对应的列车运行时间为110.05秒,运行能耗为14.12kwh,于是在对应的Pareto解集图5中寻找运行时间段[116,120]秒内的对应解。由于段内不止一个解,所以需要选出运行时间最接近118秒的解,经过计算找出的解其对应运行时间为117.96秒,对应的运行能耗为13.48kwh,其对应的运行速度曲线如图8所示,所以图8为列车早到重新挑选的速度曲线。
列车早晚点后重新挑选的速度曲线其对应能耗与原推荐速度曲线以及亦庄线目前运行数据的结果对比如表3所示:
表3实例分析结果的对比
如表3所示,原推荐速度曲线不但满足时刻表的时间要求,而且节能效果良好。假如列车晚点,重新挑选的速度曲线不仅可以弥补运行时间,而且还有一定的节能效果;而假如列车早到,那么重新选择的速度曲线的节能效果则更加明显。总的来说,仿真结果表明采用Pareto寻找最优解优的多目标优化不仅可以用于解决列车运行中的运行时间误差问题,而且节能效率甚至可以达到13.97%。
采用地铁亦庄线全线进行仿真验证。地铁亦庄线共有14个车站,全长为23.3km。其中线路两端的宋家庄、肖村、小红口、亦庄火车站、次渠、次渠南六个属于地下车站,其余八个属于高架站。
图9-10所示为采用多目标遗传算法对地铁亦庄线全线上和下行进行Pareto寻优,根据列车时刻表的要求综合实际运行情况对每个站间挑选一个满足运行时间要求的解,以此找出对应的推荐速度曲线。通过对多目标遗传算法的优化结果与实际采集的亦庄线ATO运行数据进行对比,表明本发明改进后的多目标遗传算法可以很好地满足列车运行时间的要求,有效地减少列车的运行能耗,节能率稳定在8%左右,具有较好的使用效果。
表4优化算法的节能效果与准时性
通过本次仿真可知,结合列车定时节能优化和运行时间优化节能效果更佳,对地铁实际运行调度具有一定的参考意义。
与现有技术相比,本发明的一种基于多目标优化的列车节能运行操作策略对算法进行了优化,得到的是列车实际追踪推荐速度曲线运行后的实际运行节能效果,而且还可以通过对站间选取不同的速度曲线方式对列车晚点或早到问题进行弥补,使问题得到较好的解决。
以上详细描述了本发明的较佳具体实施例,应当理解,本领域的普通方法人员无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此,凡本方法领域中方法人员依本发明构思在现有方法基础上通过逻辑分析、推理或者根据有限的实验可以得到的方法方案,均应该在由本权利要求书所确定的保护范围之中。
Claims (8)
1.一种基于多目标优化的列车节能运行方法,其特征在于,包括以下步骤:
S1、获取列车运行线路中各个运行区间的限速信息、坡度信息、线路曲线半径信息和线路隧道长度信息;
S2、根据列车的质量、所述限速信息、所述坡度信息、所述线路曲线半径信息和所述线路隧道长度信息,计算所述列车在所述各个运行区间的牵引力、制动力、基本阻力和附加阻力,以构建所述列车的单质点模型;
S3、利用所述单质点模型分别计算所述列车在各个运行区间中的牵引能耗Et、制动能耗Er和惰行能耗Eb,以构建所述列车的多目标优化模型,其中所述多目标优化模型如下:
目标函数为minY=F(x)=f(Ke,Kt,Kc)T,
约束条件为
其中,所述目标函数F(x)用于指示各个运行区间推荐速度曲线x的函数,Ke为列车的运行能耗,Kt为列车的运行时间,Kc为列车的运行舒适度,Ta为时刻表规定运行时间,T为实际运行时间,ta的值表示可以接受的时间误差,σ模型的宽度系数,ak为目标舒适度;
S4、采用多目标遗传算法对所述目标函数F(x)进行求解,得到各个运行区间的Pareto最优解集;其中,所述Pareto最优解集中的每个解用于指示一条推荐速度曲线;
S5、根据所述列车到达当前运行区间的实际运行时间和时刻表规定运行时间的时间差值,从所述下一运行区间对应的Pareto最优解集中获取最优速度曲线,以使所述列车在下一运行区间按照所述最优速度曲线运行;所述最优速度曲线用于表示时刻表规定运行时间为各个区间预设时刻表时间、预设停车时间和所述时间差值的和的推荐速度曲线。
2.根据权利要求1所述的列车节能运行方法,其特征在于,所述列车的单质点模型为:
其中,Fp表示控制每单位质量所受到的牵引力,fb表示每单位质量所受的制动力,f0(v)表示每单位质量所受基本阻力,fadd(x)表示列车受到的线路附加阻力,且fadd(x)=0,x1为列车在运行区间中的位置,v表示列车的运行速度,列车所经历的时间为t=t(x1)。
3.根据权利要求2所述的列车节能运行方法,其特征在于,还包括如下步骤:
利用PID控制结构跟踪所述当前运行区间的最优速度曲线,得到模拟速度曲线;所述模拟速度曲线用于模拟列车运行的所述最优速度曲线;
在列车完成所述当前运行区间的运行时,将所述模拟速度曲线与实际运行速度曲线进行比较,以计算所述最优速度曲线的节能率。
4.根据权利要求1所述列车节能运行方法,其特征在于,在步骤S4中通过如下步骤对所述目标函数F(x)进行求解:
S41、利用多目标遗传算法的生成所述多目标遗传算法的初始种群;所述初始种群为用于指示满足所述多目标函数的约束条件的解的集合,所述初始种群中的每一个个体对应所述多目标优化函数的一个解,所述初始种群的大小为N;
S42、对所述初始种群中的所述个体进行进化选择,并使用快速非支配排序对所述个体进行分层,以使所述Pareto最优解集逼近全局Pareto最优前沿;
S43、计算同层间个体的拥挤距离,使得Pareto最优解集均匀分布在所述全局Pareto最优前沿上;
S44、采用精英保留策略在预设的第一迭代次数内产生子代种群,得到所述多目标函数的Pareto最优前沿的非支配解集,以将所述Pareto最优前沿的非支配解集设定为所述目标函数F(x)的解。
5.根据权利要求4所述的列车节能运行方法,其特征在于,在所述步骤S42中通过如下步骤对所述个体进行分层包括如下步骤:
S421、从所述初始种群群体中获取至少两个非支配解集;
S422、针对所述两个非支配解集中的第i个非支配解集,将所述第i个非支配解集设置为第i非支配层Fi,将所述第i非支配层中的所有个体的非支配序值iranki都赋值为i,并将所述第i非支配层中所有个体从所述初始种群中除去;其中,i为整数,且1≤i≤N。
6.根据权利要求4所述的列车节能运行方法,其特征在于,在所述步骤S43中所述拥挤距离计算方法包括如下步骤:
S431、将处于同一层的个体的拥挤距离L[i]d初始化为0;
S432、对处于同一层的个体按照第m个目标函数值进行升序排列,m为整数,且,m<N;
S433、将位于排序两端的个体的拥挤距离设定为无穷大;
S434、对于排序中间的个体,采用如下公式计算其拥挤距离:
其中,L[i+1]m为第i+1个个体在第m个目标函数处的取值,和分别表示第m个目标函数的最大值和最小值;
S435、选择拥挤距离在预设数值范围内的个体组进行保留,以使所述多目标遗传算法迭代后得到的种群中的个体均匀分布于所述全局Pareto最优前沿上。
7.根据权利要求4所述的列车节能运行方法,其特征在于,在所述步骤S44中通过如下步骤得到述多目标函数的Pareto最优前沿的非支配解集:
S441、将进化选择后的初始种群和所述子代种群合并,形成大小为2N的新种群Ri;
S442、对所述新种群Ri采用快速非支配排序分层算法,获取每一个个体的非支配序值irank,并按照个体所在的非支配层Fj计算其对应的拥挤距离;
S443、在预设的第二迭代次数内根据对种群Ri计算得到的每一个个体的非支配序值irank进行择优选择,从中选取N个作为新的父代种群Pi+1,得到所述多目标函数的Pareto最优前沿的非支配解集。
8.根据权利要求3所述的列车节能运行方法,其特征在于,通过如下步骤模拟所述最优速度曲线:
对列车的状态空间模型进行线性化,得到线性化后的列车模型状态空间;
所述列车的状态空间模型为:
其中,为列车的加速度,M为列车质量,a、b、c为戴维斯参数,v为列车速度,fg为运行区间的坡道附加阻力;
所述列车模型状态空间为:
其中,为在vo处的偏导,为v在vo处的偏导,为u在vo处的偏导,vo=0;
利用PID控制结构和所述列车模型状态空间跟踪所述最优速度曲线,得到所述模拟速度曲线;
所述PID控制结构的传递函数为Gc(s)=KPKi(s)+Ki(s)/s+Kd(s)s,其中,Kp为比例控制系数,Ki为积分控制系数,Kd为微分控制系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910099114.9A CN109858154A (zh) | 2019-01-31 | 2019-01-31 | 一种基于多目标优化的列车节能运行方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910099114.9A CN109858154A (zh) | 2019-01-31 | 2019-01-31 | 一种基于多目标优化的列车节能运行方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109858154A true CN109858154A (zh) | 2019-06-07 |
Family
ID=66897239
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910099114.9A Pending CN109858154A (zh) | 2019-01-31 | 2019-01-31 | 一种基于多目标优化的列车节能运行方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109858154A (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110641523A (zh) * | 2019-10-15 | 2020-01-03 | 大连海事大学 | 一种地铁列车实时速度监测方法及*** |
CN111325462A (zh) * | 2020-02-18 | 2020-06-23 | 中国铁道科学研究院集团有限公司 | 动车组辅助驾驶方法及*** |
CN111652550A (zh) * | 2020-05-28 | 2020-09-11 | 优信数享(北京)信息技术有限公司 | 一种智能寻找最佳环线集合的方法、***及设备 |
CN112116207A (zh) * | 2020-08-24 | 2020-12-22 | 北京交通大学 | 一种多约束条件的列车运行调整的计算方法和装置 |
CN112158238A (zh) * | 2020-10-28 | 2021-01-01 | 重庆交通大学 | 基于目标和客流差异的列车自动驾驶曲线生成方法 |
CN112224244A (zh) * | 2020-10-28 | 2021-01-15 | 重庆交通大学 | 基于温度、载重条件的高速列车自动驾驶曲线生成方法 |
CN112896244A (zh) * | 2021-02-22 | 2021-06-04 | 同济大学 | 基于目标速度曲线跟踪的列车站间运行控制方法及装置 |
CN113158325A (zh) * | 2021-01-29 | 2021-07-23 | 北京交通大学 | 基于仿真的列车能耗-时间Pareto曲线生成方法 |
CN113361150A (zh) * | 2021-08-11 | 2021-09-07 | 华东交通大学 | 一种城市列车运行多目标优化方法及*** |
CN113361061A (zh) * | 2020-03-06 | 2021-09-07 | 中移智行网络科技有限公司 | 列车操纵策略优化方法、装置、存储介质和计算机设备 |
CN113591229A (zh) * | 2021-09-01 | 2021-11-02 | 北京建筑大学 | 一种高速列车制动距离计算方法及*** |
CN113665637A (zh) * | 2021-09-09 | 2021-11-19 | 中国铁道科学研究院集团有限公司通信信号研究所 | 一种基于双目标的高速铁路列车运行调整方法及装置 |
CN113971472A (zh) * | 2020-07-23 | 2022-01-25 | 南京理工大学 | 一种基于动态规划和时间分配的列车运行节能方法 |
CN113997988A (zh) * | 2020-07-28 | 2022-02-01 | 株洲中车时代电气股份有限公司 | 一种附加阻力计算方法及相关装置 |
CN114298510A (zh) * | 2021-12-19 | 2022-04-08 | 南京理工大学 | 基于nspso算法的时刻表和速度曲线优化方法 |
CN114580290A (zh) * | 2022-03-10 | 2022-06-03 | 中南大学 | 一种能耗和舒适度联合优化的高速列车组协同控制方法及*** |
CN116547188A (zh) * | 2020-11-04 | 2023-08-04 | 三菱电机株式会社 | 铁路***、运行管理装置、运行管理方法及运行管理程序 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011031697A (ja) * | 2009-07-31 | 2011-02-17 | Hitachi Ltd | 運行管理方法,運行管理装置,走行制御装置 |
CN104134378A (zh) * | 2014-06-23 | 2014-11-05 | 北京交通大学 | 一种基于驾驶经验和在线学习的城轨列车智能控制方法 |
CN106503804A (zh) * | 2016-10-11 | 2017-03-15 | 南京理工大学 | 一种基于Pareto多目标遗传算法的列车定时节能运行方法 |
CN107180282A (zh) * | 2017-07-01 | 2017-09-19 | 南京理工大学 | 一种高峰时段多列车节能优化方法 |
CN107515537A (zh) * | 2017-08-21 | 2017-12-26 | 南京理工大学 | 一种基于改进型暴力搜索法的列车单车节能运行方法 |
CN107704975A (zh) * | 2017-06-30 | 2018-02-16 | 湖南工业大学 | 基于生物地理算法的地铁列车节能运行优化方法与*** |
CN108791367A (zh) * | 2018-06-01 | 2018-11-13 | 广州地铁设计研究院有限公司 | 列车的节能操纵方法 |
-
2019
- 2019-01-31 CN CN201910099114.9A patent/CN109858154A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011031697A (ja) * | 2009-07-31 | 2011-02-17 | Hitachi Ltd | 運行管理方法,運行管理装置,走行制御装置 |
CN104134378A (zh) * | 2014-06-23 | 2014-11-05 | 北京交通大学 | 一种基于驾驶经验和在线学习的城轨列车智能控制方法 |
CN106503804A (zh) * | 2016-10-11 | 2017-03-15 | 南京理工大学 | 一种基于Pareto多目标遗传算法的列车定时节能运行方法 |
CN107704975A (zh) * | 2017-06-30 | 2018-02-16 | 湖南工业大学 | 基于生物地理算法的地铁列车节能运行优化方法与*** |
CN107180282A (zh) * | 2017-07-01 | 2017-09-19 | 南京理工大学 | 一种高峰时段多列车节能优化方法 |
CN107515537A (zh) * | 2017-08-21 | 2017-12-26 | 南京理工大学 | 一种基于改进型暴力搜索法的列车单车节能运行方法 |
CN108791367A (zh) * | 2018-06-01 | 2018-11-13 | 广州地铁设计研究院有限公司 | 列车的节能操纵方法 |
Non-Patent Citations (1)
Title |
---|
罗运真: "基于多目标优化的列车节能操纵策略研究", 《北京交通大学硕士论文集》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110641523B (zh) * | 2019-10-15 | 2021-06-29 | 大连海事大学 | 一种地铁列车实时速度监测方法及*** |
CN110641523A (zh) * | 2019-10-15 | 2020-01-03 | 大连海事大学 | 一种地铁列车实时速度监测方法及*** |
CN111325462A (zh) * | 2020-02-18 | 2020-06-23 | 中国铁道科学研究院集团有限公司 | 动车组辅助驾驶方法及*** |
CN113361061A (zh) * | 2020-03-06 | 2021-09-07 | 中移智行网络科技有限公司 | 列车操纵策略优化方法、装置、存储介质和计算机设备 |
CN111652550A (zh) * | 2020-05-28 | 2020-09-11 | 优信数享(北京)信息技术有限公司 | 一种智能寻找最佳环线集合的方法、***及设备 |
CN113971472A (zh) * | 2020-07-23 | 2022-01-25 | 南京理工大学 | 一种基于动态规划和时间分配的列车运行节能方法 |
CN113997988A (zh) * | 2020-07-28 | 2022-02-01 | 株洲中车时代电气股份有限公司 | 一种附加阻力计算方法及相关装置 |
CN112116207B (zh) * | 2020-08-24 | 2023-12-22 | 北京交通大学 | 一种多约束条件的列车运行调整的计算方法和装置 |
CN112116207A (zh) * | 2020-08-24 | 2020-12-22 | 北京交通大学 | 一种多约束条件的列车运行调整的计算方法和装置 |
CN112224244A (zh) * | 2020-10-28 | 2021-01-15 | 重庆交通大学 | 基于温度、载重条件的高速列车自动驾驶曲线生成方法 |
CN112158238A (zh) * | 2020-10-28 | 2021-01-01 | 重庆交通大学 | 基于目标和客流差异的列车自动驾驶曲线生成方法 |
CN112158238B (zh) * | 2020-10-28 | 2022-05-24 | 重庆交通大学 | 基于目标和客流差异的列车自动驾驶曲线生成方法 |
CN116547188B (zh) * | 2020-11-04 | 2024-01-30 | 三菱电机株式会社 | 铁路***、运行管理装置、运行管理方法 |
CN116547188A (zh) * | 2020-11-04 | 2023-08-04 | 三菱电机株式会社 | 铁路***、运行管理装置、运行管理方法及运行管理程序 |
CN113158325A (zh) * | 2021-01-29 | 2021-07-23 | 北京交通大学 | 基于仿真的列车能耗-时间Pareto曲线生成方法 |
CN113158325B (zh) * | 2021-01-29 | 2023-12-19 | 北京交通大学 | 基于仿真的列车能耗-时间Pareto曲线生成方法 |
CN112896244B (zh) * | 2021-02-22 | 2022-07-19 | 同济大学 | 基于目标速度曲线跟踪的列车站间运行控制方法及装置 |
CN112896244A (zh) * | 2021-02-22 | 2021-06-04 | 同济大学 | 基于目标速度曲线跟踪的列车站间运行控制方法及装置 |
CN113361150B (zh) * | 2021-08-11 | 2021-12-24 | 华东交通大学 | 一种城市列车运行多目标优化方法及*** |
US11390306B1 (en) | 2021-08-11 | 2022-07-19 | East China Jiaotong University | Method and system for multi-objective optimization of urban train operation |
CN113361150A (zh) * | 2021-08-11 | 2021-09-07 | 华东交通大学 | 一种城市列车运行多目标优化方法及*** |
CN113591229A (zh) * | 2021-09-01 | 2021-11-02 | 北京建筑大学 | 一种高速列车制动距离计算方法及*** |
CN113591229B (zh) * | 2021-09-01 | 2023-05-26 | 北京建筑大学 | 一种高速列车制动距离计算方法及*** |
CN113665637A (zh) * | 2021-09-09 | 2021-11-19 | 中国铁道科学研究院集团有限公司通信信号研究所 | 一种基于双目标的高速铁路列车运行调整方法及装置 |
CN114298510A (zh) * | 2021-12-19 | 2022-04-08 | 南京理工大学 | 基于nspso算法的时刻表和速度曲线优化方法 |
CN114580290A (zh) * | 2022-03-10 | 2022-06-03 | 中南大学 | 一种能耗和舒适度联合优化的高速列车组协同控制方法及*** |
CN114580290B (zh) * | 2022-03-10 | 2024-08-02 | 中南大学 | 一种能耗和舒适度联合优化的高速列车组协同控制方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109858154A (zh) | 一种基于多目标优化的列车节能运行方法 | |
Sun et al. | A bi-objective timetable optimization model for urban rail transit based on the time-dependent passenger volume | |
Scheepmaker et al. | Review of energy-efficient train control and timetabling | |
Fernández-Rodríguez et al. | Balancing energy consumption and risk of delay in high speed trains: A three-objective real-time eco-driving algorithm with fuzzy parameters | |
Sun et al. | Multi-objective optimization of train routing problem combined with train scheduling on a high-speed railway network | |
Xu et al. | A multi‐objective subway timetable optimization approach with minimum passenger time and energy consumption | |
CN110782079B (zh) | 基于停站方案节能的高速铁路列车运行图的调整方法 | |
CN105930937B (zh) | 一种考虑地铁快慢线的列车运行调度一体化优化方法 | |
CN109840641B (zh) | 一种列车多区间运行曲线快速优化方法 | |
CN105448082B (zh) | 一种可变间隔发车的快速公交车组合调度方法 | |
CN106652434A (zh) | 一种基于轨道交通协调的公交调度方法 | |
Luo et al. | Stop plan of express and local train for regional rail transit line | |
Zhang et al. | Data-driven train operation models based on data mining and driving experience for the diesel-electric locomotive | |
CN109398426B (zh) | 一种定时条件下基于离散蚁群算法的节能驾驶策略寻优方法 | |
Zhang et al. | Real-time energy saving optimization method for urban rail transit train timetable under delay condition | |
Li et al. | Joint optimization of delay-recovery and energy-saving in a metro system: A case study from China | |
CN114818349A (zh) | 基于时空网络客流状态估计的节能运行图优化分析方法 | |
Zhang et al. | Pareto multi-objective optimization of metro train energy-saving operation using improved NSGA-II algorithms | |
CN114298385A (zh) | 考虑客流影响和再生制动能利用的地铁列车延误调整方法 | |
CN116184907A (zh) | 一种基于多车协同优化的城轨综合***节能方法 | |
CN113592419B (zh) | 考虑客流和节能的轨道交通快慢车时刻表优化方法 | |
Fan et al. | Speed profile optimization for train operation based on ant colony algorithm | |
Chen et al. | An energy-efficient operation approach based on the utilization of regenerative braking energy among trains | |
Cao et al. | Optimisation of recommended speed profile for train operation based on ant colony algorithm | |
Yin et al. | Online adjusting subway timetable by q-learning to save energy consumption in uncertain passenger demand |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190607 |
|
RJ01 | Rejection of invention patent application after publication |