鐵路既有線復測平面曲線優(yōu)化方法*

陳 峰，辜良瑤，楊 岳，吳湘華

(1.中南大學交通運輸工程學院，湖南長沙 410075;2.中南林業(yè)科技大學機電工程學院，湖南長沙 410004)

既有鐵路的線路參數(shù)是鐵路工務部門進行軌道檢測與維保的基準。鐵路線路經(jīng)過運營與維修保養(yǎng)后，其線路參數(shù)發(fā)生改變，因此需要對線路進行復測，根據(jù)復測數(shù)據(jù)，對線路參數(shù)進行優(yōu)化［1］。根據(jù)優(yōu)化后的線路參數(shù)計算線路橫、縱段面調(diào)整量，可以用于指導大機維修施工。目前，工務管理部門均成立了專業(yè)的測量隊伍，使用全站儀對鐵路既有線復測，獲得了高精度的連續(xù)線路大地坐標。如何利用高精度的線路大地坐標數(shù)據(jù)，通過合理的優(yōu)化計算方法對線路平面曲線參數(shù)進行優(yōu)化調(diào)整就成為一項十分重要的工作。對于線路平面曲線復測優(yōu)化計算方法，文獻［2］提出基于方向加速法的平面重構(gòu)算法，該算法未實現(xiàn)緩和曲線段的調(diào)整量計算，重構(gòu)出的曲線參數(shù)不能進行人工干預，得到的曲線優(yōu)化參數(shù)過于精確，無法適應現(xiàn)場的生產(chǎn)管理。如為了測設(shè)、施工和養(yǎng)護的方便，平面曲線半徑一般取50 m和100 m的整數(shù)倍［3］;文獻［4］提出采用基于圓曲線最小二乘擬合方法的平面曲線優(yōu)化方法，在弧徑比較小的情況下，圓曲線最小二乘擬合算法無法擬合出合理的圓曲線參數(shù)，因而在進行大曲線半徑平面曲線優(yōu)化時，該方法容易產(chǎn)生病態(tài)優(yōu)化結(jié)果。為了解決上述問題，本文根據(jù)平面曲線幾何要素計算原理，提出了基于夾直線最小二乘擬合的平面曲線優(yōu)化方法。

1 平面曲線特征分界點的判別

線路平面曲線由前夾直線、前緩和曲線、圓曲線、后緩和曲線和后夾直線5個線形要素組成［2］，分別用直緩(ZH)、緩圓(HY)、圓緩(YH)、緩直(HZ)4個特征點進行分界。既有線復測后，平面曲線優(yōu)化的首要步驟是對被測曲線的4個特征分界點的點位進行判別。

文獻［4］根據(jù)測量數(shù)據(jù)進行正矢計算，由正矢的變化規(guī)律進行HY和YH分界點的判別。圖1所示為一段測量線路的超高和20 m正矢圖。從圖1可以看出:由于測量誤差的影響，正矢是波動值，很難根據(jù)正矢值判別特征分界點，而且分界點的選擇偏差對圓曲線最小二乘擬合結(jié)果敏感度較高。

線路的超高值可以直接從測量中獲取，不需要進行換算，能較為直觀的區(qū)分出平面曲線各特征段，分界點的選擇偏差對夾直線最小二乘擬合結(jié)果敏感度較低，因此線路的測量超高更適用于特征分界點的判別。

圖1 測量線路超高和正矢圖Fig.1 Superel evation and versine of measured railway line

2 基于最小二乘的圓曲線擬合方法

在判斷出既有線圓曲線、緩和曲線的特征分界點后，即可選取圓曲線上的測點進行擬合。根據(jù)最小二乘原理來逼近圓曲線的圓心與半徑［4－6］。圓曲線上各測點的殘差平方和為:

式中:R0為擬合圓半徑;(X0，Y0)為擬合圓圓心坐標;(xk，yk)為圓曲線上測點Ck坐標。根據(jù)最小二乘原理建立圓曲線擬合半徑與圓心坐標計算方程:

在弧徑比L/R很小情況下，使用最小二乘擬合得到的圓心與半徑可能會得到一病態(tài)值，以下進行分析。

以圓心在(0，0)的標準圓(R=0.5)為例，取4點(0.09，0.49)，(0.28，0.41)，(0.37，0.34)和(0.49，0.1)，弧長 L 約為0.6，采用以上算法擬合得到的圓心為(0.01，0.01)、半徑為 0.485，如圖2 左圖所示。擬合圓與標準圓基本重合，擬合結(jié)果較為合理。

