賈前，李青，劉磊

（1.西北工業(yè)大學(xué)航天學(xué)院，西安 710072；2.陜西省空天飛行器設(shè)計重點(diǎn)實驗室，西安 710072；3.清華大學(xué)航天航空學(xué)院，北京 100084）

0 引言

星載原子鐘是衛(wèi)星導(dǎo)航系統(tǒng)的核心部件，其性能決定了導(dǎo)航系統(tǒng)時間基準(zhǔn)的精度［1-2］。目前，地面光學(xué)原子鐘的穩(wěn)定度即將突破10-19量級［3-4］，JPL 的星載汞離子鐘已通過在軌驗證，比目前應(yīng)用的星載原子鐘精度提高了1 個量級以上，有望應(yīng)用于新一代衛(wèi)星導(dǎo)航系統(tǒng)［5］。然而，空間中變化的引力場以及衛(wèi)星的運(yùn)動會導(dǎo)致星載原子鐘計量的原時產(chǎn)生偏差［6］，不同衛(wèi)星的原時偏差導(dǎo)致了星間鐘差，影響對原子鐘性能極限的利用［7］。新一代導(dǎo)航衛(wèi)星均計劃或已經(jīng)通過確定星間鐘差來建立與維持星座自主時間基準(zhǔn)［8］，以實現(xiàn)導(dǎo)航系統(tǒng)的自主運(yùn)行，但衛(wèi)星原時及星間皮秒鐘差機(jī)理尚不清晰。為滿足新一代導(dǎo)航系統(tǒng)對高精度自主時間基準(zhǔn)的需求，亟需開展衛(wèi)星原時偏差以及星間皮秒鐘差模型研究。

在衛(wèi)星原時偏差方面，國內(nèi)外學(xué)者基于衛(wèi)星鐘的物理模型以及相對論效應(yīng)研究了多種模型?；谛l(wèi)星鐘的物理模型主要利用原子鐘的數(shù)據(jù)進(jìn)行模型參數(shù)求解，進(jìn)而得到鐘差數(shù)據(jù)，精度在ns 量級［9］。常用的模型有多項式模型［10］、灰色模型［11］、譜分析模型［12］、時間序列模型［13］以及組合模型［14］等。文獻(xiàn)［15］提出了二次多項式擬合預(yù)報模型，基于衛(wèi)星鐘的物理特性擬合得到衛(wèi)星鐘差，文獻(xiàn)［9］在二次多項式模型基礎(chǔ)上增加了顯著周期項，并結(jié)合BP神經(jīng)網(wǎng)絡(luò)設(shè)計了超快速鐘差預(yù)報算法實現(xiàn)流程。然而上述模型需要衛(wèi)星的精密鐘差數(shù)據(jù)與原子鐘的隨機(jī)先驗信息，且精度無法滿足新一代導(dǎo)航系統(tǒng)分米級導(dǎo)航定位需求［9］。基于相對論效應(yīng)的模型通過研究衛(wèi)星原時速率與引力場及衛(wèi)星運(yùn)動的關(guān)系得到衛(wèi)星原時偏差。文獻(xiàn)［6］與文獻(xiàn)［16］在開普勒軌道假設(shè)下研究了GPS 和北斗衛(wèi)星導(dǎo)航系統(tǒng)（BDS）衛(wèi)星原子鐘的相對論效應(yīng)，并指出周期性累積鐘差需要根據(jù)衛(wèi)星星歷實時進(jìn)行計算和改正，但未考慮攝動對星載原子鐘的影響。文獻(xiàn)［17］以地球同步軌道衛(wèi)星為研究對象分析了軌道攝動引起的鐘差，文獻(xiàn)［18］估計了地球J2 項攝動引起Galileo 衛(wèi)星的原時偏差幅值，文獻(xiàn)［19］建立了衛(wèi)星原時與坐標(biāo)時轉(zhuǎn)換模型，并指出J2 項攝動引起的原時偏差對于中軌道衛(wèi)星達(dá)到了納秒量級。根據(jù)文獻(xiàn)［6］和文獻(xiàn)［16-19］的研究可知，基于廣義相對論的鐘差模型可實現(xiàn)的精度更高，但仍需要在開普勒軌道假設(shè)的基礎(chǔ)上引入J2 項修正以達(dá)到皮秒級精度。星間鐘差主要基于星間鏈路獲得的觀測數(shù)據(jù)進(jìn)行計算，目前雙向測量工作體制的星間鏈路已經(jīng)用于高精度的鐘差測量［20］。文獻(xiàn)［21］提出了一種利用星間偽距擬合多項式和鐘差擬合多項式聯(lián)合求解星間鐘差算法，精度在5ns 以內(nèi)。文獻(xiàn)［22］針對北斗2 導(dǎo)航系統(tǒng)的星間鏈路建立了同步時分雙工（STDD）體制下星間鐘差測量模型，精度優(yōu)于1ns；文獻(xiàn)［23］研究了STDD體制下星間鐘差的測量誤差模型，并與GPS 采用的輪詢時分雙工體制（PTDD）進(jìn)行了對比，驗證了其性能。文獻(xiàn)［21-23］所采用的方法主要以激光或者微波為信號載體，通過計算兩終端發(fā)射接收信號的時間差來確定鐘差，精度取決于不同路徑中延時修正模型，且未從原時模型角度研究星間皮秒鐘差機(jī)理。文獻(xiàn)［19］在開普勒軌道的假設(shè)下利用衛(wèi)星原時模型建立了星間鐘差模型，能夠有效計算星間鐘差，但該模型并未考慮攝動加速度，誤差會隨時間增加而增大。

