包 海，方書山，張晶晶，鄢中堡，張 芯

(1.四川省第一測繪工程院，四川 成都 610100；2.中國測繪科學(xué)研究院，北京 100830)

?

局域BDS CORS網(wǎng)衛(wèi)星鐘差完備性監(jiān)測

包 海1，方書山2，張晶晶1，鄢中堡1，張 芯1

(1.四川省第一測繪工程院，四川 成都 610100；2.中國測繪科學(xué)研究院，北京 100830)

針對(duì)北斗系統(tǒng)衛(wèi)星鐘差存在一定故障(粗差)發(fā)生的概率，當(dāng)故障發(fā)生時(shí)不能保證導(dǎo)航定位的精度，而由于地面站點(diǎn)BDS觀測數(shù)據(jù)來源的限制，對(duì)應(yīng)BDS系統(tǒng)的鐘差完備性監(jiān)測方法研究論文很少的現(xiàn)狀，著重研究了局域范圍內(nèi)BDS系統(tǒng)鐘差完備性監(jiān)測方法，并以某省11個(gè)CORS站點(diǎn)實(shí)時(shí)數(shù)據(jù)進(jìn)行試驗(yàn)，利用殘差法則查找出含有粗差鐘差的BDS衛(wèi)星，同時(shí)剔除粗差衛(wèi)星并及時(shí)報(bào)警給用戶。實(shí)驗(yàn)結(jié)果表明能有效實(shí)現(xiàn)BDS衛(wèi)星的鐘差完備性監(jiān)測。

衛(wèi)星鐘差；完備性；參考鐘；北斗網(wǎng)；監(jiān)測

0 引言

隨著全球定位系統(tǒng)不斷發(fā)展，尤其是北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system，BDS)的快速組建，導(dǎo)航定位星座不斷增多，導(dǎo)航定位的精度不斷提高，特別是精密單點(diǎn)定位及網(wǎng)絡(luò)實(shí)時(shí)動(dòng)態(tài)差分法(real time kinematic，RTK)技術(shù)的發(fā)明，使得定位精度能達(dá)到動(dòng)態(tài)精度亞m級(jí)、靜態(tài)cm量級(jí)甚至mm級(jí)，已經(jīng)能夠滿足絕大多數(shù)用戶對(duì)精度的要求；同時(shí)，決定用戶安全性能的導(dǎo)航系統(tǒng)的完備性問題顯得尤其突出[1-3]。完備性是導(dǎo)航系統(tǒng)發(fā)生任何故障或者誤差超限，無法用于導(dǎo)航和定位時(shí)，系統(tǒng)向用戶及時(shí)發(fā)出報(bào)警的能力[4]。

衛(wèi)星鐘頻率漂移引起的衛(wèi)星鐘時(shí)間與BDS標(biāo)準(zhǔn)時(shí)之間的差值稱為衛(wèi)星鐘差。衛(wèi)星鐘差在精密定位中占有重要地位；高精度的BDS測量歸根到底是高精度的時(shí)間測量：衛(wèi)星鐘差質(zhì)量的好壞，直接影響到測量的結(jié)果。BDS衛(wèi)星鐘差完備性監(jiān)測就是找出含有粗差鐘差的衛(wèi)星，剔除并及時(shí)報(bào)警給用戶，給用戶提供“干凈的”、不含粗差鐘差的高質(zhì)量BDS信號(hào)，確保測量的精度，從而保證用戶的安全。

在國內(nèi)，許多省份都相繼建立了自己的BDS連續(xù)運(yùn)行參考站(continuously operating reference stations，CORS)網(wǎng)絡(luò)，可向大量用戶提供高精度、實(shí)時(shí)的定位信息；但此信息中未包含與用戶生命安全相關(guān)的完備性信息。系統(tǒng)的精度可以根據(jù)需要和可能進(jìn)行調(diào)整與控制，但系統(tǒng)的完備性信息在任何時(shí)候都是不可缺少的。從用戶的安全角度考慮，導(dǎo)航系統(tǒng)的完備性比精度有著更加重要的地位，它是保證用戶安全性的重要參數(shù)[5-7]。本文基于天寶R9提供的BDS偽距觀測值[8]，研究衛(wèi)星鐘差完備性監(jiān)測理論，實(shí)現(xiàn)BDS衛(wèi)星鐘差完備性監(jiān)測。

1 BDS衛(wèi)星鐘差完備性監(jiān)測方法

1.1 衛(wèi)星鐘差估計(jì)

