远震波场正演模拟方法及应用

地震波场数值模拟是壳幔结构成像和深部探测的重要基础.经典的远震波场数值模拟主要基于一维地球模型,包括解析法、半解析法和数值法等.这些算法能够高效地计算理论地震图,但由于将地球假设为一维层状介质,难以考虑介质的横向非均匀性对地震波场的影响.近年来随着计算机性能的不断提高,三维地震波场数值模拟方法得到快速发展,并被广泛用于局部/区域地震波场模拟及壳幔结构成像.然而由于计算成本较高,实现全球尺度模型的高频地震波场数值模拟仍存在较大挑战,因而基于远震波场的混合数值模拟方法逐渐得到关注和应用.远震波场混合数值模拟法将全球模型分解为全球区域和局部区域,在全球区域中采用一维地球模型近似,应用快速算法计算一维模型中的全球高频理论地震图;在局部目标区域内采用三维数值方法(谱元法、有限差分法等)和注入技术,模拟地震波在三维非均匀介质中的传播,从而在波场计算的效率和精度之间达到平衡.随着密集台阵观测的普及,对地下结构成像的分辨率提出了更高的要求,准确高效的远震波场混合模拟方法在高分辨率地震成像领域将发挥日益重要的作用.本文系统地总结了远震波场数值模拟的一维模拟方法,并在此基础上重点介绍远震波场混合数值模拟方法的原理及应用.