王偉才，彭成山，楊晨輝，張晨浩

（上海地礦工程勘察有限公司，上海200072）

1 引言

相對于靜力水準、電子水平尺等儀器，自動全站儀可以獲得監(jiān)測點的水平和垂直2 個方面的數據。相應的，由于在隧道內影響測角和測距精度的因素較多，其對自動全站儀監(jiān)測數據常要進行相應的處理。因此，需要對自動全站儀的測站精度進行估算與驗證。在隧道內，其氣壓、溫度、濕度與地面略有不同，在測試過程中選用了徠卡0.5s 級全站儀TM50，測距精度為0.6mm+1×10-6mm。通過對距離、高差的重復觀測，了解其測量精度的分布情況[1]。

2 精度估算

2.1 平面精度估算

試以其測角方差為0.5″，測距方差為0.8mm[2，3]。若假定距離為200m，角度為170°18′35″，那么測得其平面坐標的公式為x=x0+scosα，y=y0+ssinα（其中，x0、y0為測站坐標；s為平距；α 為測站至目標點方位角）。由于X和Y公式類似，故以下推導僅以X為例，Y方向的精度可類似推導。由于平面坐標與距離和角度是非線性關系，應用協(xié)方差傳播律求兩點間高差的方差σx2，因為觀測邊長和觀測角是相互獨立的，故其方差陣為：

從而可得其測坐標的中誤差為0.7mm，即200m 范圍內，該全站儀的估算精度在0.7mm 以內。

2.2 高差精度估算

采用強制對中觀測支架獨立測站測量的單點高差精度的推算。若測得的水平距離為D＝200m，其方差為md=0.8mm2，垂直角a=15°，其方差為0.5″，設邊長觀測值和角度觀測值是獨立觀測值，應用協(xié)方差傳播律求兩點間高差的方差，因為觀測邊長和觀測角是相互獨立的，故其方差陣為[4]：

式中，σD為平距方差；σα為角度方差，ρ 為弧度換算秒常數，取206 265。

從而可得其中誤差為0.76mm，即200m 范圍內，在不考慮大氣折光差的影響，精度是優(yōu)于0.8mm。

3 實際驗證

3.1 隧道內的大氣折光差

3.1.1 隧道內的大氣參數改正

以蘇州軌道交通1 號線為例，通過2019 年9 月（共4 次進入隧道內）對隧道內的氣壓、溫度、濕度的連續(xù)觀測，并與地面進行對比，如表1 所示。

表1 隧道內大氣參數統(tǒng)計表

TM50 全站儀，采用EDM 電子距離測量。EDM 本身是一種可見紅色激光測距，屬于1 類激光產品，其波長為785nm，最高平均輻射功率3mW，脈沖時間17ms，脈沖重復頻率小于29Hz。由TM50 可見紅色激光測距的改正公式[5，6]：

式中，ΔD1為大氣改正值，10-6m；P為氣壓，mbar（1mbar=100Pa）；T為溫度，℃；h為相對濕度，%，α=1/273.15，x=7.5T/（237.3+T）+0.7857。將溫度為24.8℃，濕度為73.6%，氣壓為1.02×105Pa，代入式（1）中進行計算得到ΔD1=10.23×10-6m。

3.1.2 幾何水準、三角高程同名點高差較差的測定

通過對隧道內大氣參數的測定，代入TM50 可見紅色激光測距的改正公式中，在200m 范圍內距離改正約2mm。具體實施過程中，由TM50 重復觀測獲得高差之差的精度較高，但通過與同點位幾何水準所得高差的比較，并不一致，如表2 所示。

表2 幾何水準與三角高程較差

表2 中，距離表示全站儀測站至各點之間的距離。水準測量按二等水準施測，嚴格控制前后視距差，采用閉合線路，觀測3 次，取平均值，并進行嚴密平差；全站儀三角高程測量，采用單向盤左、盤右觀測3 測回；通過計算，得到兩種方法較差最大2.06mm?？芍?/p>

1）在隧道內“活塞風”的作用下，設距離小于150m，垂直角小于16°，則TM50 全站儀測距精度對高差的精度影響小于0.5mm；

2）高差之差的精度影響包括垂直角和距離，測點與測站距離越小，其較差越大。因此，選擇測站時，應適當調整其測站與測點之間的距離關系；

3）因采用幾何水準測量時，嚴格控制了視距和i角等，則可認為，其較差為大氣折光差，根據表1 和表2 中數據可計算球氣差改正系數，又因監(jiān)測所獲取的成果為變形量，通過近似等距、等角的觀測，球氣差對變形成果的影響較小，對測站控制點成果有一定的影響。

3.2 全站儀精度實際驗證

測距、高差驗證。本次測試，與精度估算設置條件類似，采用1 測站盤左盤右觀測，使用ATR 自動鎖定方式，觀測14 個測回，鏡站采用大棱鏡組GPH1。經整理，其數據如表3 所示。

