A Method for Converting Triangle-Mesh-Based 3D Geological Models into Hexahedral Grids for Parallel Finite Difference Simulation.

Björn Zehner and O. Hellwig and M. Linke and I. Görz and S. Buske. ( 2015 )
in: 35th Gocad Meeting - 2015 RING Meeting, ASGA

Abstract

The 3D geological underground is often modeled using triangulated networks representing stratigraphic interfaces and structures, such as faults. If these models are to be simulated using the finite difference method, they must be converted into a representation that tessellates the full volume of the model into hexahedral cells. Often the simulations require a high grid resolution and are done using parallel computing. Then a large amount of storage space is required and it is difficult to create such a model using the standard geo-modeling packages. Since the raster representation is only required for the calculation, but not for the geometry description, we present an algorithm and concept for rasterizing geological models on the fly for the use in finite difference codes that are parallelized by domain decomposition. We implemented a C++ library and integrated it into seismic simulation software that is run as parallel code on a UNIX cluster using MPI. We can thus run the simulation with realistic and complicated surface-based geological models that are created using 3D geo-modeling software.

Download / Links

BibTeX Reference

@INPROCEEDINGS{Zehner2 GM2015,
    author = { Zehner, Björn and Hellwig, O. and Linke, M. and Görz, I. and Buske, S. },
     title = { A Method for Converting Triangle-Mesh-Based 3D Geological Models into Hexahedral Grids for Parallel Finite Difference Simulation. },
 booktitle = { 35th Gocad Meeting - 2015 RING Meeting },
      year = { 2015 },
 publisher = { ASGA },
  abstract = { The 3D geological underground is often modeled using triangulated networks representing stratigraphic interfaces and structures, such as faults. If these models are to be simulated using the finite difference method, they must be converted into a representation that tessellates the full volume of the model into hexahedral cells. Often the simulations require a high grid resolution and are done using parallel computing. Then a large amount of storage space is required and it is difficult to create such a model using the standard geo-modeling packages. Since the raster representation is only required for the calculation, but not for the geometry description, we present an algorithm and concept for rasterizing geological models on the fly for the use in finite difference codes that are parallelized by domain decomposition. We implemented a C++ library and integrated it into seismic simulation software that is run as parallel code on a UNIX cluster using MPI. We can thus run the simulation with realistic and complicated surface-based geological models that are created using 3D geo-modeling software. }
}