一種掩星探測星上反演算法

夏運兵，龔文斌，朱淑珍，劉 潔

（1.中國科學院上海微系統與信息技術研究所，上海 200050；2.上海微小衛星工程中心 上海　201203）

一種掩星探測星上反演算法

夏運兵1，2，龔文斌1，2，朱淑珍1，2，劉 潔1，2

（1.中國科學院上海微系統與信息技術研究所，上海 200050；2.上海微小衛星工程中心 上海201203）

目前掩星探測反演主要在地面進行，需要星上將掩星數據存儲，過境時回傳地面。針對星上存儲空間要求較大、數傳壓力較大、反演實時性較差的問題，本文提出一種離散化的基于TEC的電離層星上反演算法，實時進行電離層反演，直接下傳反演結果，達到減少回傳數據量，降低衛星數傳壓力、提高反演實時性的目的。

掩星探測；星上反演；星座回傳；預處理

用導航衛星無線電的掩星技術去探測地球大氣是一門上世紀80年代末發展起來的新學科[1]。它是天文學、大氣學、遙感技術等各個學科領域綜合的研究成果，是近年來空間遙測技術中最熱門的研究方向之一。與其他探測大氣的方法比較，掩星探測技術具備明顯優勢，如全天探測、準實時、高精度、高垂直分辨率、觀測點覆蓋全球并且其穩定性較好、成本較低等優點，這引起了大氣學家和大氣空間遙感及大氣環境監測甚至軍事部門的極大關注。

國內外掩星探測接收機和反演算法的研究比較深入[2]，而對于接收機輸出數據中與反演相關的數據提取技術和星上反演兩方面未見公開報道，且反演主要在地面進行。以800公里低軌衛星的STK仿真為例，其一天的可通次數為兩次，每次通信時長四五百秒，時延很大；而星座的可通信次數大大增加，為一天八十九次，極大減小時延。針對其實時性低的問題，本文提出星上反演的算法，星座回傳反演數據、接收機輸出數據的在軌預處理等技術手段，達到降低衛星數傳壓力、提高反演實時性的目的。

1　掩星探測技術

1.1技術原理

用導航衛星無線電的掩星技術去探測地球大氣的基本原理是[3]：裝載在低軌衛星上的導航衛星接收機接收導航衛星發出的電磁波信號，當電磁波信號的傳播路徑穿過大氣層時，由于受到電離層和中性大氣層的影響，電磁波信號的路徑發生延遲和彎曲，形成掩星事件，如圖1所示。

圖1　掩星探測原理圖Fig.1　Occultation diagram

掩星探測技術的主要功能：對模擬數據進行處理和分析，既可以對中性大氣的參數，如折射率、溫度、氣壓、濕度進行反演，也可以對電離層進行反演，得到電離層的參數，如電子密度和折射率隨高度變化的趨勢。

輸入和輸出為：輸入模擬電磁波的觀測信號和幅度，設置反演的類型與反演的方法；輸出使用者所需要的各種產品，如上述所說的彎曲角、折射率、中性大氣參數的溫度、氣壓、濕度和電離層的電子密度等。

1.2反演算法分析比較

掩星探測反演主要分為中性大氣和電離層兩類，它們的算法分析如下：

1）中性大氣主要的反演算法有后向傳播、滑動頻譜、正則變換、全譜反演，但是這些算法處理流程較多、計算量較大、算法復雜度較高，不利于星上在軌反演；本論文著重分析電離層反演方法。

2）電離層反演主要分為幾何光學近似Abel積分反演算法和基于TEC的Abel反演算法[4]，前者對衛星的定軌精度要求較高，且需要能夠準確計算出衛星的鐘差，處理的過程較為復雜，并不利于星上反演的實現；而基于TEC的Abel反演算法，它主要在行星間的電離層掩星探測方面應用廣泛，符合論文星上處理的需求，由于其計算的復雜度不高，并且衛星的鐘差和軌道的誤差對其不構成影響，所以將基于TEC 的Able反演的方法用于星上進行反演，直接下傳反演結果，減少衛星回傳數據量，提高反演效率。

2　星上反演算法設計

2.1TEC反演算法離散化處理

