[发明专利]基于八叉树的复杂矿体块段模型构建方法及系统

说明书

技术领域

本发明涉及地质体建模技术领域，尤其涉及基于八叉树的复杂矿体块段模型构建方法及系统。

背景技术

三维复杂矿体的建模包括三维地质结构建模和三维属性建模两个部分。结构模型，即表明模型，可以很好地表达复杂矿体的空间形态，但其内部的属性的不均质性，需要通过属性模型对其内部任意一点的属性进行表达和承载。地质体的属性建模通常是按不同的精度要求将地质体空间离散为若干个体元，这些体元可为四面体、六面体、棱柱体和多面体等多种类型，通常将六面体作为体元的属性模型称作为块段模型。

大部分矿业软件(如国外的Vulcan、Datamine、Surpac、MineSight、Micromine等，国内的Dimine、QuantyMine、3DMine等)都采用块段模型作为属性模型。属性模型必须以结构模型的空间边界(即多面体)作为约束，通过判断体元在矿体结构模型的内外关系进行构建。为了更好地拟合复杂矿体的边界，通常边界上的体元尺寸需要尽可能的小。体元尺寸越小，体元个数越多，可能达到几百万上千万甚至更多。因此块段模型构建需要高效的内位关系判断算法，方可在尽可能少点时间内处理如此多的体元。

现有技术中，通过OBB树与八叉树节点相交测试优化了采用射线法判断体元与多面体的内外关系的算法，性能有所提升，但射线与复杂模型求运算交鲁棒性不高，且该算法对多面体模型质量要求较高，不能有开口、自相交、重复面等问题，而这些问题在复杂矿体建模过程中容易出现且难以处理。采用Feito-Torres方法对算法进行了改进，虽然对结构模型的容错能力有所提升，但性能受到一定影响。采用flood-fill算法，时间复杂度降低为O(n),但需要进行三角形与立方体相交运算，且舍弃了八叉树对块段模型的高性能表达。BSP法内存消耗大，在处理非常复杂的模型时其预处理过程要消耗极大的时间和内存，甚至可能使系统崩溃。

发明内容

针对现有技术的缺陷，本发明提供基于八叉树的复杂矿体块段模型构建方法及系统，将适应于栅格扫描法的规则网格与八叉树模型进行关联，实现基于八叉树的块段模型高效构建，大大降低了算法的时间复杂度，提高了算法的鲁棒性。

第一方面，本发明提供一种基于八叉树的复杂矿体块段模型构建方法，该方法包括：

根据预设的最小体元尺寸，将目标矿体的地质空间划分为多个三维网格，所述网格行、列、层数分别为N，其中N＝2^l，l为八叉树的层级数；

根据所述三维网格对应的体元与所述地质空间边界的内外关系，通过栅格扫描法对所述体元进行筛选；

根据八叉树节点的层级关系，对筛选后的体元进行合并，形成八叉树模型。

优选地，所述根据所述三维网格对应的体元与所述地质空间边界的内外关系，通过栅格扫描法对所述体元进行筛选，包括：

针对所述多个三维网格的每一层，采用通过该层的平面切割所述矿体的表面模型，获得该层的二维轮廓线；

根据该层的网格与所述二维轮廓线的内外关系，采用栅格扫描法对所述网格进行筛选，获得在所述二维轮廓线内部的网格所对应的体元。

优选地，所述根据八叉树节点的层级关系，对筛选后的体元进行合并，形成八叉树模型，包括：

判断八叉树中每个父节点的八个子节点对应的体元是否均位于所述地质空间的内部，若是，则合并所述八个子节点所对应的体元，直至体元的尺寸等于预设的基本体元尺寸。

优选地，所述根据该层的网格与所述二维轮廓线的内外关系，采用栅格扫描法对所述网格进行筛选，包括：

初始化所有体元为“外”标记；

取第k层网格的第j列，获取x＝x_j处平行于y轴的一条射线

获取第k层对应的二维轮廓线上所有线段的两个端点的x轴坐标值(x₁,x₂)，与x_j进行比较，若x₁<x_j<x₂，则所述线段与有交点，否则所述线段与没有交点；

获得所有线段与射线的交点的y轴坐标值，并将所有y轴坐标值从小到大排序，得到序列{y₀,y₁,y₂,…y_c-1}；