基于二维三次卷积插值算法的辛几何射线追踪  

Symplectic ray-tracing based upon the bi-cubic convolution interpolation method

作　　者：李川[1] 王有学[1,2] 何晓玲[1] 刘荣平[1] 贠鹏[1] 熊彬[1] 许继峰[1,3] 

机构地区：[1]桂林理工大学地球科学学院,桂林541004 [2]桂林理工大学广西矿冶与环境科学实验中心,桂林541004 [3]中国科学院广州地球化学研究所,广州510640

出　　处：《地球物理学报》2014年第4期1235-1240,共6页Chinese Journal of Geophysics

基　　金：国家自然科学基金项目(41174077;41164004);广西壮族自治区"八桂学者"(2013);广西壮族自治区桂林市"漓江学者"(2013)专项经费联合资助

摘　　要：射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法(SAM-Symplectic Algorithm Method)及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.Ray tracing is the basic problem in seismic imaging, and the reliability of the imaging depends on the accuracies of both spatial trajectory and traveltime of the ray. Based upon the properties of Hamilton System that can keep the stability of the solution for the eikonal equation, the Sympletic algorithm method （SAM） and bi-cubic convolution algorithm are used in seismic ray tracing in this study. Due to the use of the Sympletic algorithm, the reliable seismic wavefront produces an accurate ray trajectory. Finally, the results from numerical modeling show that the use of SAM can not only keep the stability of the Hamilton System with a fast computation but also improve the accuracy of the seismic ray tracing.

关 键 词：射线追踪 辛几何算法 二维三次卷积插值 数值模拟 

分 类 号：P631[天文地球—地质矿产勘探]