求解地震波方程的一类二阶优化辛格式  

作　　者：刘少林[1] 李小凡[1] 汪文帅[1,2] 刘有山[1] 张美根[1] 张欢[1] 

机构地区：[1]中国科学院地质与地球物理研究所,中国科学院地球深部重点实验室,北京100029 [2]宁夏大学数学与计算机科学学院,银川750021

出　　处：《中国科学：地球科学》2014年第2期283-291,共9页Scientia Sinica(Terrae)

基　　金：国家自然科学基金项目(批准号:41174047,40874024,41204041)资助

摘　　要：通过将地震波位移运动方程变换至Hamilton体系,引入广义动量和广义坐标,并定义系统动能和势能的Lie算子,构造了一类新的适用于高效长时间地震波模拟的二阶辛格式,同时将此格式用Baker-Campbell-Hausdorff(BCH)公式展开,基于截断误差原理极小得到了三组优化系数.在与常见辛格式对比中,从理论上和数值实验中分析了本文构造的这类优化二阶格式高精度高效率性;在与经典Newmark格式对比中,从长时间计算角度论证了本文格式具有长时间地震波计算能力;在非均匀介质地震波模拟中,本文格式与三阶辛格式得到了一致的地面合成地震记录和单道地震记录.

关 键 词：弹性波 辛格式 截断误差 地震波场 

分 类 号：P315.31[天文地球—地震学] O241.8[天文地球—固体地球物理学]