衛(wèi)星激光測(cè)距數(shù)據(jù)處理規(guī)范中華人民共和國自然資源部發(fā)布 l I 2規(guī)范性引用文件 3術(shù)語和定義 14縮略語 25產(chǎn)品要求 25.1產(chǎn)品內(nèi)容 25.2產(chǎn)品精度 25.3產(chǎn)品格式 36數(shù)據(jù)處理 36.1一般原則 6.2數(shù)據(jù)準(zhǔn)備 6.3觀測(cè)模型 6.4動(dòng)力學(xué)模型 6.5解算策略 77精度評(píng)定 7.1衛(wèi)星激光測(cè)距觀測(cè)數(shù)據(jù)精度評(píng)定 87.2地球自轉(zhuǎn)參數(shù)精度評(píng)定 87.3地心運(yùn)動(dòng)參數(shù)精度評(píng)定 87.4測(cè)站坐標(biāo)及速度精度評(píng)定 8附錄A(規(guī)范性附錄)觀測(cè)模型 附錄B(規(guī)范性附錄)動(dòng)力學(xué)模型 本標(biāo)準(zhǔn)是“衛(wèi)星激光測(cè)距”系列標(biāo)準(zhǔn)的第二項(xiàng)標(biāo)準(zhǔn)，第一項(xiàng)、第三項(xiàng)標(biāo)準(zhǔn)分別為《衛(wèi)星激光測(cè)距數(shù)據(jù)獲取規(guī)范》和《衛(wèi)星激光測(cè)距數(shù)據(jù)庫建設(shè)規(guī)范》。本標(biāo)準(zhǔn)的起草規(guī)則依據(jù)GB/T1.1—2009。本標(biāo)準(zhǔn)由中華人民共和國自然資源部提出并歸口。本標(biāo)準(zhǔn)起草單位：中國測(cè)繪科學(xué)研究院、中國科學(xué)院上海天文臺(tái)、山東科技大學(xué)、中國礦業(yè)大學(xué)。本標(biāo)準(zhǔn)主要起草人：趙春梅、王小亞、郭金運(yùn)、王潛Ⅱ大地基準(zhǔn)的建立和維持，不僅為當(dāng)今經(jīng)濟(jì)和國防建設(shè)提供精確定位和定向服務(wù)，還為開展地球科學(xué)研究提供極其重要的參考基準(zhǔn)。衛(wèi)星激光測(cè)距是建立與維持大地基準(zhǔn)重要的觀測(cè)手段，衛(wèi)星激光測(cè)距數(shù)據(jù)處理規(guī)范關(guān)系到成果的質(zhì)量和應(yīng)用。為了指導(dǎo)衛(wèi)星激光測(cè)距數(shù)據(jù)處理，以及生產(chǎn)建立和維持地球參考框架所需的地球自轉(zhuǎn)參數(shù)、地心運(yùn)動(dòng)、測(cè)站坐標(biāo)及速度等大地基準(zhǔn)產(chǎn)品，制定本標(biāo)準(zhǔn)。1衛(wèi)星激光測(cè)距數(shù)據(jù)處理規(guī)范本標(biāo)準(zhǔn)規(guī)定了衛(wèi)星激光測(cè)距數(shù)據(jù)應(yīng)用中的產(chǎn)品要求、數(shù)據(jù)處理和精度評(píng)定，確立了相關(guān)大地基準(zhǔn)產(chǎn)本標(biāo)準(zhǔn)適用于衛(wèi)星激光測(cè)距數(shù)據(jù)處理及相關(guān)大地基準(zhǔn)產(chǎn)品的生產(chǎn)及應(yīng)用。2規(guī)范性引用文件下列文件對(duì)于本文件的應(yīng)用是必不可少的。凡是注日期的引用文件，僅所注日期的版本適用于本文件。凡是不注日期的引用文件，其最新版本(包括所有的修改單)適用于本文件。GB/T17159—2009大地測(cè)量術(shù)語IERS規(guī)范(2010)規(guī)定國際地球自轉(zhuǎn)和參考系統(tǒng)服務(wù)的規(guī)范(InternationalEarthRotationand3術(shù)語和定義下列術(shù)語和定義適用于本文件。利用激光測(cè)距儀在地面上跟蹤觀測(cè)裝有激光反射棱鏡的衛(wèi)星，以測(cè)定測(cè)站到衛(wèi)星的距離的測(cè)量技術(shù)和方法。[GB/T17159—2009,定義6.82]站心坐標(biāo)系topocentricco以測(cè)站為原點(diǎn)的坐標(biāo)系。[GB/T17159—2009,定義3.66]國際地球參考框架internationalterrestrialreferenceframe;ITRF國際地球參考系統(tǒng)的實(shí)現(xiàn)。由國際地球自轉(zhuǎn)和參考系統(tǒng)服務(wù)根據(jù)空間大地測(cè)量技術(shù)，包括甚長(zhǎng)基線干涉測(cè)量、衛(wèi)星激光測(cè)距、多里斯系統(tǒng)、全球?qū)Ш叫l(wèi)星系統(tǒng)等，所確定的相標(biāo)及速度集。衛(wèi)星激光測(cè)距參考框架satellitelaserreferenceframe;SLRF根據(jù)衛(wèi)星激光測(cè)距技術(shù)確定的相對(duì)于特定參考?xì)v元的衛(wèi)星激光測(cè)距測(cè)站坐標(biāo)及速度集，是地球參考系的一種實(shí)現(xiàn)。[GB/T17159—2009,定義3.72]2地球質(zhì)心相對(duì)地球參考框架協(xié)議原點(diǎn)的位移。地球自轉(zhuǎn)參數(shù)earthrotationparame表示地球自轉(zhuǎn)的速率、自轉(zhuǎn)軸方向及其變化的參數(shù)。地球定向參數(shù)earthorientationparameters;EOP表示地球參考系統(tǒng)相對(duì)于天球參考系統(tǒng)的定向參數(shù)。4縮略語下列縮略語適用于本文件。CSR美國空間研究中心(CenterforSpaceResearch)CSV逗號(hào)分隔值文件格式(commaseparatedvalues)ILRS國際激光測(cè)距服務(wù)(InternationalLaserRangingService)IAG國際大地測(cè)量學(xué)會(huì)(InternationalAssociationofGeodesy)LOD日長(zhǎng)(lengthofday)SINEX獨(dú)立解交換格式[solution(software/technique)independentexchangeformat]5產(chǎn)品要求衛(wèi)星激光測(cè)距大地基準(zhǔn)產(chǎn)品包括地球自轉(zhuǎn)參數(shù)、測(cè)站坐標(biāo)及速度、地心運(yùn)動(dòng)參數(shù)。地球自轉(zhuǎn)參數(shù)包括極移、日長(zhǎng)參數(shù)及其精度，以及對(duì)應(yīng)的觀測(cè)歷元。測(cè)站坐標(biāo)及速度包括測(cè)站編號(hào)、名稱、在特定參考框架下對(duì)應(yīng)觀測(cè)歷元的站坐標(biāo)和速度及其精度。地心運(yùn)動(dòng)參數(shù)包括地球瞬時(shí)質(zhì)量中心相對(duì)地球參考框架協(xié)議原點(diǎn)的位置及其精度，以及對(duì)應(yīng)的觀測(cè)歷元。地球自轉(zhuǎn)參數(shù)精度要求見表1。實(shí)際應(yīng)用中，可以表1至表3中精度的2倍為標(biāo)準(zhǔn)。3表1地球自轉(zhuǎn)參數(shù)精度要求日長(zhǎng)對(duì)于不少于5年連續(xù)觀測(cè)的衛(wèi)星激光測(cè)距測(cè)站，其測(cè)站坐標(biāo)及速度精度要求見表2。表2測(cè)站坐標(biāo)及速度精度要求地心運(yùn)動(dòng)參數(shù)精度要求見表3。表3地心運(yùn)動(dòng)參數(shù)精度要求地心運(yùn)動(dòng)5.3產(chǎn)品格式衛(wèi)星激光測(cè)距大地基準(zhǔn)產(chǎn)品應(yīng)采用CSV格式或SINEX格式。產(chǎn)品格式及內(nèi)容要求見表4。表4產(chǎn)品格式及內(nèi)容要求內(nèi)容要求SINEX格式解算模塊：觀測(cè)歷元、統(tǒng)計(jì)信息、參數(shù)先驗(yàn)值、參數(shù)估計(jì)結(jié)果、方差6數(shù)據(jù)處理6.1一般原則數(shù)據(jù)處理應(yīng)遵循以下原則：4a)采用地球動(dòng)力學(xué)衛(wèi)星觀測(cè)數(shù)據(jù)計(jì)d)數(shù)據(jù)處理模型遵循IERS規(guī)范(2010)。p1=p?+(△pm+△PRF+△PRL+△PAc+ (1)pp△pm——潮汐引起的測(cè)站位置變化而導(dǎo)致的測(cè)距誤差，是系統(tǒng)觀測(cè)誤差之一，誤差改正計(jì)算見測(cè)站坐標(biāo)應(yīng)采用IERS提供的國際地球參考框架或衛(wèi)星激光測(cè)距參考框架最新序列，如目前的5測(cè)站位置變化。其改正模型應(yīng)符合最新IERS規(guī)范，如目前的IERS規(guī)范(202)極潮改正。極潮是由地球自轉(zhuǎn)不均勻測(cè)站位置變化。其改正模型應(yīng)符合最新IERS規(guī)范，如目前的IERS規(guī)范(2010),具體模c)相對(duì)論效應(yīng)改正。廣義相對(duì)論效應(yīng)使激光在地球引力場(chǎng)中傳播時(shí)發(fā)生引力時(shí)延。對(duì)于力學(xué)衛(wèi)星，激光測(cè)距引力時(shí)延可達(dá)厘米量級(jí)，數(shù)據(jù)處理中應(yīng)考慮其影響，具體模型見公式(A.9)。d)衛(wèi)星質(zhì)心改正。由于激光反射器中心相對(duì)于衛(wèi)星質(zhì)心存在偏差，故須對(duì)激光測(cè)距進(jìn)行衛(wèi)星質(zhì)心改正。衛(wèi)星質(zhì)心改正模型取決于衛(wèi)星的外形、激光測(cè)站系統(tǒng)運(yùn)行方式，該項(xiàng)改正隨衛(wèi)星而異，常用的幾個(gè)地球動(dòng)力學(xué)衛(wèi)星質(zhì)心改正模型見e)系統(tǒng)偏差改正。衛(wèi)星激光測(cè)距系統(tǒng)偏差包括系f)測(cè)站偏心改正。衛(wèi)星激光測(cè)距儀的測(cè)量中心相對(duì)于測(cè)站標(biāo)定坐標(biāo)之間的偏差就是測(cè)站偏心改f=fm+fm+fws+fm+fgn+fsr+fer+foc+frh (2)6太陽直射到衛(wèi)星表面產(chǎn)生太陽輻射壓攝動(dòng)。對(duì)形狀復(fù)雜的衛(wèi)星，可將其分為若干個(gè)平面分別進(jìn)行計(jì)算，然后矢量求和得到所受到的太陽輻射壓攝動(dòng)加速度，具體模型見公式(B.2)。對(duì)型見公式(B.3)。大氣阻力攝動(dòng)是大氣對(duì)衛(wèi)星運(yùn)動(dòng)產(chǎn)生的阻力而引起的衛(wèi)星攝動(dòng)，其模型見公式(B.4)。對(duì)于有太陽帆板的衛(wèi)星，需要同時(shí)考慮太陽帆板產(chǎn)生的大氣阻力，其模型見公式(B.5)。地球形狀攝動(dòng)是由地球形狀不規(guī)則和質(zhì)量分布不均造成的附加攝動(dòng)。在地固坐標(biāo)系中，該攝動(dòng)對(duì)應(yīng)的引力位計(jì)算模型見公式(B.6)。固體潮使地球引力位發(fā)生變化，從而對(duì)衛(wèi)星軌道產(chǎn)生固體潮攝動(dòng)。固體潮模型可以瓦爾(Wahr)固體潮汐模型為基礎(chǔ)，考慮液核的動(dòng)力學(xué)影響，對(duì)不同的分潮波(主要是全日波)采用不同的勒夫數(shù)。勒夫數(shù)應(yīng)隨分潮波頻率的不同而不同，其模型應(yīng)符合最新IERS規(guī)范，如目前的IERS規(guī)范(2010)。具體模型見公式(B.7)、公式(B.8)、公式(B.9)。海潮使地球引力位發(fā)生變化，從而對(duì)衛(wèi)星軌道產(chǎn)生海潮攝動(dòng)。應(yīng)采用海潮模型，通過修正地球引力位系數(shù)計(jì)算海潮攝動(dòng)，其模型應(yīng)符合最新IERS規(guī)范，如目前的IERS規(guī)范(2010)。具體模型見公式大氣潮使地球引力位發(fā)生變化，從而對(duì)衛(wèi)星軌道產(chǎn)生大氣潮攝動(dòng)。應(yīng)采用大氣潮模型，通過修正地球引力位系數(shù)計(jì)算大氣潮攝動(dòng)，其模型應(yīng)符合最新IERS規(guī)范，如目前的IERS規(guī)范(2010)。具體模型見公式(B.11)、公式(B.12)。地球自轉(zhuǎn)產(chǎn)生的離心力引起地球體積和密度分布的變化，進(jìn)而使衛(wèi)星軌道攝動(dòng)產(chǎn)生一個(gè)附加攝動(dòng)。地球輻射壓攝動(dòng)包括地球反照輻射壓攝動(dòng)和紅外輻射壓攝動(dòng)，具體模型見公式(B.14)、公式月球J2項(xiàng)扁率攝動(dòng)是將月球視作扁球體時(shí)月球?qū)πl(wèi)星產(chǎn)生的附加攝動(dòng)。具體模型見公式(B.18)、公式(B.19)。76.4.8地球扁率間接攝動(dòng)地球扁率間接攝動(dòng)是將地球視作扁球體時(shí)地球?qū)πl(wèi)星產(chǎn)生的間接附加攝動(dòng)。具體模型見公式6.4.9廣義相對(duì)論效應(yīng)攝動(dòng)由于廣義相對(duì)論效應(yīng)，衛(wèi)星在地球質(zhì)心為原點(diǎn)的局部慣性坐標(biāo)系中的運(yùn)動(dòng)方程將不同于僅考慮牛頓引力場(chǎng)時(shí)的運(yùn)動(dòng)方程，這種差異可看作衛(wèi)星受到了一個(gè)附加攝動(dòng)。由于太陽引力場(chǎng)對(duì)衛(wèi)星產(chǎn)生的相對(duì)論效應(yīng)攝動(dòng)加速度較小，故可只考慮地球引力場(chǎng)引起的廣義相對(duì)論效應(yīng)攝動(dòng)。具體模型見公式6.4.10經(jīng)驗(yàn)力攝動(dòng)在衛(wèi)星激光測(cè)距數(shù)據(jù)處理中，攝動(dòng)模型不精確或者未考慮攝動(dòng)源的影響[如星際塵埃阻力、引力輻射、坡印亭-羅伯遜(Poynting-Robertson)效應(yīng)、帶電大氣阻力等]導(dǎo)致經(jīng)驗(yàn)力攝動(dòng)(或類阻力攝動(dòng))。根據(jù)實(shí)際需要，可在沿跡向、法向或徑向上引入改正，其周期通常采用衛(wèi)星軌道周期。具體模型見公式6.5解算策略6.5.1基本流程衛(wèi)星激光測(cè)距大地基準(zhǔn)產(chǎn)品解算的基本流程包括數(shù)據(jù)預(yù)處理、參數(shù)解算、精度評(píng)定。解算基本流程圖1解算基本流程6.5.2基本策略基本策略如下：a)所選的衛(wèi)星激光測(cè)距觀測(cè)數(shù)據(jù)解算應(yīng)以地球動(dòng)力學(xué)衛(wèi)星的衛(wèi)星激光測(cè)距標(biāo)準(zhǔn)點(diǎn)數(shù)據(jù)為主，宜采用多星觀測(cè)數(shù)據(jù)聯(lián)合解算；b)觀測(cè)模型應(yīng)符合6.3的規(guī)定；c)動(dòng)力學(xué)模型應(yīng)符合6.4的規(guī)定；8d)重力場(chǎng)模型采用IERS或IAG推薦的最新模型；e)大地基準(zhǔn)參數(shù)應(yīng)與其他末知參數(shù)一同進(jìn)行解算，其中動(dòng)力學(xué)參數(shù)和幾何參數(shù)須根據(jù)實(shí)際情況進(jìn)行參數(shù)個(gè)數(shù)和估計(jì)策略的調(diào)整；f)SINEX解應(yīng)具有一定的時(shí)間連續(xù)性，地球自轉(zhuǎn)參數(shù)解應(yīng)每天1組，衛(wèi)星激光測(cè)距測(cè)站坐標(biāo)解應(yīng)每7天1組，地心運(yùn)動(dòng)解應(yīng)每30天1組。7精度評(píng)定7.1衛(wèi)星激光測(cè)距觀測(cè)數(shù)據(jù)精度評(píng)定衛(wèi)星激光測(cè)距觀測(cè)數(shù)據(jù)精度按公式(3)計(jì)算，即M?=√Dπ…………(3)M?——觀測(cè)量的內(nèi)符合精度；Dπ——觀測(cè)量的方差。7.2地球自轉(zhuǎn)參數(shù)精度評(píng)定地球自轉(zhuǎn)參數(shù)精度按公式(4)計(jì)算，即 (4)MERP,y——地球自轉(zhuǎn)參數(shù)在j分量上的精度；j——地球白轉(zhuǎn)參數(shù)包括極移分量(xp,yp)和日長(zhǎng)dton,j-xp,yg,dton;n——地球自轉(zhuǎn)參數(shù)解的個(gè)數(shù)；?..,——第i個(gè)地球自轉(zhuǎn)參數(shù)解在j分量的估計(jì)值；X:,j——第i個(gè)地球自轉(zhuǎn)參數(shù)解在j分量的參考值。7.3地心運(yùn)動(dòng)參數(shù)精度評(píng)定地心運(yùn)動(dòng)(GCM)參數(shù)精度按公式(5)計(jì)算，即 (5)Mccn,/——地心運(yùn)動(dòng)位移量在j方向上的精度；n——地心運(yùn)動(dòng)位移解的個(gè)數(shù)；8.,——第i個(gè)地心運(yùn)動(dòng)位移解在j方向的估計(jì)值；X.,——第i個(gè)地心運(yùn)動(dòng)位移解在j方向的參考值。7.4測(cè)站坐標(biāo)及速度精度評(píng)定a)測(cè)站(STA)坐標(biāo)及速度參數(shù)估計(jì)精度按公式(6)計(jì)算，即 (6)9Dsr?,;——測(cè)站坐標(biāo)或速度在j方向上的方差；X!——第i個(gè)測(cè)站坐標(biāo)或速度j方向的參考值。MsrA=√MSrAx+MSrA,y+MSrA,(規(guī)范性附錄)觀測(cè)模型A.1潮汐改正A.1.1固體潮改正模型據(jù)瓦爾(Wahr)理論，當(dāng)二次勒夫(Love)下三步修正測(cè)站位置：a)計(jì)算潮汐形變引起的測(cè)站位移，即△rsTA——潮汐形變引起的測(cè)站位移；rsTA——測(cè)站在地固坐標(biāo)系中的位置矢量；rsTA——測(cè)站在地固坐標(biāo)系中的地心距；r;——月亮(j-1)或太陽(j-2)在地固坐標(biāo)系中的位置矢量；r;——月亮(j=1)或太陽(j=2)在地固坐標(biāo)系中的地心距G——萬有引力常量；Me——地球質(zhì)量；GMe——地球引力常數(shù)；GM;——月亮(j=1)或太陽(j=2)在地固坐標(biāo)系中引力常數(shù)。b)計(jì)算與頻率有關(guān)項(xiàng)，此處僅需要考慮k1頻率項(xiàng)(對(duì)應(yīng)杜德森數(shù)為165.555)引起的徑向位移，計(jì)算公式為δhfTA=-0.0253sing'cosq'sin(0,+λ)δhíTA——k?頻率引起的徑向位移；q'——測(cè)站的地心緯度；0?——格林尼治真恒星時(shí)。c)對(duì)尚未做零頻率位移改正的測(cè)站進(jìn)行徑向與北向位移改正，即δNSTA——站心坐標(biāo)系中的北向位移；φ'——測(cè)站的地心緯度。A.1.2極潮改正模型極潮改正模型為xp——極移x分量；X;(t)——t時(shí)刻主要分潮波(j)的天文學(xué)參數(shù)。ka——238.0185μm-2;b)非流體靜態(tài)力學(xué)部分的改正公式為dn——非流體靜態(tài)(nh)力學(xué)部分的對(duì)流層天頂(z)延遲；P?！獪y(cè)站水氣壓。A.2.2映射函數(shù)映射函數(shù)為m(e)——映射函數(shù)；aio——(12100.8±1.9)×10-';a13——(-1847.8±6.5ts——測(cè)站的溫度；φ——測(cè)站的大地緯度；H——測(cè)站的大地高。A.3相對(duì)論效應(yīng)改正模型廣義相對(duì)論效應(yīng)改正通常只考慮由太陽和地球引力場(chǎng)引起的相對(duì)論效應(yīng)改正，即式中：△s——太陽引力場(chǎng)引起的相對(duì)論效應(yīng)改正；△E——地球引力場(chǎng)引起的相對(duì)論效應(yīng)改正；△prn——由太陽和地球引起的相對(duì)論效應(yīng)改正；G——萬有引力常量；r?——太陽到衛(wèi)星的距離；r?——地球到衛(wèi)星的距離；R?——太陽到測(cè)站的距離；R?——地球到測(cè)站的距離；p——衛(wèi)星與測(cè)站的距離；γ——相對(duì)論效應(yīng)修正因子，可作為被估計(jì)量，正常值γ=1。A.4衛(wèi)星質(zhì)心改正表A.1給出了常用地球動(dòng)力學(xué)衛(wèi)星LAGEOS-1/2和Etalon-1/2的質(zhì)心改正模型。其中，測(cè)站不表A.1不同測(cè)站衛(wèi)星質(zhì)心改正CII/T20av——N體攝動(dòng)加速度；0:——面元i的法向與衛(wèi)星到太陽方向之間的夾角；PsR——作用在離太陽一個(gè)天文單位處黑體上的太陽輻射壓強(qiáng)；m——衛(wèi)星的質(zhì)量；CR——衛(wèi)星的表面反射系數(shù)；A——衛(wèi)星的表面積；γ——地影因子,a是太陽被蝕的視面積，b是太陽視面積，當(dāng)衛(wèi)星在日光中時(shí)γ在本影中時(shí)γ=0;G——萬有引力φ——測(cè)站的大地緯度9,——太陽或月亮在地固坐標(biāo)系中的地心緯度b)計(jì)算和頻率有關(guān)的二階項(xiàng)改正。二階項(xiàng)潮汐改正主要考慮21個(gè)長(zhǎng)周期、48個(gè)日周期和2個(gè)δk,—k在頻率?處的值k,與其標(biāo)準(zhǔn)值kzm之差；H,——頻率?項(xiàng)的振幅；△CE?=4.4228×10-?(-0.31460)k2o△Cm——永久潮汐二階改正項(xiàng)；△C,,m——海潮引起的球諧系數(shù)(余弦項(xiàng))改正值；△S,——海潮引起的球諧系數(shù)(正弦項(xiàng))改正值；C——f分潮波振幅(余弦項(xiàng)),取IERS規(guī)范(2010)表格值S,—f分潮波振幅(正弦項(xiàng)),取IERS規(guī)范(2010)表格值；0;(t)——t時(shí)刻?分潮波相位。B.4.4大氣潮攝動(dòng)由大氣潮攝動(dòng)引起的球諧系數(shù)改正為(△C,,m)Ar——大氣潮引起的球諧系數(shù)(余弦項(xiàng))改正值；(△S,m)Ar——大氣潮引起的球諧系數(shù)(正弦項(xiàng))改正值：n——杜德森系數(shù)；gE——地表平均重力加速度；G——萬有引力常數(shù)；Pw——海水的平均密度；δo,m——海潮分潮波相位；0?—μ分潮波相位；C,m——μ分潮大氣潮汐系數(shù)(余弦項(xiàng));S擊?！谭殖贝髿獬毕禂?shù)(正弦項(xiàng))。大氣潮一般只需考慮S?波的改正，n=2、m=2,具體計(jì)算時(shí)可在海潮攝動(dòng)中進(jìn)行改正，即B.4.5地球自轉(zhuǎn)形變附加攝動(dòng)由地球自轉(zhuǎn)形變附加攝動(dòng)引起的球諧系數(shù)改為(△S?,1)go——地球自轉(zhuǎn)形變攝動(dòng)引起地球引力系數(shù)S?,1的變化量；a)計(jì)算地球反照率A?和地球紅外輻射率Em,即t?——周期項(xiàng)的起始?xì)v元；an——地球反照輻射壓攝動(dòng)加速度；agw——紅外輻射壓攝動(dòng)加速度；aeR——地球反照輻射壓攝動(dòng)加速度與紅外輻別壓攝動(dòng)加速度之和；w——積分區(qū)域，即地球被衛(wèi)星可見的部分；au——天文單位；ds——面積元；rs——太陽至地球的距離；ps——衛(wèi)星受照表面的反射系數(shù)；S——地球反照和紅外輻射壓攝動(dòng)中應(yīng)考慮的衛(wèi)星截面積；m——衛(wèi)星質(zhì)量；0s——太陽入射角；α——ds的法線與p間的夾角；p——面積元ds上地球反照和紅外輻射對(duì)衛(wèi)星的壓力方向的矢量，其模為p;在實(shí)際計(jì)算中可采用近似方法，把衛(wèi)星所見表面分成若干個(gè)面積元，計(jì)算出每個(gè)面積元對(duì)衛(wèi)星的反照輻射壓攝動(dòng)加速度(daπ):和紅外輻射壓攝動(dòng)加速度