摘要：本文在传统傅里叶微分矩阵的基础上,对原始微分算子进行改进,引入了错格微分算子的傅里叶伪谱方法.尽管该方法增加了一些计算量,但却极大地提高了计算精度和稳定性.而且,该方法将微分计算过程由传统的傅里叶变换转换为一般的矩阵矢量乘积,大大降低了微分求解过程的复杂程度.在均匀介质中,将错格伪谱微分算子计算的结果和解析解进行比较,结果表明本文算子几乎达到了解析解的精度.而在分层均匀介质中的实验结果同时显示,该方法精度高、稳定性好,是一种研究层状介质中地震波传播的有效数值方法.

关键词：伪谱微分算子 波场模拟 错格傅里叶伪谱方法 

单位：中国科学院地质与地球物理研究所 北京100029 中国科学院研究生院 北京100049