利用谱元法计算牦荦犚犈犐地球的表面负荷变形-地球物理学报.pdf

第 卷 第 期 地 球 物 理 学 报 ， ６２ ７ Ｖｏｌ．６２ Ｎｏ．７ 　 　 　 　 　 年 月 ， ２０１９ ７ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｊｕｌ．２０１９ 　 　 　 廖彬彬，徐建桥，孙和平等 利用谱元法计算 地球的表面负荷变形 地球物理学报， （）： ， ： ．２０１９． ＳＮＲＥＩ ． ６２７ ２３８２２３９３ｄｏｉ１０． ／ ６０３８ｃ２０１９Ｍ００４３． ｊｇ ， ， ， ＬｉａｏＢＢ ＸｕＪＱ ＳｕｎＨ Ｐ ｅｔａｌ．２０１９．ＳｅｃｔｒａｌｅｌｅｍｅｎｔｍｅｔｈｏｄｔｏｌｏａｄｄｅｆｏｒｍａｔｉｏｎｉｎａＳＮＲＥＩＥａｒｔｈ．犆犺犻狀犲狊犲犑． ｐ （ ）， （）： ， ： ／ 犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ６２７ ２３８２２３９３ｄｏｉ１０．６０３８ｃ２０１９Ｍ００４３． 狆狔 ｊｇ 利用谱元法计算犛犖犚犈犐地球的表面负荷变形 ， ， １２ １ １２ １ 廖彬彬 ，徐建桥 ，孙和平 ，周江存 １中国科学院测量与地球物理研究所，大地测量与地球动力学国家重点实验室，武汉 ４３００７７ ２中国科学院大学，北京 １０００４９ 　 摘要 谱元法作为一种基于变分原理的数值方法，已经成为地球动力学数值模拟的重要手段 本文 旨在对 ． 　 （ ）正演模拟表面负荷问题的谱 有限元方法进行改进 简单地介绍了谱元法计算负荷边值问题的相关 Ｍａｒｔｉｎｅｃ２０００  ． 理论，并推导了中性分层液体中的平衡方程以及其对应的弱形式解 在地球半径方向上，采取高阶样条函数作为负 ． 荷解的基函数，并利用高斯 勒让德方法进行数值积分 均质地球模型数值试验表明，相较于 的方法，本文  ． Ｍａｒｔｉｎｅｃ 方法可以有效地加快数值解的收敛速度 最后，利用 模型计算了 地球的负荷勒夫数 数值计算结果 ． ＰＲＥＭ ＳＮＲＥＩ ． 表明本文方法的收敛性和稳定性都非常好 本文方法获得的负荷勒夫数与采用传统龙格 库塔积分方法获得的结 ．  果相对误差在 ０．０１％的量级． 关键词 谱元法；弱解；中性分层流体；