[发明专利]基于三角剖分模型的正演方法及计算机可读存储介质

说明书

公开了一种基于三角剖分模型的正演方法及计算机可读存储介质。该方法包括：给定离散控制点，将速度模型进行Delaunay三角剖分；在速度模型三角剖分的基础上，进行动态射线追踪；利用高斯束传播算子获得理论地震记录。本发明提供的技术方案，对复杂构造的适应性更好，减少了有限差分模拟的网格“锯齿”效应引起的散射，提高了波场模拟的精度；并且计算效率高，能直接得到去除直达波的地震记录。

技术领域

本发明涉及地震波场数值模拟领域，更具体地，涉及一种基于三角剖分模型的正演方法及计算机可读存储介质。

背景技术

波场正演模拟技术在地震勘探中应用非常广泛。通过正演模拟可以为新技术的提出、可行性分析和应用试验提供高质量的模拟数据；可以直接指导野外地震数据的采集、处理和解释。目前数值模拟主要可以归为三大类：波动方程数值解法、积分方程法和射线追踪法。高斯束正演模拟是射线追踪方法类的一种方法，计算效率较高且保留了波的动力学信息。

地震勘探中，无论是波场正演模拟还是层析反演，模型参数化是最根本的问题。当前主流技术都是基于规则网格划分(常见的有矩形网格剖分或简单的三角网格剖分)，如有限差分正演、反射波层析都是基于矩形网格的。采用矩形网格参数化时，存在很多缺点，具体如下：

(1)模型剖分灵活性差，矩形网格的均一性要求无论在构造复杂或简单区域都使用同样的精度进行划分；(2)对速度界面的描述精度差，容易产生“锯齿”状边界；(3)复杂模型时，需要精细网格描述速度模型，计算量大；(4)用于层析反演时存储量大、计算时间长、方程性态差、求解困难。

因此，有必要开发一种划分速度模型的正演模拟方法，从而克服基于规则网格模型划分灵活性差、计算量大的缺陷。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种基于三角剖分模型的正演方法及计算机可读存储介质，其能够通过Delaunay三角剖分和高斯束正演方法，实现获得理论地震记录。

根据本发明的一方面，提出了一种基于三角剖分模型的正演方法，该方法包括：给定离散控制点，将速度模型进行Delaunay三角剖分；在速度模型三角剖分的基础上，进行动态射线追踪；利用高斯束传播算子获得理论地震记录。

根据本发明的另一方面，提出了一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现以下步骤：

步骤1，给定离散控制点，将速度模型进行Delaunay三角剖分；

步骤2，在速度模型三角剖分的基础上，进行动态射线追踪；

步骤3，利用高斯束传播算子获得理论地震记录。

本发明提供的技术方案，对复杂构造的适应性更好，减少了有限差分模拟的网格“锯齿”效应引起的散射，提高了波场模拟的精度；并且计算效率高，能直接得到去除直达波的地震记录。

本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。

附图说明

通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。

图1示出了根据本发明的一种基于三角剖分模型的正演方法的步骤的流程图。

图2示出了根据本发明的一个实施例的四个控制点的两种划分方式的示意图。