如圖2所示，為電離層的掩星探測的示意圖，將L1和L2相位進行組合計算，可以得到相對總電子含量TEC。然后運用非掩星的相關輔助觀測數據來校正相對總電子含量TEC，以此來消除衛星高度以上的電子含量，得到校正以后的TEC。假設電子密度局部球對稱以及信號直線傳播，運用Abel積分變換求出電子密度的垂直廓線隨高度變化的趨勢。

圖2　掩星探測電離層示意圖Fig.2　Sketch map of Ionospheric occultation

地球的電離層一般分布在距離地球表面的100 km以上，并且其含量在300至500公里之間達到最大值。導航衛星掩星觀測信號在電離層中的彎曲角很小，基本上小于0.01o，故導航衛星信號可被近似認為是直線傳播。將導航衛星的雙頻載波相位延遲相關的觀測數據L1和L2進行計算，得出斜向總電子含量（TEC）：

其中l為GPS信號路徑，ne為電子密度（1/m3）。在此需要強調的是：電離層的雙頻反演方法需要兩個頻率信號的電磁波傳播路徑相同，故而對此進行假設，并且忽略了因為色散的影響帶來的路徑誤差。

假設電離層電子密度在地球的局部是球對稱的，低軌衛星和導航衛星之間的斜向TEC和電子密度兩者之間滿足如下關系式：

其中，r0是掩星的切點到地球質心的距離，rGPS和rLEO分別為導航衛星和低軌衛星的軌道半徑，見圖2所示。在衛星處于電離層之內并且滿足球對稱的假設下，低軌衛星的軌道高度遠低于導航衛星軌道的高度，故電離層的相關反演就存在誤差。假設低軌衛星的軌道為圓軌道并且低軌衛星的軌道面和掩星面處于同一個平面內，那么利用非掩星側的相關輔助觀測數據去校正TEC，以次抵消低軌衛星軌道高度以上的TEC0，得到最終修正的：

運用Abel積分逆變換，就可以從校正之后的電離層延遲得出電子的密度：

對積分方程進行離散化處理，進行數值積分時，在上下邊界存在反常積分的問題 （當r→rLEO時，dTEC′（r）/dr→-∞；當r→r0時，根據

假設電子的密度在相鄰的兩個采樣點之間是線性變化的，即

那么

其中系數Ak，Bk來自于對觀測序列的擬合。將方程（6）改寫成

假設在低軌衛星軌道所處高度附近的電子密度是常數，即

對低軌衛星軌道所處高度附近的校正后的TEC′可以近似表示為

利用最小二乘線的性擬合方法[5]，可以求出最高處幾公里內的電子密度。再運用方程（8）依次計算得到每個掩星切點上的電子密度。

2.2輸入輸出量分析

基于TEC的Abel反演算法改正之后的輸入量為碰撞參數、低軌衛星和導航衛星的位置信息、雙頻信號的相位信息，可以通過EGOPS軟件模擬原始數據，得到上述輸入數據量；運用3.1節所述的反演算法，可以得到電子密度廓線。

3　仿真分析

3.1EGOPS軟件模擬反演的原始數據

電離層雙頻反演方法的原始數據可以使用基于Linux系統的EGOPS[6]軟件模擬產生，得到原始數據雙頻信號的相位L1和L2，以及低軌衛星和導航衛星位置信息的離散序列。

運行EGOPS軟件，設定仿真日期為2003年3月15日，仿真時長為48個小時，地球為橢圓形，導航系統為GPS系統，反演模式為GNSS-LEO模式，反演方式為電離層反演，采樣率為10 Hz。掩星切點的高度從1公里到20公里，步長1公里，精度0.2公里；從20公里到500公里，步長20公里，精度1公里。以下仿真是基于一個掩星事件持續252.4秒，得到2 524組的數據，表1、表2為原始數據的樣例，下列坐標建立在WGS-84坐標系下。

表1　低軌衛星的數據Tab.1　Low orbit satellite data

表2　導航衛星的數據Tab.2　Navigation satellite data

3.2TEC反演算法仿真的電離層密度廓線

將上述原始數據序列對，按照3.1節離散化處理后的方法進行計算，就可以得到電離層密度和對應的高度的序列對，從而畫出電離層密度隨高度變化的廓線圖。如下圖3所示，為TEC改進算法得到的電離層廓線圖。

