北斗衛(wèi)星ECOM光壓模型參數(shù)選擇策略分析

毛 悅，宋小勇，賈小林，阮仁桂

1.地理信息工程國家重點實驗室，陜西 西安 710054；2.西安測繪研究所，陜西 西安 710054

北斗衛(wèi)星ECOM光壓模型參數(shù)選擇策略分析

毛 悅1,2，宋小勇1,2，賈小林1,2，阮仁桂1,2

1.地理信息工程國家重點實驗室，陜西 西安 710054；2.西安測繪研究所，陜西 西安 710054

以ECOM經驗光壓模型為基礎，利用北斗衛(wèi)星3年的精密星歷進行軌道擬合，獲得了ECOM光壓參數(shù)的變化規(guī)律，給出了北斗3類衛(wèi)星的ECOM光壓參數(shù)選擇策略。從北斗衛(wèi)星姿態(tài)控制模式出發(fā)，通過衛(wèi)星星體受照分析，指出在北斗衛(wèi)星地影期零偏航狀態(tài)下，由于太陽對衛(wèi)星帆板的不正照，導致與動態(tài)偏航姿態(tài)相比，光壓攝動力存在與軌道周期相關的分量，需要在ECOM 5參數(shù)的基礎上增加D向周期分量進行吸收。通過MGEX全球網數(shù)據(jù)定軌試驗，本文提出的方法可使零偏段定軌重疊段位置精度提高50%～80%。

北斗衛(wèi)星導航系統(tǒng)；太陽光壓模型；衛(wèi)星姿態(tài)控制模式；軌道確定；ECOM；地影期

隨著我國北斗區(qū)域衛(wèi)星導航系統(tǒng)的穩(wěn)定運行并逐步向全球擴展，其對導航衛(wèi)星精密定軌與軌道預報的精度需求進一步提升。而衛(wèi)星軌道確定及預報精度在很大程度上取決于衛(wèi)星動力學模型的精確程度。中高軌導航衛(wèi)星的太陽光壓攝動力是除地球引力、日月引力外最大的非保守攝動力[1]。由于光壓攝動力與衛(wèi)星姿態(tài)控制策略、星體表面材料屬性等有關，受太陽活動、衛(wèi)星姿態(tài)控制誤差以及衛(wèi)星表面材料老化等影響，使光壓攝動力成為最難以精確建模的作用力。因此，在導航衛(wèi)星精密軌道確定中，光壓建模一直是動力學定軌中一個重要研究內容。

GPS經過多年的積累，已經建立了相對較為完善的光壓模型，并不斷精化[2-8]。但GPS全部采用MEO衛(wèi)星，動態(tài)偏航姿態(tài)控制模式，而北斗衛(wèi)星導航系統(tǒng)采用GEO、IGSO、MEO 3類衛(wèi)星，并且與QZSS導航系統(tǒng)類似[9]，采用了動態(tài)偏航與零偏航兩種姿態(tài)控制模式。由此出現(xiàn)了北斗衛(wèi)星在動/零偏切換期間及零偏弧段定軌精度下降的問題[10,12,21]。目前相關文獻的處理策略是采用增加經驗力參數(shù)，或在姿態(tài)控制模式切換點重新估計光壓力[10-13]。兩種方法均可提高衛(wèi)星軌道精度，但仍然達不到動偏的水平，且不利于軌道的預報。鑒于經驗型光壓模型在應用上更具有通用性，本文以應用較為廣泛的ECOM光壓模型為基礎，在兩種姿控模式下，利用實測數(shù)據(jù)分析獲得北斗衛(wèi)星光壓參數(shù)變化規(guī)律，對ECOM光壓模型參數(shù)選擇策略進行調整，為進一步提高兩種姿控模式下的定軌及預報精度探索可能途徑，對北斗衛(wèi)星光壓模型構建有重要參考價值。

1 ECOM光壓模型參數(shù)

BERNESE ECOM光壓模型采用常數(shù)分量加周期分量的形式，在與衛(wèi)星太陽帆板指向相關的3個相互正交方向上分別使用3組參數(shù)來吸收光壓攝動力影響，使定軌精度達到厘米級[6]。其模型表達式如下

αRPR=α0+D(u)·eD+Y(u)·eY+X(u)·eX

(1)