在標準圓附近另取 4個相近點(0.360 9，0.346)，(0.360 3，0.346 7)，(0.362 2，0.344 7)和(0.362 8，0.344)，弧 長 L=0.0133，L/R =0.026 6。采用以上算法得到的擬合結(jié)果:圓心為(0.361 5，0.345 3)，半徑為0.690 7。得到的擬合結(jié)果如圖2右圖所示，擬合圓與標準圓產(chǎn)生了較大差距，擬合結(jié)果不合理。

圖2 圓曲線擬合病態(tài)結(jié)果Fig.2 Illness result of circle curve fitting

究其原因，若弧徑比過小，則測量點上的微小測量誤差，將對圓心的擬合坐標產(chǎn)生較大影響，從而導致擬合結(jié)果與實際相差巨大。

在提速線路中大量使用了大半徑平面曲線，假設(shè)圓曲線長為150 m，曲線半徑超過6 000 m的平面曲線，其弧徑比即小于0.025，使用圓曲線最小二乘擬合，就有可能導致病態(tài)的擬合結(jié)果。

相對而言，直線的最小二乘擬合不存在以上問題，而且算法簡單。由于平面曲線可由圓曲線半徑和前后緩和曲線長3個參數(shù)唯一確定，因此，采用最小二乘擬合方法確定前后夾直線后，通過建立平面曲線上各線形要素的調(diào)整量計算模型，對圓曲線半徑和前后緩和曲線長進行等步長迭代優(yōu)化計算，計算出平面曲線的最小調(diào)整量，即可獲得最優(yōu)的平面曲線參數(shù)。

3 基于夾直線最小二乘的平面曲線擬合

3.1 最小二乘的夾直線擬合

夾直線用方程y=Ax+B進行定義，則夾直線上各測點的殘差平方和為:

根據(jù)最小二乘原理，夾直線的擬合參數(shù)為:

采用最小二乘得到擬合參數(shù)A和B。夾直線的方位角β=arctan A。從而得到前夾直線方位角β1，后夾直線方位角β2，兩夾直線的交點(JD)坐標(xJD，yJD)，以及平面曲線的偏角 α =。因而，直線段上每個測點(xc，yc)的調(diào)整量為:

3.2 平面曲線的幾何要素

加設(shè)緩和曲線的平面曲線(圖3所示)的幾何要素為:偏角α，曲線半徑R，前緩和曲線長l1，后緩和曲線長l2，切線長T，曲線總長L。切線長T和曲線總長L由公式(6)計算:

式中:p為內(nèi)移距;q為切垂距;R，l1和l2為輸入?yún)?shù)。

圖3 平面曲線Fig.3 Plan curve

3.3 緩和曲線線型計算

在《線規(guī)》中明確規(guī)定:在不超過160 km/h的客貨列車共線運行的鐵路，采用三次拋物線作為緩和曲線線型。以前緩和曲線為例，建立以ZH點為原點，切線為X'軸的緩和曲線局部坐標系?？傻?次拋物線上各點的計算公式［3］:

式中:lc為緩和曲線上任一點的弧長;C=R·l1;β'c為該點的緩和曲線角。

由于平面曲線測量是在大地坐標下進行的，需要將局部坐標系轉(zhuǎn)換到大地測量坐標系。若已知ZH點大地坐標和方位角，前緩和曲線上任一點的大地坐標計算公式:

3.4 平面曲線特征點坐標計算

由于ZH和HZ點分別位于始端、終端夾直線上。根據(jù)前夾直線方位角β1，得到ZH點的大地坐標為:

根據(jù)后夾直線方位角β2，計算HZ點的大地坐標為:

將lc=l1代入式(7)和式(8)即可得到HY點坐標(xHY，yHY)及方位角 βHY。進而計算出圓心坐標:

3.5 平面曲線特征點里程計算

將整個平面曲線測量起點設(shè)置在夾直線上，將其作為優(yōu)化曲線的起點，將其測量坐標投影到所在夾直線，得到位于夾直線上的起點坐標(xQD，yQD)，設(shè)定該點的里程LQD，可由下式計算出平面曲線特征點的里程:

在后續(xù)的平面曲線中，前一平面曲線的HZ點為后一平面曲線的起點，前一平面曲線的后夾直線為后一平面曲線的前夾直線。從而計算得到整個測量曲線連續(xù)的平面曲線的特征點大地坐標及里程。

