陳振國(guó)，唐龍江

(1.江蘇省核工業(yè)二七二地質(zhì)大隊(duì)，江蘇 南京 210000；2.核工業(yè)建設(shè)集團(tuán)有限公司，江蘇 南京 210000；3.遼寧工程技術(shù)大學(xué) 測(cè)繪與地理科學(xué)學(xué)院，遼寧 阜新 123000)

BDS/GPS組合天頂對(duì)流層延遲的精度分析

陳振國(guó)1,2，唐龍江3

(1.江蘇省核工業(yè)二七二地質(zhì)大隊(duì)，江蘇 南京 210000；2.核工業(yè)建設(shè)集團(tuán)有限公司，江蘇 南京 210000；3.遼寧工程技術(shù)大學(xué) 測(cè)繪與地理科學(xué)學(xué)院，遼寧 阜新 123000)

在現(xiàn)有的精密軌道和鐘差條件下，選取8個(gè)MGEX跟蹤站2014年6—9月的觀測(cè)數(shù)據(jù)，詳細(xì)分析利用BDS/GPS組合PPP法在未固定跟蹤站坐標(biāo)和固定跟蹤站坐標(biāo)情況下估計(jì)ZTD的效果，并與IGS提供的對(duì)流層產(chǎn)品對(duì)比分析。實(shí)驗(yàn)分析表明，利用PPP法估計(jì)ZTD，BDS ZTD現(xiàn)階段的STD優(yōu)于34 mm，GPS ZTD與BDS/GPS組合現(xiàn)階段的STD相當(dāng)，均優(yōu)于14 mm。與未固定跟蹤站情形下估計(jì)的BDS ZTD相比，固定跟蹤站坐標(biāo)的方式雖然可以提高利用BDS估計(jì)ZTD的穩(wěn)定性，但不能提高精度。

BDS/GPS組合；IGS對(duì)流層產(chǎn)品；PPP法；MEGX跟蹤站；ZTD

利用GNSS估計(jì)天頂對(duì)流層延遲(zenith tropospheric delays，ZTD)，在此基礎(chǔ)上進(jìn)一步獲得大氣含水量。利用GNSS估計(jì)天頂對(duì)流層延遲的方法主要有雙差網(wǎng)解法和非差法。與非差法相比，雖然差分法是目前ZTD估計(jì)中常用的方法[1-2]，但對(duì)于長(zhǎng)距離觀測(cè)值的絕對(duì)ZTD估計(jì)需要引入與測(cè)站相距大于500 km的參考站[3-4]，而且為了取得預(yù)期的精度，必須有一定的時(shí)間間隔，一般為30 min[5]。隨著誤差改正模型的精細(xì)化和整周模糊度算法的改進(jìn)，利用非差法估計(jì)的ZTD精度與利用差分法估計(jì)的一致，且非差法更具有靈活性，因此利用非差法估計(jì)ZTD的應(yīng)用也越來(lái)越廣泛[6]。利用非差法估計(jì)ZTD又可分為非差網(wǎng)解法和精密單點(diǎn)定位法(Precise Point Positioning，PPP)。對(duì)于小區(qū)域跟蹤網(wǎng)，非差網(wǎng)解法仍然需要引入遠(yuǎn)距離測(cè)站(500～2 000 km)，以獲得絕對(duì)延遲。PPP法是指利用單站非差觀測(cè)值估計(jì)ZTD，可直接獲得絕對(duì)延遲，處理大規(guī)模數(shù)據(jù)速度快[7-8]。相關(guān)研究表明，利用PPP法估計(jì)的GPS ZTD精度優(yōu)于6 mm[9]，利用GPS/GLONASS組合估計(jì)的ZTD精度為1 cm左右[10]。

