昌平斷裂交匯區(qū)GNSS觀測網(wǎng)絡(luò)數(shù)據(jù)質(zhì)量分析

魯亞軍 馬一方 胡樂銀 李平安

(北京市地震局, 北京 100080)

0 引言

昌平-順義-海淀交界區(qū)是北東向小湯山-東北旺斷裂和黃莊-高麗營斷裂與北西向南口-孫河斷裂的交匯區(qū)(簡稱昌平斷裂交匯區(qū)),該區(qū)域是北京平原區(qū)內(nèi)活動斷裂、第四紀隱伏凹陷、現(xiàn)今地面沉降和地震活動疊加發(fā)育發(fā)生的典型地區(qū)[1-3]。利用全球?qū)Ш叫l(wèi)星系統(tǒng)(global navigation satellite system,GNSS)開展該區(qū)域高精度的地殼形變動態(tài)監(jiān)測與數(shù)據(jù)分析,能夠深化區(qū)域構(gòu)造活動、地面沉降、地震活動特征與規(guī)律及其相互關(guān)系的認識,有利于提升該區(qū)域的地震監(jiān)測能力和對地震活動趨勢的科學(xué)把握。

為了完善昌平斷裂交匯區(qū)GNSS觀測網(wǎng)絡(luò),北京市地震局于2019年8月新建了5個GNSS連續(xù)運行觀測站,與區(qū)內(nèi)4個現(xiàn)有連續(xù)運行GNSS觀測站[4]組成一個觀測網(wǎng)絡(luò),以獲取該區(qū)內(nèi)三維地殼形變信息以及3條主要活動斷裂(南口-孫河、東北旺-小湯山和黃莊高麗營斷裂)的相對運動狀態(tài)。由于GNSS原始觀測數(shù)據(jù)質(zhì)量直接關(guān)系到后續(xù)數(shù)據(jù)處理結(jié)果的精度和可靠性,對其進行質(zhì)量分析和檢查是應(yīng)用中的一項非常重要的基礎(chǔ)工作[5-8]。因此,本文選取觀測數(shù)據(jù)完整率、周跳比、多路徑效應(yīng)、信噪比等指標進行昌平斷裂交匯區(qū)GNSS連續(xù)運行基準站觀測數(shù)據(jù)的質(zhì)量評價工作。研究成果有利于后續(xù)GNSS數(shù)據(jù)處理選擇更為合適的函數(shù)模型和隨機模型,且在數(shù)據(jù)解算結(jié)果出現(xiàn)異常時能夠快速找出原因,對于GNSS監(jiān)測區(qū)域地殼運動的結(jié)果分析具有重要的指導(dǎo)作用。

1 數(shù)據(jù)情況

本次數(shù)據(jù)質(zhì)量檢查工作不僅包含昌平斷裂交匯區(qū)5個新建的GNSS觀測站,而且包含區(qū)內(nèi)4個現(xiàn)有的GNSS觀測站(北京市連續(xù)運行基準站，BJCORS),9個測站的空間分布如圖1所示,其中,五角星代表新建測站,圓點代表原有測站。對于新建測站,即板橋站(BANQ)、昌平地震臺站(CPDS)、鳳凰山陵園站(FHSL)、高崖口站(GYKC)、果莊村站(GZHC)，其觀測數(shù)據(jù)包含全球定位系統(tǒng)(global positioning system,GPS)、北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system,BDS)、格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)(GLONASS是俄語“全球衛(wèi)星導(dǎo)航系統(tǒng)global navigation satellite system”的縮寫)和伽利略(Galileo)4個系統(tǒng)。對于原有測站,即測繪院基站(CEHY)、昌平基站(CHPN)、石景山基站(SHIJ)、牛欄山基站(NLSH)，其觀測數(shù)據(jù)包含GPS和GLONASS兩個系統(tǒng)。本文選取2019年第218天和2020年第366天的觀測數(shù)據(jù)進行質(zhì)量檢查,采樣間隔為15 s,截至高度角為10°。

