王雪瑩 李 駿 安 瑋

國防科技大學(xué)電子科學(xué)與工程學(xué)院，長沙 410073

地球靜止軌道(Geostationary Orbit,GEO) 是特殊的地球同步軌道(Geosynchronous Orbit),軌道高度為35786km，軌道周期與地球自轉(zhuǎn)周期(23h56m4s)相同，因而在通信、導(dǎo)航等領(lǐng)域有其獨特的優(yōu)勢。但由于GEO上可容納的航天器數(shù)量是有限的[1]，因而這一軌道就更顯“擁擠”，而這一區(qū)域的空間碎片不但占據(jù)了寶貴的GEO 軌道資源，還嚴(yán)重影響了GEO衛(wèi)星的安全[2]，加強對GEO附近天區(qū)(軌道高度(35786±75)km，軌道傾角±15°)空間目標(biāo)尤其是碎片的監(jiān)視與編目，能夠有效利用GEO軌道，也可為在軌工作的航天器規(guī)避碰撞提供有效引導(dǎo)信息。

由于距離等約束，地基雷達對GEO目標(biāo)的觀測大多限于直徑大于1m的目標(biāo)[3]，對于尺寸較小的目標(biāo)，只能通過天基平臺觀測，而天基觀測平臺因其自身工作模式等原因，單個平臺在單周期之內(nèi)對目標(biāo)僅有數(shù)十秒到數(shù)分鐘可觀測時間[4]，如何有效利用短弧觀測成為一個關(guān)鍵問題，李駿將Milani等人[5]研究太陽系小行星運動時對短弧測角數(shù)據(jù)處理方法引入天基空間目標(biāo)監(jiān)視信息處理中，提出了一種合理應(yīng)用短弧測角數(shù)據(jù)的方法，但對新目標(biāo)的短弧測角數(shù)據(jù)而言，其處理結(jié)果仍無法直接引導(dǎo)平臺觀測。

劉光明、潘曉剛等人[6-7]參照衛(wèi)星導(dǎo)航精度衰減因子的定義[8],引入了幾何精度因子(Geometric of Precision, GOP)的概念,證明了定軌精度一定程度上依賴于幾何精度因子，對幾種極限情況做了解釋，并仿真分析了定軌精度與幾何精度因子的關(guān)系，但未對多平臺多測量弧段GOP做出定義與分析。

本文針對天基光學(xué)監(jiān)視系統(tǒng)對GEO目標(biāo)可觀測時間較短的特點，利用標(biāo)準(zhǔn)GEO高度作為約束條件，以短弧觀測內(nèi)GEO目標(biāo)做勻速運動為基本假設(shè)，解算觀測平臺位置、觀測角度與距地高度為35786km的空間球面構(gòu)成的方程，取合理解作為空間目標(biāo)的位置，再用最小二乘法擬合出目標(biāo)的位置、速度，以此為初始值利用J2的攝動模型外推，預(yù)報目標(biāo)軌道，引導(dǎo)天基平臺，對目標(biāo)進行后續(xù)觀測。并用廣義Laplace算法[9-10]對多段觀測弧段處理，得到有效初軌，為評價多弧段觀測數(shù)據(jù)定軌的觀測幾何，提出一種修正的幾何精度因子(Geometric of Precision, GOP)概念，為描述測量誤差對定軌結(jié)果的影響程度，提出幾何精度因子的誤差靈敏度(Error Sensitivity Geometric of Precision, ESGOP)概念。

仿真結(jié)果表明短弧測量條件下，利用約束條件求出的目標(biāo)位置、速度可以有效引導(dǎo)短期內(nèi)對目標(biāo)的觀測，也可作為初軌算法的初值輸入，GOP可用來表征系統(tǒng)的可觀測特性，ESGOP可以反映觀測條件對測量誤差的靈敏程度，GOP和ESGOP從本質(zhì)上反映了觀測幾何與定軌精度的關(guān)系。

1 天基光學(xué)觀測模型

1.1 觀測幾何關(guān)系

在J2000.0地心慣性坐標(biāo)系下，設(shè)地心原點為O，R為測站Δt的位置矢量，r為目標(biāo)S的位置矢量，而ρ為S相對測站的觀測矢量，如下圖所示。

圖1 觀測幾何關(guān)系

則有r=ρ+R

(1)

其中測站矢量R已知，觀測矢量ρ=ρL，L為單位觀測矢量，由目標(biāo)赤經(jīng)α和赤緯δ確定

(2)

1.2 測量模型

赤經(jīng)α和赤緯δ由J2000.0坐標(biāo)系下目標(biāo)位置矢量r=(x,y,z)和測站位置矢量R=(X,Y,Z)確定

(3)

2 幾何約束下的GEO目標(biāo)定位

,j=1,…,M

平臺的位置矢量為Rj=(Xj,Yj,Zj),j=1,…,M。

在單個天基光學(xué)平臺短弧觀測中，GEO目標(biāo)近似做勻速直線運動，像在平面上的軌跡一般近似直線[4]，通過在視場內(nèi)標(biāo)定標(biāo)準(zhǔn)同步軌道目標(biāo)軌跡的方式，可迅速判斷出同步帶目標(biāo)。