本文從廣義相對論出發(fā)，提出一種考慮J2 項攝動的衛(wèi)星原時和星間鐘差模型，根據(jù)衛(wèi)星的位置速度信息計算得到衛(wèi)星原時偏差以及星間鐘差，并通過仿真對模型進(jìn)行驗證，為星載原子鐘的校正和星間鐘差確定提供了一種可行的建模方法。

1 坐標(biāo)系與衛(wèi)星運(yùn)動方程

1.1 坐標(biāo)系描述

為建立衛(wèi)星運(yùn)動與時間的聯(lián)系，首先借助廣義相對論建立衛(wèi)星的運(yùn)動方程。運(yùn)動的描述離不開相應(yīng)的坐標(biāo)系，本文選擇地心參考坐標(biāo)系（Geocentric coordinate reference system，GCRS）建立地球引力場中衛(wèi)星的運(yùn)動模型。GCRS的原點(diǎn)在地球質(zhì)量中心，并且由四維坐標(biāo)X≡(ct，x，y，z)T所定義，其中c為真空中的光速，值為2.997 924 58×108m/s，t為坐標(biāo)系對應(yīng)的坐標(biāo)時，x，y，z為衛(wèi)星的空間位置坐標(biāo)。在GCRS中，兩個確定事件的間隔微分ds可以表示為：

式中：gμυ為GCRS 的度規(guī)張量；xμ，υ(μ，υ=0，1，2，3)為GCRS中的時空坐標(biāo)分量，μ，υ=0時，x0=ct，μ，υ=1，2，3分別對應(yīng)衛(wèi)星的三軸坐標(biāo)分量x，y，z。

度規(guī)張量可以由一組諧波勢表示，諧波勢包含標(biāo)量勢w與矢量勢wE，保留至c-4項，選取度規(guī)張量為［24］：

式中：α，β和λ為索引標(biāo)號，從1取至3；γ=diag(-1，-1，-1) ;為wE的第λ個分量。

矢量勢wE的表達(dá)式為：

式中：G為引力常數(shù)，值為6.673 5×10-11N·m2/kg2；ME為地球的質(zhì)量，為5.974 2×1024kg；r=‖r‖2為衛(wèi)星地心距；r=[x，y，z]T為衛(wèi)星在GCRS 中的位置矢量；SE=[0，0，9.8 × 108]Tm2s為地球單位質(zhì)量的角動量；“(·)×”符號表示矢量叉乘運(yùn)算的反對稱矩陣：

地球引力場中，標(biāo)量勢w由地球引起的引力勢UE和其他天體引起的潮汐勢組成，表達(dá)式為：

考慮地球引力勢至J4項，UE的表達(dá)式為：

式中：RE為地球平均赤道半徑，取值為6 378 km；J2，J3和J4為地球的二至四階帶諧系數(shù)，值分別為1.082 6×10-3，-2.532 7×10-6和-1.619 6×10-6；為簡化描述，定義UJ為地球的二至四階帶諧項引起的引力勢變化的總和。

式中：Mb為天體b 的質(zhì)量，太陽質(zhì)量MS為1.980 4×1030kg，月球質(zhì)量MM為7.336 9×1022kg；rbE為天體b與地球之間的距離；nbE=rbErbE為rbE的單位向量；符號S和M分別代表太陽和月球。

