[发明专利]一种地层界线图层的自动化生成方法

权利要求书

1.一种地层界线图层的自动化生成方法，其特征在于，该方法包括如下步骤：

(1)基于沉积、火山地层界线图层和变质地层单位图层，搜寻出图层中所有邻接地层，构建地层邻接关系图；

(2)基于地层邻接关系图提取相邻地层的公共边界，生成地层界线；

(3)基于断层图层，识别断层接触关系；

(4)基于产状图层，识别基本接触关系，完成地层界线图层的自动化生成。

2.根据权利要求1所述的一种地层界线图层的自动化生成方法，其特征在于，所述步骤(1)的具体步骤如下：

步骤11：加载沉积、火山地层界线图层和变质地层单位图层数据，并将相应的地层面要素读入集合S＝{s_i|i＝1，2，…，n}，n为地层面要素的个数；

步骤12：读出每个地层面要素s_i的顶点集v_i＝{P_j|j＝1，2，…，m}，其中P_j表示构成地层面要素s_i的顶点，m为顶点的个数；并对每个地层面要素s_i的顶点集v_i进行如下的预处理：

a)将顶点集v_i中的点P_j按顺时针排序；

b)判断每个顶点集v_i中是否存在P₁＝P_m，即首节点等于尾节点，若存在，则跳转下一步c)，对顶点集v_i进行处理；否则，跳转步骤13；

c)删除顶点集v_i的尾节点P_m，保证首尾不相接，即v_i＝{P_j|j＝1，2，…，m-1}；

步骤13：对图层中地层面要素的邻接关系进行初判：若地层面要素s_a的顶点集v_a的最小外接矩形R_a和要素s_b的顶点集v_b的最小外接矩形R_b有交点，则跳转步骤14；否则，跳转步骤16；其中，a为循环变量，a∈[1，n-1]，初值为1；b为循环变量，b∈[a+1，n]，初值为a+1；

步骤14：对两地层面要素s_a和s_b的邻接关系进行细判，具体过程如下：

a)计算顶点集v_a和v_b的公共点集T＝{p₁，p₂，...，p_k}，k为公共点的个数，并记录T中点在顶点集v_a中序号，记为N＝{n₁，n₂，...，n_k}；

b)若公共点集T中点的个数k≥2，则跳转下一步c)；否则，跳转步骤16；

c)判断公共点序号集N内是否存在连续的点号n_t，n_t+1满足∣n_t－n_t+1∣＝1，若存在，即地层面要素s₁和s₂存在邻接关系，则跳转步骤15；否则，跳转步骤16；

步骤15：构建地层邻接关系图的边d：计算顶点集v_a与v_b的重心点C_a和C_b，点C_a和C_b作为地层邻接关系图的顶点，连接点C_a和C_b作为地层邻接关系图的边d；

步骤16：若循环变量a＞n-1，则输出SARG＝{d₁，d₂，...d_w}，w为地层邻接关系图SARG的边的总个数；否则，循环变量b加1，此时若b≤n，则直接跳转步骤13；否则循环变量a加1，循环变量b＝a+1，再跳转步骤13。