本文在實(shí)時(shí)衛(wèi)星鐘差估計(jì)中，利用偽距碼為原始觀測量，雖然它們觀測噪聲大、測距精度低(衛(wèi)星鐘差估計(jì)只能達(dá)到納秒級(jí))，但是可避免周跳的探測與修復(fù)及模糊度的確定等復(fù)雜算法?；谀硞€(gè)地面CORS站及某顆BDS衛(wèi)星，組成單站單星，它們的偽距觀測方程為

Δtj(i))+ΔDw+ΔDion(i)+ΔDtro(i)+ΔDrel(i)。

(1)

對(duì)式(1)變形可得確定衛(wèi)星鐘差的誤差方程為：

(2)

(3)

當(dāng)有m個(gè)基準(zhǔn)站，同步跟蹤n顆衛(wèi)星時(shí)，設(shè)測站序列號(hào)為(1,2,…,m)，設(shè)衛(wèi)星序列號(hào)為(1,2,…,n)，省略時(shí)刻i的標(biāo)注，確定接收機(jī)和衛(wèi)星鐘差的誤差方程為

(4)

即寫成：

(5)

(6)

當(dāng)采用式(5)進(jìn)行求解可以發(fā)現(xiàn)，該方程是奇異的，無法直接解算鐘差參數(shù)；故選擇一個(gè)基準(zhǔn)站的接收機(jī)鐘或者某顆衛(wèi)星的衛(wèi)星鐘作為參考鐘，求其他接收機(jī)鐘和衛(wèi)星鐘對(duì)于參考鐘的相對(duì)鐘差。可以證明利用絕對(duì)鐘差與相對(duì)鐘差進(jìn)行定位其坐標(biāo)分量解是一致的，僅鐘差有一偏差量，只需保證基準(zhǔn)鐘的鐘差精度優(yōu)于10-6s，則相對(duì)鐘差和絕對(duì)鐘差對(duì)用戶定位結(jié)果而言是等價(jià)的，即相對(duì)鐘差的系統(tǒng)性偏差在用戶定位模型中完全被用戶接收機(jī)鐘差吸收，而不影響用戶的定位精度[9]。這就把求解絕對(duì)鐘差轉(zhuǎn)化為求解相對(duì)鐘差(對(duì)于地面測站接收機(jī)鐘和星載鐘在估計(jì)歷元時(shí)刻是否鐘差精度優(yōu)于10-6s，需要事先確認(rèn))。

1)以地面接收機(jī)鐘作為參考鐘

假設(shè)在估計(jì)歷元時(shí)刻i，第1個(gè)測站上的接收機(jī)能夠滿足要求，故選取這個(gè)基準(zhǔn)站上的接收機(jī)作為參考站，則在時(shí)刻i此測站上的觀測值誤差方程為

(7)

其他第k個(gè)地面站觀測值的誤差方程為

(8)

當(dāng)m個(gè)基準(zhǔn)站，同步跟蹤n顆衛(wèi)星時(shí)(省略時(shí)刻i的標(biāo)注)，可寫成

(9)

故為

(10)

此時(shí)變化之后的誤差方程與式(5)類似，關(guān)鍵在于未知參數(shù)為(n+m-1)個(gè)的相對(duì)鐘差參數(shù)比式(5)少1個(gè)未知數(shù)，此時(shí)方程有唯一解。設(shè)計(jì)矩陣A及未知參數(shù)X(i)都發(fā)生變化，具體表示為：

(11)

(12)

(13)

(14)

(15)

2)以星載鐘作為參考鐘

除選地面基準(zhǔn)站接收機(jī)鐘為參考鐘外，還可以選取星載鐘為參考鐘，為保持在局域網(wǎng)中對(duì)同一衛(wèi)星的可見性，一般選取局域網(wǎng)中高度角比較高的衛(wèi)星。如在局域網(wǎng)中有m個(gè)基準(zhǔn)站(序列號(hào)(1,2,…,m))，同步跟蹤到n顆衛(wèi)星(序列號(hào)為(1,2,…,n))，觀測到衛(wèi)星中，假設(shè)第1顆星(序列號(hào)為1但衛(wèi)星號(hào)不一定為1)在時(shí)刻i滿足要求，選取其作為參考站，j表示測站號(hào)，則可列出此星上觀測值誤差方程為

(16)

其他第k顆衛(wèi)星的誤差方程為

(17)

當(dāng)m個(gè)基準(zhǔn)站，同步跟蹤n顆衛(wèi)星時(shí)(省略時(shí)刻i的標(biāo)注)，可列出

(18)

故可寫成

(19)

同理此時(shí)變化之后的誤差方程與式(5)類似，未知參數(shù)也為(n+m-1)個(gè)相對(duì)鐘差參數(shù)，比式(5)少1個(gè)未知數(shù)，此時(shí)方程有唯一解。設(shè)計(jì)矩陣A及未知參數(shù)X(i)都發(fā)生變化，具體表示為:

