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

权利要求书

1.一种基于预计算插值的射线走时计算方法，其特征在于，包括以下步骤：

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

步骤2，在每一所述炮点或检波点的位置发射一组射线，射线到达所述目标层位后终止，并记录每一所述炮点或检波点发出的一组射线的射线数据；

步骤3，根据所述炮点或检波点发出的一组射线的射线数据对所述目标层位的CMP点集进行组内插值计算，获得所述目标层位的CMP点集对应所述炮点或检波点的射线数据；

步骤4，在所述炮点或检波点中选取第一炮点及第一检波点，并从所述目标层位的CMP点集对应所述炮点或检波点的射线数据中获得CMP点集对应第一炮点及第一检波点的射线数据；

步骤5，根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过迭代插值计算获得所述第一炮点及第一检波点在目标层位上的反射点；

步骤6，根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过插值计算获得所述反射点到述第一炮点及第一检波点的射线走时之和。

2.根据权利要求1所述的基于预计算插值的射线走时计算方法，其特征在于，所述步骤2中，每一所述炮点或检波点发出的一组射线的射线数据包括：

每一所述炮点或检波点发出的一组射线中的每一射线与所述目标层位的交点位置，所述每一射线在交点位置的射线方向，所述每一射线到交点位置的走时。

3.根据权利要求2所述的基于预计算插值的射线走时计算方法，其特征在于，所述步骤3，根据所述炮点或检波点发出的一组射线的射线数据对所述目标层位的CMP点集进行组内插值计算，获得所述目标层位的CMP点集对应所述炮点或检波点的射线数据包括：

步骤31，根据所述每一射线与所述目标层位的交点位置，记录所述炮点或检波点发出的一组射线中的相邻两条射线与所述目标层位的交点位置之间覆盖的CMP点；

步骤32，根据所述CMP点的位置以及所述相邻两条射线与所述目标层位的交点位置，将所述相邻两条射线在交点位置的射线方向以及所述相邻两条射线到交点位置的走时插值到所述CMP点，获得所述CMP点对应所述炮点或检波点的射线数据。

4.根据权利要求3所述的基于预计算插值的射线走时计算方法，其特征在于，所述步骤5，根据所述CMP点集对应第一炮点及第一检波点的射线数据，通过迭代插值计算获得所述第一炮点及第一检波点在目标层位上的反射点包括：

步骤51，在第一炮点及第一检波点之间的目标层位上设定一假设点；

步骤52，从所述CMP点集对应第一炮点及第一检波点的射线数据中，获取所述假设点两侧相邻的两个CMP点对应第一炮点及第一检波点的射线数据，并通过插值获得所述假设点对应第一炮点及第一检波点的射线数据；

步骤53，根据所述假设点对应第一炮点及第一检波点的射线数据，分别得到第一炮点的射线在所述假设点的射线方向与所述假设点的法向量的夹角、第一检波点的射线在所述假设点的射线方向与所述假设点的法向量的夹角；

步骤54，比较所述步骤53得到的两夹角的大小关系，并根据所述大小关系在所述假设点的相应一侧以一定步长选取另一假设点；

重复执行所述步骤52至步骤54，当选取到某一假设点，使所述某一假设点的两夹角之差小于一设定阈值时计算终止，所述某一假设点为第一炮点及第一检波点在目标层位上的反射点。