应变固体潮Matlab,应变固体潮理论值计算及其调和分析

第 卷 第 期年 月 地 球 物 理 学 报 , , 应变固体潮理论值计算及其调和分析 刘 序 俨 李 平国家地震局地震研究所 , 武汉 摘 要 本文给出了计算应变固体潮理论值的“ 封闭公式 , ,和调和分析公式 采 用 了 · 引潮位展开 , 推导了应变固体 潮的大地系数和相位改正 长周期波 、 全 日波 、 半 日波和 日波 , 可供调和分析之 用 并通过 调和分析计算 , 对理论值公式进行了验证 文 中还列 出了三个台站的首批应变固体潮观测资料的维 氏调和分析结果 引 言 近年来 , 我 国开始 了应变固体潮 的观测 和研 究工作 , 云南 、 湖北和 江苏等省有的 台站已 安装 了我国 自制的 一 型石英伸缩仪 , 河南等地也开展了体积应变仪的研 制 和观测 , 均取得了相当数量的观测 资料 , 也都分 别作了调和分析 鉴于国内还是首次进行此项 工作 , 本文将给 出有关应变固休潮理论值的最 终公式和应变固体 潮的长周期波 、 全 日波 、半 日波和 日波的全套大地系数和引角改正 二 、 计算应变 固体潮的封闭公式 根据固体 潮理论 , 所有物理过程都可 以通过 日月引潮位来表录 假定地球 内部的物 理性质具有球对称性 , 则在球坐标系 , 日, 劝 中地球内部任一点的位移分量为 〔 一 互 · 。 ,令, 口 。日 。口几 “一口一口 一亡 应变分量 为〔 口“ 。 。门门 一 下一 二砚万一 十 万 夕 口 尸 本文 年 月 日收到 , 年 月 日收到修改稿 期 刘序俨等 应变固体潮理论值计算及其调和分析 ⋯, 。盛 一卫 一 鲤只于 月 自 , 竺 日 塑 呱一一一口几 口“ 。 。 “ , 一一二 , , 万 二丁一 一 口 一 创 兄 撇一户 一一瓮 ‘ 。一合器 口 , 口鱼旦口夕 将 式代人 式 , 可得地表面的应变分量 为 、、 矗颤 · 器十 ‘· 小 众令 一 ‘· 、、蓄十 ‘·协。 艺︺︷ 一 兰 食 一些一 旦丛仁蕊 甲 甲 又 几 二 三生 三丝 、又 ︷一甲之甲甲,人之心 , 山 式 中 几, 甲 , 和 分别为观测点的地心经纬度 , 地球向径和重力加速度对于实际应用 , 月亮取至三阶位 , 太阳取至二阶位 , 且顾及到 〔, , 一 普 。 贡 ’ , 。。 ·’一 ‘, , — 行 欠 仁、 ‘ 二 一 , 甲 占 甲 占 , 式中 ‘ 为天体中心至地心的平均距离 , 为其瞬间距离 , 为观测点的地球向径 , 司 为 。 系数 , 宕 为天 沐的天顶距 , 占 为天体的黄纬 , 为天体的时角 将 式代人 式 经推导后可求得计算地表面的 应变固体潮分量的最终公式一 告 。 , 责’ 一 , 。 。一 · , 一 。一 、 一 。 万 〕 器 。 , 令’ 一 , , 器 ·, 令责‘ 一 「, 一 , ‘· “、 。 孟 。 , ‘ , , , 、 勺 , 口 门 ‘ 、 , , ‘ , ‘ 户 、 尸 了 ‘ 、 尸 , , 、 一 甲 月 一 又 ‘ 一 少 十 丁丁一 。 又尸 夕 一 屯不 气, ‘ 石 石 一 ‘ 。‘ 不 尸 亡 入 , 。 , 、 。 —行 夕夕 , 入 , 占 一 , 占 , ’ 一 , ‘ 三生 , , 、巴 去孟 —行 少 呼 , , 、 亡 、 。 , 【 , 】 尺 , , 一 二丫夕 , 占 ’ 一 , 占 〕 , , 二丫。。 二 。 。。 。 。 一 。、。 。。 。。。。 。 里丛 一 。 , 、 ‘ 令’ 一 , , 器 。 , 食‘ 二 二 。。。 , 。