付宏平，郭際明，張正祿

(1.中鐵第一勘察設(shè)計院集團(tuán)有限公司，陜西西安710043;2.武漢大學(xué)，湖北武漢430079)

一、引 言

隧道施工的特點要么是從兩個洞口相向開挖，要么是從幾個洞口(包括橫洞、斜井、豎井、平行導(dǎo)坑)開始，深入到線路中線位置后，再沿線路中線相向開挖，即目前普遍采用的“長隧短打”的施工方法。無論哪種施工方法，隧道貫通誤差均是評估隧道施工建設(shè)質(zhì)量的一個重要指標(biāo)，尤其是橫向貫通誤差和高程貫通誤差。

目前，國內(nèi)特長隧道工程已經(jīng)很多，如已建成通車的石太客專太行山隧道，全長27.839 km，是目前國內(nèi)建成通車的最長鐵路隧道。正在修建的蘭渝線西秦嶺特長隧道是目前國內(nèi)第二長大鐵路山嶺隧道，隧道全長28.2 km;西格二線關(guān)角特長隧道是目前國內(nèi)最長的鐵路山嶺隧道，隧道全長32.645 km，其貫通誤差特別值得關(guān)注。

現(xiàn)有《高速鐵路工程測量規(guī)范》(TB 10601—2009)中對特長隧道(20 km及以上)貫通誤差未作明確規(guī)定，而20 km及以上特長隧道已越來越多。如何進(jìn)行特長隧道控制網(wǎng)的設(shè)計，如何預(yù)計其貫通誤差影響值的大小，需要研究出一個切實可行、操作簡便的方法。

二、特長隧道貫通誤差預(yù)計方法研究

1.經(jīng)典方法

以前見諸各類文獻(xiàn)的隧道貫通誤差預(yù)計傳統(tǒng)方法(經(jīng)典方法)，主要包括導(dǎo)線法、誤差來源分析法、權(quán)函數(shù)法。

(1)導(dǎo)線法

導(dǎo)線法是洞外、洞內(nèi)控制測量所引起隧道橫向貫通誤差的一種近似預(yù)計方法。過去洞外平面控制采用地面三角形網(wǎng)、導(dǎo)線網(wǎng)，而洞內(nèi)只能用導(dǎo)線或?qū)Ь€網(wǎng)，此法是選擇最靠近隧道中線的一條線路(如圖1中的J-1-2-3-4-C)作為導(dǎo)線，用下述導(dǎo)線公式預(yù)計橫向貫通誤差的影響值

式中，mβ為測角中誤差為最弱邊相對中誤差;為導(dǎo)線點至貫通面垂直距離的平方和;為導(dǎo)線邊在貫通面上投影長度的平方和。

圖1 導(dǎo)線法計算貫通誤差影響值示意圖

由于洞外控制測量不再采用地面三角形網(wǎng)、導(dǎo)線網(wǎng)，該方法對洞外控制網(wǎng)不再適用，但仍可用于預(yù)計洞內(nèi)控制測量所引起的隧道橫向貫通誤差。該方法計算簡單，未考慮平差后的精度增益，預(yù)計的值偏大、偏安全。

(2)誤差來源分析法

《鐵路工程測量規(guī)范》(TB 10101—2009)及《高速鐵路工程測量規(guī)范》(TB 10601—2009)中洞外GPS平面網(wǎng)測量誤差對隧道橫向貫通的影響，采用下式預(yù)計

式中，mJ、mC為進(jìn)、出口GPS控制點的Y坐標(biāo)誤差;LJ、LC為進(jìn)、出口GPS控制點至貫通點的長度;mαJ、mαc為進(jìn)、出口GPS聯(lián)系邊的方位中誤差;θ、φ為進(jìn)、出口控制點至貫通點連線與貫通點線路切線的夾角。

后兩項也可以由下式算得

式中，mG為GPS測量定向聯(lián)系邊方向誤差，單位為(″)，為隧道設(shè)計時的先驗值，ρ=206 265。

該方法使用了隧道兩端洞口定向點點位誤差和定向邊坐標(biāo)方位角中誤差，這些值可以事先設(shè)計。因此，該方法適合在控制測量開始前對隧道貫通誤差進(jìn)行預(yù)計。

(3)權(quán)函數(shù)法

權(quán)函數(shù)法是隧道貫通誤差預(yù)計嚴(yán)密的公式，按照隧道洞外控制網(wǎng)的驗后單位權(quán)中誤差和其坐標(biāo)協(xié)因數(shù)陣來計算隧道貫通中誤差，如圖2所示。

圖2 權(quán)函數(shù)法計算貫通誤差影響值示意圖