式中，α0為先驗光壓攝動力，通常采用分析光壓模型；eD為衛(wèi)星-太陽方向的單位向量，正向指向太陽；eY為飛行器太陽能帆板軸方向的單位向量；eX由eX=eD×eY確定。u為衛(wèi)星的升交角距。D(u)、Y(u)、X(u)分別為模型三分量系數(shù)，依照對導航衛(wèi)星軌道運動特性分析，上述系數(shù)變化特性可用多項式結合周期函數(shù)模型描述

(2)

式中，D0、DC、DS、Y0、YC、YS、X0、XC和XS為ECOM的9個待估參數(shù)。在實際定軌解算中，由于參數(shù)間的相關性，并非需要解算全部9參數(shù)[3,14]，并且過度參數(shù)化，將影響軌道的預報精度。文獻[7]中通過詳細比對分析，指出解算ECOM模型中的D0、Y0、X0、XC和XS共5參數(shù)，可獲得較高軌道精度。目前CODE等多個IGS分析中心均采用ECOM 5參數(shù)進行GPS衛(wèi)星軌道確定[3]。在北斗衛(wèi)星精密定軌中，ECOM模型也是應用較為廣泛的光壓模型。文獻[15—17]采用ECOM 5參數(shù)模型，文獻[18—19]采用ECOM 9參數(shù)。目前尚沒有文獻針對北斗3類衛(wèi)星制定ECOM模型的參數(shù)選擇策略。ECOM 5參數(shù)模型是針對GPS衛(wèi)星分析獲得的，其對北斗3類衛(wèi)星的適用性需要進一步驗證。因此如何合理選擇模型參數(shù)是本文重點討論的問題。由于北斗衛(wèi)星目前沒有可靠的分析型先驗模型，因此，對北斗光壓參數(shù)的分析采用無先驗模型方式。根據(jù)文獻[20]中的研究成果，在沒有先驗模型的前提下，采用ECOM模型可取得同等的軌道確定精度，同時還可避免先驗光壓對ECOM光壓參數(shù)解算數(shù)值的影響，達到準確反映ECOM光壓參數(shù)變化規(guī)律的作用。

2 北斗衛(wèi)星ECOM參數(shù)選擇策略分析

本文以3天弧段的北斗衛(wèi)星精密星歷作為偽觀測量，利用BERNESE ECOM模型進行動力學軌道擬合，分析北斗3類衛(wèi)星在選用不同ECOM光壓參數(shù)組合情況下的變化特征及對軌道擬合精度的影響，確定北斗3類衛(wèi)星的ECOM參數(shù)選擇策略。根據(jù)文獻[3]，采用精密星歷作為偽觀測量進行軌道擬合，是開展光壓模型構建分析的通用做法。在采用L波段觀測數(shù)據(jù)定軌時，受測量誤差及其他模型化誤差的影響需要同時解算大氣、經驗力等參數(shù)，可以采用適當增加解算參數(shù)個數(shù)的方式獲得高精度的軌道。而進行光壓模型分析時，以高精度事后軌道作為偽觀測量，可將注意力主要集中在光壓參數(shù)的主項上，盡量減少其他解算參數(shù)對光壓模型參數(shù)自身的影響，通過軌道擬合分析光壓攝動模型的特征，計算效率也更高。

文中采用的數(shù)據(jù)為通過IGS網站下載的MGEX工作組武漢大學分析中心計算的北斗衛(wèi)星精密星歷，采樣間隔為15 min。時間跨度為2013-01-01—2016-03-12，共計3年3個月。根據(jù)各文獻報道，武漢大學后處理精密星歷采用3個步驟處理完成：①進行GPS/GLONASS聯(lián)合定軌；②利用定軌結果進行GPS靜態(tài)PPP解算；③進行北斗與伽利略衛(wèi)星的聯(lián)合定軌。定軌解算時將前一步PPP解算的測站坐標、鐘差以及兩小時一組的天頂對流層延遲作為已知量參與解算。光壓模型采用CODE 5參數(shù)模型，同時解算延跡方向的常數(shù)經驗加速度，并加入1.0×10-10m/s2的先驗約束[21]。軌道重疊段精度水平為：3類衛(wèi)星法向精度優(yōu)于15 cm，徑向優(yōu)于8 cm，IGSO、MEO衛(wèi)星切向精度在10～15 cm之間，GEO衛(wèi)星切向存在1～4 m量級的偏差[15-16,18,22]。GEO衛(wèi)星激光檢核視向誤差為54.5 cm，IGSO、MEO為10 cm[22]。北斗衛(wèi)星采用零偏航姿態(tài)控制模式期間的重疊軌道精度與采用動態(tài)偏航期間基本相當，但是切向和徑向精度仍舊較動偏值稍大[10]。