圖1 昌平斷裂交匯區(qū)GNSS連續(xù)運行基準站網(wǎng)的空間分布注：圖中測站,GZHC為果莊村站，F(xiàn)HSL為鳳凰山陵園站，CHPN為昌平基站，BANQ為板橋站，NLSH為牛欄山基站，GYKC為高崖口站，CPDS為昌平地震臺站，SHIJ為石景山基站，CEHY為測繪院基站。

2 參數(shù)計算

本文根據(jù)《中國地殼運動觀測技術(shù)規(guī)程》《全球定位系統(tǒng)(GPS)測量規(guī)范》(GB/T 18314—2009),選取觀測數(shù)據(jù)完整率、周跳比、多路徑效應(yīng)(MP1和MP2)、信噪比作為質(zhì)量評價指標[9-11],使用的數(shù)據(jù)質(zhì)量判定標準為數(shù)據(jù)有效率大于85%,o/slips大于200,MP1平均值小于0.5 m,MP2平均值小于0.75 m[12-14]。本文利用Fortran語言編寫程序完成支持RINEX2/3格式的GPS/BDS/GLONASS/Galileo觀測數(shù)據(jù)質(zhì)量檢查,下面對選取的相關(guān)指標的計算進行說明。

數(shù)據(jù)完整率能夠反映GNSS測站周邊環(huán)境,例如，是否存在較高建筑物遮擋,其計算公式為

(1)

式中,Na為完整觀測值個數(shù)；Ns為預(yù)期完整觀測值個數(shù),其計算方法為利用廣播文件計算相應(yīng)時刻測站觀測到的所有衛(wèi)星的高度角,如果計算高度角大于10°,則增加預(yù)期完整觀測值個數(shù)。

周跳比能夠反映數(shù)據(jù)的周跳情況,為完整觀測值個數(shù)與周跳數(shù)的比值o/slips,數(shù)值越小表明周跳越嚴重。本文采用電離層殘差組合和MW組合進行周跳的探測，如式(2)、式(3)所示。

式中,φ1和φ2為相位觀測值；P1和P2為偽距觀測值；f1和f2為載波頻率。

本文計算的MP1和MP2值分別表示L1、L2載波上的多路徑效應(yīng)對距和相位影響的綜合指標,其值越大,多路徑影響越嚴重,其計算公式為

(4)

信噪比是指有用信號和噪聲的比值,與衛(wèi)星狀況及周圍觀測環(huán)境有關(guān)。本文利用接收機給出的L1和L2相位觀測值的原始信號強度日均值給出,給出L1和L2的信噪比值SNR1和SNR2，其中，SNR是信噪比(signal-to-noise ratio)。

3 數(shù)據(jù)質(zhì)量分析結(jié)果

在數(shù)據(jù)質(zhì)量檢查之前,本文首先利用GDP軟件查看了各測站的衛(wèi)星可見性[15]。GPS和GLONASS為全星座運行,在一天內(nèi)兩衛(wèi)星系統(tǒng)均可觀測15顆衛(wèi)星以上。BDS和Galileo在2019年正處于組網(wǎng)階段,特別是Galileo可見衛(wèi)星很少且觀測質(zhì)量較差。因此本文在對新建測站的分析中,將不再包含Galileo觀測數(shù)據(jù)。對于原有測站,其觀測數(shù)據(jù)僅包含GPS和GLONASS。限于篇幅,下文隨機選取BANQ測站數(shù)據(jù)質(zhì)量檢查結(jié)果進行詳細說明,其他測站僅列出日均值結(jié)果。

