余文艷，肖志剛，李 虎

(1.中國科學(xué)院國家空間科學(xué)中心 衛(wèi)星運(yùn)控技術(shù)實(shí)驗(yàn)室，北京 100190；2.中國科學(xué)院大學(xué)，北京 100049)

0 引 言

衛(wèi)星在軌運(yùn)行時，各系統(tǒng)的監(jiān)測參數(shù)通過測控操作下達(dá)到地面站，形成大量的遙測數(shù)據(jù)。而地面支撐系統(tǒng)的工作人員會對一些關(guān)鍵性的遙測參數(shù)項(xiàng)進(jìn)行人工監(jiān)視，以確保衛(wèi)星的正常在軌運(yùn)行，如果出現(xiàn)參數(shù)項(xiàng)異常則會啟動相應(yīng)的應(yīng)急措施。隨著在軌運(yùn)行衛(wèi)星數(shù)量和應(yīng)用領(lǐng)域的不斷增多，衛(wèi)星在軌運(yùn)行實(shí)時監(jiān)視工作日益繁重，工作時間長、載荷參數(shù)多是這項(xiàng)工作的顯著特點(diǎn)。以量子科學(xué)實(shí)驗(yàn)衛(wèi)星為例，實(shí)驗(yàn)多為晚上進(jìn)行，工作人員夜間工作較疲勞，此外每顆衛(wèi)星傳回的遙測參數(shù)項(xiàng)非常之多，比如暗物質(zhì)粒子探測衛(wèi)星，其載荷參數(shù)項(xiàng)的個數(shù)多達(dá)兩千多個，人工監(jiān)視往往力不從心。因此采用計(jì)算機(jī)技術(shù)針對衛(wèi)星可能發(fā)生的異常進(jìn)行自動檢測，對于提高衛(wèi)星運(yùn)行控制的工作效率，把人從繁瑣的數(shù)據(jù)中解救出來具有重要的意義[1]。

1 衛(wèi)星遙測參數(shù)的數(shù)據(jù)特性

時間序列分析是根據(jù)系統(tǒng)觀察得到的時間序列數(shù)據(jù)，通過曲線擬合和參數(shù)估計(jì)來建立數(shù)學(xué)模型的理論和方法。時間序列模型常用于國民宏觀經(jīng)濟(jì)控制、市場潛力預(yù)測、氣象預(yù)測、農(nóng)作物害蟲預(yù)報等各個方面。而在航天領(lǐng)域中，衛(wèi)星在軌運(yùn)行期間所產(chǎn)生的大量遙測參數(shù)，如電流、電壓、溫度等，亦可以看作時間序列。

衛(wèi)星在軌運(yùn)行期間，其遙測系統(tǒng)按照一定的采樣周期對星上各部件的工作狀況和傳感器數(shù)值進(jìn)行采集，經(jīng)過A/D變換及編碼后形成相應(yīng)的衛(wèi)星遙測參數(shù)，再經(jīng)過調(diào)制、放大下傳到地面，地面接收后進(jìn)行實(shí)時逆向處理得到這些遙測參數(shù)并入庫，記為xi(t)，其中i=1,2,…,n為衛(wèi)星遙測數(shù)據(jù)的個數(shù)，t為遙測數(shù)據(jù)的采樣時間。衛(wèi)星遙測數(shù)據(jù)可以分為兩大類：數(shù)字量和模擬量[2]。數(shù)字量反映星上被測量單元的功能狀態(tài)，例如衛(wèi)星某單元的開/關(guān)狀態(tài)；模擬量是被測單元的數(shù)值測量值，通常反映被測量單元的性能狀態(tài)。目前常用的遙測參數(shù)異常檢測算法是閾值法，其算法原理為:

f{not(LiT)≤xi(t)≤HiT},then{Alarm}

(1)

其中，HiT和LiT分別為第i項(xiàng)衛(wèi)星遙測參數(shù)數(shù)據(jù)報警的上下閾值?，F(xiàn)在已有的衛(wèi)星監(jiān)視系統(tǒng)正是基于這個上下閾值進(jìn)行數(shù)據(jù)報警，一般采用將數(shù)據(jù)項(xiàng)顯示為紅色的方式，以引起值班人員的注意[3-4]。但是有一些衛(wèi)星遙測參數(shù)總是在門限內(nèi)波動，即使衛(wèi)星的某部件出現(xiàn)異常，其參數(shù)值沒有超出設(shè)定門限，從而造成故障的漏報[5]。通過對衛(wèi)星遙測參數(shù)的深入研究發(fā)現(xiàn)，當(dāng)衛(wèi)星某部件發(fā)生異常時，無論其遙測參數(shù)的數(shù)值是否超過給定的閾值，其變化規(guī)律都會發(fā)生改變，這也就是說如果不發(fā)生異常，每一個衛(wèi)星遙測參數(shù)都有其固有的變化規(guī)律。