在衛(wèi)星軌道擬合中采用的動力學模型包括，JGM3地球引力場模型，考慮日月引力，日月位置采用DE405；考慮太陽光壓力，相對論和固體潮、海潮引力位影響。地球自轉參數(shù)采用IERS B公報值。軌道擬合解算參數(shù)包括每顆衛(wèi)星6個軌道初始狀態(tài)參數(shù)以及光壓參數(shù)。

本文首先對比計算了選用ECOM 5參數(shù)(D0、Y0、X0、XC、XS)以及ECOM 9參數(shù)的計算結果，發(fā)現(xiàn)北斗3類衛(wèi)星在選用ECOM 5參數(shù)時均存在以Y0參數(shù)為代表的每半年一次的光壓參數(shù)變化異常，在時間上剛好與GEO衛(wèi)星的地影周期以及IGSO、MEO衛(wèi)星的零偏周期相同[10-13]。具體而言，對GEO衛(wèi)星由于其特殊的軌道特性，各衛(wèi)星的地影期十分接近(春分、秋分前后)[13]，由此在圖1中，各GEO衛(wèi)星的光壓參數(shù)異常期也基本相同。對IGSO、MEO衛(wèi)星，由于其軌位各不相同，各自的零偏時段也不一致，但均涵蓋了各自的光壓參數(shù)異常期。3類衛(wèi)星光壓參數(shù)的變化異常以GEO衛(wèi)星表現(xiàn)最為明顯。此現(xiàn)象也與參考文獻中的報道的定軌精度下降情況相一致。文獻[23]指出GEO衛(wèi)星在春秋分前后，即GEO衛(wèi)星的地影期，呈現(xiàn)出定軌精度下降的現(xiàn)象，UERE由正常段的1.5 m下降至7 m[23]。IGSO和MEO衛(wèi)星在姿態(tài)轉換期間和零偏期間軌道精度也大幅度降低[10,12]。分析其產生原因，北斗IGSO、MEO衛(wèi)星在太陽矢量與軌道面夾角(β角)小于一定的控制角度時，存在動態(tài)偏航與零偏航的姿態(tài)控制策略切換(與衛(wèi)星地影期接近)，而GEO衛(wèi)星則全弧段采用零偏航控制模式[10]。零偏航姿態(tài)控制模式可以避免衛(wèi)星正午及午夜機動，但不能確保太陽帆板對太陽的正照。因此這種定軌及軌道擬合精度下降應該與地影期零偏狀態(tài)下，衛(wèi)星受照情況變化所造成的光壓攝動力變化相關。而采用ECOM 9參數(shù)則可有效改善這一現(xiàn)象。北斗衛(wèi)星采用ECOM 5參數(shù)、9參數(shù)擬合結果的常數(shù)參數(shù)(D0、Y0、X0)時間序列如圖1和圖2所示。圖中豎線表示劃分年。由于IGSO、MEO衛(wèi)星的零偏弧段與地影期十分接近，因此下文在對3類衛(wèi)星綜合描述時，通稱為地影期，不再單獨強調零偏。

圖1 北斗3類衛(wèi)星ECOM 5參數(shù)常數(shù)項時間序列圖Fig.1 Time sequence diagram of the constant items using the ECOM 5 parameters model for the three types of BeiDou satellites

圖3為采用ECOM 9參數(shù)進行軌道擬合獲得的ECOM周期參數(shù)(XC、XS、DC、DS、YC、YS)時間序列圖。由圖中可得，ECOM 9參數(shù)模型擬合結果中，北斗3類衛(wèi)星的YC、YS參數(shù)較為恒定，未出現(xiàn)周期性，且解算結果接近0，推斷YC、YS參數(shù)對軌道擬合結果貢獻不大。由此本文重點分析采用ECOM 7參數(shù)(D0、Y0、X0、XC、XS、DC、DS)的軌道擬合結果。由圖4中結果，ECOM 7參數(shù)仍然可以較好地解決ECOM 5參數(shù)中地影期軌道擬合異常的問題。

圖2 北斗3類衛(wèi)星ECOM 9參數(shù)常數(shù)項時間序列圖Fig.2 Time sequence diagram of the constant items using the ECOM 9 parameters model for the three types of BeiDou satellites

圖3 北斗衛(wèi)星ECOM 9參數(shù)周期參數(shù)時間序列圖Fig.3 Time sequence diagram of the period items using the ECOM 9 parameters model for the BeiDou satellite