由平差理論，按求未知數(shù)的函數(shù)權(quán)倒數(shù)方法，在不考慮 SJG、SCG、βJ、βC的誤差時，對上兩式微分并整理得到關(guān)于ΔXG、ΔYG的權(quán)函數(shù)式(推導(dǎo)略)。寫成矩陣形式，按協(xié)因數(shù)傳播律得

對于直線隧道，當(dāng)隧道軸線為X軸時，則σΔXG為縱向貫通誤差的影響值，σΔYG為橫向貫通誤差的響值。對于曲線隧道，可以根據(jù)線路關(guān)系計算出貫通點在控制網(wǎng)坐標(biāo)系統(tǒng)中的坐標(biāo)(XG，YG)，以及該點的切線方位角(貫通點的縱向方位角)αE和橫向方位角αF=αE+90°，因此，隧道橫向貫通誤差為

該方法是一種嚴(yán)密的隧道貫通誤差預(yù)計方法。由于控制網(wǎng)的坐標(biāo)協(xié)因數(shù)陣在測量前無法構(gòu)建獲得，因此，該方法適用于控制測量施測后的隧道貫通誤差驗后預(yù)計。

2.近似方法

相關(guān)文獻(xiàn)報道過一些近似預(yù)計方法，主要包括簡化公式法、平均相對誤差估計法、設(shè)計尺寸法等，各方法分述如下。

(1)簡化公式法

簡化公式的實質(zhì)是:洞內(nèi)導(dǎo)線測量誤差所引起的隧道橫向貫通誤差完全由測角誤差所引起，測角誤差與邊長成正比，且與平均邊長有關(guān)，平均邊長越長，則導(dǎo)線測量誤差所引起的隧道橫向貫通誤差越小。洞內(nèi)導(dǎo)線測角誤差引起橫向貫通誤差可近似表示為

式中，S為導(dǎo)線的平均邊長，單位為m;n為導(dǎo)線的邊數(shù)。

該方法計算更為簡便，對于直伸隧道的洞內(nèi)導(dǎo)線來說，有參考價值。

(2)平均相對誤差估計法

按照隧道長度，用GPS控制網(wǎng)平均相對誤差1/T估計隧道貫通誤差影響值，即

式中，L為隧道長度，單位為m;1/T為GPS控制網(wǎng)平均相對誤差。

歐洲隧道公司橫穿英法海峽的歐洲隧道采用了該方法。該方法過于簡單，沒有明確區(qū)分橫向和縱向貫通誤差影響值。所謂控制網(wǎng)平均相對誤差，實際上只考慮了邊長的相對誤差，且控制網(wǎng)的邊長非常懸殊。邊長的相對誤差與邊長有關(guān)，邊長越短，邊長的相對誤差越大;邊長越長，邊長的相對誤差越小。

(3)設(shè)計尺寸法

按照隧道設(shè)計尺寸確定橫向貫通誤差允許值。有文獻(xiàn)提出:對于圓形隧道可取設(shè)計直徑的1/7～1/10作為橫向貫通誤差的允許值，對于方形隧道可取洞寬1/7～1/10作為橫向貫通誤差的允許值。這是基于當(dāng)出現(xiàn)這種橫向貫通誤差時，完全可以在貫通面附近進(jìn)行調(diào)整，既不會增加開挖量，也不會影響工程的質(zhì)量和進(jìn)度來考慮的。該方法很簡單，且與隧道的長度無關(guān)。

3.本文提出的隧道貫通誤差預(yù)計方法

(1)坐標(biāo)差統(tǒng)計法

坐標(biāo)差統(tǒng)計法原理如下:取隧道軸線為X坐標(biāo)軸，貫通面與坐標(biāo)軸垂直，貫通面與Y軸平行。對特長直伸隧道，設(shè)只有一個貫通面，且位于中間，按等邊設(shè)計并模擬從進(jìn)、出口到貫通面的洞內(nèi)導(dǎo)線網(wǎng)，如設(shè)計成六邊形的雙導(dǎo)線網(wǎng)。對進(jìn)、出口洞內(nèi)導(dǎo)線網(wǎng)分別進(jìn)行平差，計算出位于貫通面上同一點的X、Y坐標(biāo)，Y坐標(biāo)差即為洞內(nèi)導(dǎo)線測量誤差所引起的隧道橫向貫通誤差。模擬計算幾十組(大子樣)，然后對Y坐標(biāo)差值作統(tǒng)計計算，即計算Y坐標(biāo)差的均值、均值的中誤差。最后可將均值作為洞內(nèi)導(dǎo)線網(wǎng)測量誤差所引起的橫向貫通誤差。

該方法的特點是計算嚴(yán)密，使用簡便，易于編程實現(xiàn)，子樣越大，所得的結(jié)果越可靠;而且將進(jìn)、出口洞內(nèi)導(dǎo)線網(wǎng)作為一個獨立的影響因子一起預(yù)計，方法創(chuàng)新，適合各種隧道洞內(nèi)導(dǎo)線橫向和縱向貫通誤差影響值計算。

(2)坐標(biāo)中誤差法

坐標(biāo)差中誤差法原理如下:對于特長隧道，按前述坐標(biāo)系設(shè)計進(jìn)、出口端洞內(nèi)導(dǎo)線網(wǎng)，分別進(jìn)行平差，可得到不同隧道長貫通面上點的X、Y坐標(biāo)及其中誤差，可以把X、Y坐標(biāo)的中誤差視為洞內(nèi)導(dǎo)線測量誤差所引起的隧道縱、橫向貫通誤差，因為進(jìn)出口端各有一個網(wǎng)，坐標(biāo)的中誤差乘以，即為洞內(nèi)導(dǎo)線測量誤差所引起的隧道橫向貫通誤差，又因為平差是一次性的結(jié)果，在實際作業(yè)中，要對洞內(nèi)導(dǎo)線網(wǎng)作經(jīng)常性的重復(fù)測量，可以乘以一個小于1的系數(shù)，如除以，則剛好抵消，因此，可直接把進(jìn)口端洞內(nèi)導(dǎo)線網(wǎng)不同隧道長貫通面上點的Y坐標(biāo)中誤差視為洞內(nèi)導(dǎo)線測量誤差所引起的隧道橫向貫通誤差。

該方法計算非常簡單，有理論依據(jù)，可直接把進(jìn)口端洞內(nèi)導(dǎo)線網(wǎng)不同隧道長貫通面上點的Y坐標(biāo)中誤差視為洞內(nèi)導(dǎo)線測量誤差所引起的隧道橫向貫通誤差，也適合洞內(nèi)導(dǎo)線網(wǎng)邊加測陀螺方位角后的情況。

(3)微分公式法

微分公式法是在過去的隧道貫通誤差嚴(yán)密預(yù)計方法即權(quán)函數(shù)法的基礎(chǔ)上，進(jìn)行改進(jìn)和擴(kuò)展，提出的一種適合特長隧道橫向貫通誤差預(yù)計的方法，可擴(kuò)展到GPS網(wǎng)。其原理如下:設(shè)隧道洞外平面控制網(wǎng)及坐標(biāo)軸、隧道軸線如圖1所示，從進(jìn)口點J和定向點A，經(jīng)洞內(nèi)導(dǎo)線推算得貫通點PJ的坐標(biāo)，從出口點C和定向點B，經(jīng)洞內(nèi)導(dǎo)線推算得貫通點PC的坐標(biāo)，不考慮 βJ、βC、SJP、SCP的誤差，隧道洞外平面控制測量誤差所引起的橫向貫通誤差即貫通點PC、PJ的橫坐標(biāo)差 ΔyP的中誤差，對 ΔyP進(jìn)行全微分，可得

式中，方向系數(shù)aJA、bJA可按下式計算

由網(wǎng)的近似坐標(biāo)不難求得上述全微分式的系數(shù)，根據(jù)廣義誤差傳播律計算出隧道洞外平面控制測量誤差所引起的橫向貫通誤差。

由上述公式可知，隧道洞外平面控制測量誤差所引起的橫向貫通誤差與隧道兩相向開挖洞口的洞口點、定向點的位置及其精度有關(guān)，與貫通點的位置有關(guān)。對于特長隧道，無論有多少個貫通面，都可按式(9)計算出隧道洞外平面控制測量誤差所引起的橫向貫通誤差。

該方法理論嚴(yán)密，可用于洞外邊角網(wǎng)(對于不能布設(shè)GPS網(wǎng)時有意義)，也可擴(kuò)展到洞外GPS網(wǎng)。對于GPS網(wǎng)，可以看成是邊角全測的全邊角網(wǎng)，或邊長、方位角全測的網(wǎng)，因此，可用上述全微分式預(yù)計GPS網(wǎng)的貫通誤差，易于編程實現(xiàn)。

三、隧道橫向貫通誤差預(yù)計實現(xiàn)步驟

1.洞外GPS控制網(wǎng)

隧道洞外橫向貫通誤差預(yù)計采用兩種不同的形式，分邊角全測和邊長方位角全測兩種情況，邊角全測相當(dāng)于用全站儀(TPS)測量，邊長方位角全測相當(dāng)于用GPS測量。以某特長隧道(全長20.05 km)實測網(wǎng)為例，該隧道采用“長隧短打”的形式施工，其中左線有7個貫通面，右線有9個貫通面，采用2種模擬法計算出每個短隧道的洞外貫通誤差影響值(見表1)。

