汪榮峰,陳大春

(中國人民解放軍戰(zhàn)略支援部隊航天工程大學,北京 101416)

以多顆衛(wèi)星或星座對作戰(zhàn)區(qū)域進行偵察,獲取特定目標和移動目標信息,可對區(qū)域進行全時空、全方位監(jiān)控,提高戰(zhàn)場態(tài)勢感知能力和指揮作戰(zhàn)效率[1].衛(wèi)星有存儲轉發(fā)和實拍實傳2種數(shù)傳方式[2],實拍實傳數(shù)傳方式是即時數(shù)據(jù)傳輸方式,衛(wèi)星在利用星載偵察設備對地面目標進行偵察的同時,將所獲數(shù)據(jù)下傳至地面站.當采用多顆實拍實傳衛(wèi)星對目標區(qū)域進行偵察時,需進行數(shù)傳任務規(guī)劃、能力分析和仿真推演.

衛(wèi)星數(shù)傳調度問題研究的是獨立于衛(wèi)星數(shù)據(jù)獲取的數(shù)據(jù)下載決策問題,主要目的是為數(shù)傳任務分配無沖突的數(shù)傳資源(地面站數(shù)據(jù)接收天線)及時間窗口.數(shù)傳任務規(guī)劃已有大量研究成果,如劉洋等[3]提出的基于貪婪算法的衛(wèi)星地面站任務規(guī)劃算法,陳祥國等[4]提出的基于蟻群算法的規(guī)劃技術.在實拍實傳衛(wèi)星數(shù)傳調度問題方面,文獻[2]設計了一個啟發(fā)式調度算法,姜維等[5]針對存儲轉發(fā)和實拍實傳數(shù)傳模式建立統(tǒng)一模型,經飛等[6]提出的多數(shù)傳模式組合調度算法可同時處理存儲轉發(fā)和實拍實傳需求.

為衡量衛(wèi)星完成偵察任務的能力,一般以覆蓋率作為核心指標.現(xiàn)有衛(wèi)星區(qū)覆蓋分析方法主要有解析法[7]、網格點法[8]和基于幾何運算的方法[9].

王慧林等[10]設計實現(xiàn)了面向區(qū)域的電子偵察衛(wèi)星規(guī)劃系統(tǒng),采用甘特圖表現(xiàn)規(guī)劃序列,基于STK顯示規(guī)劃方案執(zhí)行過程.

本文根據(jù)實際任務所需研制規(guī)劃推演系統(tǒng),提出了新的數(shù)傳調度算法、衛(wèi)星區(qū)域覆蓋分析算法,并設計實現(xiàn)了過程推演表現(xiàn)形式和可視化方法.

1 規(guī)劃推演系統(tǒng)設計

系統(tǒng)組成如圖1所示,包括4個軟件.

指揮控制軟件是區(qū)域偵察任務的邏輯起點,根據(jù)用戶需求交互確定偵察時段、執(zhí)行偵察任務衛(wèi)星、測控站和待偵察目標區(qū)域,并將上述信息(圖1中 “任務”)寫入到共享數(shù)據(jù)庫中,供其他軟件使用.

數(shù)傳任務規(guī)劃軟件根據(jù)“任務”進行規(guī)劃,以XML文件形式將結果提供給過程推演軟件.

本研究結果表明，當合成麝香濃度較低時，合成麝香初始濃度的增加一定程度上可以促進 PAC對其的去除。當合成麝香初始濃度繼續(xù)增大時，PAC對合成麝香的去除量仍呈增加趨勢，但是去除率無明顯變化，說明合成麝香的初始濃度的提高可以提高PAC的有效利用率。

偵察能力分析軟件以給定時間段內衛(wèi)星對目標區(qū)域的覆蓋率作為核心指標,分析任務完成能力.

過程推演軟件基于OSGEarth開源引擎開發(fā),根據(jù)“任務”信息和數(shù)傳規(guī)劃結果,對偵察和數(shù)傳的過程進行仿真推演.

系統(tǒng)主要包括如下模塊:數(shù)據(jù)庫與數(shù)據(jù)文件訪問模塊,支持各軟件通過共享數(shù)據(jù)庫和XML數(shù)據(jù)文件交換信息;軌道預推模塊,支持標準二體和SGP4兩種軌道預推模型;時間窗口計算模塊,為數(shù)傳規(guī)劃算法提供支持;幾何算法,采用開源幾何引擎GEOS (Geometry Engine-Open Source),為區(qū)域覆蓋分析和三維繪制提供支持;區(qū)域覆蓋分析和數(shù)傳規(guī)劃算法,應用軟件的核心算法;圖形引擎,包括OpenGL和OpenSceneGraph(OSG),分別用于規(guī)劃、分析軟件二維繪制和過程推演軟件三維繪制.

2 面向區(qū)域偵察的數(shù)傳任務規(guī)劃算法