圖5為2013年至2016年所有非地影、地影期3天弧段軌道擬合數(shù)據(jù)點擬合殘差的綜合統(tǒng)計結果。圖中“5p”“7p”“9p”分別表示ECOM 5參數(shù)、7參數(shù)、9參數(shù)(下文同)。

在非地影弧段內GEO衛(wèi)星的擬合殘差明顯高于IGSO、MEO，這是由于GEO衛(wèi)星相對特殊的空間位置，導致地面站對GEO衛(wèi)星的觀測幾何結構不佳造成的。GEO衛(wèi)星采用ECOM 7參數(shù)的三維提升量最大，這與GEO全弧段采用零偏姿態(tài)有關，但其提升量也僅在厘米量級。IGSO、MEO的提升量均在毫米量級。建議采用ECOM 5參數(shù)模型。而地影段，與ECOM 5參數(shù)相比，3類衛(wèi)星采用ECOM 7參數(shù)進行軌道擬合的殘差提升量均可達到分米量級。ECOM 9參數(shù)擬合殘差的提升量較ECOM 7參數(shù)已較小，在厘米量級，因此建議北斗3類衛(wèi)星在地影期采用ECOM 7參數(shù)進行軌道確定。

圖4 北斗衛(wèi)星ECOM 7參數(shù)時間序列圖Fig.4 Time sequence diagram of the ECOM 7 parameters for the BeiDou satellites

圖5 軌道三維位置擬合殘差統(tǒng)計(RMS)Fig.5 Fitting residual statistics of satellite position(RMS)

為驗證前文提出的北斗3類衛(wèi)星ECOM光壓參數(shù)選取策略的正確性，本文依據(jù)衛(wèi)星在不同情況下的姿態(tài)控制策略，以衛(wèi)星受照情況為研究對象，梳理光壓參數(shù)選取策略的理論機制。

相對于衛(wèi)星本體，太陽帆板的面積最大，反射系數(shù)最低。衛(wèi)星姿態(tài)控制模式變化對太陽帆板受照的影響，是太陽光壓攝動力變化的重要影響因素。在動偏期間，衛(wèi)星通過太陽敏感器感知太陽位置，自動計算偏航角大小，使太陽帆板始終垂直于衛(wèi)星-太陽(D)方向。而零偏期間，衛(wèi)星不再跟蹤太陽位置，始終控制偏航角為零。此時，如圖6所示，帆板法向與D方向夾角隨時間變化。該入射角隨太陽與軌道面夾角β變化而變化，并非嚴格的線性趨勢，不同衛(wèi)星入射角的變化速率不同，但總體趨勢一致，對3類衛(wèi)星也沒有本質區(qū)別。

圖6 零偏帆板法向與D向夾角Fig.6 The angle between the D direction and the normal direction of the sailboard

對衛(wèi)星帆板而言，在動偏期間，衛(wèi)星太陽帆板正照，大部分光線被太陽能板吸收轉化為能量。其所產生的光壓攝動力也完全在衛(wèi)星-太陽連線方向，即ECOM模型的D方向，不存在與軌道周期相關的力。因此在動偏期間僅解算ECOM模型的5參數(shù)即可。在零偏期間，帆板太陽光線的入射角度隨衛(wèi)星位置而變化，對帆板照射所產生的光壓力即存在了與軌道周期相關的周期分量。該周期變化同時體現(xiàn)在衛(wèi)星-太陽連線方向(D方向)和太陽帆板方向(Y方向)。由于YC和XS，YS和XC存在強相關性，相關度達到0.995和0.951[24]。因此在零偏期間僅在ECOM 5參數(shù)的基礎上，增加D向周期項即可。驗證了建議北斗3類衛(wèi)星在地影期采用ECOM 7參數(shù)進行軌道確定的正確性。

3 參數(shù)選擇策略試驗驗證

為進一步驗證北斗衛(wèi)星光壓參數(shù)選擇策略在軌道確定中的應用效果，本文采用了2016年MGEX網站上提供的全球61個北斗雙頻L波段地面監(jiān)測站數(shù)據(jù)進行了單北斗系統(tǒng)的精密定軌試驗。對北斗GEO衛(wèi)星(C01-C05)，采用了2006-03-01—2016-04-08共39天的春分前后地影段數(shù)據(jù)。此外還對2016-05-19—2016-06-29共42天的非地影段數(shù)據(jù)進行了對比計算分析?？紤]到2016年6月9日至20日C11、C12(MEO)衛(wèi)星處于零偏弧段，2016年6月19日至6月29日C08(IGSO)衛(wèi)星處于零偏弧段，同樣為了對比分析，本文將C11、C12、C08衛(wèi)星的計算弧段統(tǒng)一擴大到2016-05-19—2016-06-29。