雖然國(guó)內(nèi)外學(xué)者在GNSS ZTD估計(jì)方面做了大量的研究工作，但針對(duì)BDS ZTD或BDS/GPS組合ZTD的估計(jì)研究較少。范磊等人利用中科院測(cè)地所計(jì)算的精密軌道和鐘差估計(jì)BDS/GPS組合 ZTD，以CODE提供的精密軌道和鐘差產(chǎn)品估計(jì)的GPS ZTD為參考值，BDS ZTD的估計(jì)精度為2 cm，BDS/GPS組合ZTD的估計(jì)精度與GPS相當(dāng)[11-12]。徐愛(ài)功等人以IGS精密軌道和鐘差估計(jì)的GPS ZTD為參考值，分別采用網(wǎng)解法和PPP法估計(jì)BDS ZTD，其估計(jì)的偏差和標(biāo)準(zhǔn)差分別為2 mm和5 mm[13]。已有研究均是利用自身解算的精密軌道和鐘差估計(jì)BDS ZTD。武漢大學(xué)(WHU)和德國(guó)地學(xué)中心(GFZ)在2014年3月發(fā)布了BDS系統(tǒng)的精密軌道和鐘差產(chǎn)品。為評(píng)估利用現(xiàn)階段提供的精密軌道和鐘差產(chǎn)品估計(jì)BDS/GPS ZTD的效果，分別對(duì)BDS、GPS、BDS/GPS組合在未固定跟蹤站坐標(biāo)和固定跟蹤站坐標(biāo)的情況下估計(jì)的ZTD進(jìn)行研究。

1 BDS/GPS組合PPP法估計(jì)ZTD方法

1.1 BDS/GPS組合估計(jì)ZTD的數(shù)學(xué)模型

在BDS/GPS組合PPP法估計(jì)ZTD的過(guò)程中，采用雙頻觀測(cè)值消除電離層的影響，使用精密軌道和精密鐘差來(lái)固定衛(wèi)星軌道和鐘差，同時(shí)顧及GPS系統(tǒng)和BDS系統(tǒng)之間的時(shí)空基準(zhǔn)不統(tǒng)一，引入兩類接收機(jī)鐘差，BDS/GPS組合PPP法的觀測(cè)模型表示為[14]

(1)

1.2 ZTD參數(shù)估計(jì)與誤差處理

采用WHU提供的精密軌道和鐘差產(chǎn)品，ZTD的采樣間隔為5 min，偽距先驗(yàn)方差為1 m，相位先驗(yàn)方差0.01周，數(shù)據(jù)處理流程如圖1所示。

圖1 PPP法估計(jì)ZTD流程

圖1中，數(shù)據(jù)預(yù)處理是獲得高精度參數(shù)估值的關(guān)鍵。本研究采用單站TurboEdit方法進(jìn)行數(shù)據(jù)預(yù)處理，盡可能發(fā)現(xiàn)和修復(fù)周跳，剔除異常觀測(cè)值。TurboEdit方法根據(jù)導(dǎo)航衛(wèi)星頻率自動(dòng)選擇閾值，進(jìn)行周跳探測(cè)、粗差剔除以及短弧分析。其中主要包括M-W組合法和電離層殘差法等。數(shù)據(jù)預(yù)處理后盡可能考慮多的誤差改正項(xiàng)，未能精確模型化的誤差通過(guò)參數(shù)估計(jì)吸收。采用雙頻消電離層組合觀測(cè)值消除電離層一階影響；考慮天線相位纏繞改正；接收機(jī)噪聲為白噪聲；模糊度參數(shù)對(duì)未能精確模型化的誤差有補(bǔ)償作用，所以采用實(shí)數(shù)解；采用IERS Conventions 2010標(biāo)準(zhǔn)改正地球固體潮偏差，海洋潮以及地球極潮；采用絕對(duì)相位模型改正天線相位偏差。

在利用PPP估計(jì)ZTD過(guò)程中，采用GMF投影函數(shù)加Saastamoinen模型作為對(duì)流層延遲初值，氣象數(shù)據(jù)從GPT(Global Pressure and Temperature model)獲得，剩余殘差采用隨機(jī)游走過(guò)程模擬。為了能在最小二乘估計(jì)中實(shí)現(xiàn)隨機(jī)過(guò)程估計(jì)的方法，采用同時(shí)估計(jì)狀態(tài)參數(shù)和確定性參數(shù)的最小二乘遞推算法，即首先把狀態(tài)方程及狀態(tài)參數(shù)的先驗(yàn)信息轉(zhuǎn)化為相應(yīng)的等價(jià)觀測(cè)方程；然后，在組成法方程的過(guò)程中，逐步消去一些狀態(tài)參數(shù)[15]。

2 數(shù)據(jù)處理及結(jié)果分析

為了評(píng)價(jià)利用BDS/GPS組合PPP法估計(jì)ZTD的效果，采用BDS，GPS和BDS/GPS組合PPP法在未固定跟蹤站坐標(biāo)和固定跟蹤站坐標(biāo)的情況下估計(jì)ZTD。實(shí)驗(yàn)選取8個(gè)MEGX跟蹤站2014年第152～272年積日(6—9月)的BDS/GPS觀測(cè)數(shù)據(jù)，數(shù)據(jù)采樣間隔為30°。為了評(píng)估結(jié)果的精度，以IGS提供的對(duì)流層產(chǎn)品(簡(jiǎn)稱IGS ZTD)作為參考值，將PPP法估計(jì)出的ZTD與其比較并統(tǒng)計(jì)偏差MEAN和STD。

為了分析利用上述策略估計(jì)ZTD的效果，以日本的GMSD跟蹤站為例，在未固定跟蹤站的情況下，使用BDS、GPS和BDS/GPS組合PPP法處理第218天觀測(cè)數(shù)據(jù)得到ZTD，如圖2所示。從圖2中可看出，由GPS或BDS/GPS組合估計(jì)的ZTD與IGS ZTD有較好的一致性，且ZTD估計(jì)相近，不存在明顯的系統(tǒng)性偏差；由BDS估計(jì)的ZTD與IGS ZTD變化趨勢(shì)一致，但在大部分時(shí)段存在一定的系統(tǒng)性偏差。

圖2 GMSD跟蹤站的ZTD

為了進(jìn)一步分析在未固定跟蹤站的情況下采用BDS、GPS和BDS/GPS組合PPP法估計(jì)ZTD的效果，統(tǒng)計(jì)了8個(gè)跟蹤站與IGS ZTD的4個(gè)月的差值時(shí)間序列，并按天統(tǒng)計(jì)差值的均值，如圖3所示。圖3表明，利用BDS估計(jì)的ZTD與IGS ZTD的差值浮動(dòng)范圍為-70～+50 mm， 且每個(gè)跟蹤站與其對(duì)應(yīng)IGS ZTD的差值存在規(guī)律性偏差；利用GPS或BDS/GPS組合估計(jì)的ZTD與IGS ZTD的差值浮動(dòng)范圍為-10～+10 mm，無(wú)明顯的規(guī)律性偏差。通過(guò)長(zhǎng)期的觀測(cè)數(shù)據(jù)表明，利用BDS估計(jì)的ZTD不穩(wěn)定，波動(dòng)幅度較大，這主要是受到BDS衛(wèi)星少和跟蹤站的可觀測(cè)衛(wèi)星的幾何構(gòu)型影響，中緯度跟蹤站(CUT0、GMSD等)的平均可觀測(cè)衛(wèi)星數(shù)為8顆以上，主要分布在跟蹤站一側(cè)，幾何構(gòu)型差，而其他區(qū)域的一些跟蹤站(BRST、VILL等)可觀測(cè)衛(wèi)星較少，平均可觀測(cè)衛(wèi)星數(shù)為4～6顆；利用GPS或BDS/GPS組合估計(jì)的ZTD隨時(shí)間變化幅度較小，較穩(wěn)定。因此，利用BDS/GPS組合估計(jì)ZTD可以削弱利用BDS估計(jì)ZTD帶來(lái)的差值波動(dòng)幅度大的影響。

圖3 利用BDS、GPS和BDS/GPS組合估計(jì)的ZTD與IGS ZTD差值(未固定跟蹤站坐標(biāo))

通常GNSS跟蹤站以及用來(lái)觀測(cè)水汽含量的天文參考站都是固定，坐標(biāo)變化微小，因此，可以通過(guò)固定跟蹤站坐標(biāo)減少參數(shù)估計(jì)中的未知數(shù)。針對(duì)選取的8個(gè)跟蹤站，利用PPP法解算8個(gè)跟蹤站的GPS數(shù)據(jù)得到單天解，在估計(jì)ZTD的過(guò)程中使用解算的GPS單天解固定跟蹤站坐標(biāo)，并給GPS單天解附加一個(gè)微小權(quán)約束。圖4顯示了在固定跟蹤站坐標(biāo)后，利用BDS、GPS、GPS/BDS估計(jì)的ZTD與IGS ZTD差值時(shí)間序列，并按天統(tǒng)計(jì)差值的均值。從圖4可看出，固定跟蹤站坐標(biāo)后，分別利用BDS、GPS、GPS/BDS估計(jì)的ZTD與IGS ZTD的差值浮動(dòng)范圍均為-8～+8mm，這說(shuō)明固定跟蹤站坐標(biāo)提高了BDS ZTD的穩(wěn)定性，但固定跟蹤站坐標(biāo)對(duì)GPS ZTD和GPS/BDS組合ZTD影響較小，無(wú)明顯變化。

圖4 利用BDS、GPS和BDS/GPS組合估計(jì)的ZTD與IGS ZTD差值(固定跟蹤站坐標(biāo))

以IGS ZTD為參考值，表1分別統(tǒng)計(jì)了采用BDS、GPS和BDS/GPS組合PPP法在未固定跟蹤站坐標(biāo)和固定跟蹤站坐標(biāo)的情況下估計(jì)的ZTD平均偏差MEAN和STD(4個(gè)月)。分析見(jiàn)表1。

表1 固定坐標(biāo)前后利用PPP估計(jì)ZTD的年均結(jié)果

從給出的平均偏差MEAN來(lái)看，BDS ZTD在未固定跟蹤站坐標(biāo)情況下的MEAN在23 mm以內(nèi)，固定坐標(biāo)后的MEAN在14.2 mm以內(nèi)；GPS ZTD在未固定跟蹤站坐標(biāo)情況下的MEAN在1.0 mm以內(nèi)，固定坐標(biāo)后的MEAN在1.9 mm以內(nèi)；BDS/GPS組合ZTD在未固定跟蹤站坐標(biāo)情況下的MEAN在3.0 mm以內(nèi)，固定坐標(biāo)后的MEAN在2.5 mm以內(nèi)。與未固定跟蹤站坐標(biāo)估計(jì)的BDS ZTD相比，固定跟蹤站坐標(biāo)后估計(jì)的BDS ZTD的MEAN減小了50%以上；與未固定跟蹤站坐標(biāo)估計(jì)的GPS ZTD或BDS/GPS組合ZTD相比，固定坐標(biāo)后估計(jì)的GPS ZTD或BDS/GPS組合ZTD的MEAN沒(méi)有減小，這說(shuō)明使用BDS估計(jì)ZTD時(shí)固定跟蹤站坐標(biāo)較好；而使用GPS或BDS/GPS組合估計(jì)ZTD時(shí)無(wú)需固定跟蹤站坐標(biāo)。

