CN112036030B - 矿体组合约束建模方法、装置、设备及存储介质 - Google Patents
矿体组合约束建模方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN112036030B CN112036030B CN202010892365.5A CN202010892365A CN112036030B CN 112036030 B CN112036030 B CN 112036030B CN 202010892365 A CN202010892365 A CN 202010892365A CN 112036030 B CN112036030 B CN 112036030B
- Authority
- CN
- China
- Prior art keywords
- implicit
- function
- implicit function
- ore body
- combined
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000006870 function Effects 0.000 claims abstract description 300
- 239000011435 rock Substances 0.000 claims abstract description 20
- 238000012946 outsourcing Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 12
- 210000000746 body region Anatomy 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000004441 surface measurement Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 238000012847 principal component analysis method Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 description 14
- 230000001360 synchronised effect Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000033558 biomineral tissue development Effects 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 229910006119 NiIn Inorganic materials 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 235000019800 disodium phosphate Nutrition 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000005405 multipole Effects 0.000 description 1
- 238000013138 pruning Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种矿体组合约束建模方法、装置、设备及存储介质。其中，该方法包括：基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果。本发明实施例利用组合隐式函数场，可以更灵活地构建符合地质规则的约束条件，从而提高矿体建模的精确度，且基于操作树的分层结构求取各点在所述组合隐式函数场中的函数值，可以减少计算量，从而提高建模效率。
Description
技术领域
本发明涉及矿体领域，尤其涉及一种矿体组合约束建模方法、装置、设备及存储介质。
背景技术
矿体的赋存情况本身是确定的，但是由于受地下空间、矿化特征和地勘手段的局限性，地质勘探难以获得精确刻画矿体形态及其分布的完整数据。基于隐式建模的空间插值数学模型并未反映取决于特定矿床类型的地质规则，需要根据先验信息构造约束条件来反映相应的地质规则。单隐式函数场约束方法难以反应涉及多域类型的地质规则，所重建的矿体模型容易产生不符合域间拓扑关系的问题。
发明内容
有鉴于此，本发明实施例提供了一种矿体组合约束建模方法、装置、设备及存储介质，旨在使得矿体模型尽量符合地质规则，并提高矿体建模的效率及精确度。
本发明实施例的技术方案是这样实现的：
第一方面，本发明实施例提供了一种矿体组合约束建模方法，包括：
基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果；
其中，所述求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值。
第二方面，本发明实施例还提供了一种矿体组合约束建模装置，包括：
构建模块，用于基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
曲面确定模块，用于遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果；
所述曲面确定模块求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值。
第三方面，本发明实施例还提供了一种矿体组合约束建模设备，包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器，其中，所述处理器，用于运行计算机程序时，执行本发明实施例所述方法的步骤。
第四方面，本发明实施例还提供了一种存储介质，所述存储介质上存储有计算机程序，所述计算机程序被处理器执行时，实现本发明实施例所述方法的步骤。
本发明实施例提供的技术方案，基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果。利用组合隐式函数场可以更灵活地构建符合地质规则的约束条件，从而提高矿体建模的精确度，且基于操作树的分层结构求取各点在所述组合隐式函数场中的函数值，可以减少计算量，从而提高建模效率。
附图说明
图1为本发明实施例矿体组合约束建模方法的流程示意图；
图2A至图2D为一应用示例中两个隐式函数场的轮廓线和势场示意图；
图3A至图3D示出了不同组合规则下组合隐式函数场的效果示意图；
图4A至图4D为一应用示例中矿体建模过程的示意图；
图5为本发明实施例矿体组合约束建模装置的结构示意图；
图6为本发明实施例矿体组合约束建模设备的结构示意图。
具体实施方式
下面结合附图及实施例对本发明再作进一步详细的描述。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。
本发明实施例提供了一种矿体组合约束建模方法，如图1所示，包括：
步骤101，基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
步骤102，遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果。
本发明实施例中，所述求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值。
这里，建模结果可以为基于各点的函数值得到的隐式曲面。
本发明实施例矿体组合约束建模方法，基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果。利用组合隐式函数场可以更灵活地构建符合地质规则的约束条件，从而提高矿体建模的精确度，且基于操作树的分层结构求取各点在所述组合隐式函数场中的函数值，可以减少计算量，从而提高建模效率。
在一些实施例中，所述基于矿体区域的地质规则构建多个隐式函数，包括：
所述基于矿体区域的地质规则构建第一隐式函数、第二隐式函数、第三隐式函数及第四隐式函数，其中，所述第一隐式函数表征几何边界线约束，用于控制模型的范围；所述第二隐式函数表征矿体上边界面；所述第三隐式函数表征矿体下边界面；所述第四隐式函数表征岩层界面。
这里，隐式函数f(x)可以被表示为有符号距离场。如果将矿体的表面模型看作隐式函数的零水平集S＝{x∈R3|f(x)＝0}，则矿体模型的内部和外部可以直接根据函数值的符号来区分。不失一般性地，约定模型外部的函数值的符号为正值，模型内部的函数值的符号为负值。地质域中任一点函数值的具体大小可以定义为该点到隐式曲面S上最近点的距离。
当采用单个隐式函数表示矿体模型比较复杂时，可以利用多个隐式函数组合后的组合隐式函数场来表示矿体模型。组合隐式函数场的隐式函数F(x)可以表示为
其中，fi(x)表示第i个子隐式函数，符号
表示隐式函数间的某种运算规则。不同的隐式函数fi(x)用于表示不同的组合约束。组合隐式函数场在本质上仍然是一个隐式函数。建模的目标仍然是把组合隐式函数场的隐式函数F(x)重构出来。
组合隐式函数场的各个子隐式函数场可以采用不同的隐式函数来表达，比如，fi(x)，i＝1，2，...，n可以采用不同类型的径向基类插值方法来表示。
实际应用中，为了限定模型空间中隐式函数间相互作用的范围，需要为各个子隐式函数指定模型外包。模型外包也可以看作一种隐式函数。
子隐式函数场i的模型外包的隐式函数可以表示为
当
时，子隐式函数场i的组合约束有效；当
时，子隐式函数场i的组合约束无效。因此，对于第i个子隐式函数场来说，模型外包限制了fi(x)≤0的区域。模型外包除了可以用于限制组合约束的作用域之外，还可以用于加速评估隐式函数F(x)的函数值。
在一些实施例中，可以采用主成分分析(Principal Component Analysis，PCA)方法自动计算的最佳有向包围盒来表示模型外包。如果需要调整模型外包，可以通过设置一定的比例系数来缩放最小包围盒。
示例性地，对于空间中的某个点，只要评估该点函数值的符号就可以确定该点在子隐式函数场中的位置。基于实体构造几何(CSG)的集合运算规则，定义交集、并集和差集三种最基本的组合操作。组合操作之间可以采用嵌套的方式来构造组合规则。嵌套的组合操作可以得到一颗操作树。
子隐式函数场和组合隐式函数场都可以被看作操作树中的一个节点。图2A至图2D展示了隐式函数场A和隐式函数场B两个实例。其中，图2A示出了隐式函数场A的轮廓线，图2B示出了隐式函数场A的势场；图2C示出了隐式函数场B的轮廓线，图2D示出了隐式函数场B的势场。
示例性地，图3A至图3D示出了不同组合规则下组合隐式函数场的效果示意图。
如图3A所示，当重构的模型必须限制在某个隐式函数场的内部区域时，可以采用交集操作。由隐式函数场A和隐式函数场B的交集A∩B所构成的组合约束可以表示为
F(x)＝fA∩B(x)＝max(fA(x)，fB(x))
其中，fA(x)为表示隐式函数场A的隐式函数，fB(x)为表示隐式函数场B的隐式函数。
组合隐式函数场相应的模型外包的隐式函数BF(x)可以表示为
其中，
为表示隐式函数场A的模型外包的隐式函数，
表示隐式函数场B的模型外包的隐式函数。
如图3B所示，当重构的模型应该包含某个隐式函数场的内部区域时，可以采用并集操作。由隐式函数场A和隐式函数场B的并集A∪B所构成的组合约束可以表示为
F(x)＝fA∪B(x)＝min(fA(x)，fB(x))
组合隐式函数场相应的模型外包的隐式函数BF(x)可以表示为
如图3C所示，当重构的模型必须限制在某个隐式函数场的外部区域时，可以采用差集操作。由隐式函数场A和隐式函数场B的差集
所构成的组合约束可以表示为
其中，f-B(x)为隐式函数场B的补集，满足f-B(x)＝-fB(x)。
组合隐式函数场相应的模型外包的隐式函数BF(x)可以表示为
同理，可以构建隐式函数场B与隐式函数场A的差集，如图3D所示。
在一些实施例中，所述基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值，包括：
自下而上遍历所述操作树，依次计算所述操作树上各个节点的最佳有向包围盒；其中，叶子节点的最佳有向包围盒为叶子节点的隐式函数的模型外包，非叶子节点的最佳有向包围盒基于所述非叶子节点的子节点的组合模型外包确定；
自上而下遍历所述操作树，依次判断当前点是否在各个节点的最佳有向包围盒内，若否，则忽略相应节点及所述相应节点的子节点对应的隐式函数的影响；直至遍历完毕，得到所述操作树上对当前点有影响的所有叶子节点的隐式函数；
基于所述对当前点有影响的所有叶子节点的隐式函数，求取对所述当前点有影响的所有子隐式函数场的函数值。
实际应用中，如果不对评估过程进行加速处理，则组合隐式函数场中每个点的评估需要计算多个子隐式函数场的函数值，计算量较大。
由于相比于复杂函数的评估，判断点在模型外包的内外是比较高效的。考虑到组合隐式函数场只有在包围盒重叠的区域才需要考虑多个隐式函数场的影响，因此，可以通过判断点与模型外包的位置关系来加速隐式函数场的评估。
基于此，为了提高隐式曲面重构的效率，改进了多隐式函数场组合的评估方法。本发明实施例方法的基本思路是基于操作树的分层结构采用剪枝的方式，以减少不必要的计算。嵌套的组合约束可以构成一颗操作树。操作树的一个叶子节点与一个子隐式函数场对应，一个非叶子节点与一个子组合隐式函数场对应，根节点对应最终的组合隐式函数场。
在一应用示例中，给定一个任意评估点x，F(x)的加速评估方法可以描述如下：
步骤1：求解表征各个子隐式函数场的隐式函数。
步骤2：根据用户指定的嵌套组合规则得到一颗操作树。
步骤3：自下而上遍历操作树，依次计算各个节点的最佳有向包围盒。其中，叶子节点的最佳有向包围盒为子隐式函数场的模型外包，非叶子节点的最佳有向包围盒则可以根据子节点的组合模型外包来计算。
步骤4：自上而下遍历操作树，依次判断评估点x是否在各个节点的最佳有向包围盒内。如果评估点x不在某个节点Ni的最佳有向包围盒内，则不必考虑节点Ni及其所有子节点对应的隐式函数场的影响；否则，以同样的方式继续判断评估点x是否在节点Ni的子节点的最佳有向包围盒内。
步骤5：遍历完毕之后，可以得到对评估点x有影响的在叶子节点中的子隐式函数场{fi}。
步骤6：计算对评估点x有影响的所有子隐式函数场的函数值{fi(x)}。对于径向基类型的隐式函数，可以采用快速多极方法来加速评估速度。
步骤7：自下而上遍历操作树，根据组合规则确定对评估点x有影响的各个节点的函数值。
步骤8：最终根节点的函数值即为组合隐式函数场F(x)的函数值。
实际应用中，单纯依靠采样数据的插值过程有时难以反映真实的矿化分布趋势，在一应用示例中，根据成矿规律来限制空间插值的外推区域，对于形成和分布与一定岩层相关的层控矿床，可以采用组合隐式函数场来控制矿体与岩层间的赋存关系，使重构的矿体模型限制在在某个岩层内或在某个岩层外。如此，可以先对岩层进行建模。如图4A所示，可以分别构建第一隐式函数f1(x)、第二隐式函数f2(x)、第三隐式函数f3(x)及第四隐式函数f4(x)；其中，f1(x)表示几何边界线约束，用于控制模型的范围；f2(x)表示矿体上边界面，f3(x)表示矿体下边界面，f4(x)表示岩层界面。如果要求层状薄矿体的下边界面贴合在岩层界面上，则该示例的隐式函数场组合关系为：F(x)＝f1(x)∩(f2(x)∩f4(x))。组合约束前的矿体模型如图4B所示，组合约束前的矿体和岩层模型如图4C所示，组合约束后的矿体和岩层模型如图4D所示。
为了实现本发明实施例的方法，本发明实施例还提供一种矿体组合约束建模装置，设置在矿体组合约束建模设备，如图5所示，该装置包括：构建模块501和曲面确定模块502；其中，构建模块501用于基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；曲面确定模块502用于遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果。
曲面确定模块502求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值。
在一些实施例中，曲面确定模块502基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值，包括：
自下而上遍历所述操作树，依次计算所述操作树上各个节点的最佳有向包围盒；其中，叶子节点的最佳有向包围盒为叶子节点的隐式函数的模型外包，非叶子节点的最佳有向包围盒基于所述非叶子节点的子节点的组合模型外包确定；
自上而下遍历所述操作树，依次判断当前点是否在各个节点的最佳有向包围盒内，若否，则忽略相应节点及所述相应节点的子节点对应的隐式函数的影响；直至遍历完毕，得到所述操作树上对当前点有影响的所有叶子节点的隐式函数；
基于所述对当前点有影响的所有叶子节点的隐式函数，求取对所述当前点有影响的所有子隐式函数场的函数值。
在一些实施例中，曲面确定模块502基于主成分分析方法计算所述操作树上各个节点的最佳有向包围盒。
在一些实施例中，构建模块501具体用于：
所述基于矿体区域的地质规则构建第一隐式函数、第二隐式函数、第三隐式函数及第四隐式函数，其中，所述第一隐式函数表征几何边界线约束，用于控制模型的范围；所述第二隐式函数表征矿体上边界面；所述第三隐式函数表征矿体下边界面；所述第四隐式函数表征岩层界面。
实际应用时，构建模块501及曲面确定模块502，可以由矿体组合约束建模装置中的处理器来实现。当然，处理器需要运行存储器中的计算机程序来实现它的功能。
需要说明的是：上述实施例提供的矿体组合约束建模装置在进行矿体组合约束建模时，仅以上述各程序模块的划分进行举例说明，实际应用中，可以根据需要而将上述处理分配由不同的程序模块完成，即将装置的内部结构划分成不同的程序模块，以完成以上描述的全部或者部分处理。另外，上述实施例提供的矿体组合约束建模装置与矿体组合约束建模方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
基于上述程序模块的硬件实现，且为了实现本发明实施例的方法，本发明实施例还提供一种矿体组合约束建模设备。图6仅仅示出了该设备的示例性结构而非全部结构，根据需要可以实施图6示出的部分结构或全部结构。
如图6所示，本发明实施例提供的矿体组合约束建模设备600包括：至少一个处理器601、存储器602、用户接口603和至少一个网络接口604。矿体组合约束建模设备600中的各个组件通过总线系统605耦合在一起。可以理解，总线系统605用于实现这些组件之间的连接通信。总线系统605除包括数据总线之外，还包括电源总线、控制总线和状态信号总线。但是为了清楚说明起见，在图6中将各种总线都标为总线系统605。
其中，用户接口603可以包括显示器、键盘、鼠标、轨迹球、点击轮、按键、按钮、触感板或者触摸屏等。
本发明实施例中的存储器602用于存储各种类型的数据以支持矿体组合约束建模设备的操作。这些数据的示例包括：用于在矿体组合约束建模设备上操作的任何计算机程序。
本发明实施例揭示的矿体组合约束建模方法可以应用于处理器601中，或者由处理器601实现。处理器601可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，矿体组合约束建模方法的各步骤可以通过处理器601中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器601可以是通用处理器、数字信号处理器(DSP，DigitalSignal Processor)，或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。处理器601可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者任何常规的处理器等。结合本发明实施例所公开的方法的步骤，可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于存储介质中，该存储介质位于存储器602，处理器601读取存储器602中的信息，结合其硬件完成本发明实施例提供的矿体组合约束建模方法的步骤。
在示例性实施例中，矿体组合约束建模设备可以被一个或多个应用专用集成电路(ASIC，Application Specific Integrated Circuit)、DSP、可编程逻辑器件(PLD，Programmable Logic Device)、复杂可编程逻辑器件(CPLD，Complex Programmable LogicDevice)、FPGA、通用处理器、控制器、微控制器(MCU，Micro Controller Unit)、微处理器(Microprocessor)、或者其他电子元件实现，用于执行前述方法。
可以理解，存储器602可以是易失性存储器或非易失性存储器，也可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(ROM，Read Only Memory)、可编程只读存储器(PROM，Programmable Read-Only Memory)、可擦除可编程只读存储器(EPROM，Erasable Programmable Read-Only Memory)、电可擦除可编程只读存储器(EEPROM，Electrically Erasable Programmable Read-Only Memory)、磁性随机存取存储器(FRAM，ferromagnetic random access memory)、快闪存储器(Flash Memory)、磁表面存储器、光盘、或只读光盘(CD-ROM，Compact Disc Read-Only Memory)；磁表面存储器可以是磁盘存储器或磁带存储器。易失性存储器可以是随机存取存储器(RAM，Random AccessMemory)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(SRAM，Static Random Access Memory)、同步静态随机存取存储器(SSRAM，Synchronous Static Random Access Memory)、动态随机存取存储器(DRAM，Dynamic Random Access Memory)、同步动态随机存取存储器(SDRAM，SynchronousDynamic Random Access Memory)、双倍数据速率同步动态随机存取存储器(DDRSDRAM，Double Data Rate Synchronous Dynamic Random Access Memory)、增强型同步动态随机存取存储器(ESDRAM，Enhanced Synchronous Dynamic Random Access Memory)、同步连接动态随机存取存储器(SLDRAM，SyncLink Dynamic Random Access Memory)、直接内存总线随机存取存储器(DRRAM，Direct Rambus RandomAccess Memory)。本发明实施例描述的存储器旨在包括但不限于这些和任意其它适合类型的存储器。
在示例性实施例中，本发明实施例还提供了一种存储介质，即计算机存储介质，具体可以是计算机可读存储介质，例如包括存储计算机程序的存储器602，上述计算机程序可由矿体组合约束建模设备的处理器601执行，以完成本发明实施例方法所述的步骤。计算机可读存储介质可以是ROM、PROM、EPROM、EEPROM、Flash Memory、磁表面存储器、光盘、或CD-ROM等存储器。
需要说明的是：“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。
另外，本发明实施例所记载的技术方案之间，在不冲突的情况下，可以任意组合。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
Claims (8)
1.一种矿体组合约束建模方法，其特征在于，包括：
基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果；
其中，所述求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值；
其中，所述组合隐式函数场为对所述多个隐式函数基于隐式函数间的运算规则生成的，所述运算规则包括：交集、并集和差集，所述操作树的一个叶子节点与一个子隐式函数场对应，一个非叶子节点与一个子组合隐式函数场对应，根节点对应最终的组合隐式函数场；所述基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值，包括：
自下而上遍历所述操作树，依次计算所述操作树上各个节点的最佳有向包围盒；其中，叶子节点的最佳有向包围盒为叶子节点的隐式函数的模型外包，非叶子节点的最佳有向包围盒基于所述非叶子节点的子节点的组合模型外包确定；
自上而下遍历所述操作树，依次判断当前点是否在各个节点的最佳有向包围盒内，若否，则忽略相应节点及所述相应节点的子节点对应的隐式函数的影响；直至遍历完毕，得到所述操作树上对当前点有影响的所有叶子节点的隐式函数；
基于所述对当前点有影响的所有叶子节点的隐式函数，求取对所述当前点有影响的所有子隐式函数场的函数值。
2.根据权利要求1所述的方法，其特征在于，基于主成分分析方法计算所述操作树上各个节点的最佳有向包围盒。
3.根据权利要求1所述的方法，其特征在于，所述基于矿体区域的地质规则构建多个隐式函数，包括：
所述基于矿体区域的地质规则构建第一隐式函数、第二隐式函数、第三隐式函数及第四隐式函数，其中，所述第一隐式函数表征几何边界线约束，用于控制模型的范围；所述第二隐式函数表征矿体上边界面；所述第三隐式函数表征矿体下边界面；所述第四隐式函数表征岩层界面。
4.一种矿体组合约束建模装置，其特征在于，包括：
构建模块，用于基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
曲面确定模块，用于遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果；
所述曲面确定模块求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值；
所述组合隐式函数场为对所述多个隐式函数基于隐式函数间的运算规则生成的，所述运算规则包括：交集、并集和差集，所述操作树的一个叶子节点与一个子隐式函数场对应，一个非叶子节点与一个子组合隐式函数场对应，根节点对应最终的组合隐式函数场；所述曲面确定模块基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值，包括：
自下而上遍历所述操作树，依次计算所述操作树上各个节点的最佳有向包围盒；其中，叶子节点的最佳有向包围盒为叶子节点的隐式函数的模型外包，非叶子节点的最佳有向包围盒基于所述非叶子节点的子节点的组合模型外包确定；
自上而下遍历所述操作树，依次判断当前点是否在各个节点的最佳有向包围盒内，若否，则忽略相应节点及所述相应节点的子节点对应的隐式函数的影响；直至遍历完毕，得到所述操作树上对当前点有影响的所有叶子节点的隐式函数；
基于所述对当前点有影响的所有叶子节点的隐式函数，求取对所述当前点有影响的所有子隐式函数场的函数值。
5.根据权利要求4所述的装置，其特征在于，所述曲面确定模块基于主成分分析方法计算所述操作树上各个节点的最佳有向包围盒。
6.根据权利要求4所述的装置，其特征在于，所述构建模块具体用于：
所述基于矿体区域的地质规则构建第一隐式函数、第二隐式函数、第三隐式函数及第四隐式函数，其中，所述第一隐式函数表征几何边界线约束，用于控制模型的范围；所述第二隐式函数表征矿体上边界面；所述第三隐式函数表征矿体下边界面；所述第四隐式函数表征岩层界面。
7.一种矿体组合约束建模设备，其特征在于，包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器，其中，
所述处理器，用于运行计算机程序时，执行权利要求1至3任一项所述方法的步骤。
8.一种存储介质，所述存储介质上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现权利要求1至3任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010892365.5A CN112036030B (zh) | 2020-08-31 | 2020-08-31 | 矿体组合约束建模方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010892365.5A CN112036030B (zh) | 2020-08-31 | 2020-08-31 | 矿体组合约束建模方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112036030A CN112036030A (zh) | 2020-12-04 |
| CN112036030B true CN112036030B (zh) | 2021-06-01 |
Family
ID=73586700
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010892365.5A Active CN112036030B (zh) | 2020-08-31 | 2020-08-31 | 矿体组合约束建模方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112036030B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113011040B (zh) * | 2021-04-01 | 2022-07-08 | 中南大学 | 脉状地质体建模方法、装置、设备及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109326002A (zh) * | 2018-11-27 | 2019-02-12 | 中南大学 | 基于钻孔数据的矿体建模方法、装置、系统及存储介质 |
| CN109584357A (zh) * | 2018-11-27 | 2019-04-05 | 中南大学 | 基于多轮廓线的三维建模方法、装置、系统及存储介质 |
| CN110033519A (zh) * | 2019-04-23 | 2019-07-19 | 中南大学 | 基于隐式函数的三维建模方法、装置、系统及存储介质 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1797035A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 三维复杂地质体重构的技术 |
| CN106934860B (zh) * | 2017-01-12 | 2019-12-31 | 天津大学 | 一种基于t样条的三维地质建模方法 |
| EP3531172B1 (en) * | 2018-02-23 | 2022-08-17 | Services Pétroliers Schlumberger | Meshless and mesh-based technique for modeling subterranean volumes |
-
2020
- 2020-08-31 CN CN202010892365.5A patent/CN112036030B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109326002A (zh) * | 2018-11-27 | 2019-02-12 | 中南大学 | 基于钻孔数据的矿体建模方法、装置、系统及存储介质 |
| CN109584357A (zh) * | 2018-11-27 | 2019-04-05 | 中南大学 | 基于多轮廓线的三维建模方法、装置、系统及存储介质 |
| CN110033519A (zh) * | 2019-04-23 | 2019-07-19 | 中南大学 | 基于隐式函数的三维建模方法、装置、系统及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112036030A (zh) | 2020-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11921491B2 (en) | Conversion of mesh geometry to watertight boundary representation | |
| Kim et al. | Isogeometric analysis for trimmed CAD surfaces | |
| Ramanathan et al. | Constructing medial axis transform of planar domains with curved boundaries | |
| Shephard et al. | Automatic three‐dimensional mesh generation by the finite octree technique | |
| CN109326002B (zh) | 基于钻孔数据的矿体建模方法、装置、系统及存储介质 | |
| Hu et al. | Robust interval solid modelling part II: boundary evaluation | |
| JP2010049677A (ja) | 案内曲線で定義される面の設計をなすコンピュータ実行方法 | |
| US20230177224A1 (en) | Computer aided generative design with feature thickness control to facilitate manufacturing and structural performance | |
| CN112036030B (zh) | 矿体组合约束建模方法、装置、设备及存储介质 | |
| CN110349262B (zh) | 一种基于bim的隧道超前地质预报三维建模方法及系统 | |
| Chen et al. | A graph-and feature-based building space recognition algorithm for performance simulation in the early design stage | |
| Huhnt | Reconstruction of edges in digital building models | |
| Takahashi et al. | A three‐dimensional shape optimization for transient acoustic scattering problems using the time‐domain boundary element method | |
| Asghari et al. | Advances in techniques to formulate the watertight concept for cadastre | |
| Potter et al. | Ordered line integral methods for solving the eikonal equation | |
| CN115239901A (zh) | 地质体模型构建方法、装置、设备及存储介质 | |
| US10331803B2 (en) | Computer aided modeling | |
| CN115510689A (zh) | 脉状地质体建模方法、装置、设备及存储介质 | |
| CN114357683B (zh) | 基于交叉轮廓线法向动态估计的矿体建模方法及装置 | |
| CN115267667A (zh) | 井下高精度定位修正方法、装置、设备及存储介质 | |
| Fries | Higher-order accurate integration for cut elements with Chen-Babuška nodes | |
| CN112927350B (zh) | 多域地质体模型构建方法、装置、设备及存储介质 | |
| Robinson et al. | Strategies for adding features to CAD models in order to optimize performance | |
| JP2008533614A (ja) | 整合輪郭プロファイルを生成するシステム及び方法 | |
| CN113011040B (zh) | 脉状地质体建模方法、装置、设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |