利用谱元法计算SNREI地球的表面负荷变形

• 作者: 周江存 孙和平 廖彬彬 徐建桥
• 来源期刊: 地球物理学报
• 学科: 地球科学
• 关键词: 谱元法 弱解 中性分层流体 SNREI地球 负荷勒夫数
• 年，卷（期）: 2019,（7）
• 所属期刊栏目: 空间与大气物理学·重力与大地测量学
• 页码范围: 2382-2393
• 页数: 12页
• 分类号: P312
• 字数: 8259字
• 语种: 中文
• DOI: 10.6038/cjg2019M0043

摘要

谱元法作为一种基于变分原理的数值方法,已经成为地球动力学数值模拟的重要手段.本文旨在对Martinec(2000)正演模拟表面负荷问题的谱-有限元方法进行改进.简单地介绍了谱元法计算负荷边值问题的相关理论,并推导了中性分层液体中的平衡方程以及其对应的弱形式解.在地球半径方向上,采取高阶样条函数作为负荷解的基函数,并利用高斯-勒让德方法进行数值积分.均质地球模型数值试验表明,相较于Martinec的方法,本文方法可以有效地加快数值解的收敛速度.最后,利用PREM模型计算了SNREI地球的负荷勒夫数.数值计算结果表明本文方法的收敛性和稳定性都非常好.本文方法获得的负荷勒夫数与采用传统龙格-库塔积分方法获得的结果相对误差在0.01％的量级.

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	孙和平	中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室	92	1046	19.0	28.0
5	周江存	中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室	30	267	10.0	16.0
6	徐建桥	中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室	57	487	14.0	20.0
7	廖彬彬	中国科学院测量与地球物理研究所大地测量与地球动力学国家重点实验室	1	2	1.0	1.0