定軌解算時采用SPODS軟件[25]，利用3 d弧段5 min采樣雙頻無電離層組合觀測量?？紤]衛(wèi)星天線相位中心修正、衛(wèi)星天線相位纏繞(Wind_up)、相對論效應等。對流層采用Saastamoinen天頂延遲改正及GMF映射函數(shù)。采用EGM2008 12階地球引力場，IERS2003固體潮，三體引力考慮日、月及大行星，采用JPL DE405行星星歷。光壓模型對比采用無先驗ECOM 5參數(shù)、7參數(shù)模型。解算參數(shù)包括6個衛(wèi)星初始軌道信息，12/d對流層天頂延遲，1/d對流層水平梯度，以及5～7個光壓模型參數(shù)，固定模糊度。

3類衛(wèi)星3天重疊1天的定軌重疊段時間序列及ECOM 5參數(shù)、7參數(shù)對比統(tǒng)計結果如圖7、圖8所示。圖中R、T、N分別代表徑向、切向、法向3方向。兩條豎線之間為衛(wèi)星的地影(GEO)或零偏(IGSO、MEO)弧段。由圖中結果可以明顯看出在非地影段，ECOM 5參數(shù)、7參數(shù)定軌結果精度相當，但在地影或零偏弧段內，ECOM 7參數(shù)表現(xiàn)明顯優(yōu)于5參數(shù)。

圖8中的百分比數(shù)字為相比ECOM 5參數(shù)，采用ECOM 7參數(shù)定軌重疊段精度的提升百分比。由圖中結果統(tǒng)計可得：對地影或零偏弧段而言，采用ECOM 7參數(shù)可使定軌重疊段位置精度提高50%～80%，且重點表現(xiàn)在T、N方向，提升量最高可達86.9%。對于正常弧段，采用ECOM 7參數(shù)雖然可以實現(xiàn)位置精度的13%～23%的提升，但隨之也造成了R方向38%～49%的精度下降，這應該與過度參數(shù)化有關。由此再次驗證了地影段采用ECOM 7參數(shù)、正常段采用ECOM 5參數(shù)的結論。

定軌模型參數(shù)選擇的優(yōu)劣也同樣表現(xiàn)在于預報精度上。本文中定軌結果預報1天與后處理精密軌道相比的誤差時間序列及RMS統(tǒng)計結果見圖9和圖10。由圖中結果，可以得出與定軌段相同的結論。

前文采用了重疊弧段的方式來評價定軌精度的優(yōu)劣，為了更加全面準確的評價定軌結果，本文采用了同期的激光數(shù)據(jù)，對C01、C11、C08衛(wèi)星進行精度評估。評估結果時間序列及統(tǒng)計信息見圖11、表1。對C01衛(wèi)星，由于激光數(shù)據(jù)的缺失，未能完成對比最明顯的地影中心弧段的精度評估，由此本文將ECOM 5參、7參激光檢核結果(5p、7p)與同期ECOM 5參、7參定軌R方向重疊段結果(5p-R、7p-R)繪制在一張圖上，以顯示激光數(shù)據(jù)的對應時段。由于在有激光數(shù)據(jù)的弧段內5參、7參定軌差異還不顯著，因此表1中C01衛(wèi)星統(tǒng)計結果差異也不十分顯著，但也可達到7 cm。對C11、C08衛(wèi)星均可表現(xiàn)出零偏弧段ECOM 7參數(shù)激光檢核精度的提升，最大差異為10 cm。但對于C08衛(wèi)星的正常段， ECOM 7參數(shù)的評估結果精度是下降的，這與前文重疊弧段分析結果一致。

圖7 北斗3類衛(wèi)星實測定軌重疊段時間序列Fig.7 Time sequence diagram of the BeiDou overlap orbit precision calculated using observed data

圖8 北斗衛(wèi)星實測定軌重疊段RMS統(tǒng)計結果Fig.8 RMS of the BeiDou satellite overlap orbit precision calculated using observed data