表3 外業(yè)實際驗證測試數據表

通過對觀測數據的整理及對距離、高差之方差的初步計算，并進行觀察得出初步評價：

1）在近170m 之內，平距的均方差非常小，最大僅0.25mm，其誤差的分布較為集中，且盤左盤右均值差僅為0.17mm，準確度亦較高。

2）因在自動化監(jiān)測中，常采用單向單面觀測，故而本次測試也采用該方法進行測量。在近170m 之內，盤左盤右的均方差較為接近，盤左為0.48mm，盤右為0.50mm，略低于二等水準常用的電子水準儀精度往返測0.3mm/km，但優(yōu)于估算精度。

3）在測得的高差數據中，盤左與盤右數據均值較差略大，為-13.72mm。該數據與實際采用人工照準方式測量較為接近。這表明，采用ATR 模式觀測的三角高差其精度較高，分布較為集中；不過，準確度上，存在一定的系統(tǒng)誤差。根據廠方建議，該值需要調校。

4）若采用半測回單向觀測，從監(jiān)測的重復測量的特點上看，并不影響監(jiān)測成果。但是，在控制網的測設中，若缺少對向觀測數據的條件下，進行半測回觀測，無法消除儀器自身的系統(tǒng)誤差，這不利于控制網的測設。

4 精度分析

4.1 分布情況

表3 中均方差的計算采用14 次重復觀測的數據，而在實際測量中，一周期觀測某一點，通常只觀測1 次。觀測的方式為單向，半測回。所以，實際作業(yè)過程中的觀測數據的精度不能采用重復觀測數據代替，故而需要在統(tǒng)一設計和整體布控的過程中進行合理的設計。在表2 中，某一點高差的觀測值是等精度的，通過重復觀測后，獲得一組觀測值，即其樣本各精度指標為：

若設，該組觀測值X～N（μ，σ2），則其均值Xˉ和σ2為該組觀測值的最大似然估計量。故而，在該組觀測值中其算術平均值即為點估計值。其樣本標準差與二階中心矩之差盡管很小，但樣本標準差相對來說，更貼近于實際情況，如圖1 所示[7]。

圖1 誤差分布及概率分布圖

然而，在實際作業(yè)中，為節(jié)約時間，提高觀測效率，以犧牲精度為代價，通常只進行單方向、半測回的觀測。在觀測的過程總是出現個別點變量較大的情況，這也反映出點估計值并不能完全代表實際作業(yè)過程中的誤差分布情況。

4.2 方差σ2 的包含區(qū)間

σ2 的無偏估計為s2，則有方差σ2的包含水平為1-α 的包含區(qū)間為：

在表2 高差樣本中，設標準差σ 的包含水平為0.955。則有，α/2=0.0225，n-1=13，計算得χ20.0225（13）=25.084 6，χ20.9775（13）=4.8917，S=0.51 代入式（2）中，得σ2包含水平為0.955 的包含區(qū)間為（0.13，0.69），即標準差的包含水平為0.955 的包含區(qū)間為（0.37，0.83）。

由于絕對值大于3 倍中誤差的偶然誤差出現的概率為0.3%，為小概率事件，通常以3 倍中誤差作為偶然誤差的極限值。σ 的包含水平取0.997 時，則α/2=0.001 5，計算得χ20.0015（13）=33.368 1，χ20.9985（13）=2.823 3，S=0.51 代入式（2）中，得σ2包含水平為0.997 的包含區(qū)間為（0.10，1.20），即標準差的包含水平為0.997 的包含區(qū)間為（0.32，1.09）。

由此可見，要滿足P（-3σ＜Δ＜+3σ）=99.7%的極限誤差，其絕對值的分布在區(qū)間（0.32，1.09）之內幾乎為必然事件。

5 結語

實際監(jiān)測過程中，每期數據觀測結束后，要對本次成果與上次成果的數據進行比較，判別其產生的次變量是由誤差引起的，或者是監(jiān)測點的位移量。這需要建立一個基礎數據庫，通過該基礎數據庫計算方差，并進行區(qū)間估計，在區(qū)間之內則可認為是由誤差引起的變量值；不在區(qū)間之間的則可認為是誤差產生的低概率事件，或者可認為是變形值。

這樣，在自動化監(jiān)測設計時，其高差的觀測不能僅觀測幾個測回，或者是幾個單方向的半測回，而是在初值采集過程中，有必要對其150m 范圍內的典型測量進行多測回的觀測，并計算中誤差，以作為判別累計位移量的基礎數據。

本文對自動全站儀的測站水平及垂直方向的觀測精度做了一定程度的探討，并通過實際驗證對隧道內的影響因素氣壓、濕度、溫度進行了測量計算。從實際測量模式出發(fā)，對觀測點可能出現的數值情況進行一定程度的計算，為了解自動全站儀測站精度提供了一個新的思路，供各同行參考。