重力潮汐观测数据分析方法的比较与研究.pdfVIP

重力潮汐观测数据分析方法的比较和研究 陈晓东孙和平 中国科学院测量与地球物理研究所动力大地测量学重点实验窒，武汉430077 摘要本文对目前国际上广泛使用的两种重力潮汐观测数据分析方法(Etema和Baytap—G) 进行了初步的比较和研究。文章首先讨论了相直的数学模型，然后利用武汉九峰站GWR．C032 081255 型超导重力仪的重力潮汐观测进行分析，比较结果说明，在频率从0．721499到3cpd的 范围内，所估算的主波振幅因子与理论值的差异在同一量级，经海潮负荷改正后，这一差异均小 于1％。 关键词重力潮汐观测分析方法结果比较 1 引言 无论是海潮或者是固体潮的研究之中，一个首要的任务就是把观测数据进行分析，将它 们分解成为不同周期的潮波分量，这种处理方法称为调和分析。近年来，随着计算机技术的 快速发展，利用不同的最小二乘滤波器分离干扰噪声发展了几种实用的潮汐分析方法…。 在对重力固体潮观测数据的分析中，目前国际上广泛应用的调和分析方法有二种：(1)由德 被国际地潮中心列为标准的潮汐分析方法，目前最新的版本是Etema34；(2)由日本的 Tamura研制的Baytap·G潮汐分析方法”1。 随着新型重力仪的出现，特别是超导重力仪的出现，使得重力固体潮的观测精度已达 109或更高量级”3。这样就导致了这样一个问题产生，即对于同样的观测数据，分析结果的 不同可能来源于所用分析方法的不同。本文采用了Tamura的引潮位展开表”“1．利用武汉 析方法进行了初步的比较和研究。 2数学模型和相应软件的功能 地球表面上任意一点的重力潮汐观测序列y。。(￡)不仅包含重力潮汐信号Ytlde(t)，也 含有其它信号，可以用公式(1)表示 自然科学基金委项目联合资助。 364 ￡ y。h(f)=y‰(‘)十三aIR(‘)+d(‘)+s(￡) 式中，^(t)为附加信号，为气象或水文观测序列，或者是根据IERS极坐标观测数据计算的 极潮和日长变化时间序列，a。为序列^(￡)的重力导纳值；d(￡)为仪器漂移；其它观测误 差包含在s(t)中。 重力潮汐调和分析方法对观测序列y。(t)进行修正，去掉仪器漂移项和附加信号的影 响，再利用台站坐标和引潮位展开表计算理论潮汐，通过比较理论和观测重力潮汐信号，估 计出一组台站的潮汐参数(振幅因子和相位滞后)。但是，分析方法并非能估计出引潮位展 开表中给出的每一个潮波分量的潮汐参数，根据瑞利原理，引潮位展开表中的潮波分量迭加 在一起组成波群，并假设同一波群中潮波分量的振幅因子和相位滞后相等。如此，分析方法 就可以只对每个波群的潮汐参数进行估算。 下面我们给出Etema和Baytap-G方法的数学模型和相应软件的功能，并且讨论了两个 方法中附加信号对重力潮汐观测值的影响、仪器漂移和计算结果精度的计算方法。 2．1 Eterna方法的数学模型和功能 Eterna方法采用最小二乘平差方法来估算潮汐参数、附加信号的回归系数和用来拟合 仪器漂移的切比雪夫(Tschebyscheff)多项式的系数，所用最小二乘平差的数学模型如下 Ⅳ (2) y。。(1)+r(t)=∑(丑·c。+t·s。)+∑西。·L+∑屯·Zk ^=1 m ★ 参数峨(振幅因子)和△妒。(相位滞后)的线性形式，表示为 X．=H．’cosh妒。、 (3) ■=以·sinA‘v．j C。、S。表示第n个波群的理论值，是该波群所有潮波分量的累加和，从第ol。个波开始到第 芦．个波结束，设第i潮波分量的理论振幅为A。，相位为吼，频率为l，用公式表示为