1.2 衛(wèi)星運(yùn)動方程

考慮衛(wèi)星運(yùn)動的相對論效應(yīng)時，可以在牛頓體系下衛(wèi)星的運(yùn)動方程中添加相應(yīng)的相對論修正項，也可以利用度規(guī)張量直接在廣義相對論體系下建立衛(wèi)星的運(yùn)動學(xué)方程，本文選擇第2種方式。

根據(jù)式（1）以及選擇的度規(guī)張量，可以得到兩事件的時空間隔微分為：

衛(wèi)星的運(yùn)動方程以坐標(biāo)時t為變量，因此引入拉格朗日乘子L=dsdt，根據(jù)式（8）可得拉格朗日乘子的平方表達(dá)式為：

式中：v=為衛(wèi)星運(yùn)動速度大??；vx=dxdt，vy=dydt，vz=dzdt分別為衛(wèi)星各軸向的速度大小。

對式（9）進(jìn)行開方并保留至c-4項，可得拉格朗日乘子L的表達(dá)式為：

拉格朗日乘子L應(yīng)滿足Euler-Lagrange方程：

將式（10）代入式（11），并引入非引力項以及日地相對運(yùn)動引起的加速度項［24］，推導(dǎo)可得廣義相對論體系下衛(wèi)星的運(yùn)動方程為：

式中：n=r r為衛(wèi)星位置r的單位向量；rSE為太陽相對地球的位置向量，且為非引力項產(chǎn)生的衛(wèi)星加速度。

2 考慮J2項修正的衛(wèi)星原時以及星間鐘差模型

2.1 衛(wèi)星原時模型

衛(wèi)星在軌運(yùn)行時，星載原子鐘記錄的時間稱為原時τ，原時速率會受到衛(wèi)星運(yùn)動與引力場的影響發(fā)生改變，而坐標(biāo)時僅與選定的坐標(biāo)系相關(guān)［16］，因此可通過建立坐標(biāo)時與原時的關(guān)系得到衛(wèi)星的原時模型。在衛(wèi)星附近的無窮小區(qū)域內(nèi)，度規(guī)張量為常值，且原子鐘相對衛(wèi)星靜止，即dx=dy=dz=0，因此可得：

聯(lián)立式（8）與（13）可得衛(wèi)星原時模型為：

式中：δ(r，v)為由于衛(wèi)星運(yùn)動與引力勢的存在而產(chǎn)生的原時速率偏差。假設(shè)坐標(biāo)時為導(dǎo)航系統(tǒng)的參考時間，此時原時偏差即為對應(yīng)的衛(wèi)星鐘差，如圖1所示。

圖1 引力場中原時偏差示意圖Fig.1 Schematic of proper time deviation in the gravitational field

僅考慮地球引力勢中的牛頓主項，忽略其他天體引起的潮汐勢，即標(biāo)量勢w簡化為：

在地球引力場內(nèi)，衛(wèi)星所在位置處的引力勢主要為地球產(chǎn)生的影響，并且J2 項的量級遠(yuǎn)大于高階項，因此在引力勢中引入地球非球形J2 項攝動，得到考慮J2項修正的衛(wèi)星原時模型為：

衛(wèi)星的狀態(tài)也可以通過軌道六要素來表示，將式（17）中對應(yīng)的變量轉(zhuǎn)換為軌道六要素，即可得到基于軌道六要素的衛(wèi)星原時模型。

衛(wèi)星地心距r的表達(dá)式為：

式中：a為衛(wèi)星所在軌道的半長軸；e為軌道偏心率；E為偏近點(diǎn)角。

衛(wèi)星的z軸位置分量可表示為：

式中：i為軌道傾角；f為真近點(diǎn)角；ω為近地點(diǎn)幅角。

衛(wèi)星的地心距與速度大小滿足：

聯(lián)立式（18）和（20）可得：

在軌道偏心率較小的情況下，根據(jù)軌道攝動理論，考慮J2項攝動，相關(guān)軌道參數(shù)可以表示為：

式中：a0，i0，e0，ω0和E0分別為不受攝動影響的衛(wèi)星軌道對應(yīng)的半長軸，軌道傾角，偏心率，近地點(diǎn)幅角以及偏近點(diǎn)角。

將式（18），（19），（21）以及式（22）聯(lián)立，忽略高階小量，代入式（17）可得到基于軌道六要素的衛(wèi)星原時模型為：

