用Galerkin方法研究地球本征模  

Study of the Earth's Eigen Mode by Galerkin Method

作　　者：张冕[1,2] 黄乘利[1] 

机构地区：[1]中国科学院上海天文台,上海200030 [2]中国科学院研究生院,北京100039

出　　处：《天文学进展》2009年第1期80-88,共9页Progress In Astronomy

基　　金：国家自然科学基金资助项目(10773025,10633030);中科院方向性项目(KJCX2-YW-T13);上海市科委资助课题(06DZ22101,06ZR14165)

摘　　要：地球本征模在地球内部物理结构和动力学机制研究中一直是一个重要的课题,与地球物理学、地震学联系紧密。过去采用的方法都是由一组运动方程进行积分,然后在边界上判断是否符合各物理量的约束条件,从而得到地球本征模。提出了一个新的解算此问题的方法,即用Galerkin方法计算本征模和本征函数,相应的边界条件用Tau方法处理。结果表明:在研究地球自由振荡时,Galerkin方法是一种行之有效的方法;随着所采用的地球模型的逐步完善,计算的结果也与观测更加符合。对于三层地球模型,考虑地球自转的结果比不考虑自转的有所改善,但其改善效果不如仅将地幔更细致分层对结果的改善显著。The study of the Earth＇s eigen mode is an important way to study physical properties and dynamics in the Earth＇s interior, and is closely related with geophysics and seismology closely. The Earth＇s eigen mode is usually solved by classical numerical integration of partial differential equations of motion and boundary conditions. Galerkin method is used in this paper to calculate eigen-periods and eigen-functions of free oscillations of different earth models, while boundary conditions are treated with Tau method. Results show tboldsymbol Galerkin method is an alternative tool for such study. At first Earth models are assumed to be spherical non-rotating elastic isotropic （SNREI）. The eigen modes for the free inner core is calculated separately with respect to the other two modes： the 3-layer earth mode including solid inner core, fluid outer core and solid mantle （plus crust）, and the 12-layer earth mode in which the mantle （＋crust） is separated to 10 layers. Furthermore, the eigen modes of rotating free inner core and 3-layer earth model are calculated respectively. The results are improved and approach more closely to the observation with more realized earth model. For 3-layer Earth, incorporating earth rotation will improve the result of the eigen periods comparing with non-rotating earth model, for example, the difference of the eigen period of 0S2 with observed one is improved from 3.6% to 3.1%; however, the improvement by separating the mantle to 10 layers is more remarkably than tboldsymbol by introducing earth rotation： changing the earth model from 3 layers to 12 layers （both non-rotating）, the difference of the eigen periods with observed ones arc reduced significantly from 4.3% （maximum, or 3.6% for 0S2） to less than 1.6% （maximum, or 0.78% for 0S2）.

关 键 词：天文地球动力学 地球本征模 自由振荡 GALERKIN方法 Tau方法 边值问题 

分 类 号：P542[天文地球—构造地质学]