摘要
笛卡尔坐标系中的经典程函方程在静校正、叠前偏移、走时反演、地震定位、层析成像等很多地球物理工作中都有应用,然而用其计算起伏地表的地震波走时却比较困难.本文通过把曲线坐标系中的矩形网格映射到笛卡尔坐标系的贴体网格,推导出曲线坐标中的程函方程,而后,用Lax-Friedrichs快速扫描算法求解曲线坐标系的程函方程.研究表明本文方法能有效处理地表起伏的情况,得到准确稳定的计算结果.由于地表起伏,导致与之拟合的贴体网格在空间上的展布呈各向异性,且这种各向异性的强弱对坐标变换法求解地震初至波的走时具有重要影响.本文研究表明,随着贴体网格的各向异性增强,用坐标变换法求解地表起伏区域的走时计算误差增大,且计算效率降低,这在实际应用具有指导意义.
The classical eikonal equation is commonly used in Cartesian coordinate system for problems that involve static correction,prestack migration,earthquake location and seismic tomography,but is less effective for calculating travel times in an earth model that has an irregular surface.We have presented a topography-dependent eikonal equation in a curvilinear coordinate system that makes use of the surface-fitting grid and map a rectangular grid onto a curved grid.Then,we utilized the efficient Lax-Friedrichs sweeping scheme to approximate the viscosity solutions of the eikonal equation in the curvilinear coordinate system.In this paper,we investigate the impacts due to the anisotropic stretching of the surface-fitting grid on the traveltime computation by using the topography-dependent eikonal equation,which has a significant meaning in the direction of our method in geophysical application.
出处
《地球物理学报》
SCIE
EI
CAS
CSCD
北大核心
2012年第10期3355-3369,共15页
Chinese Journal of Geophysics
基金
国家自然科学基金(41274070,41074033,40874041,41021063)资助
