相对论体系下卫星原时及星间钟差机理建模.pdf

第 45 卷第 3 期 2024 年 3 月宇 航 学 报Journal of AstronauticsNo.32024MarchVol.45相对论体系下卫星原时及星间钟差机理建模贾前1，2，李青3，刘磊1，2（1.西北工业大学航天学院，西安 710072；2.陕西省空天飞行器设计重点实验室，西安 710072；3.清华大学航天航空学院，北京 100084）摘要：面向新一代导航卫星对高精度自主时间基准的需求，针对皮秒级星间钟差机理不清晰问题，在广义相对论体系下提出一种卫星原时以及星间皮秒钟差模型，该模型可基于卫星状态计算星载原子钟相对参考时钟的偏差以及星间钟差。首先从地心参考坐标系的度规张量出发，利用广义相对论建立了卫星的运动模型；然后，根据时空间隔的不变性，提出考虑地球非球形J2摄动项的卫星原时模型，并进一步给出了包含J2项修正的星间钟差模型。最后以北斗卫星导航系统和GPS为研究对象对建立的模型进行了仿真分析，为新一代导航系统的时间校准和同步提供了一种可行的建模方法。仿真结果表明，卫星原时模型与星间钟差模型的精度都在皮秒量级，能够准确地计算星载原子钟在引力场中的原时偏差以及星间钟差。关键词：导航卫星；星间钟差；时间校准与同步；广义相对论；北斗卫星导航系统中图分类号：V412 文献标识码：A 文章编号：1000-1328（2024）03-0366-10 DOI：10.3873/j.issn.1000-1328.2024.03.004Mechanism Modeling of Satellite Proper Time and Inter-satellite Clock Difference in Relativistic SystemJIA Qian1，2，LI Qing3，LIU Lei1，2（1.School of Astronautics，Northwestern Polytechnical University，Xi an 710072，China；2.Shaanxi Aerospace Flight Vehicle Design Key Laboratory，Xi an 710072，China；3.School of Aerospace Engineering，Tsinghua University，Beijing 100084，China）Abstract:To meet the demand for high-precision autonomous time reference of new generation navigation satellites，and addressing the unclear mechanism of picosecond level inter-satellite clock difference，the proper time and inter-satellite picosecond clock difference models are proposed using general relativity theory.The models utilize the satellite state to compute the on-board atomic clock deviation relative to the reference clock and inter-satellite clock difference.Firstly，beginning with the metric tensor in geocentric coordinate reference system，the satellite motion equations are derived using general relativity theory.Then，according to the invariance of space-time interval，the proper time model considering the non-spherical J2 perturbation is established and the inter-satellite clock difference model with J2 correction is further given.Finally，the simulation is conducted by taking BeiDou Navigation Satellite System and GPS constellations for examples，providing a feasible modeling solution for time correction and synchronization in new generation navigation systems.The simulation results demonstrate the accuracy of satellite proper time and inter-satellite clock difference models is in the order of picoseconds，which can accurately calculate the proper time deviation of the on-board atomic clock in the gravitational field and inter-satellite clock difference.Key words:Navigation satellite；Inter-satellite clock difference；Time calibration and synchronization；General relativity theory；BeiDou Navigation Satellite System收稿日期：2023-07-27；修回日期：2023-09-08基金项目：国家自然科学基金（52075446，52305117）366第 3 期贾前等：相对论体系下卫星原时及星间钟差机理建模0引言星载原子钟是卫星导航系统的核心部件，其性能决定了导航系统时间基准的精度1-2。目前，地面光学原子钟的稳定度即将突破10-19量级3-4，JPL的星载汞离子钟已通过在轨验证，比目前应用的星载原子钟精度提高了1个量级以上，有望应用于新一代卫星导航系统5。然而，空间中变化的引力场以及卫星的运动会导致星载原子钟计量的原时产生偏差6，不同卫星的原时偏差导致了星间钟差，影响对原子钟性能极限的利用7。新一代导航卫星均计划或已经通过确定星间钟差来建立与维持星座自主时间基准8，以实现导航系统的自主运行，但卫星原时及星间皮秒钟差机理尚不清晰。为满足新一代导航系统对高精度自主时间基准的需求，亟需开展卫星原时偏差以及星间皮秒钟差模型研究。在卫星原时偏差方面，国内外学者基于卫星钟的物理模型以及相对论效应研究了多种模型。基于卫星钟的物理模型主要利用原子钟的数据进行模型参数求解，进而得到钟差数据，精度在 ns 量级9。常用的模型有多项式模型10、灰色模型11、谱分析模型12、时间序列模型13以及组合模型14等。文献 15 提出了二次多项式拟合预报模型，基于卫星钟的物理特性拟合得到卫星钟差，文献 9 在二次多项式模型基础上增加了显著周期项，并结合BP神经网络设计了超快速钟差预报算法实现流程。然而上述模型需要卫星的精密钟差数据与原子钟的随机先验信息，且精度无法满足新一代导航系统分米级导航定位需求9。基于相对论效应的模型通过研究卫星原时速率与引力场及卫星运动的关系得到卫星原时偏差。文献 6 与文献 16 在开普勒轨道假设下研究了GPS和北斗卫星导航系统（BDS）卫星原子钟的相对论效应，并指出周期性累积钟差需要根据卫星星历实时进行计算和改正，但未考虑摄动对星载原子钟的影响。文献 17 以地球同步轨道卫星为研究对象分析了轨道摄动引起的钟差，文献 18估计了地球J2项摄动引起Galileo卫星的原时偏差幅值，文献 19 建立了卫星原时与坐标时转换模型，并指出J2项摄动引起的原时偏差对于中轨道卫星达到了纳秒量级。根据文献 6 和文献 16-19的研究可知，基于广义相对论的钟差模型可实现的精度更高，但仍需要在开普勒轨道假设的基础上引入J2项修正以达到皮秒级精度。星间钟差主要基于星间链路获得的观测数据进行计算，目前双向测量工作体制的星间链路已经用于高精度的钟差测量20。文献 21 提出了一种利用星间伪距拟合多项式和钟差拟合多项式联合求解星间钟差算法，精度在5ns以内。文献 22 针对北斗2导航系统的星间链路建立了同步时分双工（STDD）体制下星间钟差测量模型，精度优于 1ns；文献 23 研究了 STDD体制下星间钟差的测量误差模型，并与GPS采用的轮询时分双工体制（PTDD）进行了对比，验证了其性能。文献 21-23 所采用的方法主要以激光或者微波为信号载体，通过计算两终端发射接收信号的时间差来确定钟差，精度取决于不同路径中延时修正模型，且未从原时模型角度研究星间皮秒钟差机理。文献 19 在开普勒轨道的假设下利用卫星原时模型建立了星间钟差模型，能够有效计算星间钟差，但该模型并未考虑摄动加速度，误差会随时间增加而增大。本文从广义相对论出发，提出一种考虑J2项摄动的卫星原时和星间钟差模型，根据卫星的位置速度信息计算得到卫星原时偏差以及星间钟差，并通过仿真对模型进行验证，为星载原子钟的校正和星间钟差确定提供了一种可行的建模方法。1坐标系与卫星运动方程1.1坐标系描述为建立卫星运动与时间的联系，首先借助广义相对论建立卫星的运动方程。运动的描述离不开相应的坐标系，本文选择地心参考坐标系（Geocentric coordinate reference system，GCRS）建立地球引力场中卫星的运动模型。GCRS的原点在地球质量中心，并且由四维坐标X(ct，x，y，z)T所定义，其中c为真空中的光速，值为2.997 924 58108 m/s，t为坐标系对应的坐标时，x，y，z为卫星的空间位置坐标。在GCRS中，两个确定事件的间隔微分ds可以表示为：(ds)2=gdxdx（1）式中：g为 GCRS的度规张量；x，(，=0，1，2，3)为 GCRS 中的时空坐标分量，=0时，x0=ct，=1，2，3分别对应卫星的三轴坐标分量x，y，z。度规张量可以由一组谐波势表示，谐波势包含标量势w与矢量势wE，保留至c-4项，选取度规张量为24：367宇航学报第 45 卷g00=1-2wc2+2w2c4g0=g0=-4c3wEg=()1+2wc2（2）式中：，和为索引标号，从1取至3；=diag(-1，)-1，-1；wE为wE的第个分量。矢量势wE的表达式为：wE=-GME2r3rSE（3）式中：G为引力常数，值为 6.673 510-11 Nm2/kg2；ME为地球的质量，为5.974 21024 kg；r=r2为卫星地心距；r=x，y，zT为卫星在GCRS中的位置矢量；SE=0，0，9.8 108Tm2s为地球单位质量的角动量；“()”符号表示矢量叉乘运算的反对称矩阵：r=0-zyz0-x-yx0（4）地球引力场中，标量势w由地球引起的引力势UE和其他天体引起的潮汐势utidalE组成，表达式为：w=UE+utidalE（5）考虑地球引力势至J4项，UE的表达式为：UE=GMEr 1-J2R2E2r2()3z2r2-1-J3R3E2r3(5z3r3-3zr)-J4R4E8r4()35z4r4-30z2r2+3GMEr-UJ（6）式中：RE为地球平均赤道半径，取值为 6 378 km；J2，J3和J4为地球的二至四阶带谐系数，值分别为1.082 610-3，-2.532 710-6和-1.619 610-6；为简化描述，定义UJ为地球的二至四阶带谐项引起的引力势变化的总和。其余天体只考虑太阳与月球的影响时，utidalE的表达式为：utidalE=b EGMb2r3bE3()nTbEr2-r2=GMS2r3SE3(nTSEr)2-r2+GMM2r3ME3(nTMEr)2-r2（7）式中：Mb为天体b的质量，太阳质量MS为1.980 41030 kg，月球质量MM为7.336 91022 kg；rbE为天体b与地球之间的距离；nbE=rbE rbE为rbE的单位向量；符号S和M分别代表太阳和月球。1.2卫星运动方程考虑卫星运动的相对论效应时，可以在牛顿体系下卫星的运动方程中添加相应的相对论修正项，也可以利用度规张量直接在广义相对论体系下建立卫星的运动学方程，本文选择第2种方式。根据式（1）以及选择的度规张量，可以得到两事件的时空间隔微分为：(ds)2=(1-2wc2+2w2c4)(cdt)2-(1+2wc2)(dx)2+(dy)2+(dz)2+8c2w1Edxdt+8c2w2Edydt+8c2w3Edzdt（8）卫星的运动方程以坐标时t为变量，因此引入拉格朗日乘子L=ds dt，根据式（8）可得拉格朗日乘子的平方表达式为：L2=(1-2wc2+2w2c4)c2-(1+2wc2)v2+8c2(w1Evx+w2Evy+w3Evz)（9）式中：v=v2x+v2y+v2z为卫星运动速度大小；vx=dx dt，vy=dy dt，vz=dz dt分别为卫星各轴向的速度大小。对式（9）进行开方并保留至c-4项，可得拉格朗日乘子L的表达式为：L=c1-2wc2+2w2c4-v2c2-2wv2c4+8c4()w1Evx+w2Evy+w3Evz=c 1-wc2+w22c4-v22c2-v48c4-3wv22c4+4c4()w1Evx+w2Evy+w3Evz（10）拉格朗日乘子L应满足Euler-Lagrange方程：ddt(Lx)-Lx=0，x y，z（11）将式（10）代入式（11），并引入非引力项以及日地相对运动引起的加速度项24，推导可得广义相对论体系下卫星的运动方程为：368第 3 期贾前等：相对论体系下卫星原时及星间钟差机理建模r=-GMEr3r-UJr+GMEc2r3(4GMEr-rTr)r+4(rTr)r+2GMEc2r3(r)SE-3(nTSE)n+3GMSc2r3SE(rSE)rSEr+b EGMbr3bE3()nTbEr nbE-r+ang（12）式中：n=r r为卫星位置r的单位向量；rSE为太阳相对地球的位置向量，且rSE=rSE2；ang为非引力项产生的卫星加速度。2考虑J2项修正的卫星原时以及星间钟差模型2.1卫星原时模型卫星在轨运行时，星载原子钟记录的时间称为原时，原时速率会受到卫星运动与引力场的影响发生改变，而坐标时仅与选定的坐标系相关16，因此可通过建立坐标时与原时的关系得到卫星的原时模型。在卫星附近的无穷小区域内，度规张量为常值，且原子钟相对卫星静止，即dx=dy=dz=0，因此可得：(ds)2=(cd)2（13）联立式（8）与（13）可得卫星原时模型为：d=dsc=1+()r，vdt()r，v=-wc2+w22c4-v22c2-v48c4-3wv22c4+4c4(w1vx+w2vy+w3vz)（14）式中：(r，v)为由于卫星运动与引力势的存在而产生的原时速率偏差。假设坐标时为导航系统的参考时间，此时原时偏差即为对应的卫星钟差，如图1所示。仅考虑地球引力势中的牛顿主项，忽略其他天体引起的潮汐势，即标量势w简化为：w=GMEr（15）进一步忽略矢量势分量wE，卫星的原时模型保留至c-2项，得到简化的卫星原时模型为：d=1+(r，v)dt=(1-GMEc2r-v22c2)dt（16）在地球引力场内，卫星所在位置处的引力势主要为地球产生的影响，并且J2项的量级远大于高阶项，因此在引力势中引入地球非球形J2项摄动，得到考虑J2项修正的卫星原时模型为：d=1+(r，v)dt=1-GMEc2r+J2GMER2E2c2r3(3z2r2-1)-v22c2 dt（17）卫星的状态也可以通过轨道六要素来表示，将式（17）中对应的变量转换为轨道六要素，即可得到基于轨道六要素的卫星原时模型。卫星地心距r的表达式为：r=a(1-ecosE)（18）式中：a为卫星所在轨道的半长轴；e为轨道偏心率；E为偏近点角。卫星的z轴位置分量可表示为：z=rsinisin(f+)（19）式中：i为轨道倾角；f为真近点角；为近地点幅角。卫星的地心距与速度大小满足：v22-GMEr=-GME2a（20）联立式（18）和（20）可得：v2=GMEa1+ecosE1-ecosE（21）在轨道偏心率较小的情况下，根据轨道摄动理论，考虑J2项摄动，相关轨道参数可以表示为：a=a0+3J2R2E2a0sin2i0cos2()0+fe=e0+3J2R2E2a20()1-32sin2i0cosf+14sin2i0cos()20+f 712sin2i0cos()20+3fecosE=e0cosE0+3J2R2E2a20()1-32sin2i0+5J2R2E4a20sin2i0cos2()0+f（22）图1引力场中原时偏差示意图Fig.1Schematic of proper time deviation in the gravitational field369宇航学报第 45 卷式中：a0，i0，e0，0和E0分别为不受摄动影响的卫星轨道对应的半长轴，轨道倾角，偏心率，近地点幅角以及偏近点角。将式（18），（19），（21）以及式（22）联立，忽略高阶小量，代入式（17）可得到基于轨道六要素的卫星原时模型为：d=1-3GME2a0c2-2GMEa0c2e0cosE0-7GMEJ2R2E2a30c2(1-32sin2i0)-GMEJ2R2Esin2i0a30c2cos2()0+fdt（23）根据式（23）可以看到，卫星运行过程中原时与坐标时产生的差异由四部分组成，前两项表示仅考虑牛顿主项加速度时卫星在轨运动引起的相对论效应，后两项为J2项摄动引起的相对论效应。2.2星间钟差模型根据式（17）可以看到，当卫星的状态不同时，卫星原时相对参考时间的差异不同。而导航系统通常由多颗卫星组成，这些卫星分布在不同轨道面或者同一轨道面的不同位置处，因此导航卫星的状态间存在差异，在轨运行时不同的卫星原时会产生对应的星间钟差，如图2所示。考虑导航系统中两颗卫星A和B，假设在坐标时t0时刻，两颗卫星的星载原子钟读数即原时分别为A0和B0，初始钟差为0=B0-A0，则在t1时刻两颗卫星的原时可以表示为：A=A0+t0t11+()rA，vAdtB=B0+t0t11+()rB，vBdt=A0+0+t0t11+()rB，vBdt（24）进一步对两颗卫星的原时作差，可以得到卫星A和B的星间钟差为：=B-A=0+GMEc2t0t1()1rA-1rBdt-12c2t0t1()v2A-v2Bdt+J2GMER2E2c2t0t1()3z2Br5B-3z2Ar5A+1r3A-1r3Bdt（25）式（25）即为考虑J2项修正的星间钟差模型，该模型分为四部分：（1）卫星初始钟差；（2）卫星地心距差异引起的星间钟差；（3）卫星速度差异引起的星间钟差；（4）J2项摄动引起的星间钟差。3仿真分析3.1卫星轨道误差仿真北斗卫星导航系统（BDS）是我国自主建立的全球卫星导航系统，可在全球范围内全天候、全天时为各类用户提供高精度、高可靠定位、导航、授时服务。本节以BDS中北斗三号MEO-01卫星为研究对象，对建立的卫星运动方程进行仿真分析，非引力项加速度考虑为太阳光压引起的加速度，计算方式参考文献 19。仿真中卫星的初始轨道要素如表1所示，仿真时长为两个轨道周期，仿真步长为1 s，仿真时为坐标时。考虑式（12）右端的全部加速度项时，摄动加速度引起的卫星位置误差以及速度误差分别如图 3和图 4 所示，两者均随时间的增加而增大。对于MEO-01卫星而言，在两个轨道周期时，位置误差已经达到了25.48 km，速度误差达到了3.78 m/s，根据卫星原时模型可知，这些误差均会导致原时与坐标时之间产生偏差。图2星间钟差示意图Fig.2Schematic of inter-satellite clock difference表1北斗三号-MEO-01卫星初始轨道要素Table 1Initial orbital elements of BeiDou-3 MEO-01 satellite参数半长轴a/km偏心率e轨道倾角i/()升交点赤经/()近地点幅角/()真近点角f/()数值27 9060.001 25655.76100.66296.117 50370第 3 期贾前等：相对论体系下卫星原时及星间钟差机理建模仅考虑J2项摄动加速度时，MEO-01卫星的位置误差以及速度误差如图5和图6所示。根据仿真结果可知，J2项摄动加速度引起的位置误差与速度误差与模型考虑的全部摄动加速度引起的误差变化趋势相同。两个轨道周期时，J2项摄动加速度引起的位置误差为21.11 km，速度误差为3.22 m/s。与图3和图4仿真结果对比可知，在考虑的摄动加速度中，J2项摄动加速度对MEO-01卫星的影响最大，并且对于地球引力场中的卫星，J2项的影响会随着轨道高度的减小而增大，因此，对于高精度时间基准的建立，J2项摄动加速度引起的原时偏差不可忽略。3.2卫星钟差仿真本节根据建立的卫星原时模型对卫星钟差进行仿真分析。仿真中利用卫星运动方程（12）得到卫星状态，将状态信息代入式（14）计算得到的原时偏差并将其作为实际钟差，与考虑J2项修正的卫星原时模型（17）得到的钟差进行对比，进而分析模型精度，仿真参数与3.1节相同。实际钟差、考虑J2项修正的模型（17）以及简化模型（16）的计算结果如图 7 所示。根据仿真结果可知，卫星原时与坐标时的偏差即卫星钟差随着时间的增加而增大，在偏心率接近 0时，卫星钟差与坐标时之间近似为线性关系，与文献 19 结论一致，在两个轨道周期时，MEO-01卫星的钟差已经达到了2.2110-5 s。两种模型的误差如图8所示，随着时间的增加，两种模型的计算误差都逐渐增大。在两个轨道周期内，考虑J2项修正的模型误差峰值为26.97 ps，简化模型的误均差峰值为47.95 ps，两者比实际钟差小6个数量级，即两种模型都能用于计算卫星钟差。104246810104时间/s00.51.01.52.02.5仅考虑J2项卫星位置总误差/m图5仅考虑J2项摄动加速度时位置误差Fig.5Position error only considering J2 perturbation acceleration246810104时间/s00.51.01.52.02.53.03.5仅考虑J2项卫星速度总误差/(ms1)图6仅考虑J2项摄动加速度时速度误差Fig.6Velocity error only considering J2 perturbation acceleration104104105105时间/s02468102.52.01.51.00.504.436 604.436 621.057 5201.057 515原时与坐标时偏差/s实际钟差考虑J2项修正模型简化模型图7卫星钟差计算结果Fig.7Satellite clock difference calculation results24681000.51.01.52.02.53.0考虑全部加速度卫星位置总误差/m104104时间/s图3摄动加速度引起的位置总误差Fig.3Total position error induced by perturbation acceleration246810104时间/s01234考虑全部加速度卫星速度总误差/(ms1)图4摄动加速度引起的速度总误差Fig.4Total velocity error induced by perturbation acceleration371宇航学报第 45 卷但由于简化模型未考虑J2项摄动对卫星原时速率的影响，计算得到的钟差波动大，误差的均方根为27.23 ps，而考虑J2项修正的模型计算结果波动较小，误差的均方根为16.27 ps，提高了40.26%。仿真结果表明考虑J2项修正的卫星原时模型可以有效提高卫星钟差的计算精度。模型（23）为卫星钟差的轨道要素表示形式，该形式表明J2项引起的钟差中包含常值项与周期项，轨道倾角为0时，常值项达到最大值，此最大值与轨道半长轴成反比；当轨道倾角为54.73时，常值项为0。周期项的周期为卫星轨道周期的一半，幅值与轨道半长轴及轨道倾角有关。在小偏心率情况下，真近点角可通过轨道平均角速度n近似为：f=nt=GME/a30t（26）对周期项积分可得：tJ2=path -GMEJ2R2Esin2i0a30c2cos()20+2ntdt=-GMEa30J2R2Esin2i02c2sin(20+2nt)（27）周期项幅值与轨道半长轴和轨道倾角的关系如图9所示，随着轨道倾角的增加，J2项摄动加速度引起的周期项钟差逐渐增大，并随着轨道半长轴的增加而快速减小。3.3星间钟差仿真3.2节的卫星原时仿真结果表明，卫星在轨运行期间星载原子钟会产生钟差，并且钟差与卫星轨道相关，进而导致不同卫星间产生星间钟差。本节对两卫星间的钟差进行仿真研究。文献 19 的研究表明两卫星轨道半长轴不同时，星间钟差分为线性与周期项两部分，线性部分会随着半长轴差异的增加而迅速累积，但通常卫星发射前都会根据目标轨道对原子钟进行频率校正16，因此本节只研究卫星轨道半长轴相同的情况下星间钟差的变化规律，且假设初始时刻星载原子钟均经过校准，即0=0。与卫星钟差仿真相同，仿真中两颗卫星的状态由运动方程（12）得到，卫星原时通过式（14）计算，作差后得到实际星间钟差，与式（25）计算得到的星间钟差进行对比，分析模型误差。仿真中一颗卫星轨道参数与表1相同，另一颗卫星初始真近点角为180，即两卫星轨道面相同，在轨道平面内初始位置相差180。两卫星间的钟差以及模型误差如图 10 和图 11 所示。仿真结果表明，两个轨道周期内，模型的计算结果与理论钟差基本吻合，星间钟差呈周期性变化，周期与卫星轨道周期相同，仿真中两颗卫星的钟差峰值为5.89 ns；模型误差也近似为周期变化，最大误差为8.8610-2 ps。1011104时间/s02468105432101模型误差/s考虑J2项修正模型简化模型图8原时模型误差Fig.8Proper time model error10101040123453.02.52.01.51.00.50J2项扰动引起的原时偏差/s卫星轨道半长轴/kmi=10i=30i=55i=90图9J2项摄动加速度引起的周期项钟差幅值与轨道倾角及半长轴关系Fig.9Relationship between the amplitude of periodic proper time deviation caused by J2 perturbation acceleration，orbital inclination，and semi-major axis109104时间/s02468106420246理论钟差模型计算结果星间钟差/s图10星间钟差Fig.10Inter-satellite clock difference372第 3 期贾前等：相对论体系下卫星原时及星间钟差机理建模为说明J2项修正的有效性，采用文献 19 的模型对相同参数的卫星进行星间钟差仿真，模型误差如图12所示。仿真结果表明，无J2项修正的星间钟差模型误差随时间逐渐增大，在两个轨道周期时，误差为16.37 ps，远大于本文提出的考虑J2项修正的星间钟差模型的误差。除了中轨（MEO）卫星，BDS星座还包含高轨同步（GEO）卫星及高轨倾斜（IGSO）卫星。美国的GPS星座由近圆形轨道（e=0.02）的中轨卫星组成6，选取对应轨道上的一颗卫星，另一颗卫星与其在同一轨道平面上，初始位置相差180，根据星间钟差模型计算这两颗卫星间的钟差，卫星的轨道参数以及计算结果如表2所示。根据表2的仿真结果可以看到，四种轨道平面中星间钟差峰值从小到大依次为BDS-GEO，BDS-MEO，BDS-IGSO和GPS-MEO。这是因为随着轨道偏心率的增加，卫星速度与位置的变化引起的原时偏差增大，进而引起星间钟差峰值变大。对于BDS-IGSO和GPS-MEO这两个轨道平面，由于其半长轴、偏心率和轨道倾角的影响，两个轨道平面上卫星受到的太阳及月球的影响较大，而本文建立的星间钟差模型（25）忽略了其余天体的影响，因此模型误差的峰值大于其余两个轨道平面的误差峰值，但仍在皮秒量级。这些算例的计算结果表明考虑J2项修正的星间钟差模型误差最大在皮秒量级，能够准确计算星间钟差。4结论本文针对新一代导航卫星对高精度时间基准的需求，在广义相对论体系下建立了卫星的运动方程，考虑地球引力场内卫星受到的J2项摄动影响，提出了包含J2项修正的卫星原时模型以及星间钟差模型，模型能够有效地计算出卫星在引力场中运动时产生的原时偏差以及星间钟差。相比基于开普勒轨道假设得到的模型，考虑J2项修正可以提高卫星原时以及星间钟差的计算精度。仿真结果表明，卫星运动过程中产生的原时偏差与卫星轨道要素紧密关联，并随着时间增加而增大，考虑J2项修正计算得到的原时偏差与实际钟差更为接近。同时本文提出的星间钟差模型能够在皮秒量级的精度上计算星间钟差。本文工作可为导航系统的时间基准建立与自主维持提供理论支持。参 考 文 献1 李文涛，边少锋，任青阳，等.基于粒子群优化核极限学习机的北斗超快速钟差预报 J.宇航学报，2019，40（9）：1080-1088.LI Wentao，BIAN Shaofeng，REN Qingyang，et al.Kernel extreme learning machine based on particle swarm optimization 1014104时间/s0246810108642024模型误差/s图11星间钟差模型误差Fig.11Inter-satellite clock difference model error1012104时间/s0246810201510505文献19模型误差/s图12文献 19 的星间钟差模型误差Fig.12Inter-satellite clock difference model error in Ref.19表2不同轨道的星间钟差计算结果Table 2Calculation results of inter-satellite clock difference in different orbitsNORAD 卫星编号43001453444420429486轨道平面BDS-MEOBDS-GEOBDS-IGSOGPS-MEO轨道要素(aeif)(27 906 km0.001 25655.76100.66296.120)(42 164 km0.000 2662.00310.2184.430)(42 167 km0.001 8757.0549.43206.840)(26 571 km0.010 455 154.69196.1229.180)钟差峰值5.89 ns1.53 ns10.79 ns47.87 ns模型误差峰值8.8610-2 ps6.1910-2 ps24.86 ps23.53 ps373宇航学报第 45 卷for prediction of Beidou ultra-rapid clock offsetJ.Journal of Astronautics，2019，40（9）：1080-1088.2 田婕.GPS/BDS原子钟性能分析及钟差预报模型研究 D.西安：长安大学，2015.TIAN Jie.Study on the GPS/BDS atomic clock performance and clock offset prediction model D.Xi an：Chang an University，2015.3 CAMPBELL S L，HUTSON R B，MARTI G E，et al.A Fermi-degenerate three-dimensional optical lattice clockJ.Science，2017，358（6359）：90-94.4 SCHIOPPO M，BROWN R C，MCGREW W F，et al.Ultrastable optical clock with two cold-atom ensemblesJ.Nature Photonics，2017，11：48-52.5 BURT E A，PRESTAGE J D，TJOELKER R L，et al.Demonstration of a trapped-ion atomic clock in spaceJ.Nature，2021，595：43-47.6 刘丽丽，王跃科，陈建云，等.导航星座自主时间基准的相对论效应 J.宇航学报，2015，36（4）：470-476.LIU Lili，WANG Yueke，CHEN Jianyun，et al.Effects of relativity in autonomous time reference for navigation constellationJ.Journal of Astronautics，2015，36（4）：470-476.7 戴辉.基于墨子号卫星的激光时间传递的研究 D.合肥：中国科学技术大学，2018.DAI Hui.Experimental research on laser time transter based on Micius satelliteD.Hefei：University of Science and Technology of China，2018.8 孙乐园.导航星座星间链路时间基准自主建立与维持 D.长沙：国防科技大学，2020.SUN Leyuan.Autonomously establish and keep the inter-satellite link time reference for navigation satellite constellationsD.Changsha：National University of Defense Technology，2020.9 黄观文，崔博斌，张勤，等.附加周期和神经网络补偿的实时钟差预报模型 J.宇航学报，2018，39（1）：83-88.HUANG Guanwen，CUI Bobin，ZHANG Qin，et al.Real-time clock offset prediction model with periodic and neural network corrections J.Journal of Astronautics，2018，39（1）：83-88.10 魏道坤.卫星钟差预报模型研究 D.西安：长安大学，2014.WEI Daokun.Study on the satellite clock bias forecast modelD.Xi an：Chang an University，2014.11 谭小容，许江宁，何泓洋，等.GM（1，1）模型初始条件优化及其在 GNSS 钟差预报中的应用 J.大地测量与地球动力学，2022，42（9）：919-924.TAN Xiaorong，XU Jiangning，HE Hongyang，et al.Optimized initial condition of GM（1，1）model and its application in GNSS clock offset prediction J.Journal of Geodesy and Geodynami