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

权利要求书

1.一种高效便捷的地震数据非均匀间距等值线交互追踪方法，特点是采用以下步骤实现：

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

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

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

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；

10）在种子点dp_cross处进行高程值标注，完成追踪一条等值线；将所述的等值线和追踪完的所有等值线进行合并，以避免重复的等值线；

11）重复步骤6）-10），追踪生成一系列等间距等值线；

12）重复步骤5）-11），追踪生成一系列非等间距等值线。

2.根据权利要求1的方法，特点是步骤1）所述的层位数据由一系列点d(x,y,z)构成的，其中x,y,z表示地质解释层位数据的空间三维坐标；所述的x,y为大地坐标，z为高程坐标，对于任意两个数据点，其x,y不完全相等。

3.根据权利要求1的方法，特点是步骤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坐标的最小值、最大值。

4.根据权利要求1的方法，特点是步骤9）所述的优化传播法如下：

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

采用下列公式求取等值点c(x,y,z):x=xa+(xb-xa)z-zazb-zay=ya+(yb-ya)z-zazb-zb]]>

式中：给定的高程值z，

三角形tri_cur的三条边两个顶点为a(x_a,y_a,z_a)和b(x_b,y_b,z_b)，

每个三角形在三边上有2个等值点p₁、p₂；

将与三角形1tri_cur以点p₁所在的边为相邻边的三角形2tri_ncur设为当前三角形tri_cur，点p₁所在的边为进入边；

Step2在三角形2tri_cur中搜索等值线的离去边，即找到进入相邻三角形6tri_ecur的边，假设该边为e_next，找到后内插该等值点的平面坐标，内插的方法和Step1中方法相同，只是当前边不参与计算；

Step3当内插等值点不是边的两个顶点之一，则找到e_next的相邻三角形6tri_ecur，设该三角形6为tri_cur，当前离去边为进入边。

当内插等值点是边的两个顶点之一时，需要遍历所有包含该等值点的三角形，逐个三角形进行Step2操作，找到当前点距离最短的等值点所在的离去边，设该三角形11为tri_cur，当前离去边为进入边；

Step4重复Step2-Step3，直到离去边为边界边，此时等值线为开曲线L₁

Step5追踪另一半等值线，设与Step1中的三角形1以点p₁所在的边为相邻边的三角形3为当前三角形tri_cur，以点p₂所在的边为进入边，重复Step2-Step4操作，获得另一段开曲线L₂；

Step6合并开曲线L₁和L₂，即得到一条等值线L_contour。