(20)

(21)

(22)

(23)

(24)

3)2種方式的計(jì)算

利用最小二乘法，解得式(10)、式(19)中各相對(duì)鐘差值為

(25)

式(15)、(24)中各殘差的值為

(-A(ATA)-1AT+E)L。

(26)

1.2 衛(wèi)星鐘差完備性監(jiān)測

衛(wèi)星鐘差的完備性監(jiān)測主要是指利用多個(gè)基準(zhǔn)站確定衛(wèi)星鐘差，分析同一衛(wèi)星多個(gè)站上觀測值的殘差，確定衛(wèi)星鐘差的可用性[10]。通過1.1節(jié)的衛(wèi)星鐘差估計(jì)模型得到相應(yīng)的單位權(quán)均方根(root mean square，RMS)以及相應(yīng)的觀測值殘差v，可以采用最普遍的3倍或者2倍RMS的方式進(jìn)行簡單的粗差識(shí)別和探測：當(dāng)殘差v大于2倍RMS時(shí)，認(rèn)為該觀測值含有粗差；當(dāng)某一顆衛(wèi)星的所有殘差或者大部分殘差滿足v大于2倍RMS時(shí)，則認(rèn)為該衛(wèi)星的觀測值殘差與其他衛(wèi)星殘差方差有顯著差異，且超過限值，則認(rèn)為該衛(wèi)星鐘差存在粗差，會(huì)在一定的時(shí)間限制范圍內(nèi)(延遲時(shí)間一般應(yīng)小于2 s)給出示警信息[11]，實(shí)現(xiàn)衛(wèi)星鐘差的完備性監(jiān)測。

2 實(shí)驗(yàn)與結(jié)果分析

計(jì)算采用2015-09-08某省11個(gè)CORS站(分布見圖1)的實(shí)時(shí)觀測數(shù)據(jù)(天寶R9輸出，兼容BDS系統(tǒng))，數(shù)據(jù)采樣率為1 s。前面已經(jīng)介紹，可以采取地面基準(zhǔn)站和星載鐘作為基準(zhǔn)鐘。本文采用這2種方式并進(jìn)行比較分析。

圖1 站點(diǎn)分布

1)對(duì)于地面基準(zhǔn)站的選擇：數(shù)據(jù)處理中為保持基準(zhǔn)站與其他站良好的相關(guān)性，盡量選擇處于CORS網(wǎng)中間的站點(diǎn)作為基準(zhǔn)站，因此基準(zhǔn)站點(diǎn)的接收機(jī)鐘相對(duì)穩(wěn)定，鐘差精度要優(yōu)于10-6s；經(jīng)過分析中間MC基準(zhǔn)站接收機(jī)鐘差跳變相對(duì)其他幾個(gè)站點(diǎn)更為穩(wěn)定，所以選取此站為基準(zhǔn)站。

2)對(duì)于空間的星載鐘的選擇：選擇同步軌道衛(wèi)星3號(hào)星，其位置大體不變，而且高度角比較高，且能在全部時(shí)段可見。

同時(shí)以MC站和3號(hào)同步衛(wèi)星為參考站，單個(gè)歷元計(jì)算分別得出如圖2、圖3的結(jié)果。

圖2 以MC站為參考站的殘差v

圖3 以3號(hào)星為參考站的殘差v

圖2、圖3中X軸表示當(dāng)前觀測到的衛(wèi)星，Y軸表示所采用的觀測站，Z軸(豎軸)表示相應(yīng)衛(wèi)星和觀測站的殘差值(殘差值有正有負(fù)值，為方便畫圖直觀，負(fù)值均取絕對(duì)值變成正值)。值單位為納秒(ns)。由圖中可以看出同樣10號(hào)衛(wèi)星的殘差與其他衛(wèi)星存在明顯差異，因此可基本認(rèn)定10號(hào)衛(wèi)星存在粗差。

由圖2和圖3可知，利用地面站作為參考鐘與利用星載鐘作為參考鐘計(jì)算的結(jié)果是一致的，都通過數(shù)據(jù)計(jì)算得到10號(hào)衛(wèi)星存在粗差，將粗差衛(wèi)星及時(shí)提示給用戶，保證用戶的安全，實(shí)現(xiàn)衛(wèi)星鐘差完備性監(jiān)測。

3 結(jié)束語