3.3TEC反演算法和EGOPS軟件電離層廓線圖的對比

由下圖4可得出如下結論：

1）以簡化的TEC電離層反演方法得到電離層密度隨高度變化的曲線，與EGOPS軟件仿真所得到的電離層密度廓線基本吻合，能夠用于實際應用。

2）由對比圖可知，由于星上反演方法進過了離散化處理[8]，反演曲線與地面反演相比略有差異，不夠光滑理想，但誤差很小，可以運用于星上進行反演，從而達到高實時性反演電離層密度的目的。

圖3　TEC改進算法電離層廓線圖Fig.3　TEC improved ionospheric profile graph algorithm

3）TEC的電離層反演方法運算量很小，星載計算機能夠實現。這為掩星探測星上電離層反演技術提供可行性參考，并且星上反演的高實時性需求能夠得到滿足，達到論文研究的目的。

圖4　TEC和EGOPS反演電離層密度廓線對比Fig.4　TEC and EGOPS inversion of ionospheric density profile comparison

4　結束語

本論文分析現有掩星探測與反演存在的問題，其中以低實時性問題最為突出，提出了星上實現電離層反演的算法[9]，并用EGOPS軟件仿真圖與離散化處理后的TEC星上反演算法的仿真圖對比，證明改進的TEC星上反演方法可行、可靠、高效。

星上反演和利用星座的星間鏈路回傳反演數據，可以實現高實時性反演，并且能有效降低回傳的數據量，符合預期的目標，達到了研究的目的。這對掩星探測反演技術的大規模應用來說有著重要的作用和意義。

[1]孔璐.GPS掩星探測技術[J].氣象儀器裝備，2002（2）:22-24.

[2]苗贏，孫兆偉.星載GPS測量數據預處理方法研究[J].航空學報，2010，31（3）:602-607.

[3]嚴豪健，郭鵬，張貴霞，等.上海天文臺GPS掩星技術研究現狀[J].中國科學院上海天文臺年刊，2003（24）:39-47.

[4]吳小成，胡雄，張訓械，等.電離層GPS掩星觀測改正TEC反演方法[J].地球物理學報，2006，49（2）:328-334.

[5]陳光，任志良，孫海柱.最小二乘曲線擬合及Matlab實現[J].兵工自動化，2005，24（3）:107-108.

[6]Schweitzer S，Kirchengast G，Schwaerz M，et al.Thermodynamic state retrieval from microwave occultation data and performance analysis based on end-to-end simulations[J]. Journal of Geophysical Research Atmospheres，2011，116 （D10）:2-3.

[7]王兵，殷興輝.基于GPS觀測值提取分析西南地區電離層TEC[J].電子設計工程，2013（23）:20-23.

[8]張丹.離散復雜動態網絡的一種非脆弱同步算法[J].電子科技，2014（3）:22-24，153.

[9]王兆徽，劉宇昕，宋清濤，等.星載微波輻射計的線性與非線性反演算法比較[J].航天器工程，2013（23）:130-135.

Simulation of capture algorithms of GPS signal by Matlab/Simulink

XIA Yun-bing1，2，GONG Wen-bin1，2，ZHU Shu-zhen1，2，LIU Jie1，2
（1.Shanghai Institute of Micro-system and Information Technology Chinese Academy of Science，Shanghai 200050，China；2.Shanghai Engineering Center for Micro-satellite，Shanghai 201203，China）

The occultation inversion is mainly done on the ground，and need transit star general occultation data back to the ground when crossing the back ground.For the on-board storage requirements，large data transmission pressure，poor real-time inversion problem，this paper proposes a discrete inversion algorithm based on TEC satellite ionosphere，occultation inversion in real-time，directly transfer the inversion results，and to achieve the purpose of reducing the amount of data return，low pressure，improving the inversion of satellite data transmission in real time.

radio occultation；satellite constellation inversion；constellation return；pretreatment

TN915

A

1674－6236（2016）01-0052-03

2015-05-10稿件編號：201505085

國家自然科學基金（61401278）

夏運兵（1988—），男，四川內江人，碩士研究生。研究方向：星載接收機。