• 查询稿件
  • 获取最新论文
  • 知晓行业信息
官方微信 欢迎关注

数字轨道地图线路特征分段拟合方法研究

袁祎, 陈光武

袁祎, 陈光武. 数字轨道地图线路特征分段拟合方法研究[J]. 铁路计算机应用, 2022, 31(9): 7-13. DOI: 10.3969/j.issn.1005-8451.2022.09.02
引用本文: 袁祎, 陈光武. 数字轨道地图线路特征分段拟合方法研究[J]. 铁路计算机应用, 2022, 31(9): 7-13. DOI: 10.3969/j.issn.1005-8451.2022.09.02
YUAN Yi, CHEN Guangwu. Segmented fitting method of line features of digital track map[J]. Railway Computer Application, 2022, 31(9): 7-13. DOI: 10.3969/j.issn.1005-8451.2022.09.02
Citation: YUAN Yi, CHEN Guangwu. Segmented fitting method of line features of digital track map[J]. Railway Computer Application, 2022, 31(9): 7-13. DOI: 10.3969/j.issn.1005-8451.2022.09.02

数字轨道地图线路特征分段拟合方法研究

基金项目: 国家自然科学基金资助项目(61863024);甘肃省科技引导计划项目(2020-61)
详细信息
    作者简介:

    袁 祎,在读硕士研究生

    陈光武,教授

  • 中图分类号: U212.2 : P228.4 : TP39