另外本文的激光檢核結果也與前文獻[22]中武大精密星歷檢核結果精度一致。進一步說明采用ECOM 7參數(shù)進行地影零偏弧段的定軌可以達到與武漢大學綜合定軌精度相近的水平，避免了ECOM 5參數(shù)精度下降的情況。另外本文的驗證分析均是基于全球網數(shù)據(jù)進行計算的?？紤]到區(qū)域網監(jiān)測站對北斗衛(wèi)星的監(jiān)測弧段短、觀測幾何相對不足，進而定軌精度較低，這將會淹沒光壓參數(shù)對定軌精度的影響，不利于驗證分析。但其在理論上具有一致性，因此本文不再單獨針對區(qū)域網進行驗證分析。

圖9 北斗3類衛(wèi)星預報1天軌道誤差時間序列Fig.9 Time sequence of the predicted orbit precision

圖10 北斗衛(wèi)星實測定軌預報段RMS統(tǒng)計結果Fig.10 RMS of the predicted orbit precision for BeiDou satellites

圖11 北斗3類衛(wèi)星定軌激光檢核結果時間序列Fig.11 SLR residuals for three types of BeiDou satellites

表1 北斗3類衛(wèi)星定軌激光檢核RMS統(tǒng)計結果Tab.1 RMS statistics of SLR evaluated results m

4 結 論

本文以分析BERNESE ECOM經驗型光壓模型對北斗衛(wèi)星的適用性為目標，通過以精密星歷為基礎的軌道擬合，開展了北斗3類衛(wèi)星光壓參數(shù)變化特性分析，制定了北斗衛(wèi)星光壓參數(shù)選擇策略；進行了北斗衛(wèi)星動偏、零偏姿態(tài)下衛(wèi)星星體受照情況分析，給出了光壓參數(shù)選取策略的理論依據(jù)。利用北斗MGEX全球網數(shù)據(jù)進行了定軌試驗，驗證了該策略的正確性。主要結論如下：

(1) 北斗GEO衛(wèi)星地影期，IGSO/MEO衛(wèi)星零偏弧段建議采用ECOM 7參數(shù)進行定軌解算。其他弧段建議采用ECOM 5參數(shù)。

(2) 本文提出的光壓參數(shù)選擇策略的產生機制為，在地影期零偏姿態(tài)下，太陽對衛(wèi)星帆板不完全正照，存在與衛(wèi)星軌道相關的周期項，且主要反映在D、Y方向上?？紤]參數(shù)間的相關性，應當在5參數(shù)基礎上引入D向周期參數(shù)進行吸收。

(3) 通過長達81 d的定軌試驗得出：在地影期零偏姿態(tài)下，與ECOM 5參數(shù)相比，采用ECOM 7參數(shù)可有效提高定軌位置精度，提升量達到50%～80%，且重點表現(xiàn)在T、N方向。驗證了ECOM光壓參數(shù)選擇策略的有效性。

[1] KAPLAN E D，HEGARTY C J.Understanding GPS: Principles and Applications[M].2nd ed.Norwood，MA: Artech House Inc，2005.

[2] FLIGEEL H F，GALLINI T E，SWIFT E R.Global Positioning System Radiation Force Model for Geodetic Applications[J].Journal of Geophysical Research，1992，97(B1): 559-568.

[3] SPRINGER T A.Modeling and Validating Orbits and Clocks Using the Global Positioning System[D].Berne，Switzerland: Astronomical Institute，University of Berne，1999.

[4] COLOMBO O L.The Dynamics of Global Positioning System Orbits and the Determination of Precise Ephemerides[J].Journal of Geophysical Research，1989，94(B7): 9167-9182.

[5] BAR-SEVER Y E.New and Improved Solar Radiation Models for GPS Satellites Based on Flight Data[R].Pasadena，CA: Jet Propulsion Laboratory，California Institute of Technology，1997.

[6] BEUTLER G，BROCKMANN E，GURTNER W，et al.Extended Orbit Modeling Techniques at the CODE Processing Center of the International GPS Service for Geodynamics (IGS): Theory and Initial Results[J].Manuscripta Geodaetica，1994，19(6): 367-386.

[7] SPRINGER T A，BEUTLER G，ROTHACHER M.Improving the Orbit Estimates of GPS Satellites[J].Journal of Geodesy，1999，73(3): 147-157.

[8] RODRIGUEZ-SOLANO C J，HUGENTOBLER U，STEIGENBERGER P.Adjustable Box-wing Model for Solar Radiation Pressure Impacting GPS Satellites[J].Advances in Space Research，2012，49(7): 1113-1128.

[9] 周佩元，杜蘭，方善傳，等.姿態(tài)模式切換期間QZSS衛(wèi)星軌道及其鐘差產品特性分析[J].測繪學報，2016，45(3): 274-281.DOI: 10.11947/j.AGCS.2016.20150386.