多顆實拍實傳衛(wèi)星區(qū)域偵察數(shù)傳規(guī)劃的任務是在區(qū)域約束下計算衛(wèi)星對目標的過境時間窗口,并為每個過境時間窗口分配一個或多個地面站天線作為數(shù)傳資源.

衛(wèi)星可數(shù)傳的基本條件是衛(wèi)星與地面站相互“可見”,且滿足一定角度要求.此外還需考慮衛(wèi)星數(shù)傳能力、地面站天線個數(shù)與仰角、衛(wèi)星與地面站天線頻段、地面站天線跟蹤時間等約束條件.根據(jù)任務要求做如下假設:每個地面站部署一天線,頻段與衛(wèi)星天線一致;衛(wèi)星與地面站滿足給定仰角方可進行數(shù)傳;天線對衛(wèi)星的捕獲時間、準備時間等合在一起,為“任務切換時間”.

針對任務中衛(wèi)星數(shù)量和地面站數(shù)量較少的實際,設計如下簡潔、易實現(xiàn)的規(guī)劃算法:

Step 1.地面站天線的初始時間資源為整個時間軸;

Step 2.計算衛(wèi)星對目標區(qū)域所有過境時間窗口,時間窗口計算采用文獻[11]的二次掃描算法;

Step 4.按順序取出待分配任務列表中的每個時間窗口,計算與各地面站具備數(shù)傳條件(仰角)的細分時間窗口,進行資源分配(考慮任務切換時間),同時修改地面站時間資源.

以圖2為例簡要說明上述過程的第4步:設開始有2個地面站,其時間資源均為全部時間范圍,如圖中a、b所示;有一時間窗口存在數(shù)傳需求,如圖中AD區(qū)間,根據(jù)與地面站關系,細分為3個窗口,AB段與地面站a具備數(shù)傳條件,CD段與地面站b具備數(shù)傳條件,BC段與2地面站同時具備數(shù)傳條件;AC段由地面站a進行數(shù)傳,CD段由地面站b進行數(shù)傳;修改相應地面站時間資源,減去已分配時間窗口,如圖2中a'、b'所示.

開發(fā)的數(shù)傳任務規(guī)劃軟件界面如圖3所示.其中紅色波紋扇形區(qū)表示規(guī)劃結果,對應的衛(wèi)星軌道時間范圍進行數(shù)傳.

3 基于網格的衛(wèi)星區(qū)域覆蓋分析算法

為快速計算多顆衛(wèi)星在給定時間范圍內對目標區(qū)域的覆蓋率,本文提出了一種基于剖分網格的新算法:將目標區(qū)域網格化,判斷各網格與衛(wèi)星覆蓋帶相交情況,以有交網格數(shù)除以全部網格數(shù)作為總覆蓋率,通過分類統(tǒng)計得到分類覆蓋率指標.

如圖4所示,目標區(qū)域共24個網格,衛(wèi)星覆蓋帶為abcd和ABCD.與abcd有交的網格為11、21、31、22、32、23、33、43、24、34、44,與ABCD有交的網格為41、32、42、33、43、34,被覆蓋網格共13個,網格總數(shù)為24,覆蓋率為13/24=54%.網格34、33、43、32同時被2顆衛(wèi)星覆蓋,其他被1顆衛(wèi)星覆蓋,據(jù)此可得到其他分類統(tǒng)計信息.

圖4 基于網格的衛(wèi)星區(qū)域覆蓋分析方法示意圖

根據(jù)經緯度平面網格大小,可進行均勻或非均勻網格劃分.均勻網格劃分時,網格對應的實際地理范圍大小不一致.非均勻網格劃分則按對應地理范圍近似一致的原則,不同緯度網格大小不同.

劃分網格后,需判斷各網格與衛(wèi)星覆蓋帶是否相交.逐網格判斷效率過低,但由于各行網格大小不一定相同,因此也無法整體判斷.故采用逐行判斷技術:遍歷覆蓋帶所有邊,計算邊與該行網格2水平邊所有交點;水平邊上2交點之間的網格與覆蓋帶相交.

基于覆蓋分析算法,開發(fā)偵察能力分析軟件,設計分析結果的可視化輸出方式:以柱狀圖形式輸出按衛(wèi)星、分辨率等統(tǒng)計的覆蓋結果,圖5(a)為按衛(wèi)星的柱狀統(tǒng)計圖,表現(xiàn)各衛(wèi)星及衛(wèi)星組合的覆蓋網格數(shù)與覆蓋率.

另外,以遙感影像為底圖顯示覆蓋情況:繪制目標區(qū)域的剖分網格,以不同顏色表示網格的覆蓋情況,圖5(b)為針對次數(shù)的覆蓋情況可視化,紅色網格被覆蓋1次,綠色網格被覆蓋2次,透明網格未被覆蓋.