4 平面曲線調(diào)整量計算

4.1 圓曲線調(diào)整量計算

由公式(11)計算得到了圓曲線的圓心坐標，則圓曲線上的調(diào)整量即為測量點到圓心的距離減曲線半徑:

4.2 緩和曲線調(diào)整量計算

直線和圓曲線的調(diào)整量均能通過幾何關(guān)系直接算出，在緩和曲線上則比較復雜，本文采用緩和曲線方位角搜索法進行迭代計算［7］。即在緩和曲線上找到法線通過測點的投影點，調(diào)整量即為測點與投影點的距離。

取緩和曲線段上的1個測量點C，根據(jù)C點與ZH點坐標，先計算ZH→C的方位角θZH→C，若其與ZH點的方位角的夾角＜ 0.01，則 ZH 點即為點的法線里程點。

圖4 緩和曲線調(diào)整量計算原理圖Fig.4 Adjustment calculation of transition curve

否則，計算 c點與ZH 點間距 sZH→c，令li=sZH→c× cosθsZH→c，代入式(7)和(8)，求出 E1 點處的坐標(xE1，yE1)及方位角 βE1。再計算 E1→ c的方位角θE1→c，若其與E1點的方位角的夾角＜0.01，則E1點即為c點的法線里程點。否則重復以上步驟，最后總能找到點En，使En的切線與En－c垂直。

得到En后，則該測量點的調(diào)整量為:

5 平曲曲線優(yōu)化算法

5.1 算法原理

通過建立的平面曲線調(diào)整量計算模型，在兩端夾直線不變的前提下，平面曲線的優(yōu)化問題轉(zhuǎn)換為:以R，l1和l2為優(yōu)化參數(shù)，以調(diào)整量最小為優(yōu)化目標的優(yōu)化求解問題。建立優(yōu)化方程:

Δi為每個測點的調(diào)整量。Rs，l1s和l2s為3個優(yōu)化參數(shù)的迭代步長。前、后夾直線的調(diào)整量在進行最小二乘擬合時就已確定，與R，l1和l2無關(guān)，因此，只需要對緩和曲線段和圓曲線段進行優(yōu)化。在線路設(shè)計中，平面曲線一般采用對稱基本型，即前、后緩和曲線等長，因而還可以減少一個優(yōu)化參數(shù)，提高算法效率。具體的優(yōu)化步驟如下:

(1)給定R，l1和l2初始值，根據(jù)線路養(yǎng)護管理要求設(shè)置相應步長(Rs，l1s和 l2s)，根據(jù)式(6)，計算始端切線長T1、終端的切線長T2。

(2)根據(jù)式(9)和(10)，計算ZH與 HZ點坐標。

(3)將l1代入式(7)和(8)，計算HY點坐標。

(4)將HY點坐標代入式(11)，計算得到圓曲線的圓心坐標。

(5)遍歷所有的測點，根據(jù)測點坐標值，判斷測點所在的特征段。

(6)如果測點位于緩和曲線和圓曲線區(qū)域，調(diào)用相應的調(diào)整量計算模型，計算位于緩和曲線和圓曲線區(qū)域的所有測點的調(diào)整量。

(7)按照指定步長，調(diào)整參數(shù) R，l1和l2。重復以上操作進行最小調(diào)整量搜索，直至達到各優(yōu)化參數(shù)的邊界條件，計算結(jié)束。最小調(diào)整量所處的平面曲線參數(shù)R，l1和l2即為該平面曲線的最優(yōu)參數(shù)。

5.2 算例

按照上述算法，應用C#語言開發(fā)了既有線平面曲線復測優(yōu)化軟件，并對廣深I(lǐng)線、廣深I(lǐng)I線K40～K65試驗路段進行了復測。從測量數(shù)據(jù)中取一段具備前后夾直線的完整既有線為例，對算法進行驗證。

算例一:根據(jù)工務部門提供的原有曲線參數(shù)直接進行平面曲線調(diào)整量計算，將計算結(jié)果按照測點依次繪制出調(diào)整量示意圖，如圖5(圖中上部)所示，圖中下部為測量超高。

算例二:設(shè)定曲線半徑優(yōu)化范圍為6 000～10 000 m，優(yōu)化步長取100 m，根據(jù)現(xiàn)場管理要求緩長保持不變，對測量數(shù)據(jù)進行優(yōu)化計算，優(yōu)化后的圓曲線半徑為8 300 m。圖6所示為優(yōu)化后得到的調(diào)整量圖。