構(gòu)造如下變量：

由最小二乘法的原理，可知

(4)

其中，R表示位置誤差的協(xié)方差矩陣，PX0表示估計誤差的協(xié)方差。

3 廣義Laplace初軌算法可觀測性分析

利用文獻[7-8]中算法對多段觀測數(shù)據(jù)處理，采用消元后的基本方程組

，(j=1,2,3,…，M)

(5)

(6)

上式表示線性化后系數(shù)矩陣的子塊。

對于完整的廣義Laplace方程而言，表征其可觀測性的法方程系數(shù)矩陣很難給出解析的表達式[11]，只能在一定假設(shè)條件下通過近似解釋一些現(xiàn)象，為了定量描述系統(tǒng)可觀測性，在前述的觀測模型和測量集的基礎(chǔ)上，定義2個變量。

1)幾何精度因子(Geometric of Precision, GOP)：多平臺(或單平臺)對空間目標(biāo)多個弧段的觀測數(shù)據(jù)中，任取三次角度觀測，觀測平臺指向空間目標(biāo)的單位視線矢量所構(gòu)成的四面體體積可用角度測量表示，遍歷所有的三次角度測量組合，取所有四面體體積的最大值為幾何精度因子，其表征的是某種特定觀測幾何下的系統(tǒng)可觀測性。

2)幾何精度因子的誤差靈敏度(Error Sensitivity Geometric of Precision, ESGOP)：由實際測角信息計算的幾何精度因子與平滑處理后的測角值[12-13](或理論測角值)計算的幾何精度因子的相對誤差。在特定的觀測幾何下，誤差靈敏度越小，測量誤差對定軌結(jié)果的影響越小。

|H|=2a2/3×Mm×Mn×Mk×GOPm,n,k×

vm×vn×vk

(7)

其中a為目標(biāo)軌道半長軸，ΔE為偏近點角變化量，

·r0·ΔEk-cosΔEm)-

a·sin(ΔEk-ΔEm)-(a-r0)(sinΔEk-sinΔEm)

a·sin(ΔEn-ΔEk)-(a-r0)(sinΔEk-sinΔEn)

從定義可得

不難看出其可作為衡量系統(tǒng)可觀測性的一項指標(biāo)。若將平滑處理后的測角值(或理論測角值)計算的幾何精度因子記為GOPopt，則有

×100%

4 仿真結(jié)果與分析

為仿真驗證本文算法可獲取目標(biāo)有效初軌，并驗證GOP和ESGOP可有效表征初軌算法性能，設(shè)置如表1所示的仿真場景。

表1 仿真參數(shù)

仿真過程中相機主要采用恒星跟蹤模式[14-15]，針對幾何約束下的GEO目標(biāo)定位算法和廣義Laplace初軌算法可觀測性做如下仿真。

4.1 仿真一

觀測星1的相機，在某次任務(wù)中，發(fā)現(xiàn)了GEO區(qū)的新目標(biāo)，以1s每點的頻率輸出200個角度測量值，測角誤差為100微弧度。利用(4)式估計初始時刻目標(biāo)的位置與速度，再利用J2攝動模型外推2h的軌道與真實軌道對比，得到如圖2所示結(jié)果。

從仿真結(jié)果可以看出，利用本文算法處理200s的觀測數(shù)據(jù)，可有效估計出目標(biāo)在參考時刻的位置速度，此值可以用作初軌算法的輸入，利用估計值預(yù)測目標(biāo)軌道，在1h內(nèi)位置誤差不超過80km，可有效引導(dǎo)其它平臺對目標(biāo)進行觀測。

圖2 預(yù)測軌道

4.2 仿真二

在可觀測的弧段起始處與結(jié)束處均安排平臺1對目標(biāo)進行2次100s的觀測，得到的數(shù)據(jù)稱為單平臺數(shù)據(jù)；在可觀測的弧段起始處與結(jié)束處分別安排平臺1和平臺2對目標(biāo)進行2次100s的觀測，得到的數(shù)據(jù)稱為多平臺數(shù)據(jù)(可觀測弧段不足100s時，以可觀測弧段長度為觀測弧長)。100微弧度測角誤差下，兩組數(shù)據(jù)對比得到如圖3仿真結(jié)果。

從圖3的仿真結(jié)果可看出，在觀測弧段較短的情況下(Δt小于500s)，多平臺觀測的GOP及ESGOP明顯優(yōu)于單平臺觀測，因而其定位精度更高，但是隨著可觀測弧段增長，這種優(yōu)勢也漸漸變?nèi)?。這主要是由于多平臺觀測的優(yōu)勢集中體現(xiàn)在其良好的觀測幾何中，如果單平臺有足夠長的觀測弧段，其觀測幾何也比較好，因而其定軌性能也會提升，亦即使用多平臺多觀測弧段對GEO目標(biāo)定軌，可將空間優(yōu)勢轉(zhuǎn)換為觀測弧長優(yōu)勢。

5 結(jié)束語