從給出的STD來(lái)看，在未固定跟蹤站的情況下，BDS ZTD的STD優(yōu)于34 mm，GPS ZTD的STD優(yōu)于13 mm，BDS/GPS組合ZTD優(yōu)于14 mm；在固定跟蹤站的情況下，BDS ZTD的STD優(yōu)于34 mm，GPS ZTD的STD優(yōu)于13 mm，BDS/GPS組合ZTD優(yōu)于14 mm。BDS ZTD的STD明顯高于GPS ZTD或BDS/GPS組合ZTD的STD且固定跟蹤站坐標(biāo)對(duì)STD幾乎沒(méi)有影響，這是因?yàn)楸倍范l(wèi)星導(dǎo)航系統(tǒng)沒(méi)有準(zhǔn)確的姿態(tài)模型以及天線相位中心模型，利用BDS估計(jì)的ZTD存在系統(tǒng)誤差；利用BDS/GPS組合與GPS估計(jì)的ZTD精度近乎一致，均在13 mm內(nèi)，這說(shuō)明使用BDS/GPS組合與GPS估計(jì)的ZTD與IGS產(chǎn)品符合的比較好。與GPS ZTD相比，雖然BDS/GPS組合估計(jì)的ZTD沒(méi)有獲得較高的精度，但BDS/GPS組合觀測(cè)增加了可視衛(wèi)星數(shù)，這對(duì)觀測(cè)條件惡劣的地方(比如城市、山區(qū))，具有重要意義。另外，比較固定跟蹤站坐標(biāo)估計(jì)的BDS ZTD和未固定跟蹤站坐標(biāo)估計(jì)的BDS/GPS組合ZTD發(fā)現(xiàn)，雖然二者M(jìn)EAN的互差僅為幾個(gè)mm，但RMS互差為1～2 cm，由此可見(jiàn)，利用BDS/GPS組合估計(jì)ZTD的方式優(yōu)于在固定跟蹤站坐標(biāo)后利用BDS估計(jì)ZTD的方式。值得說(shuō)明的是，這六種方法估計(jì)的ZTD與IGS ZTD差值的STD均超過(guò)了6.4 mm，這種差異是利用PPP法估計(jì)ZTD過(guò)程中具體的數(shù)據(jù)處理策略不一致，如參數(shù)設(shè)置、投影函數(shù)等。

3 結(jié) 論

為分析利用BDS/GPS組合PPP法估計(jì)ZTD的效果，采用WHU提供的精密軌道和鐘差產(chǎn)品，對(duì)8個(gè)MEGX跟蹤站2014年6月—9月的數(shù)據(jù)進(jìn)行了測(cè)試，系統(tǒng)比較了采用BDS,GPS和BDS/GPS組合PPP法在未固定跟蹤站坐標(biāo)和固定跟蹤站坐標(biāo)的情況下估計(jì)ZTD的效果。結(jié)果表明，利用PPP法估計(jì)ZTD，BDS ZTD現(xiàn)階段的STD優(yōu)于34 mm，GPS ZTD現(xiàn)階段的STD優(yōu)于13 mm，BDS/GPS組合現(xiàn)階段的STD優(yōu)于14 mm。與利用GPS估計(jì)ZTD的方式相比，利用BDS/GPS組合估計(jì)ZTD的STD與GPS相當(dāng)。另外，通過(guò)固定跟蹤站坐標(biāo)的方式雖然可以提高利用BDS估計(jì)ZTD的穩(wěn)定性，但不能提高精度。

[1] 曹云昌, 方宗義, 夏青. 軌道誤差對(duì)近實(shí)時(shí) GPS遙感水汽的影響研究[J]. 氣象科技, 2004, 32(4): 229-232.

[2] SHI J, XU C, GUO J, et al. Real-time GPS precise point positioning-based precipitable water vapor estimation for rainfall monitoring and forecasting[J]. Geoscience and Remote Sensing, IEEE Transactions on, 2015, 53(6): 3452-3459.

[3] 盛崢, 黃思訓(xùn), 張流青, 等. 利用 GPS 技術(shù)探測(cè)大氣中的水汽含量[J]. 沈陽(yáng)師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2005, 23(1): 33-37.

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

