[发明专利]一种地质体变属性插值方法及系统

权利要求书

1.一种地质体变属性插值方法，其特征在于，包括：

步骤1：加载在设定工区范围内已解释的地质体成果数据；

步骤2：根据地质体的控制点数据，构建地质体的表面网格模型；

步骤3：以块体几何表面三角形网格数据作为限定条件，对限定条件进行规则化处理，获得规则的限定点、限定线段、限定面，获得规则化表面网格模型；

步骤4：采用三维限定DELAUNAY剖分算法对所述规则化表面网格模型进行剖分，获得具有拓扑一致性的四面体网格；

步骤5：利用基于梯度的变属性插值算法对步骤4得到的四面体网格进行属性插值，获得插值结果；

利用几何造型方法，构建复杂地质体的表面网格模型，几何造型方法包括：二维限定Delaunay剖分、曲面求交、曲面缝合和复杂地质体快速追踪法；

步骤3具体包括：

步骤301：根据所述三角形网格数据，设置约束线段集和离散点集；

步骤302：将所述约束线段集包含的所有线段的端点记入所述离散点集；

步骤303：选择所述约束线段集中的任意两条线段，判断所述两条线段是否相交，若两条线段相交且共线，则将共线的部分记入所述约束线段集，且在所述约束线段集中删除所述两条线段，若两条线段仅相交不共线，则将交点记入所述离散点集；

步骤304：若三角形网格不满足规则化的要求，则增加新点，判断所述新点是否属于所述离散点集，若不属于所述离散点集，则将所述新点记入所述离散点集；判断所述约束线段集中的线段是否包含所述新点，若包含所述新点，则以所述新点作为端点，将包含所述新点的线段划分为两条新线段，将所述两条新线段记入所述约束线段集，且删除包含所述新点的线段；

步骤305：重复执行步骤303-304，直到所有三角形网格满足规则化的要求为止，获得规则的限定点、限定线段、限定面，获得由规则的限定点、限定线段、限定面组成的规则化表面网格模型。

2.根据权利要求1所述的变属性插值方法，其特征在于，根据三角形的面积和三角形的最小角判定是否满足规则化的要求。

3.根据权利要求1所述的变属性插值方法，其特征在于，步骤4具体包括：

步骤401：采用三维限定DELAUNAY剖分算法对所述规则化表面网格模型进行剖分获得四面体网格；

步骤402：检查限定条件三角形网格数据中每个点、每条边、每个三角形是否在剖分结果中，若不在剖分结果中，则将交点分别加入到当前三角形网格和四面体网格；

步骤403：采用变网格方法对所述四面体网格进行优化，获得具有拓扑一致性的四面体网格。

4.根据权利要求3所述的变属性插值方法，其特征在于，所述变网格方法为根据地质体构造形态增加辅助点来构造四面体网格。

5.根据权利要求1所述的变属性插值方法，其特征在于，采用基于梯度的变属性插值算法对复杂地质体进行属性插值，所述变属性插值算法使用的公式为：

SampleValue[i] ＝ mBaseValue +mGX*gx +mGY*gy + mGZ*gz (1)

其中，SampleValue[i]代表第i个四面体的属性插值结果，mBaseValue代表待插四面体所在块体的基础属性值，mGX、mGY、mGZ分别代表X、Y、Z坐标方向的属性变化梯度，gx、gy、gz分别代表X、Y、Z坐标方向的坐标增量。