根據(jù)式（23）可以看到，衛(wèi)星運(yùn)行過程中原時與坐標(biāo)時產(chǎn)生的差異由四部分組成，前兩項表示僅考慮牛頓主項加速度時衛(wèi)星在軌運(yùn)動引起的相對論效應(yīng)，后兩項為J2項攝動引起的相對論效應(yīng)。

2.2 星間鐘差模型

根據(jù)式（17）可以看到，當(dāng)衛(wèi)星的狀態(tài)不同時，衛(wèi)星原時相對參考時間的差異不同。而導(dǎo)航系統(tǒng)通常由多顆衛(wèi)星組成，這些衛(wèi)星分布在不同軌道面或者同一軌道面的不同位置處，因此導(dǎo)航衛(wèi)星的狀態(tài)間存在差異，在軌運(yùn)行時不同的衛(wèi)星原時會產(chǎn)生對應(yīng)的星間鐘差，如圖2所示。

圖2 星間鐘差示意圖Fig.2 Schematic of inter-satellite clock difference

考慮導(dǎo)航系統(tǒng)中兩顆衛(wèi)星A和B，假設(shè)在坐標(biāo)時t0時刻，兩顆衛(wèi)星的星載原子鐘讀數(shù)即原時分別為τA0和τB0，初始鐘差為Δτ0=τB0-τA0，則在t1時刻兩顆衛(wèi)星的原時可以表示為：

進(jìn)一步對兩顆衛(wèi)星的原時作差，可以得到衛(wèi)星A和B的星間鐘差為：

式（25）即為考慮J2 項修正的星間鐘差模型，該模型分為四部分：（1）衛(wèi)星初始鐘差；（2）衛(wèi)星地心距差異引起的星間鐘差；（3）衛(wèi)星速度差異引起的星間鐘差；（4）J2項攝動引起的星間鐘差。

3 仿真分析

3.1 衛(wèi)星軌道誤差仿真

北斗衛(wèi)星導(dǎo)航系統(tǒng)（BDS）是我國自主建立的全球衛(wèi)星導(dǎo)航系統(tǒng)，可在全球范圍內(nèi)全天候、全天時為各類用戶提供高精度、高可靠定位、導(dǎo)航、授時服務(wù)。本節(jié)以BDS 中北斗三號MEO-01 衛(wèi)星為研究對象，對建立的衛(wèi)星運(yùn)動方程進(jìn)行仿真分析，非引力項加速度考慮為太陽光壓引起的加速度，計算方式參考文獻(xiàn)［19］。仿真中衛(wèi)星的初始軌道要素如表1所示，仿真時長為兩個軌道周期，仿真步長為1 s，仿真時為坐標(biāo)時。

表1 北斗三號-MEO-01衛(wèi)星初始軌道要素Table 1 Initial orbital elements of BeiDou-3 MEO-01 satellite

考慮式（12）右端的全部加速度項時，攝動加速度引起的衛(wèi)星位置誤差以及速度誤差分別如圖3和圖4 所示，兩者均隨時間的增加而增大。對于MEO-01 衛(wèi)星而言，在兩個軌道周期時，位置誤差已經(jīng)達(dá)到了25.48 km，速度誤差達(dá)到了3.78 m/s，根據(jù)衛(wèi)星原時模型可知，這些誤差均會導(dǎo)致原時與坐標(biāo)時之間產(chǎn)生偏差。

圖4 攝動加速度引起的速度總誤差Fig.4 Total velocity error induced by perturbation acceleration

僅考慮J2 項攝動加速度時，MEO-01 衛(wèi)星的位置誤差以及速度誤差如圖5 和圖6 所示。根據(jù)仿真結(jié)果可知，J2 項攝動加速度引起的位置誤差與速度誤差與模型考慮的全部攝動加速度引起的誤差變化趨勢相同。兩個軌道周期時，J2 項攝動加速度引起的位置誤差為21.11 km，速度誤差為3.22 m/s。與圖3和圖4仿真結(jié)果對比可知，在考慮的攝動加速度中，J2 項攝動加速度對MEO-01 衛(wèi)星的影響最大，并且對于地球引力場中的衛(wèi)星，J2 項的影響會隨著軌道高度的減小而增大，因此，對于高精度時間基準(zhǔn)的建立，J2項攝動加速度引起的原時偏差不可忽略。

圖5 僅考慮J2項攝動加速度時位置誤差Fig.5 Position error only considering J2 perturbation acceleration