ZHOU Peiyuan，DU Lan，F(xiàn)ANG Shanchuan，et al.Analysis of Characteristics of QZSS Satellite Orbit and Clock Products during Yaw Attitude Model Switching[J].Acta Geodaetica et Cartographica Sinica，2016，45(3): 274-281.DOI: 10.11947/j.AGCS.2016.20150386.

[10] 郭靖.姿態(tài)、光壓和函數(shù)模型對導航衛(wèi)星精密定軌影響的研究[D].武漢: 武漢大學，2014.

GUO Jing.The Impacts of Attitude，Solar Radiation and Function Model on Precise Orbit Determination for GNSS Satellites[D].Wuhan: Wuhan University，2014.

[11] 毛悅，宋小勇，王維，等.北斗IGSO/MEO衛(wèi)星姿態(tài)控制及光壓差異分析[J].測繪科學，2015，40(8): 129-134.

MAO Yue，SONG Xiaoyong，WANG Wei，et al.BeiDou IGSO and MEO Navigation Satellites’ Yaw-Steering and Orbit-normal Attitude Control Modes and Solar Radiation Pressure Difference Analysis[J].Science of Surveying and Mapping，2015，40(8): 129-134.

[12] 毛悅，宋小勇，王維，等.IGSO姿態(tài)控制模式切換期間定軌策略研究[J].武漢大學學報(信息科學版)，2014，39(11): 1352-1356.

MAO Yue，SONG Xiaoyong，WANG Wei，et al.IGSO Satellite Orbit Determining Strategy Analysis with the Yaw-Steering and Orbit-normal Attitude Control Mode Switching[J].Geomatics and Information Science of Wuhan University，2014，39(11): 1352-1356.

[13] 毛悅，宋小勇，賈小林，等.北斗導航衛(wèi)星地影狀態(tài)分析[J].測繪學報，2014，43(4): 353-359.DOI: 10.13485/j.cnki.11-2089.2014.0053.

MAO Yue，SONG Xiaoyong，JIA Xiaolin.Earth Eclipse Status Analysis of BeiDou Navigation Satellites[J].Acta Geodaetica et Cartographica Sinica，2014，43(4): 353-359.DOI: 10.13485/j.cnki.11-2089.2014.0053.

[14] 宋小勇，毛悅，賈小林.BERNESE光壓模型參數(shù)的統(tǒng)計分析[J].測繪科學，2009，34(3): 25-27.

SONG Xiaoyong，MAO Yue，JIA Xiaolin.Analysis of the Parameter of BERNESE Radiation Force Model by Statistic Method[J].Science of Surveying and Mapping，2009，34(3): 25-27.

[15] 李敏.多模GNSS融合精密定軌理論及其應用研究[D].武漢: 武漢大學，2011.

LI Min.Research on Multi-GNSS Precise Orbit Determination Theory and Application[D].Wuhan: Wuhan University，2011.

[16] ZHAO Qile，GUO Jing，LI Min，et al.Initial Results of Precise Orbit and Clock Determination for COMPASS Navigation Satellite System[J].Journal of Geodesy，2013，87(5): 475-586.

[17] 劉偉平，郝金明，李建文，等.多GNSS融合的北斗衛(wèi)星精密定軌[J].測繪學報，2014，43(11): 1132-1138.DOI: 10.13485/j.cnki.11-2089.2014.0186.

LIU Weiping，HAO Jinming，LI Jianwen，et al.Multi-GNSS Joint Precise Orbit Determination of BeiDou Navigation Satellites System[J].Acta Geodaetica et Cartographica Sinica，2014，43(11): 1132-1138.DOI: 10.13485/j.cnki.11-2089.2014.0186.

[18] GE Maorong，ZHANG Hongping，JIA Xiaolin，et al.What is Achievable with Current COMPASS Constellation[DB/OL].[2012-11-01].http:∥www.gpsworld.com/what-is.-achievable-with-the-current-compass-constellation.

[19] 劉偉平，郝金明，魏偉，等.兩種北斗衛(wèi)星精密定軌方法的對比[J].大地測量與地球動力學，2016，36(3): 211-215.

LIU Weiping，HAO Jinming，WEI Wei，et al.Comparison of Two Methods of Precise Orbit Determination of BeiDou Satellites[J].Journal of Geodesy and Geodynamics，2016，36(3): 211-215.