通過查閱文獻(xiàn)并結(jié)合量子科學(xué)實(shí)驗(yàn)衛(wèi)星的實(shí)際遙測參數(shù)數(shù)據(jù)，可以把遙測參數(shù)大致分成三大類：第一類參數(shù)是平穩(wěn)型遙測參數(shù)，這一類遙測參數(shù)其數(shù)值隨著時間的變化比較平穩(wěn)緩慢，例如衛(wèi)星的星敏計(jì)數(shù)。第二類參數(shù)是突變型遙測參數(shù)，這一類遙測參數(shù)其數(shù)值隨著時間的變化呈跳躍性變化，例如衛(wèi)星上各部件的開/關(guān)狀態(tài)。第三類參數(shù)是周期性遙測參數(shù)，這一類遙測參數(shù)其數(shù)值隨著時間的變化呈周期性變化，且具有高噪聲，難預(yù)測的特點(diǎn)，例如衛(wèi)星的太陽能帆板外板的溫度、太陽矢量與帆板的夾角[6]。

通過對以上三類遙測參數(shù)數(shù)據(jù)的研究發(fā)現(xiàn)，前兩類遙測參數(shù)的變化規(guī)律比較容易掌握，通過一些時間序列算法甚至一些簡單的數(shù)學(xué)函數(shù)擬合就能很好地預(yù)測它們未來的變化趨勢。根據(jù)這個變化趨勢，就能夠發(fā)現(xiàn)遙測參數(shù)的細(xì)小變化，從而提供一種有效判斷衛(wèi)星異常的新方法。而第三類參數(shù)的變化規(guī)律由于其高噪聲、周期性的特點(diǎn)，相對比較難預(yù)測，這也使得衛(wèi)星周期性遙測參數(shù)預(yù)測成為了當(dāng)前衛(wèi)星遙測數(shù)據(jù)分析領(lǐng)域一大研究熱點(diǎn)。文獻(xiàn)[7]針對衛(wèi)星周期性遙測參數(shù)難以預(yù)測的問題，提出了一種時間序列分解的衛(wèi)星周期性參數(shù)預(yù)測方法。

2 衛(wèi)星遙測參數(shù)的預(yù)測方法

當(dāng)前，國內(nèi)外對于不同環(huán)境中的電子機(jī)械設(shè)備的趨勢分析及其參數(shù)預(yù)測展開了較深入的研究，但是受實(shí)驗(yàn)條件和成本的限制，對于衛(wèi)星這類復(fù)雜環(huán)境中的參數(shù)預(yù)測研究非常有限[8-10]。文獻(xiàn)[11]利用超1 000小時的陀螺儀監(jiān)測數(shù)據(jù)，擬合得到一個線性高階多項(xiàng)式方程來預(yù)測性能趨勢；文獻(xiàn)[12]將預(yù)測區(qū)間技術(shù)用于衛(wèi)星遙測數(shù)據(jù)趨勢預(yù)測；文獻(xiàn)[13]進(jìn)一步使用灰度理論進(jìn)行設(shè)備故障預(yù)測；Gebraeel和Lawley利用神經(jīng)網(wǎng)絡(luò)預(yù)測軸承性能退化及壽命預(yù)測[14]；文獻(xiàn)[15]將支持向量機(jī)算法用于設(shè)備性能退化軌跡建模。對于衛(wèi)星遙測數(shù)據(jù)而言，時間序列法中的ARMA預(yù)測方法是一種常用的方法，既可以用于衛(wèi)星故障預(yù)測[16-17]，也可以針對衛(wèi)星遙測數(shù)據(jù)做短期預(yù)測[18]。除此之外多項(xiàng)式擬合外推，非參數(shù)回歸對于衛(wèi)星遙測參數(shù)的預(yù)測也有較好的效果。針對量子科學(xué)實(shí)驗(yàn)衛(wèi)星的實(shí)際遙測數(shù)據(jù)特點(diǎn)，從已有的算法中選擇了兩種具有代表性的數(shù)據(jù)分析處理方法進(jìn)行研究。