[5] 畢研盟, 毛節(jié)泰, 楊光林, 等. 地基 GPS 遙感觀測(cè)安徽地區(qū)水汽特征[J]. 氣象科技, 2004, 32(4): 225-228.

[6] CHEN J, ZHANG Y, WANG J, et al. A simplified and unified model of multi-GNSS precise point positioning[J]. Advances in Space Research, 2015, 55(1): 125-134.

[7] DUAN J, BEVIS M, FANG P, et al. GPS meteorology: Direct estimation of the absolute value of precipitable water[J]. Journal of Applied Meteorology, 1996, 35(6): 830-838.

[8] 葉世榕, 張雙成, 劉經(jīng)南. 精密單點(diǎn)定位方法估計(jì)對(duì)流層延遲精度分析[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2008, 33(8): 788-791.

[9] 何錫揚(yáng), 張小紅, 李星星, 等. PPP 估計(jì)天頂對(duì)流層延遲方法與結(jié)果分析[J]. 測(cè)繪信息與工程, 2010 (1): 3-5.

[10] 蔡昌盛, 夏朋飛, 史俊波,等.利用 GPS/GLONASS 組合精密單點(diǎn)定位方法估計(jì)天頂對(duì)流層延遲[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2013, 33(2): 54-57.

[11] 范磊,鐘世明,潭冰峰,等.基于BDS與GPS精密單點(diǎn)定位的天頂對(duì)流層延遲估計(jì)比較[J]. 大地測(cè)量與地球動(dòng)力學(xué),2015(1):67-71.

[12] 王珍.網(wǎng)絡(luò)RTK動(dòng)態(tài)中長(zhǎng)基線模糊度解算方法比較研究[J].測(cè)繪與空間地理信息，2016,39(2)：71-74.

[13] XU A, XU Z, GE M, et al. Estimating zenith tropospheric delays from BeiDou navigation satellite system observations[J]. Sensors, 2013, 13(4): 4514-4526.

[14] 李星星. GNSS精密單點(diǎn)定位及非差模糊度快速確定方法研究[D].武漢：武漢大學(xué),2013.

[15] 葛茂榮, 劉經(jīng)南. GPS定位中對(duì)流層折射估計(jì)研究[J].測(cè)繪學(xué)報(bào), 1996, 25(4): 285-291.

[責(zé)任編輯：李銘娜]

Precision analysis on the zenith tropospheric delays using combined BDS/GPS Data

CHEN Zhenguo1,2, TANG Longjiang3

(1. Geological Brigade 272 of Jiangsu Nuclear Industry, Nanjing 210000, China; 2. The Nuclear Industry Construction Group Co.,Ltd., Nanjing 210000, China; 3. School of Geomatics of Liaoning Technical University, Fuxin 123000, China)

In order to analyze the precision of the zenith tropospheric delays (ZTD) using combined BDS/GPS data under the conditions of the existing final products, the PPP method is used to calculate observation data from June to September in 2014 collected from 8 MGEX stations fixed and unfixed. ZTD result from different conditions is compared with ZTD products of IGS. The result shows that the standard deviation (STD) of BDS ZTD is better than 34mm and in addition.A similar result is achieved by combined BDS/GPS compared to that of GPS only by using PPP method. Compared with BDS ZTD under the unknown of MGEX station, the BDS ZTD under the known of MGEX station is more stable. However, the precision can’t be improved.

combined BDS/GPS; ZTD products of IGS; Precision Point Positioning method; MGEX station; ZTD

2016-08-09

國(guó)家青年科學(xué)基金資助項(xiàng)目(41504030)

陳振國(guó)(1991-)，男，助理工程師.

著錄：陳振國(guó)，唐龍江.BDS/GPS組合天頂對(duì)流層延遲的精度分析[J].測(cè)繪工程，2017，26(9)：5-9.

10.19349/j.cnki.issn1006-7949.2017.09.002

P228

A

1006-7949(2017)09-0005-05