王博1，韩峰2
(1. 兰州交通大学 测绘与地理信息学院，甘肃 兰州 730070；2. 兰州交通大学 土木工程学院，甘肃 兰州 730070)
摘 要：既有曲线整正需要与最新的勘测技术相结合，进行简洁高效的线路坐标测量并设计相应的曲线整正算法是研究的关键所在。以实测线路坐标数据为基础，采用几何重心法求算曲线中心区域点位对应的初始圆心坐标及对应半径，利用圆心到切线垂距与初始半径的差值估算初始缓和曲线长度。以拨道量平方和最小为目标函数，设计相应的半径及缓和曲线长度优化算法，并结合现场实测数据，采用计算机程序实现设计的算法，取得较为满意的计算结果。
关键词：坐标测量；曲线整正；优化算法；几何重心法
轨道结构具有边运营边维修的特点，曲线位置变化一直是影响列车运营舒适度、安全性的重要因素。铁路工作者一直伴随测绘技术的不断发展研究相应算法。从绳正法、偏角法、坐标法直到轨检小车测量，测量的手段和效率不断提升，计算结果的可靠性也更为准确[1−4]。对于大量的既有线路而言，采用坐标法进行线路曲线测量具有易选取置镜点位置、整体测量对线路运营干扰小且工作效率比较高的特点而得以广泛应用[5−7]。目前基于坐标法进行曲线测量和整正的方法主要存在2方面的问题：一是对初始半径及缓和曲线长的计算，刘文涛等[8−9]采用圆曲线段最小二乘法求得初始半径及缓和曲线长，但在数据处理过程中圆曲线段长度难以确定且计算数据较大。候茂盛[10]采用圆曲线段3点数据计算初始半径及缓和曲线长，在圆曲线段较长时缺乏代表性。二是在拨道量优化算法的设计，任少伟等[11−12]通过最小拨道量求出最优半径及缓和曲线长，覃乃轩等[13]采用半径和缓和曲线长的不同组合求出最小拨道量来确定最优半径及缓长，算法虽然减少了对曲线的扰动但在计算过程中拨正前后曲线长度变化引起的无缝线路实际锁定轨温变化，不利于养护维修工作。本文根据实测线路坐标，采用几何重心法计算初始半径及缓和曲线长。结合最小曲线扰动和轨道长度不变性原则，利用初始曲线长，计算最优半径、缓长及拨道量。
曲线整正问题的实质就是如何利用现有的实测数据获取最为可靠的初始半径及缓和曲线长度，并借助优化算法实现其最优的半径及缓和曲线长度组合，获取满意的线形。基于坐标测量的曲线整正算法基本思路如图1所示。
从流程图中可以看出，曲线整正计算的关键环节主要包括曲线偏角、初始半径、初始缓和曲线长度计算，优化算法设计等方面。
DONG Jun-feng, TENG Fei, LI Pei-lei, SUN Ke-yan, FU Hong, GUO Wen-yuan, DING Guo-shan
1) 铁路既有线整正时整正前后必须满足：
① 曲线转角值不变；
② 整正前后轨道长度保持不变；
③ 坐标测量的起点和终点在曲线外的直线上，起点和终点的整正拨距为0，保证线路两端的位置是固定不变以确定直线的方向不变；
陈洋十指冰冷，白皙的皮肤透出暗黑的青紫色，罗恬觉得自己的视线变得有些模糊不清了，她隐约看见楼下的角落里，有个熟悉的身影，她呢喃地问：“是杜朗吗？”
2) 拨道量计算结果要求：在满足设计的各条件要求下各测点整正的拨道量平方和为最小值。
3) 特殊位置既有线整正要求：在实际的铁路既有线上整正过程中，当遇到整正既有曲线的控制点时，例如桥梁、隧道、涵洞、等大型建筑物，这些控制点的整正拨道量值要小于规定的数值。
4) 半径与缓和曲线的要求：优化拨道量的过程中半径的取值要符合《铁路线路设计规范》(TB11009—2017)表5.4.1中规定的最小曲线半径值。缓和曲线要符合《铁路线路设计规范》表5.4.3-1中规定的缓和曲线最小长度[14]。采用的曲线半径在困难条件下，可采用规定范围内10 m的整倍数。特殊困难条件下按曲线半径进整值进行取值。
图1 基于坐标测量的曲线整正算法流程图
Fig. 1 Flow chart of the modification algorithm of existing railway curves based on coordinate measurement
确定一条完整的铁路曲线位置需要确定线路的曲线转角α，交点坐标(XJD,YJD)，曲线半径R，缓和曲线长度l01和l02。
目前常用的初始半径与圆心坐标计算方法主要有3点法与最小二乘法。3点法是取圆曲线范围内的中间3个测点坐标求出初始半径与圆心坐标，当圆曲线段线段较长时，缺乏代表性。并且当中间3点的测点受列车影响较大时，计算结果缺乏可靠性。最小二乘法求初始半径与圆心坐标是通过圆曲线段所有测点经过最小化误差的平方和寻找数据的最佳值。首先通过n个测点3点定圆的方式求出n−2个圆心与半径，做出半径变化图或者做出曲率变化图，通过变化图中半径或曲率变化较大的点近似确定圆缓点与缓圆点确定圆曲线范围。在确定半径或曲率变化较大的点时，变化量没有统一的定义，圆曲线范围难以确定。并且最小二乘法求初始半径与圆心坐标计算量较大。
按照3点定圆的方法，从曲线中心第n点向两侧分别选取相邻的2个点位，共计5个点。利用n−2，n−1和n 3个点位确定可行半径Rn−1及圆心坐标(OXn−1,OYn−1)，利用n−1，n和n+1 3个点位确定可行半径Rn及圆心坐标(OXn,OYn)，利用n，n+1和n+2 3个点位确定可行半径Rn+1及圆心坐标(OXn+1,OYn+1)，如图2所示。
图2 几何重心求圆心示意图
Fig. 2 Schematic diagram of calculating the coordinates of the center based on geometrical-centering method
则初始半径Rs：
可将相应的(OXn−1,OYn−1)，(OXn,OYn)，(OXn+1, OYn+1)3点所构成的三角形的重心坐标作为初始圆心坐标，计算公式为：
对于较长的曲线可以向两侧选取较多的点位参与初始半径计算，利用多边形几何重心法求解相应的初始圆心坐标及半径。
多边形几何重心计算公式：较长的曲线向两侧选取较多的点位，根据3点成圆的方法，确定n个圆心坐标，由这n个坐标可以构成n边形。如图3所示，以A1，A2，…，An为顶点各坐标为(X1,Y1)，(X2,Y2)，…，(Xn,Yn)，可以把n边形分成n−2个三角形，各三角形重心为G1(x1,y1)，G2(x2,y2)，…，Gn−2(xn−2,yn−2)，面积分别为σ1，σ2，…，σn−2。则多边形几何重心如下所示。
式中：
图3 多边形几何重心示意图
Fig. 3 Schematic diagram of the geometrical-center of a polygon
通过几何重心法与传统的3点法和最小二乘法求初始的半径结合实例对比可以得出，几何重心法计算简单，计算结果符合实际。如表1所示。
根据上述数据可知，路面竣工后空隙率、稳定度和流值均满足规范要求，面层结构为灌入式半柔性面层，且与基层结合很好。说明该半柔性路面施工技术的应用，能有效保证路面施工质量要求。
表1 半径计算结果对比
Table 1 Comparison of radius calculation results
测区3点法几何重心法最小二乘法 所需测点半径/m所需测点半径/m所需测点半径/m 13605.65593.28605.6 23583.36602.714601.2 33810.95806.414805.4 43589.65603.111606.4 531 038.361 002.1251 002.4 63764.75795.615802.6 73778.65806.619806.9 83806.56801.418797.9 93609.36599.719603.3
由图4可以看出，初始圆心到2条切线的距离D1和D2理论上是初始半径与前后缓和曲线对应的内移距之和：
根据求解的P1和P2，选择初始缓和曲线长度：
当l01和l02相差不超过10 m时，取二者平均值作为初始缓和曲线长度，按照等长缓和曲线进行计算。
当l01和l02相差超过10 m时，各自取整，按照不等长缓和曲线进行计算。
目前常用的优化方法是组合法，以不同半径及缓长组合成多组(R, l01, l02)组合，通过计算机辅助计算求出最优的拨道量。此类优化方法虽能求出拨道量最优值但拨正前后曲线长度变化引起的无缝线路实际锁定轨温变化，不利于养护维修工作[15]，基于这种情况以既有线初始半径及缓和曲线长下的曲线长L为约束条件对拨道量进行优化，按缓和曲线长的设置情况其优化可分为以下2种情况。
1) 对称优化：即前后缓和曲线长度相等，此时优化参数有2个：缓长l0和半径R。半径步长ΔR以曲线转角和步长关系进行取值，按正负2个方向取n组值。由式(10)可推导出等长缓和曲线长下半径与缓长关系式(11)，计算各R±ΔR下对应的l0，并将l0取整为10 m的整倍数,得到不同R与l0的组合，计算各组拨道量。以初始R与l0所对应的初始拨道量为参考值，以目标函数最小值，比选出最优的半径、缓长及拨道量。
2) 非对称优化：即前后缓和曲线长度不相等，此时优化参数有3个：缓长l01，l02和半径R。半径步长ΔR以曲线转角和步长关系进行取值，按正负2个方向取n组值。由式(10)可推导出不等长缓和曲线长下半径与缓长关系式(12)，固定l01，计算各R±ΔR下对应的l02，并将l02取整为10 m的整倍数,改变l01重复上述步骤，得到不同R与l01和l02的组合，计算各组拨道量。以初始R与l0所对应的初始拨道量为参考值，以目标函数最小值，比选出最优的半径、缓长及拨道量。
曲线拨道量目标函数：
图4 曲线偏角及交点坐标计算示意图
Fig. 4 Schematic diagram of calculating the deflection angle and intersection coordinate
曲线拨道量目标函数可以验证曲线整正结果是否良好，以各测点(即各里程点)的拨道量的平方和为计算结果。
通过组合法优化拨道量与曲线长约束优化拨道量结合实例得到的最优值看出，组合法拨道量最优值略好于曲线约束法，但曲线约束法优化拨道量对轨道曲线段长度的影响较小，可以减少无缝线路的养护维修工作，更好的满足曲线整正要求，如表2所示。
表2 优化方法计算对比
Table 2 Comparison of optimization methods
测区优化方法目标函数/m2曲线段长度变化/m 1曲线约束法0.0221.77 组合法0.011 5−15.51 2曲线约束法0.059 33.08 组合法0.034 1−8.30 3曲线约束法0.026 70 组合法0.023−4.63
由于铁路既有线整正是在铁路轨道丈量完成的条件下进行的，所以本算法首先要求完成轨道里程的丈量工作，曲线地段每20 m，曲线外侧明显的直线地段测量至少2个点，如图5中的K0，K1以及Kn−1，Kn以确定直线方向。剩余直线段每20 m为一个测点。测量过程中，置镜在轨道外，建立相对坐标系，得到各测点坐标Xi,Yi，及确定始点里程K0和计算步距ΔK，如图5所示。
图5 既有曲线坐标测量示意图
Fig. 5 Diagram of the existing curve of coordinate measurement
根据本文所设计的算法，结合铁路既有线实测数据(表5中原X与Y值)，通过计算机程序语言编制相应的程序，使计算更加快速准确。
为了培养学生的独立自主的精神和良好的自学品质，提高学生自我探索和发展新知识的兴趣和能力，我是从以下几方面入手：①从当代高科技社会对人才的要求出发，向学生阐明要学会生存、竞争，就得学会怎样学习的道理，使学生树立学会自学方法的使命感和紧迫感；②向学生介绍高尔基、毛泽东、周恩来、童第周等伟人、名人刻苦学习、立志成才的故事，鼓励学生树立信心和斗志，提高自学水平；③开展自学竞赛活动，交流自学方法、自学成果，鼓励创造新的自学方法，使学生获得掌握自学方法的成就感、自豪感，以此激发学生的自学热情。
1) 交点坐标及曲线转角计算
如图4所示，按照实测直线段点位A和B的坐标计算直线AB的方位角A1及直线方程，按照实测直线段点位C，D的坐标计算直线CD的方位角A2及直线方程，则曲线转角：
利用推求的直线方程可以求解曲线交点坐标(XJD,YJD)，如表3所示。
表3 曲线转角及交点坐标
Table 3 Deflection angle and intersection coordinate
名称曲线转角A交点坐标X交点坐标Y 值−0.3428 000.7175 004.177
2) 初始圆心、半径及缓长
在风景园林建设的多种要素中，植物和植物景观是与碳循环有着密切联系的核心要素[4]。因为在绿道建设中其他要素低碳效应在于减少碳排放，而植物和植物景观既与减少碳排放有关，又与碳固定有联系，因此在绿道规划设计中，需要找到合适的植物景观规划设计策略，才能最大化的减少碳排放。
根据本文所设计的算法，求出初始半径与缓和曲线长，如表4所示。
表4 既有线初始圆心、半径及缓和曲线长
Table 4 Initial center, radius and easement curve of the existing railway
名称初始圆心X初始圆心Y初始半径R初始l01初始l02 值/m7 421.4624 806.1126005050
3) 设计坐标及拨道量计算
根据选取的初始半径Rs，缓和曲线长度l01，l02，考虑交点坐标(XJD,YJD)，偏角α，初始里程LC等因素，采用文献16的算法可以简洁计算出曲线上各点理论坐标。
可以利用实测坐标和理论坐标之间的距离作为拨道量，并根据其与圆心之间的距离判断拨道方向。
4) 曲线整正优化
根据本文所设计的算法，结合铁路既有线实测数据确定最优半径为600 m缓和曲线长为50 m时拨道量值为最优值，曲线拨后坐标和拨道量如表5所示，最优目标函数值为0.029 5。
表5 拨道量及整正前后坐标值
Table 5 Correction distance and coordinate values before and after curve realignmen
里程/m原X值/m原Y值/m拨后X值/m拨后Y值/m拨道量/m 63 6407 923.5765 145.2727 923.5765 145.2720.000 63 6607 933.1535 127.7557 933.1535 127.7550.000 63 6807 942.7615 110.2047 942.7625 110.174−0.030 63 7007 952.2395 092.6277 952.2535 092.5690.060 63 7207 961.3945 074.8437 961.4145 074.7910.056 63 7407 970.0155 056.7427 970.0325 056.7440.017 63 7607 978.0425 038.3727 978.0445 038.4200.048 63 7807 985.4205 019.8247 985.4415 019.8390.026 63 8007 992.1925 001.0007 992.2145 001.0220.031 63 8207 998.3434 981.9807 998.3564 981.9900.016 63 8408 003.8604 962.7098 003.8614 962.7630.054 63 8608 008.7224 943.3268 008.7224 943.3640.038 63 8808 012.9424 923.7638 012.9344 923.8130.051 63 9008 016.5734 904.0658 016.5544 904.1140.061 63 9208 019.8414 884.3638 019.8154 884.4120.055 63 9408 022.9904 864.6128 022.9724 864.629−0.024 63 9608 026.1174 844.8948 026.1224 844.878−0.016 63 9808 029.2724 825.1288 029.2724 825.1280.000 64 1008 048.1654 706.6658 048.1654 706.6650.000
本文关于既有线曲线整正的算法在《铁路线路设计规范》(TB11009—2017)中客货共线铁路、重载铁路所要求的半径及缓和曲线长的情况下进行设计，所设计的算法也即适用于客货共线铁路、重载铁路的整正设计。
1) 通过实测数据计算采用几何重心法所求得的初始半径、缓长及圆心计算简单准确，计算结果符合工程实际。采用的优化算法可以减少无缝线路的养护维修工作，更好的满足曲线整正要求。
2) 坐标法在对既有线进行曲线整正计算中操作灵活，理论严密，计算简洁。并且能直观体现半径及缓和曲线长改变对坐标值与拨道量的影响。结合计算机语言辅助计算可以使计算过程更加快速，准确。
参考文献：
[1] 全志强, 杨育林. 既有铁路曲线整正测量[J]. 测绘科学, 2015, 40(4): 147−150. QUAN Zhiqiang, YANG Yulin. Realignment and survey of exiting line[J]. Science of Surveying and Mapping, 2015, 40(4): 147−150.
[2] 张金良. 偏角法和坐标法在既有铁路线路整正中的拓展研究及应用[D]. 南昌: 华东交通大学, 2017.ZHANG Jinliang. Study and application of deflection angle method and coordinate method based on existing railway lines[D]. Nanchang: East China Jiaotong University, 2017.
[3] YANG H, LI Long.Existing railway curve realignment constrained optimization algorithm research based on coordinates[J].Mathematics in Practice & Theory, 2009, 39(24): 166−171.
[4] 段晓峰, 韩峰, 李建强. 基于点云信息的既有线整正方法研究[J]. 铁道标准设计, 2013(8): 15−18. DUAN Xiaofeng, HAN Feng, LI Jianqiang.Method study on curve realignment of existing railway based on point cloud data[J]. Railway Standard Design, 2013(8): 15−18.
[5] Birkholz A B, Schlegel H B.Path optimization by a variational reaction coordinate method. I. Development of formalism and algorithms[J].Journal of Chemical Physics, 2015, 143(24): 363−368
[6] 于庆乐. 全站仪在铁路曲线测量上的应用[J]. 煤炭技术, 2012, 31(4): 152−153. YU Qingle. Aplication of whole station railway curve[J].Coal Technology, 2012, 31(4): 152−153.
[7] 孔晨辉, 翁呷. 线路曲线的测设[J]. 测绘地理信息, 2017, 42(6): 69−72. KONG Chenhui, WENG Xia. Surveying of the line curve[J]. Journal of Geomatics, 2017, 42(6): 69−72.
[8] 刘文涛. 基于坐标算法的铁路线路整正系统研究[D].兰州: 兰州交通大学, 2013.LIU Wentao.Research the railway realignment system based on the coordinate algorithm[D].Lanzhou: Lanzhou Jiaotong University, 2013.
[9] HOU M S, LIU Y X. Research on constrained optimization of coordinate method in railway curve realignment[J].Journal of Lanzhou Institute of Technology, 2014, 21(4): 38−41.
[10] 候茂盛. 基于最优化理论的既有铁路曲线整正方法研究及应用[D]. 兰州: 兰州交通大学, 2015.HOU Maosheng. Research on calculation method in existing railway curve realignment by optimization theory and its application[D].Lanzhou:Lanzhou Jiaotong University, 2015.
[11] 全志强, 杨育林. 既有铁路曲线整正测量[J]. 测绘科学, 2015, 40(4): 147−150. QUAN Zhiqiang, YANG Yulin.Realignment and survey of exiting line[J]. SurveyScience of Surveying and Mapping, 2015, 40(4): 147−150.
[12] 任少伟, 胡欣明. 运用坐标法进行既有复曲线拨距计算[J]. 铁道运营技术, 2011, 17(3): 58−60. REN Shaowei, HU Xinming. Coordinate method is used to calculate existing complex curve adjustment[J].Railway Operation Technology, 2011, 17(3): 58−60.
[13] 覃乃轩, 孙和金. 坐标法曲线整正计算方法探讨[J]. 铁道运营技术, 2013, 19(3): 49−51, 54. TAN Naixuan, SUN Hejin. Discussion on the method of existing railway curve correction by coordinate method [J]. Railway Operation Technology, 2013, 19(3): 49−51, 54.
[14] TB10098—2017, 铁路线路设计规范[S].TB10098—2017, Code for design of railway line[S].
[15] WANG B C, HAN F.Research on optimization method for the existing railway curve realignment based on coordinate measurement[J].Journal of Lanzhou Jiaotong University, 2008, 1(1): 11−13.
[16] 赵长石, 韩峰. 铁路中桩坐标统一计算模型[J]. 兰州交通大学学报, 2009, 28(1): 49−51. ZHAO Changshi, HAN Feng. Unified calculation model of railway curve coordinate[J].Journal of Lanzhou Jiaotong University, 2009, 28(1): 49−51.
Design and practice of existing railway curve realignment algorithm based on coordinate measurement
WANG Bo1, HAN Feng2
(1. Faculty of Geomatics, Lanzhou Jiaotong University, Lanzhou 730070, China; 2. College of Civil Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)
Abstract:Track geometry misalignment inevitably occurs, especially in curves. Modification of existing railway curves shall incorporate the latest surveying technology, the key to research lies in measuring coordinates with efficiency and in designing corresponding curve correction algorithm. The initial coordinate of the center and the radius of the curve with geometrical-centering method were calculated based on the measured coordinate data, the author also estimated the initial length of the transition curve using the difference between the initial radius and the vertical distance from curve center to tangent. Taking the quadratic sum of the modification volume of the track as objective function, an algorithm was designed to optimize the design of the length of transition curves and the corresponding radii. The algorithm was realized via computer program associated with the field measured data, and the satisfactory results was obtained.
Key words:coordinate measurement; curve realignment; optimization algorithm; geometric center-of-gravity method
中图分类号：U216.42+6
文献标志码：A
文章编号：1672 − 7029(2019)06− 1398 − 07
DOI: 10.19713/j.cnki.43−1423/u.2019.06.006
收稿日期：2018−07−26
基金项目：国家自然科学基金资助项目(51568037)
通信作者：韩峰(1975−)，男，陕西蒲城人，教授，博士，从事道路与铁路现代勘测与选线设计的研究及教学工作；E−mail：153025377@qq.com
(编辑 涂鹏)
联系客服