本文提出一種利用天基光學(xué)短弧觀測對GEO區(qū)目標(biāo)定位的方法，在GEO區(qū)半徑、偏心率等約束條件下，僅使用數(shù)百秒觀測即可估計出目標(biāo)的位置與速度，此結(jié)果可以用作初軌算法的初始值，也可以其為初值進行軌道預(yù)測，數(shù)小時內(nèi)預(yù)測結(jié)果可有效引導(dǎo)其它平臺對目標(biāo)觀測。為表征觀測幾何對廣義Laplace初軌算法定軌性能的影響，將幾何精度因子(GOP)的概念擴展至多平臺多觀測弧段，并引入幾何精度因子的誤差靈敏度(ESGOP)表征測量誤差對觀測幾何的影響，仿真結(jié)果表明：對GEO目標(biāo)而言，當(dāng)GOP不小于0.03、ESGOP不大于0.005時，對目標(biāo)的初軌精度可以控制在50km之內(nèi)。

圖3 不同觀測弧長下單平臺與多平臺數(shù)據(jù)定軌性能比較

參 考 文 獻

[1] 李恒年, 高益軍, 余培軍, 等.地球靜止軌道共位控制策略研究[J] .宇航學(xué)報, 2009,30 (3) : 967-973.(LI Heng-nian,GAO Yi-jun, YU Pei-jun, et al.The Strategies and Algorithms Study for Multi-GEO Satellites Collocation [J].Journal of Astronautics,2009, 30(3): 967-973.)

[2] 梁斌, 徐文福, 李成, 等.地球靜止軌道在軌服務(wù)技術(shù)研究現(xiàn)狀與發(fā)展趨勢[J].宇航學(xué)報,2010,31(1):1-13.(LIANG Bin, XU Wen-fu, LI Cheng, et al.The Status and Prospect of Orbital Servicing in the Geostationary Orbit[J].Journal of Astronautics, 2010,31(1):1-13.)

[3] 唐軼峻, 姜曉軍, 魏建彥, 等.高軌空間碎片光電觀測技術(shù)綜述[J].宇航學(xué)報, 2008, 29 (4) : 1094-1098.(TANG Yi-jun,JIANG Xiao-jun, WEI Jian-yan, et al .Review of Optical Observations of High Apogee Space Debris[J].Journal of Astronautics, 2008,29 (4):1094-1098.)

[4] 李駿，安瑋，周一宇.天基光學(xué)短弧初軌的約束微分修正方法[J].宇航學(xué)報,2009,30(2):669-774.(Li Jun, AN Wei, ZHOU Yi-yu. Constrained Differential Correction in Initial Orbit Determination with Short-arcs in Optical Space-based Space Surveillance[J].Journal of Astronautics, 2009,30 (2) :669-774.)

[5] A.Milani and G.Gronchi .Theory of Orbit Determination[M].Cambridge: Cambridge University Press,2010.

[6] 潘曉剛.空間目標(biāo)定軌的模型與參數(shù)估計方法研究及應(yīng)用[D].長沙: 國防科技大學(xué), 2009.(Pan Xiao-gang.Research on Space Target Orbit Determination Model and Paraments Estimation Algorithms[D].Changsha: National University of Defense Technology, 2009.)

[7] 劉光明.基于天基測角信息的空間非合作目標(biāo)跟蹤算法及相關(guān)技術(shù)研究[D].長沙: 國防科技大學(xué), 2011.(Liu Guangming.Research on non-cooperative Target Tracking Algorithms and Related Technologies Using Space-based Bearings-only Measurements [D].Changsha: National University of Defense Technology, 2011.)

[8] Cai Z W.Research on Autonomous Orbit Determination of Navigation Satellite Based on Crosslink Range and Orientation Parameters Constraining[J].Geo-spatial Information of Science, 2006, 9(1):18-23.

[9] 劉林.航天器軌道理論[M].北京:國防工業(yè)出版社, 2000.

[10] 劉林.關(guān)于初軌計算[J].飛行器測控學(xué)報，2004, 23(3): 41-50.(Liu Lin.On Initial Orbit Determination[J].Journal of Spacecraft TT&C Technology,2004, 23(3): 41-50.)

[11] 李駿.空間目標(biāo)天基光學(xué)監(jiān)視跟蹤關(guān)鍵技術(shù)研究[D].長沙: 國防科技大學(xué), 2009.(Li Jun.Research on Technologies of Space Objects Surveillance and Tracking in Space-based Optical Surveillance [D].Changsha: National University of Defense Technology, 2009.)

[12] 張玉祥.人造衛(wèi)星測軌方法[M].北京:國防工業(yè)出版社, 2007.

[13] M.Ju.Sokolaskaya.On the Laplacian Orbit Determination of Asteroids[J].Planet.Space Sci., 1998, 45(12): 1575-1580.

[14] Grant H, Stokes, Curt von Braun, et al.Space-based Visible Program[J].Lincoln Laboratory Journal,1998,11(2):205-229.

[15] Jayant Sharma, Grant H, Stokes, et al.Toward Operational Space-based Space Surveillance[J].Lincoln Laboratory Journal,2002,13 (2):309-334.