圖5 R=8 000調(diào)整量圖Fig.5 Adjustment when radius is equal to 8 000

圖6 R=8 300調(diào)整量圖Fig.6 Adjustment when radius is equal to 8 300

從兩圖中可以看出，優(yōu)化后的平面曲線，從ZH點到HZ點的調(diào)整量明顯減少。優(yōu)化前后的平面曲線參數(shù)如表1所示。

表1 優(yōu)化結(jié)果對比Table 1 Optimization result comparison m

6 結(jié)論

(1)基于復測線路的大地測量坐標，建立了基于R，l1和l23個參數(shù)的平面曲線特征點大地坐標及里程的計算模型。

(2)采用緩和曲線方位角搜索法，實現(xiàn)了緩和曲線段調(diào)整量的計算，從而實現(xiàn)了整個平面曲線任一測點的調(diào)整量計算。

(3)根據(jù)以上計算模型，提出了基于夾直線最小二乘的平面曲線優(yōu)化方法，算法簡單，易于編程實現(xiàn)。平面曲線優(yōu)化參數(shù)個數(shù)、初始值、步長可以根據(jù)現(xiàn)場養(yǎng)護需求進行靈活控制。該算法不僅避免了基于圓曲線最小二乘的優(yōu)化方法出現(xiàn)的病態(tài)擬合結(jié)果，而且在直線段可加長測點間隔，提高測量效率。

根據(jù)廣深I(lǐng)線、廣深I(lǐng)I線K40～K65試驗路段的復測數(shù)據(jù)，采用本文提出的基于夾直線擬合的優(yōu)化方法，進行了試驗區(qū)段的平面曲線優(yōu)化，并根據(jù)優(yōu)化后的參數(shù)指導大機精確搗固，搗固后試驗路段總體質(zhì)量逐步提高。廣深Ⅰ線搗固后(2011－10－08)質(zhì)量指數(shù)TQI為3.96，比搗固前(7月22日)降低了0.48，廣深Ⅱ線質(zhì)量指數(shù) TQI為3.95，比搗固前降低了 0.61。

［1］李家穩(wěn)，陳 峰，張海燕，等.坐標法優(yōu)化撥距計算方法的研究［J］.北方交通大學學報，2004，28(4):34－36.LI Jia-wen，CHEN Feng，ZHANG Hai-yan，et al.Study of the move distance calculation method［J］.Journal of Northern Jiaotong University，2004，28(4):34 －36.

［2］李 偉，蒲 浩，彭先寶.基于方向加速法的鐵路既有線平面重構(gòu)優(yōu)化算法［J］.鐵道科學與工程學報，2009，6(3):47－51.LI Wei，PU Hao，PENGXian-bao.Existing railway plane Line reconstruction algorithm based on direction acceleration method［J］.Journal of Railway Science and Engineering，2009，6(3):47－51.

［3］李遠富.線路勘測設(shè)計［M］.北京:高等教育出版社，2009.LI Yuan-fu.Route survey and design［M］.Beijing:Higher Education Press，2009.

［4］李紅艷，陳治亞，邢 誠，等.鐵路既有線曲線復測計算方法［J］.中國鐵道科學，2009，30(2):18－22.LI Hong-yan，CHEN Zhi-ya，XING Cheng，et al.Calculation method for curve re－surveying of the existing railway line［J］.China Railway Science，2009，30(2):18 －22.

［5］魯永辰.坐標法復測既有鐵路曲線［J］.鐵道勘察，2006，32(3):15 －17.LU Yong-chen.Re－surveying the curves of the existing railways with the method of coordinate［J］.Railway Investigation and Surveying，2006，32(3):15 －17.

［6］王 磊.既有鐵路平面改建CAD系統(tǒng)理論與方法的研究［D］.長沙:中南大學土木建筑學院，2008.WANG Lei.Theory and method research on reconstruction CAD system of existing railway［D］.Changsha:School of Civil Engineering，Central South University，2008.

［7］李富強，陳 斐，王 強.已知點坐標換算路線里程的方法及程序［J］.山西交通科技，2005(4):65－67.LI Fu-qiang，CHEN Fei，WANG Qiang.A method and program of using point coordinate to convert route mileage［J］.Shanxi Science ＆ Technology of Communications，2005(4):65－67.