CN112036030A - 矿体组合约束建模方法、装置、设备及存储介质 - Google Patents

矿体组合约束建模方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN112036030A
CN112036030A CN202010892365.5A CN202010892365A CN112036030A CN 112036030 A CN112036030 A CN 112036030A CN 202010892365 A CN202010892365 A CN 202010892365A CN 112036030 A CN112036030 A CN 112036030A
Authority
CN
China
Prior art keywords
implicit
function
ore body
implicit function
operation tree
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.)
Granted
Application number
CN202010892365.5A
Other languages
English (en)
Other versions
CN112036030B (zh
Inventor
钟德云
王李管
贾明滔
毕林
王晋淼
陈鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN202010892365.5A priority Critical patent/CN112036030B/zh
Publication of CN112036030A publication Critical patent/CN112036030A/zh
Application granted granted Critical
Publication of CN112036030B publication Critical patent/CN112036030B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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)可以表示为
Figure BDA0002657226630000041
其中，fi(x)表示第i个子隐式函数，符号
Figure BDA0002657226630000051
表示隐式函数间的某种运算规则。不同的隐式函数fi(x)用于表示不同的组合约束。组合隐式函数场在本质上仍然是一个隐式函数。建模的目标仍然是把组合隐式函数场的隐式函数F(x)重构出来。
组合隐式函数场的各个子隐式函数场可以采用不同的隐式函数来表达，比如，fi(x)，i＝1，2，...，n可以采用不同类型的径向基类插值方法来表示。
实际应用中，为了限定模型空间中隐式函数间相互作用的范围，需要为各个子隐式函数指定模型外包。模型外包也可以看作一种隐式函数。
子隐式函数场i的模型外包的隐式函数可以表示为
Figure BDA0002657226630000055
当
Figure BDA0002657226630000053
时，子隐式函数场i的组合约束有效；当
Figure BDA0002657226630000054
时，子隐式函数场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)可以表示为
Figure BDA0002657226630000061
其中，
Figure BDA0002657226630000065
为表示隐式函数场A的模型外包的隐式函数，
Figure BDA0002657226630000066
表示隐式函数场B的模型外包的隐式函数。
如图3B所示，当重构的模型应该包含某个隐式函数场的内部区域时，可以采用并集操作。由隐式函数场A和隐式函数场B的并集A∪B所构成的组合约束可以表示为
F(x)＝fA∪B(x)＝min(fA(x)，fB(x))
组合隐式函数场相应的模型外包的隐式函数BF(x)可以表示为
Figure BDA0002657226630000062
如图3C所示，当重构的模型必须限制在某个隐式函数场的外部区域时，可以采用差集操作。由隐式函数场A和隐式函数场B的差集
Figure BDA0002657226630000067
所构成的组合约束可以表示为
Figure BDA0002657226630000063
其中，f-B(x)为隐式函数场B的补集，满足f-B(x)＝-fB(x)。
组合隐式函数场相应的模型外包的隐式函数BF(x)可以表示为
Figure BDA0002657226630000064
同理，可以构建隐式函数场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 (10)