2.1 基于多項(xiàng)式擬合外推的預(yù)測技術(shù)

基于多項(xiàng)式擬合外推的預(yù)測方法是根據(jù)衛(wèi)星在軌運(yùn)行的歷史數(shù)據(jù)和實(shí)時數(shù)據(jù)，尋求在軌衛(wèi)星隨時間推移而發(fā)生變化的規(guī)律，通過對參數(shù)變化曲線進(jìn)行多項(xiàng)式擬合，得到擬合公式，用解析的擬合公式再對數(shù)據(jù)進(jìn)行變化趨勢的預(yù)測。它包括全局和分段多項(xiàng)式擬合兩種不同的處理方式。其中分段多項(xiàng)式擬合預(yù)測，依據(jù)參數(shù)數(shù)據(jù)的周期性特點(diǎn)(如圖1中的量子科學(xué)實(shí)驗(yàn)衛(wèi)星太陽帆板外板溫度)，將原始數(shù)據(jù)分成兩段，將第一段進(jìn)行多項(xiàng)式擬合，得出擬合公式，再將第二段作為函數(shù)關(guān)系的輸入，從而預(yù)測未來第三段數(shù)據(jù)的結(jié)果。全局多項(xiàng)式擬合預(yù)測則是根據(jù)最小二乘法，通過尋找選取數(shù)據(jù)的最佳匹配函數(shù)(即最佳擬合多項(xiàng)式)，并將最佳匹配函數(shù)應(yīng)用到數(shù)據(jù)中進(jìn)行預(yù)測，得到預(yù)測值。實(shí)驗(yàn)發(fā)現(xiàn)，這種情況只適合變化緩慢的衛(wèi)星遙測數(shù)據(jù)，無法擬合數(shù)據(jù)趨勢，因此預(yù)測準(zhǔn)確率較低。圖1是對量子科學(xué)實(shí)驗(yàn)衛(wèi)星太陽帆板外板溫度遙測數(shù)據(jù)進(jìn)行多項(xiàng)式擬合預(yù)測的結(jié)果，表明其短期預(yù)測(半個周期內(nèi))的效果最好，時間過長則會產(chǎn)生過擬合現(xiàn)象，適合周期性隨時間變化緩慢的遙測參數(shù)預(yù)測。

圖1 多項(xiàng)式擬合預(yù)測實(shí)驗(yàn)

2.2 基于自回歸滑動平均模型的預(yù)測技術(shù)

時間序列預(yù)測是預(yù)測方法體系中的重要組成部分，其中ARMA模型是迄今為止理論最為完善的時間序列預(yù)測方法,它是解決非平穩(wěn)趨勢數(shù)據(jù)項(xiàng)預(yù)測問題最常用的模型[19]。自回歸滑動平均模型(auto-regressive and moving average model，ARMA)是一種非平穩(wěn)的時間序列模型。現(xiàn)實(shí)中的序列很多時候并不是平穩(wěn)的，數(shù)學(xué)中一般采用單位根校驗(yàn)法來判斷某序列是否平穩(wěn)，檢驗(yàn)序列存在單位根則表示該序列不平穩(wěn)。研究發(fā)現(xiàn)，通過差分操作可以將絕大多數(shù)非平穩(wěn)序列轉(zhuǎn)變成平穩(wěn)序列，使得非平穩(wěn)序列經(jīng)差分操作之后表現(xiàn)出平穩(wěn)序列的一些性質(zhì)，這樣的序列稱之為差分平穩(wěn)序列[20]。文中后面的研究就是采用ARMA模型預(yù)測方法對一些符合差分平穩(wěn)性質(zhì)的衛(wèi)星遙測參數(shù)進(jìn)行預(yù)測。

ARMA模型預(yù)測方法又稱博克斯-詹金斯法，比較適合處理復(fù)雜、具有多種模式的時間序列。模式中可包含趨勢變化、季節(jié)變化、隨機(jī)變化、循環(huán)變化等因素的綜合或單一影響。ARMA模型預(yù)測的過程可描述如下：

(1)獲取被觀測系統(tǒng)時間序列數(shù)據(jù)；

(2)利用檢驗(yàn)統(tǒng)計(jì)量(如Q統(tǒng)計(jì)量、LB統(tǒng)計(jì)量)對序列進(jìn)行白噪聲檢驗(yàn)；