[20] MONTENBRUCK O，STEIGENBERGER P，HUGEN TOBLER U.Enhanced Solar Radiation Pressure Modeling for Galileo Satellites[J].Journal of Geodesy，2015，89(2): 283-297.

[21] GUO Jing，XU Xiaolong，ZHAO Qile，et al.Precise Orbit Determination for Quad-constellation Satellites at Wuhan University: Strategy，Result Validation，and Comparison[J].Journal of Geodesy，2016，90(2): 143-159.DOI: 10.1007/s00190-015-0862-9.

[22] 耿濤，蘇醒，許小龍，等.北斗衛(wèi)星導航系統(tǒng)精密定軌和廣播星歷軌道精度分析[J].中國科技論文，2015，10(9): 1023-1026，1032.

GENG Tao，SU Xing，XU Xiaolong，et al.Research on Precise and Broadcast Ephemeris Orbit of BeiDou Navigation Satellite System[J].China Science Paper，2015，10(9): 1023-1026，1032.

[23] LI X J，ZHOU J H，HU X G，et al.Orbit Determination and Prediction for BeiDou GEO Satellites at the Time of the Spring/Autumn Equinox[J].Science China Physics，Mechanics & Astronomy，2015，58(8): 089501.DOI: 10.1007/s11433-015-5675-6.

[24] 劉永，張傳定，李珊珊，等.基于IGS精密星歷的BERNESE光壓模型參數(shù)分析[J].測繪工程，2016，25(1): 29-33，38.

LIU Yong，ZHANG Chuanding，LI Shanshan，et al.Analysis of BERN SRP Parameters Based on IGS Ephemeris[J].Engineering of Surveying and Mapping，2016，25(1): 29-33，38.

[25] 魏子卿，阮仁桂，賈小林，等.衛(wèi)星定位定軌系統(tǒng)SPODS: 理論與測試[J].測繪學報，2014，43(1): 1-4.DOI: 10.13485/j.cnki.11-2089.2014.0001.

WEI Ziqing，RUAN Rengui，JIA Xiaolin，et al.Satellite Positioning and Orbit Determination System SPODS: Theory and Test[J].Acta Geodaetica et Cartographica Sinica，2014，43(1): 1-4.DOI: 10.13485/j.cnki.11-2089.2014.0001.

Analysis about Parameters Selection Strategy of ECOM Solar Radiation Pressure Model for BeiDou Satellites

MAO Yue1,2,SONG Xiaoyong1,2,JIA Xiaolin1,2,RUAN Rengui1,2

1.State Key Laboratory of Geo-information Engineering，Xi’an 710054，China; 2.Xi’an Research Institute of Surveying and Mapping，Xi’an 710054， China

In this paper，the variations of ECOM solar radiation pressure model parameters were analyzed for BeiDou satellites through orbit fitting using 3 years precise ephemeris data.The ECOM parameter selection strategies for the three types of BeiDou satellites were confirmed.Based on the satellite attitude control modes，satellite illumination was theoretically analyzed.We pointed out that solar panels are perpendicular illuminated，and then solar radiation pressure has the periodic component associated with the orbital period when satellite is at orbit-normal attitude control mode.Periodic parameters of D direction need to be added to enhance the ECOM 5 parameters.Orbit determination tests were carried out using MGEX global network data.Tests prove that using ECOM 7 parameters the orbit overlap accuracy improvement could reach up to 50%～80%.

BeiDou navigation satellite system；solar radiation pressure model；satellite attitude control mode；orbit determination；ECOM；earth eclipse segment

The National Natural Science Foundation of China (No.41774012)

MAO Yue(1981—)，female，engineer，PhD，majors in satellite navigation researches and its applications.

毛悅，宋小勇，賈小林，等.北斗衛(wèi)星ECOM光壓模型參數(shù)選擇策略分析[J].測繪學報，2017,46(11)：1812-1821.

10.11947/j.AGCS.2017.20160485.

MAO Yue,SONG Xiaoyong,JIA Xiaolin,et al.Analysis about Parameters Selection Strategy of ECOM Solar Radiation Pressure Model for BeiDou Satellites[J].Acta Geodaetica et Cartographica Sinica,2017,46(11):1812-1821.DOI:10.11947/j.AGCS.2017.20160485.

P237

A

1001-1595(2017)11-1812-10

國家自然科學基金(41774012)

(責任編輯：宋啟凡)

2016-09-30

修回日期：2017-09-18

毛悅(1981—)，女，工程師，博士，主要從事衛(wèi)星導航及其應用研究。

E-mail：maoyue0810@163.com