王樹東，萬(wàn) 軍，2

(1.建設(shè)綜合勘察研究設(shè)計(jì)院有限公司，北京 100007；2.中國(guó)測(cè)繪科學(xué)研究院，北京 100830)

0 引言

全球定位系統(tǒng)(global positioning system，GPS)測(cè)量數(shù)據(jù)處理是研究全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system，GNSS)定位技術(shù)的一項(xiàng)重要內(nèi)容，其中尤為重要的是GPS基線解算。目前，高精度GPS數(shù)據(jù)處理軟件國(guó)際上比較知名的是GAMIT數(shù)據(jù)處理軟件，采用雙差觀測(cè)量解算模式，通過最小二乘算法進(jìn)行參數(shù)估計(jì)，其優(yōu)點(diǎn)是可以完全消除衛(wèi)星鐘差與接收機(jī)鐘差產(chǎn)生的影響，同時(shí)也能大幅度減弱大氣折射誤差、軌道誤差等的影響。GAMIT基線解算過程中通常要求精密的軌道信息文件作為位置基準(zhǔn)，但是在實(shí)際工作中由于最終精密星歷的播發(fā)是延后的，所以有時(shí)在一些需要實(shí)時(shí)進(jìn)行解算基線的工程中，不能滿足GAMIT基線解算的要求[1-5]。

本文通過分析星歷精度對(duì)于定位精度的影響，探討基于不同的星歷進(jìn)行GAMIT基線解算與平差計(jì)算的方法。

1 星歷誤差

1.1 星歷類型

按照精度劃分，GPS衛(wèi)星星歷可以分為精密星歷和廣播星歷。國(guó)際GNSS服務(wù)組織(International GNSS Service，IGS)負(fù)責(zé)提供各種導(dǎo)航定位系統(tǒng)的星歷，由7個(gè)IGS處理分析中心的結(jié)果進(jìn)行加權(quán)平均之后得到，主要分為最終星歷、快速衛(wèi)星歷和超快速星歷3種。廣播星歷較精密星歷的精度稍差，目前主要用于導(dǎo)航定位、低等級(jí)或短基線的工程測(cè)量[6-8]。

星歷的更新時(shí)刻為世界協(xié)調(diào)時(shí)(coordinated universal time，UTC)。各種星歷的具體情況如表1所示。

表1 星歷類型

1.2 星歷誤差對(duì)定位精度的影響

衛(wèi)星星歷誤差指衛(wèi)星星歷給出的衛(wèi)星位置與衛(wèi)星實(shí)際位置之差。衛(wèi)星在空中運(yùn)行時(shí)，會(huì)受到很多力的影響，如天體引力攝動(dòng)、太陽(yáng)輻射壓力、衛(wèi)星鐘差等。經(jīng)過理論的分析與實(shí)驗(yàn)研究，衛(wèi)星星歷誤差對(duì)相對(duì)定位的影響結(jié)果為

(1)

表2 不同星歷的相對(duì)基線誤差

2 計(jì)算流程

2.1 基線解算原理

假設(shè)在2點(diǎn)同步觀測(cè)s、k衛(wèi)星，可以組成2個(gè)單差觀測(cè)方程為：

(2)

(3)

式中：A、B表示2臺(tái)接收機(jī)；s、k表示衛(wèi)星；f為相位頻率；τ為波長(zhǎng)；ρAB為t時(shí)刻測(cè)站到衛(wèi)星的幾何距離；δAB為接收機(jī)鐘差；NAB為t時(shí)刻的相位整周未知數(shù)；ΔAB為大氣誤差。

僅考慮頻率相同的情況，即τs、τk相同時(shí)用τ表示，對(duì)式(2)及式(3)求差化簡(jiǎn)可得

(4)

由式(4)可見通過雙差觀測(cè)方程消除了衛(wèi)星鐘差與接收機(jī)鐘差。對(duì)于短距離基線來說，大氣誤差影響較小，可忽略不計(jì)。此時(shí)的未知參數(shù)為2測(cè)站的相對(duì)矢量和整周模糊度[9-10]。

2.2 解算策略設(shè)置

本文中的解算軟件為GAMIT 10.5，其需要在UNIX或LINUX系統(tǒng)中運(yùn)行，同時(shí)需要GCC編譯器的支持。軟件的具體安裝步驟可參考相關(guān)文獻(xiàn)，本文不再對(duì)軟件的安裝進(jìn)行詳細(xì)講述。

比特幣交易是相互關(guān)聯(lián)的輸入輸出交易事務(wù)，使用未經(jīng)使用的交易的輸出（UTXO）作為輸入，并生成新的輸出。同一筆交易可以有多個(gè)輸入和多個(gè)輸出，具體交易格式如表1所示。