(3)判斷序列是否存在單位根，對序列進(jìn)行平穩(wěn)性檢驗(yàn)，若不平穩(wěn)，則進(jìn)行d階差分直至平穩(wěn)；

(4)繪制自相關(guān)圖(ACF)和偏自相關(guān)圖(PACF)創(chuàng)建ARMA預(yù)測模型；

(5)對得到的ARMA預(yù)測模型進(jìn)行殘差分析，確認(rèn)模型的有效性；

(6)利用得到的有效ARMA模型對遙測數(shù)據(jù)進(jìn)行預(yù)測。

在對衛(wèi)星實(shí)際的遙測參數(shù)數(shù)據(jù)進(jìn)行ARMA建模的過程中發(fā)現(xiàn)，該方法適合處理周期變化的參數(shù)和遙測參數(shù)平穩(wěn)變化或維持在某恒定值附近的參數(shù)預(yù)測。

3 基于時間序列模型的衛(wèi)星異常檢測方法

前面已經(jīng)提到在實(shí)際的衛(wèi)星在軌運(yùn)行任務(wù)中，大量的工作是通過人工完成的，比如衛(wèi)星在軌運(yùn)行實(shí)時監(jiān)視工作，由于工作時間長、工作強(qiáng)度大以及衛(wèi)星遙測參數(shù)項(xiàng)過多，人工操作往往有點(diǎn)力不從心。目前普遍采用閾值報警法，但是許多異常產(chǎn)生并不會超出閾值的上下界限，導(dǎo)致異常的漏報。如何使用計(jì)算機(jī)技術(shù)盡可能地對異常進(jìn)行自動檢測，減少人工成本、減少異常漏報的概率，已經(jīng)成為當(dāng)下衛(wèi)星在軌運(yùn)行任務(wù)的亟待解決的問題。

文中針對衛(wèi)星遙測數(shù)據(jù)的時間序列性、多樣性、高噪聲等特點(diǎn)，提出了一種基于時間序列模型的衛(wèi)星異常檢測方法。其框架如圖2所示。

圖2 基于時間序列模型的衛(wèi)星異常檢測模型

如圖2所示，為了解決衛(wèi)星遙測參數(shù)項(xiàng)過多的問題，采用分治的思想，先將經(jīng)過預(yù)處理后的衛(wèi)星遙測參數(shù)數(shù)據(jù)進(jìn)行分類，然后針對不同類別的衛(wèi)星遙測參數(shù)數(shù)據(jù)建立相應(yīng)的預(yù)測模型，每個模型輸出一個預(yù)測結(jié)果，將預(yù)測結(jié)果同新產(chǎn)生的遙測數(shù)據(jù)進(jìn)行比較，若當(dāng)前時刻某參數(shù)項(xiàng)的數(shù)值與預(yù)測結(jié)果不在合理范圍內(nèi)波動，則表示衛(wèi)星當(dāng)前狀態(tài)可能發(fā)生異常，系統(tǒng)則會給值班人員發(fā)出預(yù)警消息。當(dāng)新的數(shù)據(jù)入庫后，模型會實(shí)現(xiàn)增量的更新預(yù)測模型，以適應(yīng)新的遙測參數(shù)變化規(guī)律。為了不影響正常的衛(wèi)星運(yùn)行任務(wù)，模型更新的時間，一般會選擇沒有數(shù)傳和遙測事件的時間。

4 模型實(shí)現(xiàn)與驗(yàn)證

模型采用Python語言進(jìn)行衛(wèi)星遙測數(shù)據(jù)的處理，并將擬采用的時間序列算法(ARMA、多項(xiàng)式擬合外推)進(jìn)行編程實(shí)現(xiàn)。以暗物質(zhì)粒子探測衛(wèi)星和量子科學(xué)實(shí)驗(yàn)衛(wèi)星真實(shí)的遙測數(shù)據(jù)為數(shù)據(jù)源，進(jìn)行模型訓(xùn)練和實(shí)驗(yàn)驗(yàn)證。通過對特定遙測參數(shù)的隨時間變化規(guī)律的挖掘，來實(shí)現(xiàn)衛(wèi)星異常自動檢測，并在第一時間發(fā)出報警信息。采用對比實(shí)驗(yàn)的方法，同一遙測參數(shù)用不同的算法進(jìn)行實(shí)驗(yàn)，同一算法對不同的遙測參數(shù)進(jìn)行實(shí)現(xiàn)，最后對算法進(jìn)行組合改進(jìn)，構(gòu)建最優(yōu)的預(yù)測模型。

