IMPLICIT MODELLING OF GEOLOGICAL STRUCTURES: A CARTESIAN GRID METHOD HANDLING DISCONTINUITIES WITH GHOST POINTS

Julien Renaudeau and Modeste Irakarama and Gautier Laurent and Frantz Maerten and Guillaume Caumon. ( 2018 )
in: BEM/MRM 41, pages 189 - 199

Abstract

In geology, implicit structural modelling constructs the geometry of geological structures (e.g. layers) by interpolating between sparse field data. A model is represented by a volumetric scalar field which is discontinuous on structural discontinuities such as faults or stratigraphic unconformities. The management of such discontinuities may involve boolean operations on several scalar fields or the creation of conformal meshes. Instead, we propose a ghost cell technique for the cartesian grid together with a set of relations between the ghost points, the regular nodes, and the discontinuities. Consequently, poor quality meshes are avoided and only the resolution of the grid has an impact on the solution. The modelling problem is posed as a least square’s minimization of a bending energy penalization on data mismatch functions and approximated by finite difference. As all relationships in the grid are implicitly defined, except close to the discontinuities, this algorithm is computationally efficient. We provide some benchmarks of the method on two-dimensional examples with folds, faults, and erosions.

Download / Links

BibTeX Reference

@inproceedings{renaudeau:hal-02136480,
 abstract = {In geology, implicit structural modelling constructs the geometry of geological structures (e.g. layers) by interpolating between sparse field data. A model is represented by a volumetric scalar field which is discontinuous on structural discontinuities such as faults or stratigraphic unconformities. The management of such discontinuities may involve boolean operations on several scalar fields or the creation of conformal meshes. Instead, we propose a ghost cell technique for the cartesian grid together with a set of relations between the ghost points, the regular nodes, and the discontinuities. Consequently, poor quality meshes are avoided and only the resolution of the grid has an impact on the solution. The modelling problem is posed as a least square’s minimization of a bending energy penalization on data mismatch functions and approximated by finite difference. As all relationships in the grid are implicitly defined, except close to the discontinuities, this algorithm is computationally efficient. We provide some benchmarks of the method on two-dimensional examples with folds, faults, and erosions.},
 address = {New Forest, United Kingdom},
 author = {Renaudeau, Julien and Irakarama, Modeste and Laurent, Gautier and Maerten, Frantz and Caumon, Guillaume},
 booktitle = {{BEM/MRM 41}},
 doi = {10.2495/BE410171},
 hal_id = {hal-02136480},
 hal_version = {v1},
 keywords = {energy penalization ; ghost points ; mesh reduction methods ; 2D structural modelling},
 month = {September},
 number = {11},
 pages = {189 - 199},
 series = {Boundary Elements and other Mesh Reduction Methods XXXXI},
 title = {{IMPLICIT MODELLING OF GEOLOGICAL STRUCTURES: A CARTESIAN GRID METHOD HANDLING DISCONTINUITIES WITH GHOST POINTS}},
 url = {https://hal.univ-lorraine.fr/hal-02136480},
 volume = {122},
 year = {2018}
}