本文基于完備性監(jiān)測的理念研究了BDS衛(wèi)星鐘差完備性監(jiān)測理論。基于省局域網(wǎng)的11個(gè)BDS CORS站，以地面站和星載鐘作為參考鐘，同時(shí)實(shí)時(shí)進(jìn)行了衛(wèi)星鐘差完備性監(jiān)測。由結(jié)果可知，2種方式計(jì)算的結(jié)果是一致的，能夠辨識(shí)出粗差衛(wèi)星。

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

[2] 秘金鐘.GNSS完備性監(jiān)測理論與應(yīng)用[M].北京：測繪出版社，2012：1-3.

[3] 中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室.北斗衛(wèi)星導(dǎo)航系統(tǒng)空間信號(hào)接口控制文件公開服務(wù)信號(hào)B1I(1.0版)[EB/OL].(2012-12-27)[2016-02-24].http：//www.beidou.gov.cn/attach/2012/12/27/2012122755318f7eabbe451aa6d052f829f92e50.pdf.

[4] 秘金鐘，李毓麟.衛(wèi)星導(dǎo)航完備性監(jiān)測的最新進(jìn)展[J].測繪科學(xué)，2004，29(1)：64-67.

[5] 秘金鐘.GNSS完備性監(jiān)測方法、技術(shù)與應(yīng)用[D].武漢：武漢大學(xué)，2010：1-4.

[6] 郭英.GALILEO系統(tǒng)完備性理論及其應(yīng)用研究[D].青島：山東科技大學(xué)，2006：2-4.

[7] 陳金平.GPS完善性增強(qiáng)研究[D].鄭州：信息工程大學(xué)測繪學(xué)院，2001：1-6.

[8] ANONYMOUS.Trimble NetR9 GNSS reference receiver[EB/OL].(2010-5)[2016-02-24].http：//trl.trimble.com/docushare/dsweb/Get/Document-495804/NetR9-UserGuide-13506.pdf.

[9]李星星，徐運(yùn)，王磊.非差導(dǎo)航衛(wèi)星實(shí)時(shí)/事后精密鐘差估計(jì)[J].武漢大學(xué)學(xué)報(bào)?信息科學(xué)版，2010，35(1)：661-664.

[10]GROVER R B.A baseline GPS RAIM scheme and a note on the equivalence of three RAIM methods navigation[J].Journal of Institute of Navigation，1992，39(3)：301-316.

[11]劉經(jīng)南，陳俊勇，張燕平，等.GPS廣域差分定位原理與方法[M].北京：測繪出版社，1999：57-63.

Integrity monitoring of satellite clock offset based on local area BDS CORS network

BAO Hai1，F(xiàn)ANG Shushan2，ZHANG Jinjin1，YAN Zhongbao1，ZHANG Xin1

(1.Sichuan First Institute of Surveying and Mapping Engineering，Chengdu，Sichuan 610100，China；2.Chinese Academy of Surveying and Mapping，Beijing 100830，China)

Aiming at the current situation that there is some probability that BDS satellite clock would have trouble，and the accuracy of navigation is not guaranteed when failures occur，while it is lack of related research on integrity monitoring of satellite clock offset for BDS due to the limits of data source，the paper focused on the methods of integrity monitoring of BDS satellite clock in local area.Through the experiment of 11 CORS-site data，the satellites with clock error were found by the residual error rule，meanwhile the error satellites were removed and then users were informed in time.Result showed that the method could implement integrity monitoring of BDS satellite clock error effectively.

satellite clock error；integrity；reference clock；BDS network；monitoring

2016-02-26

國家自然科學(xué)基金項(xiàng)目(41304030)；國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFB0502105)；國家863計(jì)劃項(xiàng)目(2015AA124001)；四川省測繪地理信息局科技支撐項(xiàng)目(J2015ZC01)；四川省科技廳項(xiàng)目(2015SZ0046)；中國測繪科學(xué)研究院科研業(yè)務(wù)費(fèi)支持項(xiàng)目(7771519，7771503)。

包海(1982—)，男，四川達(dá)州人，碩士，工程師，研究方向?yàn)樾l(wèi)星導(dǎo)航與定位。

方書山(1983—)，男，湖北咸寧人，博士，研究方向?yàn)樾l(wèi)星導(dǎo)航與定位。

包海，方書山，張晶晶，等.局域BDS CORS網(wǎng)衛(wèi)星鐘差完備性監(jiān)測[J].導(dǎo)航定位學(xué)報(bào)，2016，4(4)：65-68,116.(BAO Hai,FANG Shushan,ZHANG Jinjin，et al.Integrity monitoring of satellite clock offset based on local area BDS CORS network[J].Journal of Navigation and Positioning，2016，4(4)：65-68,116.)

10.16547/j.cnki.10-1096.20160413.

P228

A

2095-4999(2016)04-0065-05