最新地震模拟与逆时偏移中的无单元法及其效率分析.pdf

摘要 有限差分法和有限元法是两种常见的也是最成熟的数值计算方法。这两种方 法需要事先划分大量的网格，前期预处理工作比较繁琐，随着研究问题的深入和 研究领域的扩展，有限元法由于单元的限制而无法很好的解决某些问题。无单元 法则从根本上抛弃了单元划分这一思想，只需要研究区域中节点和边界的信息， 具有局部拟合的特点，因此使得问题的处理相对简单、灵活。 无单元法的缺陷是计算量大，成本高，目前主要用于二维的线性问题上， 而实际问题都是三维的，成本问题严重限制了其应用和发展。本文通过压缩稀疏 矩阵和优化稀疏矩阵的操作，在一定程度上节省了无单元法计算的成本，提高了 其在地震模拟和成像中应用的可行性。 无单元法用于地震逆时偏移成像的问题主要集中在成本方面。在地震波场 迭代过程中，涉及到多个大型稀疏矩阵相加、相乘和求逆。为了节省成本，提高 计算效率，我们对稀疏矩阵进行压缩存储，只保存非零元，矩阵操作也只针对非 零元：同时，稀疏矩阵的求逆非常困难，既消耗时间又消耗内存，我们采用求解 线性方程组的方法避免求逆。这样在很大程度上节省了内存要求，提高了计算效 率；另外，我们对迭代方程进行了简单的变形，使计算成本进一步降低。通过这 些手段，我们在保证成像精度的前提下，提高了计算效率，节省了计算成本，为 无单元法在地震波场模拟与偏移成像中的推广应用提供了可行性。 本文首先详细介绍了地震波场模拟与逆时偏移成像技术的一些概念和方法， 紧接着介绍了无单元法的历史发展情况，无单元法的核心优势和基本理论，随后 将无单元法引入到地震波动方程的求解中，并通过压缩稀疏矩阵和优化稀疏矩阵 操作提高了计算效率，降低了计算成本，然后将高效率的无单元法用于求解地震 波波动方程，并对几种比较有代表性的速度模型进行了模拟与成像(包括叠后与 叠前)，采用了成像时刻和互相关两种成像条件。计算结果表明无单元法对于各 种简单和复杂模型都有很好的适应性，结果精度高，成像清晰。 关键词：地震模拟逆时偏移无单元法稀疏矩阵 Abstract Thefinitedif琵rence thefiniteelement memod(FDM)and numerical allkindsof andmaturememodsin themostcommon excluslVe element．free many ments， withthesetwomethods，the method(EFM)shows ot istheabsence whichis inrecent ofthismethod keypoint developedyears．The method andmore dataare makesthis cheaper nodal elements，only required，which and moretlme flexible． saV