一种基于BRep模型的多面域地层有限元网格生成方法

技术特征：
1.一种基于brep模型的多面域地层有限元网格生成方法，其特征在于：包括以下步骤：步骤1：建立露天矿采场、地层界面和断层界面的三角形网格；步骤2：建立以采场为表面的模型范围网格，具体方法为：步骤2.1：按照给定的建模范围的长、宽和高尺寸，建立模型范围的四边形网格；步骤2.2：采场界面约束模型范围网格，形成以采场为表面的模型范围网格集m
′
；步骤3：地层界面与断层界面之间进行相互约束，具体方法为：步骤3.1：对地层界面进行约束，得到有效的地层面三角形集l'；步骤3.2：将约束边p按delaunay准则加入到断层界面三角形网格f，形成断层面f'；步骤4：建立多面域地层空间brep模型，具体方法为：步骤4.1：使用由露天矿采场面构建的范围网格m'对地层界面l'进行多面域brep模型建立；步骤4.2：使用由露天矿采场面构建的外围网格m”对断层界面f'进行多面域brep模型建立；步骤5：多面域地层空间brep模型四面体剖分生成有限元网格。2.根据权利要求1所述的一种基于brep模型的多面域地层有限元网格生成方法，其特征在于：所述步骤1的具体方法为：对露天矿采场台阶线进行约束delaunay三角形剖分，生成露天矿采场面；以钻孔岩性分布数据为样本，使用距离平方反比法对所需要的岩性数据进行插值，分别生成各地层界面和断层界面的三角形网格。得到露天矿采场面的三角形集t＝{t1,t2,...,t
i
,...,t
n
}，其中，t
i
为露天矿采场面的第i个三角形，i∈[1,n]，n为露天矿采场面的三角形总数；得到地层界面的三角形集l＝{l
11
,l
12
,...,l
re
,...,l
oq
}，其中，l
re
为第r层地层界面上的第e个三角形，q为第o层地层界面上的三角形总数；得到断层界面的三角形集f＝{f
11
,f
12
,...,f
hg
,...,f
vb
}，其中，f
hg
为第h层断层界面上的第g个三角形，b为第v层断层界面上的三角形总数。3.根据权利要求1所述的一种基于brep模型的多面域地层有限元网格生成方法，其特征在于：所述步骤2的具体方法为：步骤2.1：按照给定的建模范围的长、宽和高尺寸，建立模型范围的四边形网格；进一步地，得到四边形集m＝{m1,m2,...,m
x
,...,m
z
}，其中，m
x
为模型范围四边形网格上的第x个四边形，x∈[1,z]，z为模型范围四边形的总数。步骤2.2：采场界面约束模型范围网格，形成以采场为表面的模型范围网格集m
′
，具体方法为：步骤2.2.1：用露天矿采场界面与模型范围网格求交，得到交线p＝{p1,p2,...,p
y
,...,p
u
}，作为约束边，其中，p
y
为约束线上的第y个点，y∈[1,u]，u为约束线上点的总数；步骤2.2.2：将约束边p按delaunay准则分别加入到采场界面三角形集t和模型范围网格集m，利用约束边为界删除多余的采场界面三角形和模型范围网格四边形，保留有效的采场三角形和模型范围四边形，并进行组合，形成以采场为表面的模型范围网格集m
′
。4.根据权利要求1所述的一种基于brep模型的多面域地层有限元网格生成方法，其特征在于：所述步骤3的具体方法为：步骤3.1：对地层界面进行约束，得到有效的地层面三角形集l'；
进一步地，两个相互约束的地层界面三角形网格l和断层界面三角形网格f求交，得地层界面与断层界面之间的交线，将该交线作为地层界面的约束边p，将约束边p按delaunay准则加入到地层界面三角形集l，并以约束边p为界删除地层界面多余三角形，保留有效地层面三角形集l'。步骤3.2：将约束边p按delaunay准则加入到断层界面三角形网格f，形成断层面f'；进一步地，按照各个地层的生成机理来确定地层间的拓扑关系，根据地层界面顺序对以上方法进行重复操作，直到完成所有地层界面与断层界面之间的求交并约束断层界面三角形网格，即实现了实现地层界面与断层界面之间的约束。5.根据权利要求1所述的一种基于brep模型的多面域地层有限元网格生成方法，其特征在于：所述步骤4的具体方法为：步骤4.1：使用由露天矿采场面构建的范围网格m'对地层界面l'进行多面域brep模型建立，具体建立方法为：步骤4.1.1：根据建立两个相互约束的地层界面l'和模型范围网格m
′
求得地层界面与模型范围网格之间的交线，将该交线作为地层界面的约束线p按delaunay准则加入到地层界面三角形网格，以约束线p为界删除地层界面多余三角形网格，保留有效地层面l”；步骤4.1.2：用地层界面l'对外围网格m'进行约束操作，将约束后的外围网格进行组合，形成外围网格m'
″
，根据地层界面顺序对以上方法进行重复操作，直到完成所有地层界面与范围网格之间的约束操作，即实现实现地层界面与范围网格之间的约束。步骤4.2：使用由露天矿采场面构建的外围网格m”对断层界面f'进行多面域brep模型建立，具体建立方法为：步骤4.2.1：求得断层界面f'与外围网格m”之间的交线，将该交线作为断层界面的约束线p按delaunay准则加入到断层界面三角形网格，以约束线p为界删除断层界面多余三角形网格，保留有效断层面f”；步骤4.2.2：用断层界面f'对外围网格m”进行约束操作，将约束后的外围网格进行组合，形成外围网格m
”′
，根据顺序对以上方法进行重复操作，直到完成所有断层界面与外围网格之间的约束操作，即完成建立多面域地层空间brep模型。6.根据权利要求1所述的一种基于brep模型的多面域地层有限元网格生成方法，其特征在于：所述步骤5的具体方法为：进一步地，对建立的多面域地层空间brep模型进行delaunay不规则四面体剖分，实现基于brep模型的多面域地层空间有限元网格生成。

技术总结
本发明提供一种基于BRep模型的多面域地层有限元网格生成方法，涉及三维地质建模技术领域。主要研究内容包括：对建模区域的钻孔数据进行处理，以Delaunay三角剖分为基础，生成露天矿采场面，结合距离幂次反比法对地层界面以及断层界面等进行插值生成；根据建模所需范围建立模型范围网格；对地层界面与断层界面之间进行处理，建立地层界面与断层界面之间的BRep模型；然后，建立模型范围网格与地层界面以及断层界面之间的BRep模型，即为多面域地层的BRep模型；对多面域地层空间进行不规则四面体剖分，实现基于BRep模型的多面域地层空间有限元网格生成。结果表明，该方法能精确地生成多面域地层空间的有限元网格，有效提高工作效率。率。率。


技术研发人员：陈应显 王鹏飞 李广贺 王东 贾兰 姜聚宇 张立国
受保护的技术使用者：辽宁工程技术大学
技术研发日：2022.08.22
技术公布日：2022/11/18