Finite {{Difference Implicit Structural Modeling}} of {{Geological Structures}}

Modeste Irakarama and Gautier Laurent and Julien Renaudeau and Guillaume Caumon. ( 2020 )
in: Mathematical Geosciences

Abstract

We introduce a new method for implicit structural modeling. The main developments in this paper are the new regularization operators we propose by extending inherent properties of the classic one-dimensional discrete second derivative operator to higher dimensions. The proposed regularization operators discretize naturally on the Cartesian grid using finite differences, owing to the highly symmetric nature of the Cartesian grid. Furthermore, the proposed regularization operators do not require any special treatment on boundary nodes, and their generalization to higher dimensions is straightforward. As a result, the proposed method has the advantage of being simple to implement. Numerical examples show that the proposed method is robust and numerically efficient.

Download / Links

BibTeX Reference

@ARTICLE{Irakarama2020MG,
    author = { Irakarama, Modeste and Laurent, Gautier and Renaudeau, Julien and Caumon, Guillaume },
     title = { Finite {{Difference Implicit Structural Modeling}} of {{Geological Structures}} },
     month = { "sep" },
   journal = { Mathematical Geosciences },
      year = { 2020 },
      issn = { 1874-8953 },
       doi = { 10.1007/s11004-020-09887-w },
  abstract = { We introduce a new method for implicit structural modeling. The main developments in this paper are the new regularization operators we propose by extending inherent properties of the classic one-dimensional discrete second derivative operator to higher dimensions. The proposed regularization operators discretize naturally on the Cartesian grid using finite differences, owing to the highly symmetric nature of the Cartesian grid. Furthermore, the proposed regularization operators do not require any special treatment on boundary nodes, and their generalization to higher dimensions is straightforward. As a result, the proposed method has the advantage of being simple to implement. Numerical examples show that the proposed method is robust and numerically efficient. }
}