Segmented fitting method of line features of digital track map

  • 摘要: 基于全球卫星导航系统(GNSS ,Global Navigation Satellite System)的列车定位技术降低了铁路建设与运营的成本,数字轨道地图作为地图匹配方法的基础,通常作为辅助定位手段对卫星定位偏差进行校正,对轨道线形进行准确描述及约简是轨道电子地图构建存储的前提。文章采用方位角近似估计曲率的方法对轨迹平面线形进行初步识别分类,为提高线形的拟合精度,引入自适应权值因子,对最小二乘拟合法进行优化,以消除大偏差观测数据对拟合趋势的影响,从而保证整体轨迹数据的平滑拟合及精度。实测车辆轨迹数据验证表明,该算法可对车辆轨迹线路特征进行有效提取及准确描述。
    Abstract: The train positioning technology based on the Global Navigation Satellite System (GNSS) reduces the cost of railway construction and operation. As the basis of map matching method, digital orbit map is usually used as an auxiliary positioning means to correct satellite positioning deviation. Accurate description and reduction of track alignment are the premise of the construction and storage of track electronic map. This paper used the azimuth approximate estimation of curvature to preliminarily identify and classify the track plane alignment. In order to improve the fitting accuracy of the line shape, the paper introduced an adaptive weight factor to optimize the least squares fitting method, and eliminate the influence of large deviation observation data, thereby ensure the smooth fitting and accuracy of the overall trajectory data. The verification of measured vehicle trajectory data shows that the algorithm can effectively extract and accurately describe the characteristics of vehicle trajectory.
  • 随着我国铁路建设的不断推进,高速列车在运行中能耗问题的重要性愈发彰显。列车在运行过程中有4种工况,能源消耗主要集中在牵引和制动这两种工况。列车运行曲线对应不同工况组合下的列车能耗,因此研究能耗最小的列车运行曲线具有重要的现实意义。

    目前,众多学者对能耗优化问题进行了研究,文献[1]将列车节能问题转换为一定约束条件下的寻优问题;文献[2]利用极大值原理求解列车转换点;文献[3]提出通过划分坡道求解列车转换原则和运行位置的关系;文献[4]引入遗传算法,优化列车运行曲线的惰行距离,通过调整惰行点实现节能效果;文献[5]基于遗传算法,讨论列车运行曲线优化及晚点恢复等问题;文献[6]利用遗传算法解决列车运行调整问题;文献[7]利用遗传算法对列车运行过程中最优的组合序列进行求解。由于遗传算法的固有缺陷,一些学者开始研究针对列车节能问题的算法改进。文献[8]利用改进后变长算子的遗传算法,对城轨列车节能运行问题进行优化;文献[9]通过在算法中加入准点调整和局部搜索的引导机制对传统遗传算法进行改进。

    本文在前人研究的基础上,通过建立约束条件下的列车运行模型和对遗传算法中算子的改进,增强算法的全局搜索能力和收敛能力,同时,引入自适应机制,并通过Matlab仿真,验证了改进后算法的有效性。

    列车在运行过程中受到牵引力、制动力、运行阻力的共同作用[10],根据牛顿第二定律,对运行列车进行受力分析可得,列车在正常行驶过程中的动力学方程为:

    $$\frac{{dv(x)}}{{dt}} = \frac{{F(v) - B(v)}}{M} - {w_0}(v) - {w_f}(x)$$ (1)

    其中,$v$是列车当前运行速度;$t$是列车行驶时间;$F(v)$为列车在当前速度下的牵引力;$B(v)$为列车在当前速度下的制动力;${w_0}(v)$为列车在当前速度下的单位运行基本阻力,计算方式如式(2)所示;${w_f}(x)$为列车当前位置运行单位附加阻力;$M$为列车总质量;$x$为列车当前运行里程数。

    $${w_0}(v) = a + bv + c{v^2}$$ (2)

    其中,$a、b、c$均为常数,一般根据运行车型来取值。

    列车运行单位附加阻力${w_f}(x)$ 包含单位坡道附加阻力${w_p}(x)$、单位曲线附加阻力${w_q}(x)$、单位隧道附加阻力${w_s}({x})$。由于列车长度较长,计算阻力时将其视为单一质点并不可靠,因此考虑将列车的每节车厢视为单质点,来计算列车单位附加阻力,而列车所受的总附加阻力为各个车厢所受附加阻力之和,如公式(3)所示:

    $$\begin{array}{l} {w_f}(x) = \displaystyle\sum\limits_{n = 1}^m {[{w_p}({x_n}) + {w_q}({x_n}) + {w_s}({x_n})]} \\ \left\{ \begin{array}{l} {w_p}\left( {{x_n}} \right) = p\left( {{x_n}} \right)\\ {w_q}\left( {{x_n}} \right) = \dfrac{{600}}{{R\left( {{x_n}} \right)}}\\ {w_s}\left( {{x_n}} \right) = 0.000\;13 \cdot {L_s}\left( {{x_n}} \right) \end{array} \right. \end{array}$$ (3)

    其中,${w_p}({x_n})$为第$n$节列车在$x$处的单位坡道附加阻力;$p({x_n})$是该位置的坡道千分数;${w_q}({x_n})$是第$n$节列车在x处所受到的单位曲线附加阻力;$R({x_n})$是第$n$节列车在$x$位置处的曲线半径;${w_s}({x_n})$是第$n$节列车在$x$处的单位隧道附加阻力;${L_s}({x_n})$是第$n$节列车在$x$位置处的轨道长度。

    线路附加阻力在同一位置是固定的,为方便计算,考虑提前处理线路附加阻力,并将其存储起来,后续使用时根据车头与列车各车厢相对距离判断所处区间,从而得到附加阻力。

    将列车运行能耗最小作为优化目标,同时需要满足准时、定点等安全约束条件,建立列车节能运行数学模型[11]如下:

    $$\begin{array}{l} \min E = \min \displaystyle\int_0^T {F(t) \cdot v} (t)dt\\ s.t.\left\{ \begin{array}{l} v = \dfrac{{ds}}{{dt}},\\ v(0) = v(S) = 0,\\ v(x) \leqslant {v_{\lim }}(x),\\ \displaystyle\int_0^S {v{{(x)}^{ - 1}}dx} = T,\\ t(0) = 0,t(S) = T \end{array} \right. \end{array}$$ (4)

    其中,$E$是列车从出发到终点所消耗的能量,本文只考虑列车牵引能耗;$F(t)$$t$时刻的牵引力;$v(t)$$t$时刻列车的速度;$S$是出发站到终点站的距离;$T$是列车总运行时间;${v_{\lim }}(x)$是当前路段的最大速度限制。

    列车在同一路线中采取不同的运行策略,运行能耗也不同。本文将(xi , Mqi)称为工况转换点,其中,xi为列车工况切换点的位置;Mqi为对应工况。根据列车节能操纵原则[12],列车起动时以最大牵引力起动,制动时以最大制动力停车,有利于节能。如图1所示,列车运行工况可分为牵引(${{\rm{SQ}}}$)、匀速(${{\rm{SY}}}$)、惰性(${{\rm{SD}}}$)和制动(${{\rm{SZ}}}$ ),列车以${\rm{SQ}}$工况从${{\rm{s_1}}}$处加速行驶至${{\rm{s_2}}}$,切换成${\rm{SY}}$工况匀速行驶至${{\rm{s_3}}}$,再切换为${\rm{SD}}$工况至${{\rm{s_4}}}$处,接着继续以${\rm{SY}}$工况行驶,到${{\rm{s_5}}}$转换为${\rm{SZ}}$工况直至终点。$[({{\rm{s}}_{\rm{1}}},{\rm{SQ}}),({{\rm{s}}_2},{\rm{SY}}),$$({{\rm{s}}_{\rm{3}}},{\rm{SD}}),({{\rm{s}}_{\rm{4}}},{\rm{SY}}), $$ ({{\rm{s}}_{\rm{5}}},{\rm{SZ}})]$ 称为列车的一个操纵序列可行解。实际运行中存在着无数种工况转换点的组合,因此需要找出所有可行解中能耗最低的运行策略。

    由本文设置的列车操纵序列推导出的列车运行曲线未必全符合运行条件,当出现停车位置或运行时间超过一定范围的不可行解时,则不予保留。本文停车精确度取10 m,运行时间精确度取60 s。

    图  1  列车工况组合

    列车制动时按最大减速度${a_{\max }}$制动,制动曲线从停车点开始反算。如图2所示,当列车运行曲线与停车制动曲线相交时,可切换为制动工况,因此只要制动曲线落于${S_{QMNL}}$内且与$x$轴相交,则可满足停车精确度要求。

    图  2  列车停车点设置

    数值法求解列车节能问题模型较复杂,而遗传算法适用于解决复杂情况的寻优问题[13],但传统的遗传算法存在易早熟和收敛性差的问题,因此本文考虑改进算法中的遗传算子,并引入进化逆转算子。

    (1)编码

    列车运行策略是由一系列上文讨论过的工况组成,将工况转换点$({{{x_i}}},{{{M_{qi}}}})$做为基因,运行区间内所有工况转换点的集合做为染色体。本文采用实数编码。染色体结构如图3所示,其中$ l $为染色体长度。

    图  3  染色体结构

    (2)种群初始化

    根据线路情况和列车运行工况设置原则,随机生成原始工况转换点位置,初始种群规模为n。本文设置初始种群规模为100。

    (3)适应值计算

    列车节能优化的目标函数为$\min E$,为方便计算,取目标函数的倒数为适应度值$fitness = \dfrac{1}{{\min E}}$。适应度越高的染色体越接近理想求解。

    (4)收敛条件

    以是否达到最大进化次数作为收敛条件,如果没有达到条件,则继续进化,否则退出。本文设置最大进化代数为200。

    采用轮盘赌法选择染色体以保证在进化过程中适应度较大的个体有较大的概率,进入下一代种群。

    通常,在算法初期普遍种群适应度值较低,染色体质量较差,后期染色体适应度值和质量都会相应提升,因此,考虑将普通遗传算法中交叉算子的概率改变为随着迭代次数变化的动态值:

    $$P_c^m = \left\{ \begin{array}{l} {P_{c,\max }} \cdot {e^{ - \frac{m}{M}}},{\text{若}}{P_{c,\max }} \cdot {e^{ - \frac{m}{M}}} < {P_{c,\min }}\\ {P_{c,\min }},{\text{其他}} \end{array} \right.$$ (5)

    其中,$m$表示当前进化代数;$M$为总进化代数;${P_{c,\max }}$${P_{c,\min }}$分别是预先设定的最大交叉概率和最小交叉概率。交叉概率随着进化代数的增加而变小,以便初期能够得到更多的新染色体,增大解空间,提高算法全局搜索能力,后期交叉概率相应减小,保证算法稳定收敛。

    实际运行中,一般加速转匀速的工况转换点不会轻易变化,只有${{({{{x}}}_{2},\;{{{x}}}_{3},\;{{{x}}}_{7},\;{{{x}}}_{8})}}$这类工况转换位置可以改变,而传统的变异算子随机性过强,无法保留上一代进化过程中的重要信息,因此针对这类关键基因点,考虑引入3代评价函数来逼近搜索,得到最佳值。如图4所示。

    图  4  列车转换点示意

    具体操作为,记录下每个关键转换点的更新方向,称为标志位,以及更新后的适应度值。标志位前进为1,后退为−1,若当前标志位连续3代适应度值都增加,则继续在该方向上搜索;若3代中连续2代适应度值减小,则更新标志位方向,即从第1代开始的位置反方向变异搜索。其他变异规则如表1所示。

    表  1  3代逼近搜索规则
    第1代++++
    第2代++++
    第3代+++
    下次变异方向保持当前标志位方向回朔保持当前标志位方向更新标志位方向回朔保持当前标志位方向更新标志位方向
    下载: 导出CSV 
    | 显示表格

    变异算子通过选择染色体中某个基因进行变异,以增强算法的局部搜索能力,提高种群的多样性,防止出现早熟现象。变异概率应随着迭代次数增多而趋小,且适应度较好的个体变异概率也应较小,以保留优良特性,适应度较低个体的变异概率应较大。

    $$\qquad\qquad\qquad\qquad\qquad P_b^m = \left\{ \begin{array}{l} {P_{b,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \alpha \right\} }}, {\text{若}}{P_{b,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \alpha \right\} }} > {P_{b,\min }}\\ {P_{b,\min }},{\text{其他}} \end{array} \right.$$ (6)

    其中,$fitnes{s_{\max }}$为当前种群最大适应度;$fitness({X_i})$为当前个体的适应度值;${P_{b,\max }}$为预设的最大变异概率;${P_{b,\min }}$为预设的最小变异概率;$\alpha $为系数调节因子。

    传统遗传算法易陷入局部搜索,且传统交叉算子的随机性易使父代的优秀基因在进化过程中消失。本文结合上述引导机制,为增大解空间,考虑引入逆转算子[14],经逆转操作后,适应度值变大的种群参与下一代进化,适应度趋近或更低的种群则舍弃。算法初期,整体种群适应度较差,多样性较低,逆转概率可以稍大,以增大解空间;算法后期,整体适应度较高,此时逆转操作的概率应当较低以保留优良个体。逆转概率设置如下:

    $$\qquad\qquad\qquad\qquad\qquad P_r^m = \left\{ \begin{array}{l} {P_{r,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \frac{m}{M} \cdot \beta \right\} }}, {\text{若}}{P_{r,\max }} \cdot {e^{ - \left\{ \left(\left| {1 - \frac{{fitness({X_i})}}{{fitnes{s_{\max }}}}} \right|\right) \cdot \frac{m}{M} \cdot \beta \right\} }} > {P_{r,\min }}\\ {P_{r,\min }},{\text{其他}} \end{array} \right.$$ (7)

    其中,${P_{r,\max }}$为预设的最大逆转概率;${P_{r,\min }}$为预设的最小逆转概率;$\beta $为系数调节因子。

    考虑本文所设置的染色体情况,进行逆转操作时,需要对染色体进行如下处理,工况转换点中运行工况应与逆转前和其位置最相近点的工况保持一致,以防过多不可行解的出现。

    为验证本文算法的有效性,选定运行线路为合福高铁(合肥—福州)中,合肥至无为的线路,仿真车型为CRH380B。列车共8节编组,线路总长度100.59 km,限速300 km/h,在66~68 km路段限速250 km/h,图定运行时间为1424 s。列车参数如表2所示。列车牵引制动曲线如图5所示。

    表  2  列车参数
    参数名称参数值
    列车总重495 t
    列车长度200.3 m
    最大运行速度300 km/h
    单位阻力(0.0257+0.050 7 v+0.000 505 57 v2)kN
    v为列车当前运行时速
    定员510 人
    下载: 导出CSV 
    | 显示表格

    将列车参数及线路信息输入列车运行仿真模型,对算法进行初始化。经过仿真分析,比较正常行驶下的列车运行曲线及根据本文提出的算法优化后的列车运行曲线,如图6所示,可以看出,优化后的列车运行曲线明显增加了惰性区段长度。计算可得优化前后的列车运行数据变化,如表3所示。

    图  5  列车特性曲线
    图  6  优化前后列车运行曲线仿真对比
    表  3  改进前后列车运行数据对比分析
    评价指标优化前指标优化后指标性能比较
    运行时间1430 s1459 s+2%
    能耗6.87×10 9 J6.13×10 9 J−10.7%
    下载: 导出CSV 
    | 显示表格

    表3可知,对高速列车运行曲线进行优化后,列车运行时间变为1459 s,满足列车准点性要求,而能耗由6.87×109 J降为6.13×109 J,降低了10.7%的能源消耗,本论文提出的算法在满足运行时间指标的情况下有效降低了列车运行中的能耗。

    (1)本文基于既有的列车节能操纵策略,建立了以列车安全性、准时性为指标的列车优化模型,讨论了列车停车制动点的设置,为后续列车能耗计算提供基础;

    (2)在遗传算法的迭代过程中,提出一种3代逼近搜索的寻优引导机制,解决传统遗传算法求解速度较慢的问题。同时为了扩大解空间,引入逆转算子,进一步提高了算法搜索能力;

    (3)以合福高铁区间线路数据和CRH380B型列车为基础进行仿真,结果表明,采用本文提出的方法,列车运行能耗降低10.7%,同时满足列车停车的安全性和准点性,具有较好的理论参考价值。

  • 图  1   铁路线形平面示意

    图  2   圆弧切线法确定方位角示意

    图  3   平面线形识别算法流程

    图  4   车辆原始轨迹平面线形

    图  5   车辆轨迹点曲率变化示意

    图  6   直线、曲线初始分段属性

    图  7   轨迹拟合后的曲率示意

    图  8   车辆轨迹线形分段拟合示意

    表  1   铁路线形平面参数特点

    平面线形直线缓和曲线圆曲线
    半径/$ R \to \infty $$ R $
    方位角定值二次抛物线变化线性递增/递减
    曲率恒为0线性变化$ 1/R $
    下载: 导出CSV

    表  2   算法性能指标对比

    性能指标标准LS改进LS
    AVG1.04210.8213
    MAX1.65420.9356
    VAR0.35240.0572
    MSE0.54460.3741
    下载: 导出CSV
  • [1]

    OTEGUI J, BAHILLO A, LOPETEGI I, et al. A survey of train positioning solutions [J]. IEEE Sensors Journal, 2017, 17(20): 6788-6797. DOI: 10.1109/JSEN.2017.2747137

    [2] 夏 青,杨剑锋. 数字轨道地图线路分段表示方法研究 [J]. 铁路计算机应用,2012,21(4):9-11. DOI: 10.3969/j.issn.1005-8451.2012.04.003
    [3] 倪晓林,杨 妮,张 东. 列控数字轨道地图自动生成算法研究 [J]. 铁道运营技术,2010,16(1):24-27. DOI: 10.3969/j.issn.1006-8686.2010.01.009
    [4] 左自辉,王开锋,许 聪,等. 生成高精度轨道电子地图的数据处理方法 [J]. 中国铁道科学,2016,37(4):134-138. DOI: 10.3969/j.issn.1001-4632.2016.04.21
    [5] 龚 利,沙建领,赵 超,等. 基于北斗定位技术的铁路机车轨道电子地图生成算法 [J]. 计算机系统应用,2019,28(6):172-177. DOI: 10.15888/j.cnki.csa.006945
    [6] 刘 江,蔡伯根,唐 涛,等. 基于GPS的列控轨道地图数据生成方法研究 [J]. 测绘学报,2011,40(1):111-117.
    [7] 陈德旺,蔡伯根,王 剑,等. 轨道交通GPS数据约简的数学模型与算法研究 [J]. 铁道学报,2008(4):116-119. DOI: 10.3321/j.issn:1001-8360.2008.04.023
    [8] 曾 强,陈德旺,王丽娟,等. 基于主曲线和自适应半径的多GPS轨迹数据融合算法 [J]. 铁道学报,2015,37(2):46-51. DOI: 10.3969/j.issn.1001-8360.2015.02.007
    [9] 陈德旺. 铁路卫星定位数据的智能处理[M]. 北京: 北京交通大学出版社, 2014.
    [10]

    Camacho-Torregrosa F J, AM Pérez-Zuriaga, JM Campoy-Ungría, et al. Use of Heading Direction for Recreating the Horizontal Alignment of an Existing Road [J]. Computer-Aided Civil and Infrastructure Engineering, 2015, 30(4): 282-299. DOI: 10.1111/mice.12094

    [11]

    Winter H , Luthardt S , Willert V , et al. Generating Compact Geometric Track-Maps for Train Positioning Applications[C]//2019 IEEE Intelligent Vehicles Symposium (IV), June 09-12, 2019, Paris, France. New York, USA: IEEE, 2019.

    [12] 陶维杰,蔡伯根,王 剑,等. 数字轨道地图平面线形特征提取方法研究 [J]. 铁道学报,2019,41(9):81-87. DOI: 10.3969/j.issn.1001-8360.2019.09.011
    [13] 陈明晶,方源敏,陈 杰. 最小二乘法和迭代法圆曲线拟合 [J]. 测绘科学,2016,41(1):194-197,202. DOI: 10.16251/j.cnki.1009-2307.2016.01.038
  • 期刊类型引用(3)

    1. 刘文斌,马子彧,张春来,徐晓磊. 基于GIS技术的铁路网能力利用分析系统设计与应用. 铁路计算机应用. 2025(03): 38-43 . 本站查看
    2. 黄方圆,徐泽儒,郭禅. 基于GIS技术的多模式航运运输里程精细化分析. 中国航务周刊. 2024(48): 59-61 . 百度学术
    3. 杨文澈,刘茂朕,史晓磊,董卓皇,陈亚勋. 铁路货车里程管理及辅助分析应用系统设计及应用. 铁路计算机应用. 2024(12): 59-64 . 本站查看

    其他类型引用(0)

图(8)  /  表(2)
计量
  • 文章访问数:  112
  • HTML全文浏览量:  283
  • PDF下载量:  21
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-03-08
  • 刊出日期:  2022-09-28

目录

/

返回文章
返回