1.一种矿体组合约束建模方法，其特征在于，包括：
基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果；
其中，所述求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值。
2.根据权利要求1所述的方法，其特征在于，所述基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值，包括：
自下而上遍历所述操作树，依次计算所述操作树上各个节点的最佳有向包围盒；其中，叶子节点的最佳有向包围盒为叶子节点的隐式函数的模型外包，非叶子节点的最佳有向包围盒基于所述非叶子节点的子节点的组合模型外包确定；
自上而下遍历所述操作树，依次判断当前点是否在各个节点的最佳有向包围盒内，若否，则忽略相应节点及所述相应节点的子节点对应的隐式函数的影响；直至遍历完毕，得到所述操作树上对当前点有影响的所有叶子节点的隐式函数；
基于所述对当前点有影响的所有叶子节点的隐式函数，求取对所述当前点有影响的所有子隐式函数场的函数值。
3.根据权利要求2所述的方法，其特征在于，基于主成分分析方法计算所述操作树上各个节点的最佳有向包围盒。
4.根据权利要求1所述的方法，其特征在于，所述基于矿体区域的地质规则构建多个隐式函数，包括：
所述基于矿体区域的地质规则构建第一隐式函数、第二隐式函数、第三隐式函数及第四隐式函数，其中，所述第一隐式函数表征几何边界线约束，用于控制模型的范围；所述第二隐式函数表征矿体上边界面；所述第三隐式函数表征矿体下边界面；所述第四隐式函数表征岩层界面。
5.一种矿体组合约束建模装置，其特征在于，包括：
构建模块，用于基于矿体区域的地质规则构建多个隐式函数，并对所述多个隐式函数进行组合，得到表征所述矿体区域的矿体与岩层间的赋存关系的组合隐式函数场；
曲面确定模块，用于遍历所述矿体区域内各点，求取各点在所述组合隐式函数场中的函数值，得到所述矿体区域的建模结果；
所述曲面确定模块求取各点在所述组合隐式函数场中的函数值，包括：
基于所述组合隐式函数场生成操作树；
基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值；
自下而上遍历所述操作树，确定所述操作树上对所述当前点有影响的各个节点的函数值；
基于所述对所述当前点有影响的各个节点的函数值，求取所述操作树的根节点的函数值，作为所述当前点在所述组合隐式函数场中的函数值。
6.根据权利要求5所述的装置，其特征在于，所述曲面确定模块基于所述操作树求取对当前点有影响的所有子隐式函数场的函数值，包括：
自下而上遍历所述操作树，依次计算所述操作树上各个节点的最佳有向包围盒；其中，叶子节点的最佳有向包围盒为叶子节点的隐式函数的模型外包，非叶子节点的最佳有向包围盒基于所述非叶子节点的子节点的组合模型外包确定；
自上而下遍历所述操作树，依次判断当前点是否在各个节点的最佳有向包围盒内，若否，则忽略相应节点及所述相应节点的子节点对应的隐式函数的影响；直至遍历完毕，得到所述操作树上对当前点有影响的所有叶子节点的隐式函数；
基于所述对当前点有影响的所有叶子节点的隐式函数，求取对所述当前点有影响的所有子隐式函数场的函数值。
7.根据权利要求6所述的装置，其特征在于，所述曲面确定模块基于主成分分析方法计算所述操作树上各个节点的最佳有向包围盒。
8.根据权利要求5所述的装置，其特征在于，所述构建模块具体用于：
所述基于矿体区域的地质规则构建第一隐式函数、第二隐式函数、第三隐式函数及第四隐式函数，其中，所述第一隐式函数表征几何边界线约束，用于控制模型的范围；所述第二隐式函数表征矿体上边界面；所述第三隐式函数表征矿体下边界面；所述第四隐式函数表征岩层界面。
9.一种矿体组合约束建模设备，其特征在于，包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器，其中，
所述处理器，用于运行计算机程序时，执行权利要求1至4任一项所述方法的步骤。
10.一种存储介质，所述存储介质上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现权利要求1至4任一项所述方法的步骤。
CN202010892365.5A 2020-08-31 2020-08-31 矿体组合约束建模方法、装置、设备及存储介质 Active CN112036030B (zh)

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 true CN112036030A (zh) 2020-12-04
CN112036030B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113011040A (zh) * 2021-04-01 2021-06-22 中南大学 脉状地质体建模方法、装置、设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1797035A (zh) * 2004-12-29 2006-07-05 中国石油天然气集团公司 三维复杂地质体重构的技术
CN106934860A (zh) * 2017-01-12 2017-07-07 天津大学 一种基于t样条的三维地质建模方法
CN109326002A (zh) * 2018-11-27 2019-02-12 中南大学 基于钻孔数据的矿体建模方法、装置、系统及存储介质
CN109584357A (zh) * 2018-11-27 2019-04-05 中南大学 基于多轮廓线的三维建模方法、装置、系统及存储介质
CN110033519A (zh) * 2019-04-23 2019-07-19 中南大学 基于隐式函数的三维建模方法、装置、系统及存储介质
US20190265386A1 (en) * 2018-02-23 2019-08-29 Schlumberger Technology Corporation Meshless and Mesh-Based Technique for Modeling Subterranean Volumes

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1797035A (zh) * 2004-12-29 2006-07-05 中国石油天然气集团公司 三维复杂地质体重构的技术
CN106934860A (zh) * 2017-01-12 2017-07-07 天津大学 一种基于t样条的三维地质建模方法
US20190265386A1 (en) * 2018-02-23 2019-08-29 Schlumberger Technology Corporation Meshless and Mesh-Based Technique for Modeling Subterranean Volumes
CN109326002A (zh) * 2018-11-27 2019-02-12 中南大学 基于钻孔数据的矿体建模方法、装置、系统及存储介质
CN109584357A (zh) * 2018-11-27 2019-04-05 中南大学 基于多轮廓线的三维建模方法、装置、系统及存储介质
CN110033519A (zh) * 2019-04-23 2019-07-19 中南大学 基于隐式函数的三维建模方法、装置、系统及存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113011040A (zh) * 2021-04-01 2021-06-22 中南大学 脉状地质体建模方法、装置、设备及存储介质
CN113011040B (zh) * 2021-04-01 2022-07-08 中南大学 脉状地质体建模方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN112036030B (zh) 2021-06-01

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
CN109326002B (zh) 基于钻孔数据的矿体建模方法、装置、系统及存储介质
Ramanathan et al. Constructing medial axis transform of planar domains with curved boundaries
Hu et al. Robust interval solid modelling part II: boundary evaluation
JP2010049677A (ja) 案内曲線で定義される面の設計をなすコンピュータ実行方法
CN105469446A (zh) 点云网格简化系统及方法
Razafindrazaka et al. Perfect matching quad layouts for manifold meshes
CN112036030B (zh) 矿体组合约束建模方法、装置、设备及存储介质
Zhu et al. Analytical solutions for sketch-based convolution surface modeling on the GPU
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
Chen et al. A graph-and feature-based building space recognition algorithm for performance simulation in the early design stage
Canelas et al. A new method for reliability analysis and reliability-based design optimization
CN115267667B (zh) 井下高精度定位修正方法、装置、设备及存储介质
CN115239901A (zh) 地质体模型构建方法、装置、设备及存储介质
US20230177224A1 (en) Computer aided generative design with feature thickness control to facilitate manufacturing and structural performance
US10331803B2 (en) Computer aided modeling
Potter et al. Ordered line integral methods for solving the eikonal equation
Fries Higher-order accurate integration for cut elements with Chen-Babuška nodes
CN112927350B (zh) 多域地质体模型构建方法、装置、设备及存储介质
JP2008533614A (ja) 整合輪郭プロファイルを生成するシステム及び方法
CN113011040B (zh) 脉状地质体建模方法、装置、设备及存储介质
CN114357683B (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