圖2給出了BANQ站的數(shù)據(jù)質(zhì)量檢查結(jié)果,由圖可知,在數(shù)據(jù)完整率方面,GPS數(shù)據(jù)完整率均超過85%,而BDS數(shù)據(jù)完整率較低(<80%),在觀測時段內(nèi)(共計514天),GLONASS數(shù)據(jù)僅有107天的數(shù)據(jù)完整率超過85%;在周跳比方面,3個系統(tǒng)數(shù)據(jù)均有較好的結(jié)果,表明數(shù)據(jù)接收穩(wěn)定,受衛(wèi)星失鎖、信號中斷等影響較小;在多路徑效應(yīng)方面,除了個別觀測天數(shù)MP1值超過0.5 m,其他時間內(nèi)MP1和MP2值均符合數(shù)據(jù)判斷標準要求;在信噪比方面,3個系統(tǒng)接收信號的強度均較為良好。表1～表3給出了所有測站在觀測時段內(nèi)各項數(shù)據(jù)質(zhì)量檢查指標的日均值。結(jié)果表明,新建測站在數(shù)據(jù)有效率方面較原有測站較差,新建測站GPS觀測的數(shù)據(jù)完整率比較穩(wěn)定,其他兩個系統(tǒng)數(shù)據(jù)完整率較低,后續(xù)有待進行相應(yīng)改進。原有測站的o/slips較新建測站離散程度小,尤其是新建測站BDS和GLONASS兩個系統(tǒng)的周跳相對嚴重。在信噪比方面,除了GPS和GLONASS的L2相位觀測值信號強度相對較弱以外,其他GNSS信號強度均較為良好。所有測站受到多路徑的影響均較小,特別是新建測站的BDS觀測數(shù)據(jù)。因此,昌平斷裂交匯區(qū)GNSS觀測網(wǎng)絡(luò)觀測質(zhì)量較好,后續(xù)可繼續(xù)升級硬件系統(tǒng)以滿足BDS應(yīng)用的需要。

圖2 BANQ站數(shù)據(jù)質(zhì)量檢查結(jié)果注：圖中，BDS為北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system)；GPS為全球定位系統(tǒng)(global positioning system)；GLONAS為格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)(是俄語“全球衛(wèi)星導(dǎo)航系統(tǒng)global navigation satellite system”的縮寫)；SNR為信噪比(signal-to-noise ratio); MP1 和 MP2 值分別表示L1、L2 載波上的多路徑效應(yīng)對距和相位影響的綜合指標。

表1 觀測時段內(nèi)各測站數(shù)據(jù)質(zhì)量檢查指標平均值(GPS)

表2 觀測時段內(nèi)各測站數(shù)據(jù)質(zhì)量檢查指標平均值(GLONASS)

表3 觀測時段內(nèi)各測站數(shù)據(jù)質(zhì)量檢查指標平均值(BDS)

4 結(jié)束語

本文通過選取數(shù)據(jù)完整率、多路徑效應(yīng)、信噪比以及周跳比等數(shù)據(jù)質(zhì)量評價指標,對昌平斷裂交匯區(qū)9個GNSS連續(xù)運行基準站的觀測數(shù)據(jù)質(zhì)量進行了分析與評估,得出以下結(jié)論：(1)所有測站的GPS觀測數(shù)據(jù)完整率均符合規(guī)范要求,新建測站的BDS和GLONASS數(shù)據(jù)完整率較低,不符合規(guī)范要求;(2)所有測站受到多路徑效應(yīng)的影響均較小;(3)所有測站接收到信號的強度良好,尤其是新建測站的BDS觀測數(shù)據(jù);(4)所有測站的周跳現(xiàn)象均較少,觀測數(shù)據(jù)穩(wěn)定可靠。昌平斷裂交匯區(qū)9個GNSS連續(xù)運行基準站的GPS觀測數(shù)據(jù)質(zhì)量較好,可為后續(xù)地殼形變的監(jiān)測提供完整、可靠、高精度的觀測數(shù)據(jù),新建測站的硬件設(shè)備有待升級,以滿足后續(xù)BDS的推廣應(yīng)用。