(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111554016.3 (22)申请日 2021.12.17 (71)申请人 中国地质大 学 （武汉） 地址 430074 湖北省武汉市洪山区鲁磨路 388号 (72)发明人 李章林　王平　张夏林　汪宁　 廉令坤　许颖颖　翁正平　张志庭　 张军强　 (74)专利代理 机构 北京东方盛凡知识产权代理 事务所(普通 合伙) 11562 代理人 李娜 (51)Int.Cl. G06T 17/00(2006.01) G06Q 10/04(2012.01) G06Q 50/02(2012.01) (54)发明名称 一种基于八叉树的局部异向性搜索椭球体 动态建模方法 (57)摘要 本发明公开了一种基于八叉树的局部异向 性搜索椭球体动态建模方法， 包括： 设定约束条 件并基于待测矿体构建矿体八叉树模 型； 基于约 束条件对矿体 八叉树模型进行剖分， 获得若干个 子空间并在每个子空间内生成局部异向性搜索 椭球体模型。 本发明在研究区范围较大或矿体形 态复杂的情况下， 通过八叉树模 型对研究空间进 行合理划分， 实现了对海量局部异向性搜索椭球 体进行合理 高效的组织和管 理， 减少了建模的过 程中存储空间浪费和效率低下的问题， 从而提高 地质统计学储量估算方法的适用性和实用性。 权利要求书2页 说明书5页 附图4页 CN 114170388 A 2022.03.11 CN 114170388 A 1.一种基于八叉树的局部异向性搜索 椭球体动态建模方法， 其特 征在于， 包括： 设定约束条件并基于待测矿体构建矿体八叉树模型； 基于所述约束条件对所述矿体八叉树模型进行剖分， 获得若干个子空间并在每个所述 子空间内生成局部异向性搜索 椭球体模型。 2.根据权利要求1所述的方法， 其特 征在于， 所述约束条件 包括： 对结点的约束： 结点的尺度大于最小允许值； 对子空间的约束： 每 个子空间的尺度小于等于最大允许值； 对差异度的约束： 差异度小于预设的差异度阈值。 3.根据权利要求1所述的方法， 其特征在于， 基于待测矿体构建矿体八叉树模型的过程 包括： 计算待测矿体整体对应的长方体包围盒， 将所述长方体包围盒作为所述矿体八叉树 模型的根结点， 并分别计算所述矿体八叉树模型中每个结点的结点坐标的最大值和最小 值。 4.根据权利要求2所述的方法， 其特征在于， 基于所述约束条件对所述矿体八叉树模型 进行剖分， 获得若干个子空间并在每个所述子空间内生成局部异向性搜索椭球体模型的过 程包括： 判断结点是否满足所述约束条件， 若满足则对所述结点进行剖分， 获得若干个子空间； 判断若干个所述子空间是否满足所述约束条件， 若满足则分别在每个所述子空间内生 成局部异向性搜索 椭球体模型； 计算所述若干个子空间对应的局部异向性搜索椭球体之间的差异度， 若满足所述约束 条件， 则结束剖分。 5.根据权利要求4所述的方法， 其特征在于， 所述判断结点是否满足所述约束条件， 若 满足则对所述结点进行剖分， 获得若干个子空间的过程中包括： 若所述结点的结点尺度大于最小允许值， 则将结点均匀剖分成八个所述子空间， 所述 结点的尺度若 小于等于最小允许值则不必进行剖分。 6.根据权利要求4所述的方法， 其特征在于， 判断若干个所述子空间是否满足所述约束 条件， 若满足则分别在每 个所述子空间内生成局部异向性搜索 椭球体模型的过程中包括： 若存在任意一个所述子空间的尺度大于最大允许值， 则再次进行剖分， 直到所有所述 子空间的尺度均小于等于最大允许值时， 分别在每个所述子空间内生成局部异向性搜索椭 球体模型。 7.根据权利要求4所述的方法， 其特征在于， 计算所述若干个子空间对应的局部异向性 搜索椭球体之间的差异度， 若满足所述约束条件， 则结束剖分的过程中包括： 预设所述差异度阈值， 计算每个所述子空间对应的局部异向性搜索椭球体之间的差异 度， 并与所述差异度阈值比较， 若所述差异度大于所述差异度阈值， 则需继续剖分； 若所述差异度小于所述差异度阈值则结束剖分。 8.根据权利要求7所述的方法， 其特征在于， 计算每个所述子空间对应的椭球体之间的 差异度的过程中包括： 分别基于角度方差和变异系数计算每个所述局部异向性搜索椭球体的差异度， 其中所 述差异度包括角度差异度和轴长 差异度。 9.根据权利要求8所述的方法， 其特 征在于， 计算所述角度差异度的公式Sθ为：权　利　要　求　书 1/2 页 2 CN 114170388 A 2上式中， θ表示任意 一个角度参数。 10.根据权利要求8所述的方法， 其特 征在于， 计算所述轴长 差异度的公式为： 上式中， x表示任意 一个轴长参数， CVx表示轴长向量x的变异性， μ表示8个轴长的均值。权　利　要　求　书 2/2 页 3 CN 114170388 A 3