用Galerkin方法研究地球本征模.pdfVIP

第27卷第1期 天文学进展 Vbl．27．No．1 2009年3月 PROGRESSINASTRONOMYMar．．2009 文章编号：100m8349(2009)01-008m09 张 冕1，一，黄乘利1 (1．中国科学院上海天文台，上海200030；2．中国科学院研究生院，北京100039) 摘要t地球本征模在地球内部物理结构和动力学机制研究中一直是一个重要的课题，与地球物理 学、地震学联系紧密。过去采用的方法都是由一组运动方程进行积分，然后在边界上判断是否符合 各物理量的约束条件，从而得到地球本征模．提出了一个新的解算此问题的方法，即用Galerkin 方法计算本征模和本征函数，相应的边界条件用rIhu方法处理。结果表明：在研究地球自由振荡 时，Galerkin方法是一种行之有效的方法；随着所采用的地球模型的逐步完善，计算的结果也与 观测更加符合．对于三层地球模型，考虑地球自转的结果比不考虑自转的有所改善，但其改善效果 不如仅将地幔更细致分层对结果的改善显著． 关键词：天文地球动力学；地球本征模；自由振荡； Galerkin方法；rIau方法；边值问题 中图分类号tP542 文献标识码：A —1上．—j一 日fJ 吾 地球在经历地震时也会产生振动。这种振动只与地球本身的内部属性有关，是地球的各 种本征模的叠加。通常人们将周期短于1h的本征模称为自由振荡，另外还有周日频段的周 究地球的本征模可以让人类更多地了解地壳及地球内部物质密度分布及流变性质，是地震学 与大时间尺度地球动力学之间联系的窗口。 荡周期是60min，这与今天实际观测的54 2008_0410 收稿日期：2008-03-14；修回日期t 市科委资助课题(06Dz22101，06zRl4165) 1期 张冕，等人。用Galerkin方法研究地球本征模 81 利用了一个两层的地球模型，假设地球由一个液核和地幔组成，计算得到环形场最长的周期 为42．5 后Alterman等人【4】系统研究了本征模的计算问题，把二阶偏微分方程组转换为一阶常微分 方程组，奠定了现代计算机数值求解的基础。 地球模型的简正模的数值方法。方明(1991)【6】在他的博士论文中更详细研究了地球自由振荡 问题。但自20世纪90年代后，地球自由振荡的研究就渐渐趋少了． 现今应用的方法，都是建立在Alterman等人【4】方法的基础上：首先把二阶的偏微分方 程组转换成一阶的常微分方程组，然后给定初值，对常微分方程组进行积分，然后在边界上 (从一层到十二层、从不考虑自转到考虑自转)的本征模进行了计算，并对结果进行了讨论。 2用Galerkin方法研究地球的自由振荡 2．1 基本方程和边界条件 以流体静平衡态作为初始应力状态，在考虑自转、自引力及预应力等条件下，弹性地球的 运动方程为 (1) 其中，u是位移场，Q是地球自转角速度矢量，Ⅵ是由于质点弹性位移造成质量重新分布 产生的附加势。90，加分别是流体静平衡态时地球的重力加速度和密度．对各向同性的弹性 体，其中的应力T为 T=入(V·u)J+p【Vu+(Vu)T】 (2) 其中，A、“为拉梅参量。泊松方程为 V2Ⅵ=47rGV·(肋t正) (3) 由于扁率的影响只有1／300，与其他量相比是个微小量，所以本文研究的是一个简化 的球对称模型，所有讨论均在球坐标中进行，把球面上的矢量场较方便地表示成球形场和 环形场。 将球面上的标量场Ⅵ表示为 Ⅵ=∑∑Ⅵ。。(r)可。。(p，咖) (4) n=O