长安大学设计综述.doc

课 程 设 计 题 目: 有限频射线层析成像 专 业：地球物理学 姓 名： 学 号： 指导教师：白超英 教授 完成时间：2012年03月01日 有限频射线层析成像 1.引言 与探索太空一样，长久以来，了解地下世界的面貌也是人类的美好愿景，这也是地震学研究的主要任务之一。上天容易入地难！一直以来，人们主要通过地震资料，特别是天然地震资料来认识地球深部构造。地震层析成像方法广泛应用于探测地壳和地幔结构，其发展相对较早，技术方法也比较成熟。我们所要介绍的有限频射线层析成像方法就是其中的一种。 众所周知，在地震层析成像中，基于射线理论的正演解决方案具有局限性。为正确地表达地震层析成像中的正演问题，需要应用三维地震波理论。但就地震研究的规模而言，目前为止，这种方法仍是不可行的。使用地震波菲涅尔体的概念，产生了一套更加完整、更符合实际情况、更加精确的正演解决方案。一种在地震层析成像中使用有限频射线的方法，称为有限频射线层析成像。有限频理论更加符合实际情况，有限频射线层析成像产生的反演结果一定程度上优于射线层析成像产生的反演结果。 2.现有方法的缺点 自上世纪年代以来，地震层析成像主要是根据射线理论，假设地震波为无限高频（或无限频宽），按照费马原理，即地震波沿着耗时最少的路径传播，在不同台站上观测到的到时的差异完全取决于射线路径上的速度结构。然而，实际地震波并不是无限频宽，其频率成分是有限的，影响地震波传播的主要构造也不限于中心射线，而是集中于射线邻近区域，射线理论只是无限高频条件下的一种近似。该理论只适用于地球内部速度变化平缓且速度异常体变化的尺度远大于地震波本身波长的情形，这在很大程度上限制了地震层析成像对细微速度结构的解析能力。 世纪年代就已经提出了不需要高频近似的波动方程层析成像方法，在理论上比走时层析成像具有更高的反演分辨率。由于波动方程层析成像计算效率非常低，其反演的目标函数和速度摄动之间表现为强烈的非线性关系，对初始模型要求很高，再加上地震子波反演困难、地震信号的信噪比较低、实际地震波射线追踪的方法种类较多。传统方法有基于初值问题的打靶法（Shooting method, i.e., Julian and Gubbins，1977和基于边值问题的弯曲法（Bending method, i.e., Um and Thurber, 1987），但是他们不能处理介质中较强的速度变化，有时无法求出全局最小走时，计算效率较低，阴影区内无射线。随着射线追踪方法的发展，出现了直接从Huygens原理或Fermat原理出发，采用等价波前来描述地震波场特征的方法。Vidale（1988，1990）和Podvin（1990）等人则从程函方程出发，首先求出走时场分布，再用计算走时场的最速下降方向的办法，得到每一条接收点到震源的射线路径。随后，Qin（1992）等人对Vidale的方法作了改进，提出了波前扩展方法。黄联捷、李幼铭、吴如山(1992)基于Huygens 原理提出WFRT 方法，根据计算精度的要求，逐次细化网格，进而求得最小走时。Sava 和Fomel（1998）提出了HWT（huygencs wavefront tracing）法。Moser（1991）提出了根据费马原理及图论的最短路径法。 地震层析成像反演方法可以分为两类：第一类是基于算子的线性或拟线性反演方法，又称为“确定性反演方法”；另一类是基于模型的完全非线性反演方法，又称为“随机反演方法”。 ①线性反演方法，主要包括代数重建法（ART）、同时迭代重建法（SIRT）、及其改进型；②拟线性反演方法，主要包括高斯-牛顿法（Gauss-Newton method or Quasi-Newton method）、阻尼最小二乘法（Damped-least squares）、最速下降法（Steepest descent method）、共扼梯度法（Conjugate gradients）、亚空间法（Subspace method）、及其改进型；③完全非线性反演方法，主要包括蒙特卡罗类方法（Monto Carlo-type method），如：遗传算法（Genetic algorithm）、神经网络法（Neural network method）、模拟退火法（Simulated annealing）、及其改进型。上述反演算法又可与空间平滑（smoothing）、滤波（filtering）、约束（constrained）、规则化（regulation）等组合形成不同的反演方法，有关局域解的反演算法及其相互关系，相关文献有Zhao and Kayal，2000；Rawlinson and Sambirdge，2003；Aster et al，2005；