基于程函方程快速行进法地震走时成像方法研究-地球探测与信息技术专业论文.docxVIP

Seismic Traveltime Tomography by Using a Fast Marching Method Based on Solving Eikonal Equation A Dissertation Submitted for the Degree of Master Candidate：Li Qiang Supervisor：Prof. Bai Chaoying Chang?an University, Xi?an, China I 摘要 摘 要 地震层析成像技术已成为研究地球内部的结构及各向异性最为行之有效的方法之 一，其中走时层析成像技术是一种研究地下速度结构的方法。一般来讲，地震射线追踪 的正演是走时层析成像技术中非常之重要的一个环节。地震射线追踪的方法有很多，例 如基于射线理论的打靶法、弯曲法(伪弯曲法)、高斯射线束算法等；基于网格单元扩展 的有限差分解程函方程法、最短路径算法; 以及结合射线和网格单元扩展的波前构造法 等。目前还出现了多次波射线追踪技术、以及多值波前追踪技术(如水平集算法、向空 间算法)。 一般来讲，正演的问题比较完善，其因果关系相对于反演则较为明确。而反演的问 题在实际上是相当困难的，实际中反演问题的多解性一直困扰着地球物理学家们，导致 一些现象比较难以解释。进行地震反演的方法也有很多，例如最小二乘法、蒙特卡洛法、 最速下降法、牛顿法、共轭梯度法、以及子空间法等等。 本文内容主要包括三部分：(1) 在地震正演算法中重点研究了快速行进法；(2) 在 地震反演算法中重点讲述了子空间反演技术；(3) 完成了快速行进法结合子空间反演方 法的成像技术。最后成功的实现了初至波和反射波的反演成像，并取得了良好的图像效 果。本文的成果对多次波的应用和成像精度的改善都具有着重要的应用价值和研究意 义。同时，也对复杂构造的反演与解释有着广泛的应用价值。 关键词：地震走时成像，地震射线追踪，快速行进法,共轭梯度法，子空间反演法 Abstr Abstract 目录第一章 目 录 第一章 绪论 1 1.1 本课题研究的重要性及意义 1 1.2 地震层析成像的发展及现状 1 1.3 本文研究的主要内容 3 1.3.1 正演——基于程函方程的快速行进法 3 1.3.2 反演——子空间方法 4 PAGE PAGE IV III III Seismic tomography technology has become one of the most effective methods on studying the earths internal structure and the anisotropy, and seismic traveltime tomography is a method of studying underground velocity structure. Generally speaking, seismic ray tracing is an important part of seismic traveltime tomography methods. There are many ray tracing methods, such as the shooting method, the bending method (or pseudo bending method) and the Gaussian beam method; the grid-cell based method, such as the finite-difference Eikonal equation solver (FD) and the shortest-path method (SPM); and the combined method, such as the wave front construction method.Meanwhile there are some current algorithms, such as multiple arrival tracking and multi-valued wave front tracing (e g, the level set method and the phase space method). Generally speaking, the forward problem is perfect, and its causality is clearly, but the inversion problem is very difficult essentially, its multi-solutions ar