宋 艷,應斌杰
(國網浙江省電力有限公司 麗水供電公司,浙江 麗水 323000)
GPS等衛星導航系統對軍事和民用領域都產生了革命性的影響。基于衛星的定位、導航和授時(Positioning,Navigation and Timing,PNT)設備廣泛應用于各個國民經濟行業[1]。近年來,GNSS/GPS干擾已經逐漸出現在實際應用過程中,衛星導航存在先天的脆弱性,信號強度弱,頻點和格式信息公開,易受阻擋、干擾,甚至是欺騙。有報道指出,一個價值29美元的設備就可以對GPS信號進行阻塞甚至欺騙[2]。欺騙是一種干擾形式,欺騙干擾源使用類似GNSS信號來影響導航接收機的正常工作,使其計算出錯誤的位置、速度或時間。2008年,英國政府測試使用兩個低成本干擾器阻塞了北海30 km范圍內的GPS信號[3]。2015年,阿里巴巴移動通信安全部門的研究人員利用簡易的欺騙干擾源成功實施了對Apple Watch和華為手機的時間和位置欺騙[4]。虛假GNSS/GPS信號被播發到接收機上會造成位置、時間、導航出現誤差,從而影響人們的交通、導航以及日常生活,體現了GNSS系統脆弱的一面。
經過多年的研究與發展,衛星授時抗欺騙干擾研究取得了一定的成果,隨著技術的進步,有以下幾個研究方向和發展趨勢。
欺騙干擾技術不斷發展,欺騙干擾功率越來越低,尤其入侵式欺騙干擾隱蔽性很強[5],早期在跟蹤環路很難發現,因此需要研究多手段跟蹤環路技術有效監測到入侵式干擾信號。例如,通過對PLL環鑒相器輸出信號進行譜分析就能檢測出信號中是否含有欺騙信號,利用載波相位跟蹤譜分析實現入侵式欺騙信號檢測。
自動增益控制(Automatic Gain Control,AGC)是一種在輸入信號變化很大的情況下,通過調節可控增益放大器,使輸出信號幅值基本恒定或僅在較小范圍內變化的系統[6]。傳統的AGC系統根據提取的反饋信號可大致分為基于功率檢測和基于峰值檢測兩類。但隨著干擾信號樣式的不同,其期望功率及期望電平不斷變化,傳統AGC系統不能滿足模數轉換器(Analog-to-Digital Converter,ADC)接近最小信噪比損耗的要求,因此需要研究自適應AGC調整算法,使得ADC之后損耗接近最小值,從而實現智能干擾檢測[7]。
繼美國GPS衛星導航系統投入使用,并取得重大軍事和經濟效益之后,俄羅斯的GLONASS系統、我國的北斗系統以及歐洲的GALILEO系統先后開始建設,衛星導航接收機的多模式兼容已成為衛星導航的發展趨勢[8,9]。每種衛星導航模式中又有多個工作頻點,充分利用多模式多頻的衛星導航信號,可大幅提高抗欺騙干擾接收機的性能。
傳統欺騙檢測技術應軍事背景需求而產生,基本是應用在強對抗、瞄準式以及持續的干擾場景下,使得復雜度和成本等都比較大,而民用授時環境是靜態的,有高穩時鐘源,因此其對抗技術有特殊性[10]。在分析研究民用授時環境下欺騙干擾特征及作用機理的基礎上,結合欺騙干擾檢測技術和守時技術,可大大提高授時接收機的抗欺騙干擾性能,同時降低目標實現的復雜度、施工難度、成本、體積以及功耗等。
欺騙干擾的實現基于偽距定位原理,對于一個正常工作的接收機而言,必須得到4顆及以上有效衛星觀測數據,然后求解定位方程組來獲取接收機的位置信息。接收機位置坐標(x,y,z)和接收機鐘差Δtu組成的4個未知量的獲取則依賴于衛星位置、衛星鐘差以及觀測偽距ρ。對應的PVT解算方程為:

篡改導航電文或影響偽距觀測結果,將導致PVT解算結果出錯。
這里分析通過拉偏碼偽距觀測影響輸出PVT結果的方法,接收機的碼偽距觀測信息包含了偽距、偽距的1階量(碼速率偏移量)、2階量(碼加速率偏移量)以及3階量(碼加加速率偏移量)。開展碼偽距拉偏時,1階量的變化范圍可以很大,但2階量和3階量要做拉偏時,變化范圍是有限的,尤其對定時型接收機而言,其應用場景是靜態的,一般環路動態可承受變化范圍都比較小。針對不同干擾場景情況下的PVT結算結果,檢測情況如表1所示。

表1 不同干擾場景下的檢測情況
時鐘模型建立包括了時鐘模型描述以及參數優化計算兩個部分。導航接收機正常解算后,就可以利用頻差變化來建立二項式時鐘偏差模型。對于參數的優化計算可以采用最小二乘法、Kalman濾波法以及自適應Kalman濾波法等多種手段,這些手段可以有效克服測量噪聲以及測量野值的影響,獲得可靠的模型參數估計值。
目前,基于線性預報的多項式模型通過若干階多項式完成對鐘差Clk(t)的描述為:

式中,a0,…,am表示多項式系數,m為階數,t0為參考時刻。
在建立模型過程中,首先要確定模型的階數,時鐘為銣鐘時通常用二階模型就可以加以精確描述,即:

式中,x0,y0,a分別表示銣鐘的鐘差、頻差以及頻差變化率,n(t)為銣鐘的隨機誤差。
在對接收機鐘差進行預測時,本項目采用Kalman濾波完成對模型參數的優化估計。Kalman濾波算法是最小方差意義上的最優估計算法。采用動力學方程(狀態方程)描述被估計量的動態變化規律,被估計量的動態統計信息由激勵白噪聲的統計信息和動力學方程確定,由于考慮了銣鐘隨機誤差部分的影響,因而其預測精度較高且跟蹤性能較強。狀態變量包括銣鐘的鐘差、頻差以及頻差變化率,連續形式的狀態方程為:

式中,n1(t)和n2(t)為相互獨立的白噪聲,分別代表隨機噪聲部分的頻率白噪聲和頻率隨機游走噪聲。若測量間隔為τ,對式(4)進行離散化可得狀態轉移矩陣Φk+1,k和噪聲方差陣Q為:

式中,方差σ12和σ22可由和σ22=π2h-2來確定。
觀測量為鐘差值,則觀測方程為:

式中,n3(k)代表隨機噪聲部分的相位白噪聲,n4(k)為相位量測白噪聲。
觀測方程和狀態轉移方程建立后,依據Kalman濾波的過程即可完成參數模型的實施更新。
衛星導航信號弱、頻點和格式公開,因此較易被干擾和欺騙,而假的GNSS/GPS信號被播發到接收機上會造成位置、時間以及導航誤差。抗干擾技術是實現衛星精準授時的基礎,必將在電網發展的將來起到重要作用,同時也將為電網的發展做出貢獻。