[发明专利]一种基于预计算插值的射线走时计算方法及系统

说明书

技术领域

本发明是关于石油勘探技术，主要涉及地震数据处理技术，尤指一种基于预计算插值的射线走时计算方法。

背景技术

在地震数据处理技术中，介质速度的计算可以利用相干反演技术，根据给定的若干个扫描速度，对其中每一个扫描速度生成一个相应的介质模型，在此模型上计算射线走时并与道集数据比较得出在当前速度模型下的相似度，最终获得介质速度。因此，射线走时的计算很大程度影响到最终结果的精确程度和计算效率。现有的射线走时计算方法一般是采用射线追踪方法，这种方法具有很高的精度，但是计算效率低下，而且在大规模数据处理中往往存在严重的重复计算问题。

发明内容

本发明为了克服现有的射线走时计算方法的不足，提出了一种基于预计算插值的射线走时计算方法及系统，以达到高效率的走时计算并能保证计算精度。

为达到上述目的，本发明提出的基于预计算插值的射线走时计算方法包括以下步骤：步骤1，采集实验地区数据建立所述实验地区的速度模型，在所述速度模型上设定炮点、检波点，并设定所述速度模型的目标层位；步骤2，在每一所述炮点或检波点的位置发射一组射线，射线到达所述目标层位后终止，并记录每一所述炮点或检波点发出的一组射线的射线数据；步骤3，根据所述炮点或检波点发出的一组射线的射线数据对所述目标层位的CMP点集进行组内插值计算，获得所述目标层位的CMP点集对应所述炮点或检波点的射线数据；步骤4，在所述炮点或检波点中选取第一炮点及第一检波点，并从所述目标层位的CMP点集对应所述炮点或检波点的射线数据中获得CMP点集对应第一炮点及第一检波点的射线数据；步骤5，根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过迭代插值计算获得所述第一炮点及第一检波点在目标层位上的反射点；步骤6，根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过插值计算获得所述反射点到述第一炮点及第一检波点的射线走时之和。

为达到上述目的，本发明还提出了一种基于预计算插值的射线走时计算系统，包括：模型建立模块，用于根据实验地区的数据建立所述实验地区的速度模型，在所述速度模型上设定炮点、检波点，并设定所述速度模型的目标层位；射线数据采集模块，用于在每一所述炮点或检波点的位置发射一组射线，射线到达所述目标层位后终止，并记录每一所述炮点或检波点发出的一组射线的射线数据；插值计算模块，用于根据所述炮点或检波点发出的一组射线的射线数据对所述目标层位的CMP点集进行组内插值计算，获得所述目标层位的CMP点集对应所述炮点或检波点的射线数据；射线数据选取模块，用于在所述炮点或检波点中选取第一炮点及第一检波点，并从所述目标层位的CMP点集对应所述炮点或检波点的射线数据中获得CMP点集对应第一炮点及第一检波点的射线数据；反射点获取模块，用于根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过迭代插值计算获得所述第一炮点及第一检波点在目标层位上的反射点；射线走时计算模块，用于根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过插值计算获得所述反射点到述第一炮点及第一检波点的射线走时之和。

本发明的基于预计算插值的射线走时计算方法及系统能够高效地计算相干反演算法中所需要的射线走时信息，并且也能达到较高的精度，同时该算法可以并行运行，提高了地震勘探数据处理的效率。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1为本发明一实施例的基于预计算插值的射线走时计算方法流程图。

图2为本发明一实施例的基于预计算插值的射线走时计算方法中步骤3的详细流程图。

图3为本发明一实施例的基于预计算插值的射线走时计算方法中步骤5的详细流程图。

图4为本发明一具体实施例的射线方向与目标层位的法向量的夹角示意图。

图5为本发明一实施例的基于预计算插值的射线走时计算系统的结构示意图。

图6为本发明一实施例的基于预计算插值的射线走时计算系统中插值计算模块的结构示意图。

图7为本发明一实施例的基于预计算插值的射线走时计算系统中反射点获取模块的结构示意图。

具体实施方式

以下配合图式及本发明的较佳实施例，进一步阐述本发明为达成预定发明目的所采取的技术手段。

图1为本发明一实施例的基于预计算插值的射线走时计算方法流程图。如图1所示，该方法包括以下步骤：

步骤1，采集实验地区数据建立实验地区的速度模型，在速度模型上设定炮点、检波点，并设定速度模型的目标层位。