CN114911153B - 基于星-地超长基线干涉测量的世界时ut1测定方法 - Google Patents
基于星-地超长基线干涉测量的世界时ut1测定方法 Download PDFInfo
- Publication number
- CN114911153B CN114911153B CN202210333431.4A CN202210333431A CN114911153B CN 114911153 B CN114911153 B CN 114911153B CN 202210333431 A CN202210333431 A CN 202210333431A CN 114911153 B CN114911153 B CN 114911153B
- Authority
- CN
- China
- Prior art keywords
- vlbi
- satellite
- ground
- observation station
- observation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005305 interferometry Methods 0.000 title description 6
- 239000005433 ionosphere Substances 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 8
- 239000005436 troposphere Substances 0.000 claims abstract description 8
- 238000012937 correction Methods 0.000 claims description 15
- 230000000694 effects Effects 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000033001 locomotion Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000001934 delay Effects 0.000 claims description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 4
- 229910052739 hydrogen Inorganic materials 0.000 claims description 4
- 239000001257 hydrogen Substances 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于星‑地超长基线VLBI观测世界时UT1测定方法及***,通过地面VLBI观测站与卫星VLBI观测站组成星地超长基线,计算目标射电源信号到地面VLBI观测站与卫星VLBI观测站的观测时延;最后将误差修正过的观测时延进行处理,解算UT1。本发明可以突破国土限制,摆脱地球尺寸对地面VLBI基线长度限制,提高VLBI角分辨率,从而提高UT1的解算精度;同时由于卫星的高度较高,使其受到的对流层及电离层影响较小,可以在更高的频率观测。
Description
技术领域
本发明属于卫星授时技术领域,具体涉及一种基于星-地超长基线干涉测量的世界时UT1测定方法。
背景技术
世界时(Universal Time,UT1)是以地球自转定义的时间尺度,是构成标准时间,即协调世界时(Universal Time Coordinated,UTC)的重要组成部分。在UTC的产生与保持过程中,通过UTC闰秒机制来保证UTC与UT1的偏差在±0.9s以内。世界时UT1与极移、章动改正项组成了地球定向参数(Earth Orientation Parameter,EOP),UT1还是实现天球与地球坐标变换的核心参数,是深空探测、卫星导航、航天器测控等领域中不可或缺的重要参量。
当前世界时UT1测量主要使用VLBI技术,VLBI技术无需考虑动力学因素效应,是唯一能够直接测定UT1的空间大地测量技术,也是目前唯一能够解算全部地球定向参数的空间技术。VLBI技术因具有极高的测角精度和长期的稳定性,在解算EOP、维持地球、天球参考框架方面发挥着不可替代的作用。根据电磁波的干涉可知,望远镜的角分辨率约等于观测波长除以有效口径。对于VLBI,其空间分辨率约等于观测的波长除以基线长度(即两个VLBI射电望远镜之间的距离)。在观测波长一定的情况下,基线长度是影响UT1精度的主要因素。
当前国际基线(1万公里)的UT1精度约10微秒,不能满足未来深空探测对UT1的精度需求。为了提高UT1的测量精度,现有技术中采用地面多台站VLBI组网观测进行观测以增加VLBI基线长度,在该方法中面多台站VLBI组网=是由全球分布的多个VLBI台站参加,一般是8-10个站。该方法虽然极大的利用了地球的尺寸,使VLBI的测量精度有所提升,但最长基线也只能接近地球直径,无法对UT1的测量精度有质的提升;且海外不易建站,需要多个国家共同协商观测,极易受到影响,时效性较差,大约滞后2周。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于星-地超长基线干涉测量的世界时UT1测定方法。本发明要解决的技术问题通过以下技术方案实现:
本发明提供的一种基于星-地超长基线VLBI观测世界时UT1测定方法,应用于星-地VLBI***,星-地VLBI***包括地面VLBI观测站以及卫星VLBI观测站,地面VLBI观测站与卫星VLBI观测站组成星地超长基线,卫星VLBI观测站以及地面VLBI观测站上设置有VLBI射电望远镜,基于星-地超长基线VLBI观测世界时UT1测定方法包括:
卫星VLBI观测站以及地面VLBI观测站,用于基于相同的观测纲要,在同一观测时刻观测多组已知精确位置的目标射电源,并接收同一观测时刻同一目标射电源发送的同一目标射电源信号;
卫星VLBI观测站,用于将接收的同一目标射电源信号的第一接收时刻采用RDF数据格式进行数据采集,并以特定频率实时发送至地面观测站;
地面VLBI观测站,用于将接收的RDF数据格式进行解析以及记录,得到第一接收时刻,并根据自身接收同一目标射电源信号的第二接收时刻,计算同一目标射电源信号到达卫星VLBI观测站以及地面VLBI观测站的几何时延;
地面VLBI观测站,用于对几何时延进行误差修正,得到实际时延,并建立实际时延与UT1的依赖关系表达式;
地面VLBI观测站,用于对依赖关系表达式进行求解,得到地球自转角,并根据地球自转角测定世界时UT1。
可选的,卫星VLBI观测站包括:高轨卫星、高精度测定轨设备、高精度可展开VLBI射电望远镜、信号接收处理***、高稳定频率基准***、高速RDF数据记录***以及高码率传输***;高精度可展开VLBI射电望远镜,用于观测目标射电源;高精度测定轨设备,用于对高轨卫星进行精密定轨;高稳定频率基准***采用星载氢钟,用于在高轨卫星上产生卫星所使用的基准信号;高速RDF数据记录***,用于在高轨卫星上采用RDF数据格式进行记录数据;高码率传输***,用于采用RDF数据格式进行记录的数据的进行以特定频率实时发送至地面VLBI观测站。
其中,几何时延表示为:
其中,τg为地面站到卫星观测站的时间延迟,T1和T2分别为同一目标射电源信号到达卫星VLBI观测站和地面VLBI观测站的时间,下标1表示最先接收到目标射电源信号的VLBI观测站,基线向量是地面VLBI观测站G到卫星VLBI观测站S的距离,方向为地面观测站G到卫星观测站S,卫星空间坐标矢量S又是卫星轨道参数Ei(i=1,2,3,4,5,6)和重力场与非重力场Pk(k=1,2,..,m,..,n)等参数的函数,这样的几何关系就变成了动力学关系,为目标射电源单位方向向量,可由目标射电源在天球参考架下的赤经赤纬得到。
其中,地面VLBI观测站,进一步用于将观测时刻的地面VLBI观测站的坐标从国际地球参考系转换到国际天球参考系;
其中,转化坐标的转化矩阵表示为:
[ICRS]=Q(t)·R(t)·W(t)·[ITRS]
其中,Q(t)、R(t)、W(t)分别由天极在ICRS中运动的岁差章动旋转矩阵、地球绕极轴旋转的地球自转角旋转矩阵、极点在ITRS中运动形成的极移转换矩阵;
国际天球参考系下的几何时延;
其中,国际天球参考系下的几何时延表示为:
其中,地面VLBI观测站,进一步用于消除国际天球参考系下的几何时延中因行星引力、电离层、对流层、地面VLBI观测站与卫星VLBI观测站之间的钟差、仪器设备、地球潮汐以及广义及狭义相对论效应影响造成的误差,得到国际天球参考系下的实际时延。
其中,国际天球参考系下的实际时延表示为:
τobs=τg+τion+τtro+τclk+τinstr+τdef+ε
其中,τion为电离层延迟;τtro为对流层延迟;τclk为两VLBI观测站的钟差;τinstr为设备时延;τdef为天线形变时延;ε为其他时延。
其中,依赖关系表达式表示为:
其中,
其中,地面VLBI观测站,进一步用于根据地球自转角与UT1之间的线性关系,得到测定的世界时UT1;
线性关系表示为:
ERA=2π(0.7790572732640+1.00273781191135448Tu)
其中,Tu为J2000.0起算的UT1的儒略日计数,Tu=Julian UT1 date-2451545.0。
可选的,地面VLBI观测站,在进行观测目标射电源之前,用于:
确定卫星观测站的卫星轨道;
根据卫星轨道确定卫星VLBI观测站和地面VLBI观测站的观测频段;
在观测频段内观测多个射电源;
根据观测到的射电源的分布、源流量、方位以及俯仰角选择卫星地面VLBI观测站以及地面VLBI观测站能够在相同时段接收到射电源信号的目标射电源;
根据选择的目标射电源的观测时间、观测频率以及观测带宽,生成UT1的观测纲要。
第二方面,本发明提供的一种基于星-地超长基线VLBI观测世界时UT1测定***,包括地面VLBI观测站以及卫星VLBI观测站,地面VLBI观测站与卫星VLBI观测站组成星地超长基线,卫星VLBI观测站以及地面VLBI观测站上设置有VLBI射电望远镜;
卫星VLBI观测站以及地面VLBI观测站,用于基于相同的观测纲要,在同一观测时刻观测多组已知精确位置的目标射电源,并接收同一观测时刻同一目标射电源发送的同一目标射电源信号;
卫星VLBI观测站,用于将接收的同一目标射电源信号的第一接收时刻采用RDF数据格式进行采集,并以特定频率实时发送至地面观测站;
地面VLBI观测站,用于将接收的RDF数据格式的数据进行解析以及记录,得到第一接收时刻,并根据自身接收同一目标射电源信号的第二接收时刻,计算同一目标射电源信号到达卫星VLBI观测站以及地面VLBI观测站的几何时延;
地面VLBI观测站,用于对几何时延进行误差修正,得到实际时延,并建立实际时延与UT1的依赖关系表达式;
地面VLBI观测站,用于对依赖关系表达式进行求解,得到地球自转角,并根据地球自转角测定世界时UT1。
本发明提供了一种基于星-地超长基线VLBI观测世界时UT1测定方法及***,通过地面VLBI观测站与卫星VLBI观测站组成星地超长基线,计算目标射电源信号到地面VLBI观测站与卫星VLBI观测站的观测时延;最后将误差修正过的观测时延进行处理,解算UT1。本发明可以突破国土限制,摆脱地球尺寸对地面VLBI基线长度限制,提高VLBI角分辨率,从而提高UT1的解算精度;同时由于卫星的高度较高,使其受到的对流层及电离层影响较小,可以在更高的频率观测。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的星-地VLB***示意图;
图2是本发明实施例提供的一种基于星-地超长基线VLBI观测世界时UT1测定方法的流程示意图;
图3是本发明实施例提供的星-地UT1观测解算仿真流程图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
如图1所示,本发明提供的一种基于星-地超长基线VLBI观测世界时UT1测定方法应用于星-地VLBI***,星-地VLBI***包括地面VLBI观测站以及卫星VLBI观测站,地面VLBI观测站与卫星VLBI观测站组成星地超长基线,卫星VLBI观测站以及地面VLBI观测站上设置有VLBI射电望远镜。
其中,卫星VLBI观测站包括:高轨卫星、高精度测定轨设备、高精度可展开VLBI射电望远镜、信号接收处理***、高稳定频率基准***、高速RDF数据记录***以及高码率传输***;高精度可展开VLBI射电望远镜,用于观测目标射电源;高精度测定轨设备,用于对高轨卫星进行精密定轨;高稳定频率基准***采用星载氢钟,用于在高轨卫星上产生卫星所使用的基准信号;高速RDF数据记录***,用于在高轨卫星上采用RDF数据格式进行记录数据;高码率传输***,用于采用RDF数据格式进行记录的数据的进行以特定频率实时发送至地面VLBI观测站。
RDF数据格式是一种由地面VLBI观测站记录并输出的一种基于磁盘存储的数据格式。RDF格式由于只有一个文件头的额外数据记录,其余的记录都是原始甚长基线干涉测量数据,相较而言RDF具有更高的记录效率和较好的可靠性。其中高精度测定轨设备提供精密轨道,优于1cm的定轨精度,使卫星可作为移动测站,另外通过提高频段来减小卫星天线口径,使得星上搭载可行。
星-地基线VLBI观测时,最主要的是卫星精密轨道的确定,观测时刻卫星位置的任何误差将直接影响时间延迟及延迟率等VLBI观测量,导致由观测量解算的UT1精度恶化,因此卫星位置的确定必须达到很高精度。卫星轨道的误差源主要为太阳光压以及相对论效应。目前VLBI的观测精度已达到1-3cm,SLR测距精度已达到2-5cm,联合VLBI及SLR以及高精度测定轨技术等多技术融合定轨,有望使定轨精度优于1cm。
本发明的卫星搭载VLBI望远镜,与地面VLBI望远镜组成星-地VLBI,提高基线长度,从而提高UT1解算精度,保障深空等领域的进一步发展。
值得说明的是:利用甚长基线干涉测量时,其角分辨率θ由工作波长λ和两面垂直于目标射电源方向的投影距离D决定,公式如下:
从上式可以明显看出,当波长一定时,角分辨率与两天线距离(基线长度)呈反比,即提高基线长度,可以提高角分辨率。由于卫星VLBI与地面VLBI之间的基线距离可以做到比地基VLBI基线距离更长,故星-地VLBI可通过增大基线长度,提高角分辨率,从而达到更高的UT1观测解算精度。
参考图2,本发明提供的一种基于星-地超长基线VLBI观测世界时UT1测定方法包括:
S1,卫星VLBI观测站以及地面VLBI观测站,用于基于相同的观测纲要,在同一观测时刻观测多组已知精确位置的目标射电源,并接收同一观测时刻同一目标射电源发送的同一目标射电源信号;
本发明在观测之间,地面VLBI观测站,在进行观测目标射电源之前,用于:
确定卫星观测站的卫星轨道;
根据卫星轨道确定卫星VLBI观测站和地面VLBI观测站的观测频段;
在观测频段内观测多个射电源;
根据观测到的射电源的分布、源流量、方位以及俯仰角选择卫星地面VLBI观测站以及地面VLBI观测站能够在相同时段接收到射电源信号的目标射电源;
根据选择的目标射电源的观测时间、观测频率以及观测带宽,生成UT1的观测纲要。
值得说明的是:在实际时一般是观测50个射电源左右;不同的射电源观测的次数时长不同;根据射电源的流量、方位以及俯仰、分布来判断出不同的射电源,在实际过程中一次观测时长为2到3个小时,一共观测50到60个目标射电源左右。地面VLBI观测站会将观测纲要上传至卫星VLBI观测站。
S2,卫星VLBI观测站,用于将接收的同一目标射电源信号的第一接收时刻采用RDF数据格式进行数据采集,并以特定频率实时发送至地面观测站;
S3,地面VLBI观测站,用于将接收的RDF数据格式的数据进行解析以及记录,得到第一接收时刻,并根据自身接收同一目标射电源信号的第二接收时刻,计算同一目标射电源信号到达卫星VLBI观测站以及地面VLBI观测站的几何时延;
星-地VLBI观测时,需要地面VLBI观测站和卫星VLBI观测站共同观测同一目标射电源。同一目标射电源的信号到达两VLBI观测站时间不同。由于目标射电源离两VLBI观测站太远,故信号从目标射电源到两VLBI观测站的传播路径可视为平行线,根据几何关系,同一信号到达两观测站的几何时延可以表示为:
其中,τg为地面站到卫星观测站的时间延迟,T1和T2分别为同一目标射电源信号到达卫星VLBI观测站和地面VLBI观测站的时间,下标1表示最先接收到目标射电源信号的VLBI观测站,基线向量是地面VLBI观测站G到卫星VLBI观测站S的距离,方向为地面观测站G到卫星观测站S,卫星空间坐标矢量S又是卫星轨道参数Ei(i=1,2,3,4,5,6)和重力场与非重力场Pk(k=1,2,..,m,..,n)等参数的函数,这样的几何关系就变成了动力学关系,为目标射电源单位方向向量,可由目标射电源在天球参考架下的赤经赤纬得到。
由于计算两观测站之间时延差,需对观测时刻的观测站坐标进行转换,本发明的地面VLBI观测站,进一步用于将观测时刻的地面VLBI观测站的坐标从国际地球参考系转换到国际天球参考系;
其中,转化坐标的转化矩阵表示为:
[ICRS]=Q(t)·R(t)·W(t)·[ITRS]
其中,Q(t)、R(t)、W(t)分别由天极在ICRS中运动的岁差章动旋转矩阵、地球绕极轴旋转的地球自转角旋转矩阵、极点在ITRS中运动形成的极移转换矩阵。
同样的,几何时延也需要转化,本发明的地面VLBI观测站,进一步用于将位于国际地球参考系下的几何时延转化为国际天球参考系下的几何时延;
其中,国际天球参考系下的几何时延表示为:
S4,地面VLBI观测站,用于对几何时延进行误差修正,得到实际时延,并建立实际时延与UT1的依赖关系表达式;
值得说明的是:由于信号从目标射电源到观测站的过程中会受到行星引力、电离层、对流层等影响,以及两观测站之间的钟差、仪器设备的影响,地球潮汐、广义及狭义相对论效应等影响。因此本发明需要对几何延迟进行误差修正。
地面VLBI观测站,进一步用于消除国际天球参考系下的几何时延中因行星引力、电离层、对流层、地面VLBI观测站与卫星VLBI观测站之间的钟差、仪器设备、地球潮汐以及广义及狭义相对论效应影响造成的误差,得到国际天球参考系下的实际时延。
国际天球参考系下的实际时延表示为:
τobs=τg+τion+τtro+τclk+τinstr+τdef+ε
其中τg为几何时延,主要与VLBI观测站的相对位置、目标射电源的位置及地球的方向有关。τion为电离层延迟,主要是电磁波在介质中传播的速度与真空中不同,且信号在介质中传播时会因为弯曲效应增加传播距离,可通过双频观测减轻电离层干扰。τtro为对流层延迟,主要与传播路径的温度、压力和湿度有关,并随观测地点的方向和时间变化。由于卫星高度原因,卫星VLBI受到的大气影响要比地面VLBI小一些。τclk为两VLBI观测站的钟差,主要与观测站间不同的氢钟振荡速率有关。τdef为天线形变时延,主要与温度、风荷载及重力等造成的结构形变有关。τinstr为设备时延,主要是信号从馈源接收到数据记录之间的设备所引入的时延,一般情况下,需要进行相应的信号标定。ε为其他时延,主要是由于相对论效应等引入的时延。在星-地超长基线VLBI观测时,相对论效应对***的固有影响在本质上比对地基VLBI的影响要复杂,原因在于卫星VLBI观测天线相对于地面VLBI天线在高速运动。在星-地VLBI观测时,卫星-地面之间基线两个端点的相对运动是地球自转和卫星运动的结果。由于两观测天线之间的相对速度而引起对卫星VLBI与地面VLBI之间时间延迟观测量的相对论影响,称作下垂弯曲改正,它依赖于卫星相对于地面的速度,所以需要进行额外的相对论效应修正。
若源矢量K→、基线矢量B→、极移、岁差章动等均为已知量,则VLBI几何时延与UT1的依赖关系为:
其中,
S5,地面VLBI观测站,用于对依赖关系表达式进行求解,得到地球自转角,并根据地球自转角测定世界时UT1。
本步骤中地面VLBI观测站,进一步用于根据地球自转角与UT1之间的线性关系,得到测定的世界时UT1;
线性关系表示为:
ERA=2π(0.7790572732640+1.00273781191135448Tu)
其中,Tu为J2000.0起算的UT1的儒略日计数,Tu=Julian UT1 date-2451545.0。
值得说明的是,本发明通过测定多组已知精确位置目标射电源,得到测量时延,然后通过建立时延模型,误差方程,再通过平差算法即可精确测定地球自转角ERA。ERA和世界时则是简单的线性关系,进行最小二乘拟合,即可得到世界时UT1。
本发明提供一种基于星-地超长基线VLBI观测世界时UT1测定方法及***,通过地面VLBI观测站与卫星VLBI观测站组成星地超长基线,计算目标射电源信号到地面VLBI观测站与卫星VLBI观测站的观测时延;最后将误差修正过的观测时延进行处理,解算UT1。本发明可以突破国土限制,摆脱地球尺寸对地面VLBI基线长度限制,提高VLBI角分辨率,从而提高UT1的解算精度;同时由于卫星的高度较高,使其受到的对流层及电离层影响较小,可以在更高的频率观测。
本发明提供的一种基于星-地超长基线VLBI观测世界时UT1测定***,包括地面VLBI观测站以及卫星VLBI观测站,地面VLBI观测站与卫星VLBI观测站组成星地超长基线,卫星VLBI观测站以及地面VLBI观测站上设置有VLBI射电望远镜;
卫星VLBI观测站以及地面VLBI观测站,用于基于相同的观测纲要,在同一观测时刻观测多组已知精确位置的目标射电源,并接收同一观测时刻同一目标射电源发送的同一目标射电源信号;
卫星VLBI观测站,用于将接收的同一目标射电源信号的第一接收时刻采用RDF数据格式进行采集,并以特定频率实时发送至地面;
地面VLBI观测站,用于将接收的RDF数据格式的数据进行解析以及记录,得到第一接收时刻,并根据自身接收同一目标射电源信号的第二接收时刻,计算同一目标射电源信号到达卫星VLBI观测站以及地面VLBI观测站的几何时延;
地面VLBI观测站,用于对几何时延进行误差修正,得到实际时延,并建立实际时延与UT1的依赖关系表达式;
地面VLBI观测站,用于对依赖关系表达式进行求解,得到地球自转角,并根据地球自转角测定世界时UT1。
参考图3,结合软件仿真对本发明过程进行说明。
1、.由于该方法现处于预研阶段,不能利用搭载VLBI望远镜的卫星进行实际的UT1测量解算实验,因此前期采用部分仿真数据进行计算,选取一颗搭测定轨设备的卫星,通过地面测轨网获取该卫星的精密轨道,得到其精密的坐标,利用仿真软件模拟该卫星VLBI望远镜与地面VLBI望远镜观测同一目标射电源的观测时延;
2、对该模拟观测时延进行相关处理,得到各通道的频谱值,再经过后处理从相位信息中提取时延值,对该时延进行钟差校准、设备校准、对流层、电离层时延校准,相对论效应等修正,并进行带宽综合处理,提取出经过校准的目标射电源观测群时延;
3、再结合温度、湿度、气压数据,生成分析软件所需的NGS数据格式文件;
4、将生成的NGS文件通过VIEVS软件利用最小二乘算法进行解算,得到UT1。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (10)
1.一种基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,应用于星-地VLBI***,所述星-地VLBI***包括地面VLBI观测站以及卫星VLBI观测站,所述地面VLBI观测站与所述卫星VLBI观测站组成星地超长基线,所述卫星VLBI观测站以及地面VLBI观测站上设置有VLBI射电望远镜,所述基于星-地超长基线VLBI观测世界时UT1测定方法包括:
所述卫星VLBI观测站以及地面VLBI观测站,用于基于相同的观测纲要,在同一观测时刻观测多组已知精确位置的目标射电源,并接收同一观测时刻同一目标射电源发送的同一目标射电源信号;
所述卫星VLBI观测站,用于将接收的同一目标射电源信号的第一接收时刻采用RDF数据格式进行数据采集,并以特定频率实时发送至地面观测站;
所述地面VLBI观测站,用于将接收的RDF数据格式进行解析以及记录,得到第一接收时刻,并根据自身接收同一目标射电源信号的第二接收时刻,计算同一目标射电源信号到达卫星VLBI观测站以及地面VLBI观测站的几何时延;
所述地面VLBI观测站,用于对所述几何时延进行误差修正,得到实际时延,并建立实际时延与UT1的依赖关系表达式;
所述地面VLBI观测站,用于对所述依赖关系表达式进行求解,得到地球自转角,并根据地球自转角测定世界时UT1。
2.根据权利要求1所述的基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,所述卫星VLBI观测站包括:高轨卫星、高精度测定轨设备、高精度可展开VLBI射电望远镜、信号接收处理***、高稳定频率基准***、高速RDF数据记录***以及高码率传输***;所述高精度可展开VLBI射电望远镜,用于观测目标射电源;所述高精度测定轨设备,用于对高轨卫星进行精密定轨;所述高稳定频率基准***采用星载氢钟,用于在高轨卫星上产生卫星所使用的基准信号;所述高速RDF数据记录***,用于在高轨卫星上采用RDF数据格式进行记录数据;所述高码率传输***,用于采用RDF数据格式进行记录的数据的进行以特定频率实时发送至地面VLBI观测站。
3.根据权利要求1所述的基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,所述几何时延表示为:
5.根据权利要求4所述的基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,所述地面VLBI观测站,进一步用于消除国际天球参考系下的几何时延中因行星引力、电离层、对流层、地面VLBI观测站与卫星VLBI观测站之间的钟差、仪器设备、地球潮汐以及广义及狭义相对论效应影响造成的误差,得到国际天球参考系下的实际时延。
6.根据权利要求5所述的基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,所述国际天球参考系下的实际时延表示为:
τobs=τg+τion+τtro+τclk+τinstr+τdef+ε
其中,τion为电离层延迟;τtro为对流层延迟;τclk为两VLBI观测站的钟差;τinstr为设备时延;τdef为天线形变时延;ε为其他时延。
8.根据权利要求7所述的基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,所述地面VLBI观测站,进一步用于根据地球自转角与UT1之间的线性关系,得到测定的世界时UT1;
所述线性关系表示为:
ERA=2π(0.7790572732640+1.00273781191135448Tu)
其中,Tu为J2000.0起算的UT1的儒略日计数,Tu=Julian UT1 date-2451545.0。
9.根据权利要求8所述的基于星-地超长基线VLBI观测世界时UT1测定方法,其特征在于,所述地面VLBI观测站,在进行观测目标射电源之前,用于:
确定卫星观测站的卫星轨道;
根据卫星轨道确定卫星VLBI观测站和地面VLBI观测站的观测频段;
在所述观测频段内观测多个射电源;
根据观测到的射电源的分布、源流量、方位以及俯仰角选择卫星地面VLBI观测站以及地面VLBI观测站能够在相同时段接收到射电源信号的目标射电源;
根据选择的目标射电源的观测时间、观测频率以及观测带宽,生成UT1的观测纲要。
10.一种基于星-地超长基线VLBI观测世界时UT1测定***,其特征在于,包括地面VLBI观测站以及卫星VLBI观测站,所述地面VLBI观测站与所述卫星VLBI观测站组成星地超长基线,所述卫星VLBI观测站以及地面VLBI观测站上设置有VLBI射电望远镜;
所述卫星VLBI观测站以及地面VLBI观测站,用于基于相同的观测纲要,在同一观测时刻观测多组已知精确位置的目标射电源,并接收同一观测时刻同一目标射电源发送的同一目标射电源信号;
所述卫星VLBI观测站,用于将接收的同一目标射电源信号的第一接收时刻采用RDF数据格式进行采集,并以特定频率实时发送至地面观测站;
所述地面VLBI观测站,用于将接收的RDF数据格式的数据进行解析以及记录,得到第一接收时刻,并根据自身接收同一目标射电源信号的第二接收时刻,计算同一目标射电源信号到达卫星VLBI观测站以及地面VLBI观测站的几何时延;
所述地面VLBI观测站,用于对所述几何时延进行误差修正,得到实际时延,并建立实际时延与UT1的依赖关系表达式;
所述地面VLBI观测站,用于对所述依赖关系表达式进行求解,得到地球自转角,并根据地球自转角测定世界时UT1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210333431.4A CN114911153B (zh) | 2022-03-31 | 2022-03-31 | 基于星-地超长基线干涉测量的世界时ut1测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210333431.4A CN114911153B (zh) | 2022-03-31 | 2022-03-31 | 基于星-地超长基线干涉测量的世界时ut1测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114911153A CN114911153A (zh) | 2022-08-16 |
CN114911153B true CN114911153B (zh) | 2023-02-21 |
Family
ID=82762551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210333431.4A Active CN114911153B (zh) | 2022-03-31 | 2022-03-31 | 基于星-地超长基线干涉测量的世界时ut1测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114911153B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118091718B (zh) * | 2024-04-17 | 2024-07-05 | 中国科学院国家授时中心 | 通过低轨卫星下行导航信号改善ut1解算精度的方法 |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006093436A1 (fr) * | 2005-03-04 | 2006-09-08 | Valeri Grigirievich Grachev | Systeme espace - temps fonctionnant en temps reel |
CN101915926A (zh) * | 2010-08-05 | 2010-12-15 | 中国科学院云南天文台 | 一种卫星激光测距中单向距离的确定方法 |
CN103078688A (zh) * | 2012-12-28 | 2013-05-01 | 中国电子科技集团公司第五十四研究所 | 干涉测量时射电星信号和航天器信号时延不一致的校准方法 |
CN105659869B (zh) * | 2011-09-27 | 2014-06-25 | 北京航天飞行控制中心 | 一种基于数传信号的甚长基线干涉测量处理方法 |
CN106643702A (zh) * | 2016-11-09 | 2017-05-10 | 中国科学院西安光学精密机械研究所 | 基于x射线的vlbi测量方法和***以及地面验证装置 |
RU2624912C1 (ru) * | 2016-06-14 | 2017-07-10 | ОО Международная академия наук экологии, безопасности человека и природы | Радиоинтерферометрический способ исследования объектов ближнего и дальнего космоса и система для его реализации |
CN107643527A (zh) * | 2017-09-08 | 2018-01-30 | 桂林电子科技大学 | Gps共视信号模拟生成方法及gps共视信号模拟器 |
CN107782317A (zh) * | 2017-09-20 | 2018-03-09 | 西安空间无线电技术研究所 | 一种倒vlbi深空单向导航定位方法 |
CN108919634A (zh) * | 2018-08-13 | 2018-11-30 | 中国科学院国家授时中心 | 一种北斗三频非差非组合观测值时间传递***及方法 |
CN109031198A (zh) * | 2018-07-30 | 2018-12-18 | 中国科学院国家授时中心 | 一种利用非合作卫星信号实现定位的方法 |
CN109298392A (zh) * | 2018-09-06 | 2019-02-01 | 中国人民解放军63921部队 | 一种完全空间基线vlbi***的稀疏测量方法 |
CN110018505A (zh) * | 2019-04-22 | 2019-07-16 | 中国科学院国家授时中心 | 一种利用通信卫星干涉测量实现舰船定位的方法 |
RU2702098C1 (ru) * | 2018-07-25 | 2019-10-04 | федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации | Способ определения параметров орбиты искусственного спутника земли с использованием приемных опорных реперных станций |
CN110554373A (zh) * | 2019-08-25 | 2019-12-10 | 中国科学院国家授时中心 | 干涉时间测量与测距方法 |
CN110865200A (zh) * | 2019-11-25 | 2020-03-06 | 北京无线电计量测试研究所 | 地球自转角速率测量方法 |
CN111308413A (zh) * | 2019-12-03 | 2020-06-19 | 中国人民解放军63921部队 | 一种大地测量vlbi射电源天区覆盖评价方法 |
CN112964229A (zh) * | 2021-02-08 | 2021-06-15 | 西安电子科技大学 | 目标天区覆盖性的星地联合观测确定方法 |
CN113108755A (zh) * | 2021-04-14 | 2021-07-13 | 南京中探海洋物联网有限公司 | 一种基于vlbi的高精度太空飞行器轨道测量方法 |
CN113589679A (zh) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | 一种基于干涉测量的卫星精密时间传递方法 |
CN113885056A (zh) * | 2021-09-07 | 2022-01-04 | 中国科学院上海天文台 | 星地双向时间同步设备时延的动力学标定***及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2784745T3 (es) * | 2010-03-19 | 2020-09-30 | Sky Perfect Jsat Corp | Sistema de posicionamiento para satélite artificial geoestacionario |
-
2022
- 2022-03-31 CN CN202210333431.4A patent/CN114911153B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006093436A1 (fr) * | 2005-03-04 | 2006-09-08 | Valeri Grigirievich Grachev | Systeme espace - temps fonctionnant en temps reel |
CN101915926A (zh) * | 2010-08-05 | 2010-12-15 | 中国科学院云南天文台 | 一种卫星激光测距中单向距离的确定方法 |
CN105659869B (zh) * | 2011-09-27 | 2014-06-25 | 北京航天飞行控制中心 | 一种基于数传信号的甚长基线干涉测量处理方法 |
CN103078688A (zh) * | 2012-12-28 | 2013-05-01 | 中国电子科技集团公司第五十四研究所 | 干涉测量时射电星信号和航天器信号时延不一致的校准方法 |
RU2624912C1 (ru) * | 2016-06-14 | 2017-07-10 | ОО Международная академия наук экологии, безопасности человека и природы | Радиоинтерферометрический способ исследования объектов ближнего и дальнего космоса и система для его реализации |
CN106643702A (zh) * | 2016-11-09 | 2017-05-10 | 中国科学院西安光学精密机械研究所 | 基于x射线的vlbi测量方法和***以及地面验证装置 |
CN107643527A (zh) * | 2017-09-08 | 2018-01-30 | 桂林电子科技大学 | Gps共视信号模拟生成方法及gps共视信号模拟器 |
CN107782317A (zh) * | 2017-09-20 | 2018-03-09 | 西安空间无线电技术研究所 | 一种倒vlbi深空单向导航定位方法 |
RU2702098C1 (ru) * | 2018-07-25 | 2019-10-04 | федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации | Способ определения параметров орбиты искусственного спутника земли с использованием приемных опорных реперных станций |
CN109031198A (zh) * | 2018-07-30 | 2018-12-18 | 中国科学院国家授时中心 | 一种利用非合作卫星信号实现定位的方法 |
CN108919634A (zh) * | 2018-08-13 | 2018-11-30 | 中国科学院国家授时中心 | 一种北斗三频非差非组合观测值时间传递***及方法 |
CN109298392A (zh) * | 2018-09-06 | 2019-02-01 | 中国人民解放军63921部队 | 一种完全空间基线vlbi***的稀疏测量方法 |
CN110018505A (zh) * | 2019-04-22 | 2019-07-16 | 中国科学院国家授时中心 | 一种利用通信卫星干涉测量实现舰船定位的方法 |
CN110554373A (zh) * | 2019-08-25 | 2019-12-10 | 中国科学院国家授时中心 | 干涉时间测量与测距方法 |
CN110865200A (zh) * | 2019-11-25 | 2020-03-06 | 北京无线电计量测试研究所 | 地球自转角速率测量方法 |
CN111308413A (zh) * | 2019-12-03 | 2020-06-19 | 中国人民解放军63921部队 | 一种大地测量vlbi射电源天区覆盖评价方法 |
CN112964229A (zh) * | 2021-02-08 | 2021-06-15 | 西安电子科技大学 | 目标天区覆盖性的星地联合观测确定方法 |
CN113108755A (zh) * | 2021-04-14 | 2021-07-13 | 南京中探海洋物联网有限公司 | 一种基于vlbi的高精度太空飞行器轨道测量方法 |
CN113589679A (zh) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | 一种基于干涉测量的卫星精密时间传递方法 |
CN113885056A (zh) * | 2021-09-07 | 2022-01-04 | 中国科学院上海天文台 | 星地双向时间同步设备时延的动力学标定***及方法 |
Non-Patent Citations (3)
Title |
---|
Geostationary orbit determination using SATRE;Lei Hui etc;《Advances in Space Research》;20110508;全文 * |
The NTSC VLBI System and its application in UT1 measurement;Dang Yao etc;《Research in A stronomy and Astrophysics》;20201231;第20卷(第6期);全文 * |
基于13米宽带VLBI***的UTI自主测定研究;姚当;《中国博士学位论文全文数据库 基础科学辑》;20220215;1-100 * |
Also Published As
Publication number | Publication date |
---|---|
CN114911153A (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schuh et al. | Very long baseline interferometry for geodesy and astrometry | |
Schuh et al. | VLBI: A fascinating technique for geodesy and astrometry | |
Huang et al. | Improvement of orbit determination for geostationary satellites with VLBI tracking | |
Petrachenko et al. | VLBI2010: next generation VLBI system for geodesy and astrometry | |
Delva et al. | Chronometric geodesy: methods and applications | |
Bi et al. | An introduction to China FY3 radio occultation mission and its measurement simulation | |
Schreiner et al. | Quality assessment of COSMIC/FORMOSAT-3 GPS radio occultation data derived from single-and double-difference atmospheric excess phase processing | |
EP1862810A1 (en) | Space-time real-time system | |
Anderson et al. | Simulations of VLBI observations of a geodetic satellite providing co-location in space | |
Ma | Very long baseline interferometry applied to polar motion, relativity and geodesy | |
CN114911153B (zh) | 基于星-地超长基线干涉测量的世界时ut1测定方法 | |
Wang | Integrated processing of GNSS and VLBI on the observation level | |
Resch et al. | The media calibration system for Cassini radio science: Part II | |
CN114035418B (zh) | 基于脉冲星的共视时间比对方法 | |
Wu et al. | Simulation experiments on high-precision VGOS time transfer for future geopotential difference determination | |
Heinkelmann | Vlbi geodesy: Observations, analysis and results | |
Marz et al. | Geosynchronous satellites expanding a future GNSS satellite constellation: A precise orbit determination study | |
MacMillan | Comparison of EOP and scale parameters estimated from the three simultaneous CONT17 VLBI observing networks | |
Koyama | The use of very long baseline interferometry for time and frequency metrology | |
Montenbruck et al. | GNSS orbit determination and time synchronization | |
Ryan et al. | Crustal dynamics project data analysis fixed station VLBI geodetic results | |
Yao et al. | The NTSC VLBI System and its application in UT1 measurement | |
Renzetti et al. | The Deep Space Network. An instrument for radio navigation of deep space probes | |
Kareinen | Observation strategies for current and future geodetic very long baseline interferometry | |
CN118151189B (zh) | 基于实时信息的电离层掩星信号测试方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |