Implicit {Tetrahedral} {Mesh} {Updating} for {Subsurface} {Modeling}.
Ayoub Belhachmi and Guillaume Caumon and Charles Dapogny. ( 2026 )
in: 2026 {RING} meeting, pages 293--315, ASGA
Abstract
Generating a valid, high-quality subsurface mesh, that conforms to geological interfaces is a challenging task because of the complex geometrical configurations encountered in geosciences. Moreover, uncertainties regarding location, geometry, or even the existence of geological discontinuities raise the need for fast, robust local mesh updating techniques. Such methods allow the mesh to be updated efficiently when new discontinuities or interfaces are introduced, or when existing ones are modified, without remeshing the entire domain from scratch. The aim of this article is to introduce a novel mesh-updating method that enables the insertion of finite surfaces such as faults in existing multi-material tetrahedral meshes. In this approach, the surface to be inserted is implicitly represented using two level set functions: it is defined as the intersection of the zero-level set of one function and the negative subdomain of the other. These functions can be obtained by interpolation, subject to a smoothness criterion. From this input, we perform an implicit surface insertion strategy based on the open-source MMG remeshing library. The proposed level set representation also allows fault interactions to be controlled, so that one fault can branch onto another and form a Y-type configuration rather than an undesired X-type crossing. After insertion, the same representation can be used to duplicate nodes along the fault surfaces, introducing displacement discontinuities and enabling horizon updates with fault offsets, which open the way for incremental structural modeling.
Download / Links
BibTeX Reference
@inproceedings{Belhachmi2026RM,
abstract = {Generating a valid, high-quality subsurface mesh, that conforms to geological interfaces is a challenging task because of the complex geometrical configurations encountered in geosciences. Moreover, uncertainties regarding location, geometry, or even the existence of geological discontinuities raise the need for fast, robust local mesh updating techniques. Such methods allow the mesh to be updated efficiently when new discontinuities or interfaces are introduced, or when existing ones are modified, without remeshing the entire domain from scratch. The aim of this article is to introduce a novel mesh-updating method that enables the insertion of finite surfaces such as faults in existing multi-material tetrahedral meshes. In this approach, the surface to be inserted is implicitly represented using two level set functions: it is defined as the intersection of the zero-level set of one function and the negative subdomain of the other. These functions can be obtained by interpolation, subject to a smoothness criterion. From this input, we perform an implicit surface insertion strategy based on the open-source MMG remeshing library. The proposed level set representation also allows fault interactions to be controlled, so that one fault can branch onto another and form a Y-type configuration rather than an undesired X-type crossing. After insertion, the same representation can be used to duplicate nodes along the fault surfaces, introducing displacement discontinuities and enabling horizon updates with fault offsets, which open the way for incremental structural modeling.},
author = {Belhachmi, Ayoub and Caumon, Guillaume and Dapogny, Charles},
booktitle = {2026 {RING} meeting},
language = {en},
pages = {293--315},
publisher = {ASGA},
title = {Implicit {Tetrahedral} {Mesh} {Updating} for {Subsurface} {Modeling}.},
year = {2026}
}
