[发明专利]一种地震数据非均匀间距等值线交互追踪方法

说明书

技术领域

本发明涉及地球物理勘探数据等值线成图技术，属于地球物理勘探数据处理中数据成图范畴，具体是一种地震数据非均匀间距等值线交互追踪方法。

背景技术

地球物理勘探数据处理时，需要将数据生成各种图形，在地质层位数据成图中，等值线成图是一个重要的组成部分，其中的等值线追踪又在等值线成图中起到举足轻重的作用。常规等值线追踪都是采取的自动追踪方法，这类方法针对给定的离散数据，采用四角网格或者三角网将整个数据划分为一个个四边形或三角形。然后，从数据边界出发，按照线性插值等方式逐个四边形或三角形进行插值，并在插值基础上追踪，一次性生成所有的等值线，这种全局性等值线追踪的速度慢，等值线显示疏密难以控制，在数据变化复杂的地方很难完全满足用户的需求。

采用全局追踪方法，算法要求等值线间距、标注的设定全都全局统一。在此前提下极易发生在数据变化平滑地区等值线过于稀疏，而在数据变化剧烈地区等值线过于密集，甚至发生不同值等值线重叠等现象。而且，由于等值线标注信息位置是自动计算生成的，常常发生等值线信息覆盖其他有效信息的问题。

发明内容

本发明目的是提供一种高效便捷的地震数据非均匀间距等值线交互追踪方法。

本发明的具体实施步骤如下：

1）在工区采集任意一层地质解释层位数据；

所述的层位数据由一系列点d(x,y,z)构成的，其中x,y,z表示地质解释层位数据的空间三维坐标。

所述的x,y为大地坐标，z为高程坐标，对于任意两个数据点，其x,y不完全相等。

2）采用逐点比较法求出层位数据的三维数据包围盒，并以此作为数据的有效范围；

所述的三维数据包围盒是(x_min,x_max)、(y_min,y_max)、(z_min,z_max)，其中(x_min,x_max)分别是所有数据点x坐标的最小值、最大值，(y_min,y_max)分别是所有数据点y坐标的最小值、最大值，(z_min,z_max)分别是所有数据点z坐标的最小值、最大值；

3）在大地坐标x,y平面上，采用带约束条件的Delaunay(狄洛尼)不规则三角网（TIN）生成算法构建层位数据TIN网；

4）采用逐点比较法统计不规则三角网中所有点的z值范围，采用尾数截断法自动计算等值线起始值z_start和步长z_step；

5）根据数据范围确定等值线追踪的起始值z_start和步长z_step；

6）任意连续点选在层位数据生成并显示在图上的不规则三角网上，并根据二次曲线法依据这些点定义一条曲线L；

7）采用两段直线求焦法求取曲线L与不规则三角网的交点系列SP_cross，并根据双线性法计算交点系列SP_cross中每个交点p_cross的z值；

8）连接交点系列SP_cross，针对交点系列SP_cross中的任意两两相邻点p_i、p_i+1采用线性平均插值方法获得满足定义的等值线点系列DP_cross；

9）以等值线点系列DP_cross为种子点系列，针对种子点系列中的每个种子点dp_cross进行如下处理：

采用最小面积法计算出种子点dp_cross在不规则三角网中所在的三角形tri_cur，以三角形tri_cur为起点，采用优化传播法生成一条等值线L_contour；

所述的优化传播法如下：

Step1将采用最小面积法计算所得到的三角形1tri_cur设为当前三角形tri_cur；