以衛(wèi)星星敏復(fù)位這個異常為例，暗物質(zhì)粒子探測衛(wèi)星自2015年12月發(fā)射升空以來，發(fā)生過好幾次星敏復(fù)位的事件，盡管事后都得到了及時的處理，但是人工監(jiān)視有它的弊端，容易造成疏忽漏報的情況出現(xiàn)，因此采用多項(xiàng)式外推的方式能很好地自動檢測出星敏感器復(fù)位事件，并在第一時間向相關(guān)負(fù)責(zé)人發(fā)出警報信息。其時間序列趨勢如圖3所示。

圖3 暗物質(zhì)粒子探測衛(wèi)星星敏感器APPS計(jì)數(shù)

此外對于量子科學(xué)實(shí)驗(yàn)衛(wèi)星Y帆板外板溫度，重點(diǎn)采用了ARMA算法進(jìn)行實(shí)驗(yàn)研究。選取85個時間單位(約5小時)的衛(wèi)星歷史遙測數(shù)據(jù)參照ARMA模型訓(xùn)練過程，預(yù)測未來15個時間單位的數(shù)據(jù)(約1小時)，經(jīng)過前期的數(shù)據(jù)預(yù)處理、序列平穩(wěn)性校驗(yàn)、序列白噪聲校驗(yàn)，然后根據(jù)ACF與PACF圖確認(rèn)ARMA模型的p、q參數(shù)，再通過殘差分析確認(rèn)ARMA模型的有效性，最后使用最終確定的ARMA(3,2)模型對未來15個時間單位的溫度數(shù)據(jù)進(jìn)行預(yù)測，取得了較好的實(shí)驗(yàn)結(jié)果。

預(yù)測效果如圖4所示，其中實(shí)線表示真實(shí)數(shù)據(jù)，虛線表示預(yù)測數(shù)據(jù)。

圖4 基于ARMA模型對衛(wèi)星遙測數(shù)據(jù)的預(yù)測實(shí)驗(yàn)

圖3中的星敏感器計(jì)數(shù)，是衛(wèi)星在軌運(yùn)行實(shí)時監(jiān)視工作重點(diǎn)關(guān)注的參數(shù)，利用時間序列模型對其進(jìn)行趨勢預(yù)測，實(shí)現(xiàn)異常的自動檢測，可以有效地輔助人工判斷，提高衛(wèi)星在軌任務(wù)運(yùn)行的工作效率，同時也降低了異常的漏報率。而圖4中的Y帆板溫度數(shù)據(jù)，在衛(wèi)星在軌運(yùn)行實(shí)時監(jiān)視工作期間工作人員并沒有特別的關(guān)注，如果出現(xiàn)閾值內(nèi)的溫度異常，人工是很難發(fā)現(xiàn)的，因此通過時間序列模型對其進(jìn)行趨勢預(yù)測，可以及時檢測到異常，并發(fā)出報警信息。

5 結(jié)束語

主要討論了基于衛(wèi)星遙測數(shù)據(jù)的時間序列特性，將時間序列模型應(yīng)用于衛(wèi)星異常檢測。以暗物質(zhì)粒子探測衛(wèi)星和量子科學(xué)實(shí)驗(yàn)衛(wèi)星在軌運(yùn)行一年多的時間里所產(chǎn)生的遙測數(shù)據(jù)為基礎(chǔ)，結(jié)合中科院空間科學(xué)先導(dǎo)專項(xiàng)實(shí)際空間科學(xué)衛(wèi)星運(yùn)控任務(wù)背景，對衛(wèi)星有效載荷的異常檢測進(jìn)行實(shí)驗(yàn)研究，提出了一種基于時間序列模型的衛(wèi)星有效載荷異常檢測方法。列出了兩個典型的應(yīng)用案例，分析了該方法在實(shí)際的衛(wèi)星在軌運(yùn)行任務(wù)中的應(yīng)用價值。

然而，衛(wèi)星遙測數(shù)據(jù)的高噪聲、難預(yù)測等特點(diǎn)使得該方法仍存在一定的局限性，不能夠?qū)崿F(xiàn)高準(zhǔn)確率的預(yù)測，因此在實(shí)際的任務(wù)場景中的應(yīng)用會受到一定的限制。未來還需要開展更多領(lǐng)域更深層次的研究，比如結(jié)合測控事件構(gòu)建物理模型進(jìn)行參數(shù)預(yù)測或者引入新的預(yù)測算法。