Finite Element Implicit 3D Subsurface Structural Modeling

Modeste Irakarama and Morgan Thierry-Coudon and Mustapha Zakari and Guillaume Caumon. ( 2022 )
in: Computer-Aided Design, 149 (103267)

Abstract

We introduce a method for 3D implicit geological structural modeling from sparse sample points, where several conformable geological surfaces are represented by one single scalar field. Laplacian and Hessian regularization energies are discretized on a tetrahedral mesh using finite elements. This scheme is believed to offer some geometrical flexibility as it is readily implemented on both structured and unstructured grids. While implicit modeling on unstructured grids is not new, methods based on finite elements have received little attention. The finite element method is routinely used to solve boundary value problems. However, because boundary conditions are typically unknown in implicit subsurface structural modeling, the traditional finite element method requires some adjustments. To this end, we present boundary free discretizations of the Laplacian and Hessian energies that do not assume vanishing Neumann boundary conditions, thereby eliminating the boundary artifacts usually associated with that assumption. Furthermore, we argue that while an appropriate discretization of the Laplacian can be used to minimize the curvature of a function on triangulated meshes, it may fail to do so on tetrahedral meshes.

Download / Links

BibTeX Reference

@article{irakarama:hal-03654609,
 abstract = {We introduce a method for 3D implicit geological structural modeling from sparse sample points, where several conformable geological surfaces are represented by one single scalar field. Laplacian and Hessian regularization energies are discretized on a tetrahedral mesh using finite elements. This scheme is believed to offer some geometrical flexibility as it is readily implemented on both structured and unstructured grids. While implicit modeling on unstructured grids is not new, methods based on finite elements have received little attention. The finite element method is routinely used to solve boundary value problems. However, because boundary conditions are typically unknown in implicit subsurface structural modeling, the traditional finite element method requires some adjustments. To this end, we present boundary free discretizations of the Laplacian and Hessian energies that do not assume vanishing Neumann boundary conditions, thereby eliminating the boundary artifacts usually associated with that assumption. Furthermore, we argue that while an appropriate discretization of the Laplacian can be used to minimize the curvature of a function on triangulated meshes, it may fail to do so on tetrahedral meshes.},
 author = {Irakarama, Modeste and Thierry-Coudon, Morgan and Zakari, Mustapha and Caumon, Guillaume},
 doi = {10.1016/j.cad.2022.103267},
 hal_id = {hal-03654609},
 hal_version = {v1},
 journal = {{Computer-Aided Design}},
 keywords = {Subsurface modeling ; Implicit modeling ; Data interpolation ; Finite elements ; Laplacian energy ; Hessian energy},
 month = {August},
 pages = {103267},
 pdf = {https://hal.univ-lorraine.fr/hal-03654609/file/Irakarama_et_al.pdf},
 publisher = {{Elsevier}},
 title = {{Finite Element Implicit 3D Subsurface Structural Modeling}},
 url = {https://hal.univ-lorraine.fr/hal-03654609},
 volume = {149},
 year = {2022}
}