一种基于预计算插值的射线走时计算方法及系统技术方案

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

【技术实现步骤摘要】
【专利摘要】本专利技术公开了一种基于预计算插值的射线走时计算方法及系统，其中，该方法包括：建立实验地区的速度模型，设定炮点、检波点、目标层位；在每一炮点或检波点的位置发射一组射线，射线到达目标层位后终止，记录每一炮点或检波点发出的一组射线的射线数据；根据射线数据对目标层位的CMP点集进行组内插值计算，获得目标层位的CMP点集对应炮点或检波点的射线数据；选取第一炮点及第一检波点，获得CMP点集对应第一炮点及第一检波点的射线数据；根据CMP点集对应第一炮点及第一检波点的射线数据，通过迭代插值计算获得第一炮点及第一检波点在目标层位上的反射点；通过插值计算获得反射点到述第一炮点及第一检波点的射线走时之和。【专利说明】一种基于预计算插值的射线走时计算方法及系统
本专利技术是关于石油勘探技术，主要涉及地震数据处理技术，尤指一种基于预计算插值的射线走时计算方法。
技术介绍
在地震数据处理技术中，介质速度的计算可以利用相干反演技术，根据给定的若干个扫描速度，对其中每一个扫描速度生成一个相应的介质模型，在此模型上计算射线走时并与道集数据比较得出在当前速度模型下的相似度，最终获得介质速度。因此，射线走时的计算很大程度影响到最终结果的精确程度和计算效率。现有的射线走时计算方法一般是采用射线追踪方法，这种方法具有很高的精度，但是计算效率低下，而且在大规模数据处理中往往存在严重的重复计算问题。
技术实现思路
本专利技术为了克服现有的射线走时计算方法的不足，提出了一种基于预计算插值的射线走时计算方法及系统，以达到高效率的走时计算并能保证计算精度。 为达到上述目的，本专利技术提出的基于预计算插值的射线走时计算方法包括以下步骤:步骤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，采集实验地区数据建立实验地区的速度模型，在速度模型上设定炮点、检波点，并设定速度模型的目标层位。 步骤2，在每一炮点或检波点的位置发射一组射线，射线到达目标层位后终止，并记录每一炮点或检波点发出的一组射线的射线数据。 在本步骤中，每一炮点或检波点发出的一组射线的射线数据包括:每一炮点或检波点发出的一组射线中的每一射线与目标层位的交点位置，每一射线在交点位置的射线方向，每一射线到交点位置的走时。 步骤3，根据炮点或检波点发出的一组射线的射线数据对目标层位的CMP点集进行组内插值计算，获得目标层位的CMP点集对应炮点或检波点的射线数据。 步骤4，在炮点或检波点中选取第一炮点及第一检波点，并从目标层位的CMP点集对应炮点或检波点的射线数据中获得CMP点集对应第一炮点及第一检波点的射线数据。 步骤5，根据CMP点集对应第一炮点及第一检波点的射线数据，通过迭代插值计算获得第一炮点及第一检波点在目标层位上的反射点。 步骤6，根据CMP点集对应第一炮点及第一检波点的射线数据，通过插值计算获得反射点到述第一炮点及第一检波点的射线走时之和。 在本实施例中，对于步骤3有以下详细流程，结合图2所示，步骤3具体包括: 步骤31，根据每一射线与目标层位的交点位置，记录炮点或检波点发出的一组射线中的相邻两条射线与目标层位的交点位置之间覆盖的CMP点。 步骤32，根据CMP点的位置以及相邻两条射线与目标层位的交点位置，将相邻两条射线在交点位置的射线方向以及相邻两条射线到交点位置的走时插值到CMP点，获得CMP点对应炮点或检波点的射线数据。 在本实施例中，对于步骤5有以下详细流程，结合图3所示，步骤5具体包括: 步骤51，在第一炮点及第一检波点之间的目标层位上设定一假设点。 步骤52，从CMP点集对应第一炮点及第一检波点的射线数据中，获取假设点两侧相邻的两个CMP点对应第一炮点及第一检波点的射线数据，并通过插值获得假设点对应第一炮点及第一检波点的射线数据。 步骤53，根本文档来自技高网...

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

【技术特征摘要】

【专利技术属性】
技术研发人员：梁兼栋，郝晓光，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11