GAMIT數(shù)據(jù)處理主要的解算策略設(shè)置，采用軌道松弛模式；LC_AUTCLN解算類型；光壓模型采用BERNE；對(duì)流層天頂距延遲參數(shù)采取2 h估計(jì)1個(gè)，采用GMF映射函數(shù)；海潮改正模型采用FES2004；衛(wèi)星截止高度角設(shè)置為10°。

2.3 解算流程

2.3.1 基線解算處理流程

1)GAMIT的參數(shù)設(shè)置。采用RELAX解模式，解的數(shù)據(jù)類型為L(zhǎng)C-AUTCLN，衛(wèi)星截止高度角為10°，對(duì)流層延遲每2 h設(shè)置1個(gè)，對(duì)基準(zhǔn)站N、E、U方向的約束都為5 cm。

2)GAMIT數(shù)據(jù)解算分步進(jìn)行處理，編制調(diào)用GAMIT有關(guān)模塊的命令。通過創(chuàng)建一個(gè)年積日文件夾，將觀測(cè)文件、廣播星歷文件以及tables中的文件復(fù)制到年積日的文件夾下[11]：

①執(zhí)行makexp，檢查要解算的RINEX格式或GAMIT認(rèn)可的X開頭的觀測(cè)文件，正常運(yùn)行后，按照提示依次輸入有關(guān)工程的各類信息和參數(shù)；

②執(zhí)行sh_sp3fit，生成g文件以及t文件；

③執(zhí)行sh_check_sess，對(duì)生成的g文件進(jìn)行檢查，以查看session.info中有對(duì)應(yīng)的衛(wèi)星；

④執(zhí)行makej，由廣播星歷文件生成鐘差文件，即參考時(shí)刻衛(wèi)星鐘的鐘偏以及鐘漂；

⑤執(zhí)行sh_check_sess，檢查J文件，以查看session.info中有對(duì)應(yīng)的衛(wèi)星；

⑦執(zhí)行fixdrv，生成接收機(jī)鐘差多項(xiàng)式系數(shù)I文件以及批處理文件(b****y.bat)和相應(yīng)的輸入文件；

⑧執(zhí)行csh b****y.bat，得到每天的計(jì)算結(jié)果文件。

2.3.2 平差處理流程

1)利用GLOBK的htoglb模塊將H文件轉(zhuǎn)換為二進(jìn)制文件。

2)將所有*.glx文件做一個(gè)列表文件，命名為*.gdl。

3)執(zhí)行g(shù)lobk單天平差計(jì)算，輸入“glred 6 glred.prt glred.log *.gdl globk_comb.cmd”，得到單天平差結(jié)果文件。

2.4 解算結(jié)果分析

單時(shí)段解算的標(biāo)準(zhǔn)化均方根殘差(normalized root mean square value，NRMS)是判斷GAMIT基線解算結(jié)果質(zhì)量好壞的一個(gè)重要指標(biāo)。通常來說，NRMS值越小，基線解算精度越高；反之，精度越低。依據(jù)Q文件中的NRMS值對(duì)其結(jié)果進(jìn)行衡量分析，比較理想的NRMS值應(yīng)小于0.25；若大于0.5，說明處理過程中有可能存在未消除的大周跳或某一參數(shù)的解存在較大偏差；查看O文件中的基線解算結(jié)果和精度，從而進(jìn)一步對(duì)基線解算結(jié)果進(jìn)行分析[12-15]。

3 算例及精度分析

3.1 算例

選用2012-10-02(年積日第289天)福建省7個(gè)連續(xù)運(yùn)行參考站(continuously operating reference stations，CORS)的原始觀測(cè)數(shù)據(jù)進(jìn)行解算，觀測(cè)時(shí)間為24 h，歷元間隔30 s。7個(gè)CORS站點(diǎn)分布如圖1所示。

用7種星歷對(duì)CORS站的數(shù)據(jù)進(jìn)行了解算，并從解算完的基線中選擇了14條，具體情況如表3所示。

表3 GAMIT解算實(shí)驗(yàn)基線

3.2 精度分析

通常只要NRMS值符合要求，則可認(rèn)為基線解算結(jié)果比較可靠。表4為不同的星歷解算得到的NRMS值。從表中可以看到，幾種星歷解算的NRMS值都小于0.25，并且與最終精密星歷解算的NRMS值之差僅為0.012，說明由快速星歷和超快速星歷得到的解算結(jié)果是合理的。

表4 7種星歷解算驗(yàn)后NRMS值

為進(jìn)一步比較解算結(jié)果的精度，對(duì)不同星歷平差完成的基線坐標(biāo)進(jìn)行比較分析。解算結(jié)果以利用最終精密星歷進(jìn)行解算的結(jié)果為基準(zhǔn)，將其他星歷解算得到的基線坐標(biāo)分量與之求差，得到基線解分別在北方向(N)、東方向(E)、天向(U)3個(gè)方向的差值，如圖2、圖3、圖4所示。

