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