(1)按GPS網(wǎng)進(jìn)行估算

1)人工生成隧道洞外GPS網(wǎng)方案文件(* .GFA)。

2)人工生成隧道洞外GPS網(wǎng)貫通誤差引導(dǎo)文件(*GPS.GTI)。

3)自動生成洞外GPS網(wǎng)的平差結(jié)果文件。

4)估算的貫通誤差影響值。

表1 隧道各貫通面的貫通誤差影響值(GPS網(wǎng)法)

(2)按邊角網(wǎng)進(jìn)行估算

1)人工生成邊角網(wǎng)方案文件(*.FA2)。

2)自動生成邊角網(wǎng)觀測值文件(*.IN2)。

3)人工生成的隧道洞外GPS網(wǎng)貫通誤差引導(dǎo)文件(* .GTI)。

4)自動生成洞外GPS網(wǎng)的平差結(jié)果文件平差結(jié)果文件。

5)估算的貫通誤差影響值(見表2)。

表2 隧道各貫通面貫通誤差影響值(TPS網(wǎng)法-全邊角網(wǎng)法)

GPS網(wǎng)法和TPS網(wǎng)法-全邊角網(wǎng)法兩種方法的計算結(jié)果基本是一致的，相差都在2 mm以內(nèi)。

2.洞內(nèi)導(dǎo)線網(wǎng)

設(shè)隧道長50 km，平均邊長500 m，按方向中誤差為1.0″(測角中誤差1.4″)、邊長固定誤差為2 mm、比例誤差為2 mm進(jìn)行導(dǎo)線網(wǎng)模擬，并對進(jìn)、出口洞內(nèi)導(dǎo)線網(wǎng)分別進(jìn)行平差，計算出位于貫通面上同一點的X、Y坐標(biāo)。估算步驟如下:

1)人工生成50 km隧道洞內(nèi)進(jìn)(或出)口導(dǎo)線網(wǎng)方案文件(*.FA2)。

2)自動生成50 km隧道洞內(nèi)進(jìn)(或出)口導(dǎo)線網(wǎng)觀測值文件(*.IN2)。

3)平差得50 km隧道洞內(nèi)進(jìn)(或出)口導(dǎo)線網(wǎng)平差結(jié)果文件(*.OU2)。

4)由平差坐標(biāo)及其精度可計算進(jìn)出口洞內(nèi)導(dǎo)線網(wǎng)在貫通面上點的坐標(biāo)差，同樣可獲得洞內(nèi)導(dǎo)線網(wǎng)點的Y坐標(biāo)中誤差，把其視為洞內(nèi)導(dǎo)線測量誤差所引起的隧道橫向貫通誤差(見表3)。

表3 洞內(nèi)導(dǎo)線網(wǎng)測量誤差引起的橫向貫通誤差估值及中誤差 mm

［1］付宏平.烏鞘嶺特長隧道控制網(wǎng)布設(shè)以及貫通誤差預(yù)計［J］.蘭州交通大學(xué)學(xué)報，2008，27(5):54-57.

［2］張正祿，張松林，伍志剛，等.20～50 km超長隧道(洞)橫向貫通誤差允許值研究［J］.測繪學(xué)報，33(1):83-88.

［3］張項鐸，張正祿.隧道工程測量［M］.北京:測繪出版社，1998:44-66.

［4］TB 10601—2009，高速鐵路工程測量規(guī)范［S］.北京:中國鐵道出版社，2010.

［5］張項鐸.淺論隧道橫向貫通誤差估計方法［J］.隧道建設(shè)，2003(5):5-7，10.

［6］劉成龍，王鵬.基于坐標(biāo)協(xié)因數(shù)陣的GPS網(wǎng)隧道橫向貫通中誤差嚴(yán)密計算方法研究［J］.鐵道學(xué)報，2009(1):74-77.

［7］張懷良.隧道貫通誤差預(yù)計方法探討［J］.山西建筑，2004(20):230-231.

［8］張項鐸.隧道GPS網(wǎng)布測中有關(guān)問題的探討［J］.中國鐵道科學(xué)，2001，22(6):73-78.

［9］孫家兵，薛利華.GPS控制網(wǎng)隧道貫通影響值計算編程實現(xiàn)要點［J］.大壩與安全，2004(S0):28-30.

［10］吳連雄，徐一鳴，劉成龍.羅峰隧道施工GPS控制網(wǎng)建網(wǎng)及其貫通誤差預(yù)計［J］.黑龍江工程學(xué)院學(xué)報，2002，16(4):36-38.