4 目標區(qū)域和偵察范圍組合可視化方法

實拍實傳衛(wèi)星區(qū)域偵察過程推演軟件的任務是:根據(jù)指揮控制軟件所確定的地面站、衛(wèi)星、地面目標區(qū)域、偵察時段等信息,以及數(shù)傳規(guī)劃軟件給出的數(shù)傳資源與時間窗口,在三維場景中將衛(wèi)星區(qū)域偵察的要素以及偵察、數(shù)傳過程表現(xiàn)出來.采用OSGEarth為開發(fā)引擎,OSGEarth是基于OSG開發(fā)的地形模型加載和渲染工具.圖6為推演軟件顯示效果圖.

圖6中,三維地球、A處衛(wèi)星、B處地面站、C處空間鏈路、D處圓錐傳感器覆蓋、E處矩形傳感器覆蓋等可視化要素,OSGEarth可直接支持.

推演系統(tǒng)中最核心的2項可視化要素是F處大范圍半透明多邊形所代表的目標區(qū)域和H、G處所表現(xiàn)的隨時間動態(tài)變化的已偵察范圍.OSGEarth引擎無法直接支持,需設計專門的算法來實現(xiàn).

目標區(qū)域是定義在經緯度平面的多邊形,已偵察范圍是衛(wèi)星覆蓋帶與目標區(qū)域的交.三維場景中進行可視化時,一是需將其轉化為大地坐標系中的曲面,二是避免繪制共面多邊形產生的顯示混亂.

因此設計目標區(qū)域與已偵察范圍的組合繪制算法:首先在經緯度平面將目標區(qū)域離散為可作為曲面表達的小多邊形;然后當衛(wèi)星位置改變時,計算覆蓋帶與每個離散多邊形的交與差,交表示已偵察范圍,差表示目標區(qū)域;最后進行繪制.

圖5 實拍實傳衛(wèi)星區(qū)域偵察能力分析軟件輸出結果

圖6 衛(wèi)星區(qū)域偵察過程推演顯示效果

(1)目標區(qū)域離散

目標區(qū)域一般采用經緯度平面的多個點定義,但其實際形狀為貼合于地球表面的曲面,必須離散為多邊形表示,采用如下技術:在經緯度平面計算目標區(qū)域的矩形包圍盒,并根據(jù)預設離散精度將包圍盒范圍離散為均勻網格,如圖7所示,目標區(qū)域為ABCD,經緯度平面包圍盒離散為24個均勻網格;目標區(qū)域與各網格的交作為離散結果,如網格32位于目標區(qū)域內部,離散結果為矩形defg,網格51與目標區(qū)域部分相交,離散結果為三角形abc.

圖7 目標區(qū)域與偵察范圍組合可視化方法

(2)生成覆蓋多邊形

根據(jù)衛(wèi)星位置、載荷類型、視場角等,每個衛(wèi)星采樣點生成2個覆蓋帶邊界點,將邊界點按順序組織在一起,構成覆蓋帶多邊形.為優(yōu)化效率,采用增量式計算方法:如新形成2點所形成區(qū)域與目標區(qū)域無交,不必繼續(xù)處理;否則,或創(chuàng)建新的覆蓋帶多邊形(原來無覆蓋多邊形),或更新原覆蓋多邊形.

(3)計算覆蓋多邊形與區(qū)域離散結果的交與差

交多邊形作為已偵察區(qū)域繪制用多邊形,差多邊形作為目標區(qū)域繪制用多邊形.以圖7中多邊形defg為例說明,覆蓋帶多邊形與該多邊形相交的部分是v0v1v2v3v4p4p3p2p1p0,交點為h、i、j、k,以多邊形hiv1v2v3jkp3p2p1p0表示已偵察區(qū)域,以多邊形iv1v2v3jfe和kp3p2p1hdg表示目標區(qū)域.

(4)繪制

上述交、差多邊形需轉換為大地坐標進行繪制.轉化后多邊形與地表并不完全貼合,即使完全貼合,采用多邊形偏移繪制技術也會導致消隱不正確、時現(xiàn)時隱.采用解決方法為:根據(jù)經緯度計算大地坐標時增加固定高程,使其高于地表;為避免視點貼近地表時觀察到懸空效果,交互操作限制視點高度.

5 結束語

研制的實拍實傳衛(wèi)星規(guī)劃推演系統(tǒng)是更大規(guī)模信息系統(tǒng)的組成部分,成功應用于某任務的指揮決策、規(guī)劃分析和仿真驗證.本文提出的面向區(qū)域偵察的實拍實傳衛(wèi)星數(shù)傳規(guī)劃算法、基于網格的衛(wèi)星區(qū)域覆蓋分析算法、基于OSGEarth的目標區(qū)域和偵察范圍組合可視化方法均為首次提出,對相關研究具有借鑒和參考意義.