数值微分法在潮汐资料的Nakai处理方法中的应用.pdfVIP

维普资讯 第 9卷第4期 地 壳 形 变 与 地 震 Vo1．9．]qo．4 1989年工1月 CRUSTAL DEFORMATION ^NDEARTHQUAKE N 0v． ， 1989 数值微分法在潮汐资料的Nakai处理 方 法 中 的 应 用 张 雁 滨 (国家地震局地 震研究所 ) 摘 要 本 文在潮汐资料~Nakal处理方法计算 中使用了数值微分法 申的插值微分公式 ，加上 插值 点的治 当选取，使计算量及计算时 阏大幅度 的减少，并使Nakai法能在Pc—l5O0机上进行，为 NakM 法的广泛应用提供 了方便 。 一 、 弓『 言 目前，在对潮汐资料进行调和分析时，包括国际固体潮中心在内的研究人员通常都惯用 Nakai法对潮汐资料进行预处理，Nakai处理方法在对潮汐资料的评价与分析处理 中有着 重 要的意义 、‘，由于其中要计算潮{夕理论值对时间的导数矗 ()，若用封 闭公式来进行复 台 求 导而推导 (z)的表达式则太繁，所 以传统的方法是在程序里采用分波叠加的办法来计 算 。 这样却增大了计算量， 也无法在 PC一1500这样的机子上进行， 而 目前大多数 台站只 有 Pc一 1500机。即使在 IBM机上也需较长的计算时间，因而使得 Nakai法的应用受到限制 。针 对 以上情况，1987年底在 《地应变连续观测程序集 》中推出了采用数值微分法中的插值微分公 式来计算 R ()的Nakai法，加上插值点的巧妙选取，在精度得到完全保证的基础上，计算 量及计算时间都大大减少 (在不加任何协处理器的 IBM286型机 上，处理一个月的应变资料 不到 8分钟 )，在 PC一1500机上也能很方便地进行，从而使台站的同志能简单、方 便地应 用Nakai法进行潮汐资料处理 。 二 、Nakai法简介 [1] 对某一 时刻 的潮汐观洲值可写成： y l： f+D (1) 本文1988年3且1j日艘到． I)奎相束．张雁滨，N h i法在处理台站删 蹙料中的用途-形变台站通讯，lg84年增刊 维普资讯 80 地 壳 形 变 与 地 震 9卷 一 l _ 一 一 一 一 一 一 一 — 一 一 ——= 一 其中T 是潮汐部分，D。是漂移部分。Nakai采用多项式扭合法对潮汐部分用潮汐理论值进 行拟台，对漂移部分用时间的多项式来近似 ，并只取到二次项 ，于是 (1)式为 y = ． (￡)+卢 + 。+ I￡+ 。 (2) 其中 ()是潮汐理论值， 是比例因子，它与潮汐因子以及格值有关， =，／8．(f是平均 潮汐因子， 是灵敏度 )。卢=一 ·6t，而 是平均的时间滞后，它和相位滞后有关 。 。 是常数项 ， 反映线性漂移 ， ：反 映二次漂移 (非线性 )。 取 个时刻的观测值为一组，当计算出各时刻的潮汐理论值 ()及理论值的导数 () 时，就可 由 (2)式得到 个方程。用最小二乘法解这 个方程就可得到 (2)式 中 的 5 个未知系数。若把 5个系数再代 回 (2)式，则得到各 t时刻的拟台值 l， 。用它可对观测 值进行补缺，并且还能得到反映拟台程度的均方误差； ： - _k-