圖6 僅考慮J2項攝動加速度時速度誤差Fig.6 Velocity error only considering J2 perturbation acceleration

3.2 衛(wèi)星鐘差仿真

本節(jié)根據(jù)建立的衛(wèi)星原時模型對衛(wèi)星鐘差進(jìn)行仿真分析。仿真中利用衛(wèi)星運(yùn)動方程（12）得到衛(wèi)星狀態(tài)，將狀態(tài)信息代入式（14）計算得到的原時偏差并將其作為實際鐘差，與考慮J2 項修正的衛(wèi)星原時模型（17）得到的鐘差進(jìn)行對比，進(jìn)而分析模型精度，仿真參數(shù)與3.1節(jié)相同。

實際鐘差、考慮J2 項修正的模型（17）以及簡化模型（16）的計算結(jié)果如圖7 所示。根據(jù)仿真結(jié)果可知，衛(wèi)星原時與坐標(biāo)時的偏差即衛(wèi)星鐘差隨著時間的增加而增大，在偏心率接近0時，衛(wèi)星鐘差與坐標(biāo)時之間近似為線性關(guān)系，與文獻(xiàn)［19］結(jié)論一致，在兩個軌道周期時，MEO-01 衛(wèi)星的鐘差已經(jīng)達(dá)到了2.21×10-5s。

圖7 衛(wèi)星鐘差計算結(jié)果Fig.7 Satellite clock difference calculation results

兩種模型的誤差如圖8所示，隨著時間的增加，兩種模型的計算誤差都逐漸增大。在兩個軌道周期內(nèi)，考慮J2項修正的模型誤差峰值為26.97 ps，簡化模型的誤均差峰值為47.95 ps，兩者比實際鐘差小6個數(shù)量級，即兩種模型都能用于計算衛(wèi)星鐘差。但由于簡化模型未考慮J2 項攝動對衛(wèi)星原時速率的影響，計算得到的鐘差波動大，誤差的均方根為27.23 ps，而考慮J2 項修正的模型計算結(jié)果波動較小，誤差的均方根為16.27 ps，提高了40.26%。仿真結(jié)果表明考慮J2 項修正的衛(wèi)星原時模型可以有效提高衛(wèi)星鐘差的計算精度。

圖8 原時模型誤差Fig.8 Proper time model error

模型（23）為衛(wèi)星鐘差的軌道要素表示形式，該形式表明J2 項引起的鐘差中包含常值項與周期項，軌道傾角為0°時，常值項達(dá)到最大值，此最大值與軌道半長軸成反比；當(dāng)軌道傾角為54.73°時，常值項為0。周期項的周期為衛(wèi)星軌道周期的一半，幅值與軌道半長軸及軌道傾角有關(guān)。在小偏心率情況下，真近點(diǎn)角可通過軌道平均角速度n近似為：

對周期項積分可得：

周期項幅值與軌道半長軸和軌道傾角的關(guān)系如圖9所示，隨著軌道傾角的增加，J2項攝動加速度引起的周期項鐘差逐漸增大，并隨著軌道半長軸的增加而快速減小。

圖9 J2項攝動加速度引起的周期項鐘差幅值與軌道傾角及半長軸關(guān)系Fig.9 Relationship between the amplitude of periodic proper time deviation caused by J2 perturbation acceleration，orbital inclination，and semi-major axis

3.3 星間鐘差仿真

3.2 節(jié)的衛(wèi)星原時仿真結(jié)果表明，衛(wèi)星在軌運(yùn)行期間星載原子鐘會產(chǎn)生鐘差，并且鐘差與衛(wèi)星軌道相關(guān)，進(jìn)而導(dǎo)致不同衛(wèi)星間產(chǎn)生星間鐘差。本節(jié)對兩衛(wèi)星間的鐘差進(jìn)行仿真研究。文獻(xiàn)［19］的研究表明兩衛(wèi)星軌道半長軸不同時，星間鐘差分為線性與周期項兩部分，線性部分會隨著半長軸差異的增加而迅速累積，但通常衛(wèi)星發(fā)射前都會根據(jù)目標(biāo)軌道對原子鐘進(jìn)行頻率校正［16］，因此本節(jié)只研究衛(wèi)星軌道半長軸相同的情況下星間鐘差的變化規(guī)律，且假設(shè)初始時刻星載原子鐘均經(jīng)過校準(zhǔn)，即Δτ0=0。與衛(wèi)星鐘差仿真相同，仿真中兩顆衛(wèi)星的狀態(tài)由運(yùn)動方程（12）得到，衛(wèi)星原時通過式（14）計算，作差后得到實際星間鐘差，與式（25）計算得到的星間鐘差進(jìn)行對比，分析模型誤差。

仿真中一顆衛(wèi)星軌道參數(shù)與表1 相同，另一顆衛(wèi)星初始真近點(diǎn)角為180°，即兩衛(wèi)星軌道面相同，在軌道平面內(nèi)初始位置相差180°。兩衛(wèi)星間的鐘差以及模型誤差如圖10 和圖11 所示。仿真結(jié)果表明，兩個軌道周期內(nèi)，模型的計算結(jié)果與理論鐘差基本吻合，星間鐘差呈周期性變化，周期與衛(wèi)星軌道周期相同，仿真中兩顆衛(wèi)星的鐘差峰值為5.89 ns；模型誤差也近似為周期變化，最大誤差為8.86×10-2ps。

圖10 星間鐘差Fig.10 Inter-satellite clock difference

圖11 星間鐘差模型誤差Fig.11 Inter-satellite clock difference model error

為說明J2 項修正的有效性，采用文獻(xiàn)［19］的模型對相同參數(shù)的衛(wèi)星進(jìn)行星間鐘差仿真，模型誤差如圖12所示。仿真結(jié)果表明，無J2項修正的星間鐘差模型誤差隨時間逐漸增大，在兩個軌道周期時，誤差為16.37 ps，遠(yuǎn)大于本文提出的考慮J2 項修正的星間鐘差模型的誤差。

除了中軌（MEO）衛(wèi)星，BDS 星座還包含高軌同步（GEO）衛(wèi)星及高軌傾斜（IGSO）衛(wèi)星。美國的GPS星座由近圓形軌道（e=0.02）的中軌衛(wèi)星組成［6］，選取對應(yīng)軌道上的一顆衛(wèi)星，另一顆衛(wèi)星與其在同一軌道平面上，初始位置相差180°，根據(jù)星間鐘差模型計算這兩顆衛(wèi)星間的鐘差，衛(wèi)星的軌道參數(shù)以及計算結(jié)果如表2 所示。根據(jù)表2 的仿真結(jié)果可以看到，四種軌道平面中星間鐘差峰值從小到大依次為BDS-GEO，BDS-MEO，BDS-IGSO 和GPS-MEO。這是因為隨著軌道偏心率的增加，衛(wèi)星速度與位置的變化引起的原時偏差增大，進(jìn)而引起星間鐘差峰值變大。對于BDS-IGSO 和GPS-MEO 這兩個軌道平面，由于其半長軸、偏心率和軌道傾角的影響，兩個軌道平面上衛(wèi)星受到的太陽及月球的影響較大，而本文建立的星間鐘差模型（25）忽略了其余天體的影響，因此模型誤差的峰值大于其余兩個軌道平面的誤差峰值，但仍在皮秒量級。這些算例的計算結(jié)果表明考慮J2 項修正的星間鐘差模型誤差最大在皮秒量級，能夠準(zhǔn)確計算星間鐘差。

表2 不同軌道的星間鐘差計算結(jié)果Table 2 Calculation results of inter-satellite clock difference in different orbits

4 結(jié)論

本文針對新一代導(dǎo)航衛(wèi)星對高精度時間基準(zhǔn)的需求，在廣義相對論體系下建立了衛(wèi)星的運(yùn)動方程，考慮地球引力場內(nèi)衛(wèi)星受到的J2 項攝動影響，提出了包含J2 項修正的衛(wèi)星原時模型以及星間鐘差模型，模型能夠有效地計算出衛(wèi)星在引力場中運(yùn)動時產(chǎn)生的原時偏差以及星間鐘差。相比基于開普勒軌道假設(shè)得到的模型，考慮J2 項修正可以提高衛(wèi)星原時以及星間鐘差的計算精度。仿真結(jié)果表明，衛(wèi)星運(yùn)動過程中產(chǎn)生的原時偏差與衛(wèi)星軌道要素緊密關(guān)聯(lián)，并隨著時間增加而增大，考慮J2 項修正計算得到的原時偏差與實際鐘差更為接近。同時本文提出的星間鐘差模型能夠在皮秒量級的精度上計算星間鐘差。本文工作可為導(dǎo)航系統(tǒng)的時間基準(zhǔn)建立與自主維持提供理論支持。