從圖2～圖4中可以看出由快速星歷解算得到的結(jié)果和最終精密星歷解算得到的結(jié)果在3個(gè)方向上的差值都為零。這說明在本次基線解算實(shí)驗(yàn)中，快速精密星歷和最終精密星歷基線解算的測(cè)站坐標(biāo)結(jié)果一致。超快速精密星歷解算得到的結(jié)果與最終精密星歷解算得的結(jié)果在3個(gè)方向上存在差異；當(dāng)基線長(zhǎng)度在370 km以內(nèi)時(shí)，在N、E方向上差值基本在0.5 mm以內(nèi)；當(dāng)基線長(zhǎng)度達(dá)到432和494 km時(shí)，差值在1～2 mm之間，且呈現(xiàn)隨著基線增長(zhǎng)差值逐漸增大的趨勢(shì)。

4 結(jié)束語(yǔ)

本文使用不同星歷對(duì)實(shí)測(cè)CORS站數(shù)據(jù)進(jìn)行基線解算，通過這次基線解算實(shí)驗(yàn)以及工作中大量的基線解算經(jīng)驗(yàn)，驗(yàn)證了在特定情況下，使用其他2種星歷代替最終精密星歷進(jìn)行基線解算的可行性與準(zhǔn)確性，且得出：當(dāng)在無(wú)法獲取最終精密星歷但可以獲取快速精密星歷的的特殊情況下，可以使用快速精密星歷代替最終精密星歷進(jìn)行高精度的基線解算與平差處理；當(dāng)在無(wú)法獲取上述2種星歷但可獲取超快速星歷的特殊情況下，在處理小于400 km以內(nèi)的中短基線數(shù)據(jù)時(shí)，可以使用超快速精密星歷代替精密星歷進(jìn)行高精度的基線解算與平差處理。

[1] HERRING T A，KING R W，MCCLUSKY S C.GAMIT reference manual GPS analysis at MIT(Release 10 3)[EB/OL].(2015-06-16)[2017-01-19].http：//www.doc88.com/p-817688430780.html.

[2] 李征航，張小紅.衛(wèi)星導(dǎo)航定位技術(shù)及高精度數(shù)據(jù)處理方法[M].武漢：武漢大學(xué)出版社，2010.

[3] 黨亞民，秘金鐘，成英燕.全球?qū)Ш叫l(wèi)星系統(tǒng)原理與應(yīng)用[M].北京：測(cè)繪出版社，2007.

[4] 劉大杰，施一民，過靜珺.全球定位系統(tǒng)(GPS)的原理與數(shù)據(jù)處理[M].武漢：武漢大學(xué)出版社，2005.

[5] 張勤，李加權(quán).GPS測(cè)量原理及應(yīng)用[M].北京：科學(xué)出版社，2005.

[6] 李征航，黃勁松.GPS測(cè)量與數(shù)據(jù)處理[M].武漢：武漢大學(xué)出版社，2010.

[7] 高旺，高成發(fā)，潘樹國(guó)，等.基于廣播星歷的GAMIT基線解算方法及精度分析[J].測(cè)繪科學(xué)，2014，39(8)：54-57.

[8] 曹炳強(qiáng)，成英燕，許長(zhǎng)輝.海潮模型對(duì)連續(xù)跟蹤參考站數(shù)據(jù)解算的影響：以福建省為例[J].測(cè)繪科學(xué)，2015，40(12)：108-111.

[9] 黃功文，王斌，王延偉.高精度GPS控制網(wǎng)基線重復(fù)性檢驗(yàn)與質(zhì)量分析[J].測(cè)繪通報(bào)，2011(7)：9-11.

[10] 王霞迎.BDS-GPS基線解算算法研究[D].北京：中國(guó)測(cè)繪科學(xué)研究院，2014.

[11] 王武星，顧國(guó)華.利用基準(zhǔn)站觀測(cè)資料確定GPS衛(wèi)星軌道[J].大地測(cè)量與地球動(dòng)力學(xué)，2003，23(3)：112-115.

[12] 張雙成，曹海洋，高涵，等.基于GAMIT的GPS短基線解類型分析及應(yīng)用[J].測(cè)繪通報(bào)，2011(10)：27-29.

[13] 趙建三，楊創(chuàng)，聞德保.利用GAMIT進(jìn)行高精度GPS基線解算的方法及精度分析[J].測(cè)繪通報(bào)，2011(8)：5-8.

[14] 曹炳強(qiáng)，成英燕，許長(zhǎng)輝，等.間距分區(qū)法在解算衛(wèi)星連續(xù)運(yùn)行站數(shù)據(jù)中的應(yīng)用[J].測(cè)繪通報(bào)，2016(11)：15-17.

[15] 萬(wàn)軍，成英燕，黨亞民，等.陸態(tài)網(wǎng)測(cè)站高精度解算分區(qū)方案[J].測(cè)繪科學(xué)，2016，40(4)：1-4.