Interactive editing of 3D geological structures and tectonic history sketching via a rigid element method

Numerical models of geological structures are generally built with a geometrical approach, which lacks an explicit representation of the deformation history and may lead to incompatible structures. We advocate that the deformation history should be investigated and represented from the very first steps of the modelling process, provided that a series of rapid, interactive or automated, deformation tools are available for local editing, forward modelling and restoration. In this paper, we define the specifications of such tools and emphasise the need for rapidity and robustness. We briefly review the different applications of deformation tools in geomodelling and the existing deformation algorithms. We select a deformation algorithm based on rigid elements, first presented in the Computer Graphics community, which we refer to as Reed. It is able to rapidly deform any kind of geometrical object, including points, lines or volumes, with an approximated mechanical behaviour. The objects to be deformed are embedded in rigid cells whose displacement is optimised by minimising a global cost function with respect to displacement boundary conditions. This cost function measures the difference in displacement between neighbouring elements. The embedded objects are then deformed based on their original position with respect to the rigid elements. We present the basis of our implementation of this algorithm and highlight its ability to fulfil the specifications we defined. Its application to geomodelling specific problems is illustrated through the construction of a synthetic structural model of multiply deformed layers with a forward modelling approach. A special boundary condition adapted to restore large folds is also presented and applied to the large anticline of Han-sur-Lesse, Belgium, which demonstrates the ability of this method to efficiently perform a volumetric restoration without global projections. Graphical abstract(A) Initial structural model (Annot, France) adapted from Salles 2010, (b) Rigid Element based deformation tool defining a large folding event and (c) Applied deformation.Display Omitted HighlightsWe propose specifications for deformation tools for 3D geological structure modelling.We select a deformation algorithm based on rigid elements and deformation energy.We propose a process for modelling structures with several folding events.We implement boundary conditions for geological structure restoration.We implement anisotropy and local volume conservation rules.

[1]  N. Cardozo Trishear modeling of fold bedding data along a topographic profile , 2005 .

[2]  Hans-Peter Seidel,et al.  Vector field based shape deformations , 2006, SIGGRAPH 2006.

[3]  Markus H. Gross,et al.  Adaptive Space Deformations Based on Rigid Cells , 2007, Comput. Graph. Forum.

[4]  F. Maerten,et al.  Chronologic modeling of faulted and fractured reservoirs using geomechanically based restoration: Technique and industry applications , 2006 .

[5]  Ming C. Lin,et al.  Fast volume-preserving free-form deformation using multi-level optimization , 2000, Comput. Aided Des..

[6]  Ivan Viola,et al.  Modeling Terrains and Subsurface Geology , 2013, Eurographics.

[7]  Guillaume Caumon,et al.  Handling natural complexity in three-dimensional geomechanical restoration, with application to the recent evolution of the outer fold and thrust belt, deep-water Niger Delta , 2013 .

[8]  M. Jessell,et al.  Structural geophysics: Integrated structural and geophysical modelling , 1996 .

[9]  F. Donze,et al.  A DEM model for soft and hard rocks: Role of grain interlocking on strength , 2013 .

[10]  Oskar Vidal-Royo,et al.  Multiple mechanisms driving detachment folding as deduced from 3D reconstruction and geomechanical restoration: the Pico del Águila anticline (External Sierras, Southern Pyrenees) , 2012 .

[11]  Maximilien Siavelis Modélisation numérique X-FEM de grands glissements avec frottement le long d'un réseau de discontinuités , 2011 .

[12]  G. Caumon,et al.  Surface-Based 3D Modeling of Geological Structures , 2009 .

[13]  Andrew Nealen,et al.  Physically Based Deformable Models in Computer Graphics , 2006, Comput. Graph. Forum.

[14]  Bertrand Guillier,et al.  Restoration and balance of a folded and faulted surface by best-fitting of finite elements: principle and applications , 1991 .

[15]  Guillaume Caumon,et al.  Balanced restoration of geological volumes with relaxed meshing constraints , 2010, Comput. Geosci..

[16]  J. Suppe,et al.  3-D Restoration of Complexly Folded and Faulted Surfaces Using Multiple Unfolding Mechanisms , 2000 .

[17]  A. R. Syversveen,et al.  Fault displacement modelling using 3D vector fields , 2010, Computational Geosciences.

[18]  I. Moretti,et al.  KINE3D: a New 3d Restoration Method Based on a Mixed Approach Linking Geometry and Geomechanics , 2006 .

[19]  J. Mallet,et al.  Building and Editing a Sealed Geological Model , 2004 .

[20]  Rollin T. Chamberlin,et al.  The Appalachian Folds of Central Pennsylvania , 1910, The Journal of Geology.

[21]  Dominique Bechmann,et al.  A survey of spatial deformation from a user-centered perspective , 2008, TOGS.

[22]  Vaneeda Allken,et al.  Three‐dimensional numerical modeling of upper crustal extensional systems , 2011 .

[23]  Isabelle Bonniver,et al.  Etude hydrogéologique et dimensionnement par modelisation du "système-traçage" du réseau karstique de Han-sur-Lesse (Massif de Boine - Belgique). , 2011 .

[24]  S. Aanonsen,et al.  Optimized trishear inverse modeling , 2009 .

[25]  Geir Evensen,et al.  Structural Surface Uncertainty Modeling and Updating Using the Ensemble Kalman Filter , 2010 .

[26]  C. D. A. Dahlstrom Balanced cross sections , 1969 .

[27]  Marcin Dabrowski,et al.  Evolution of large amplitude 3D fold patterns: A FEM study , 2008 .

[28]  Jean Braun,et al.  Subducting slabs: Jellyfishes in the Earth's mantle , 2010 .

[29]  Thomas W. Sederberg,et al.  Free-form deformation of solid geometric models , 1986, SIGGRAPH.

[30]  J. Braun,et al.  Numerical models of the evolution of accretionary wedges and fold-and-thrust belts using the distinct-element method , 2002 .

[31]  Matthias Müller,et al.  A versatile and robust model for geometrically complex deformable solids , 2004, Proceedings Computer Graphics International, 2004..

[32]  M. Jessell,et al.  A parametric method to model 3D displacements around faults with volumetric vector fields , 2013 .

[33]  A. Richards,et al.  RESTORATION AND BALANCE OF COMPLEX FOLDED AND FAULTED ROCK VOLUMES : FLEXURAL FLATTENING, JIGSAW FITTING AND DECOMPACTION IN THREE DIMENSIONS , 1997 .

[34]  Serena Jones,et al.  A new technique for 3-D flexural-slip restoration , 2002 .

[35]  I. Moretti,et al.  Working in complex areas: New restoration workflow based on quality control, 2D and 3D restorations , 2008 .

[36]  Ari Rappoport,et al.  Simple constrained deformations for geometric modeling and interactive design , 1994, TOGS.

[37]  Lise Salles Contrôles structuraux en 3 dimensions de la sédimentation turbiditique dans les chaînes plissées : exemple des Grès d’Annot (SE de la France) , 2010 .

[38]  Oskar Vidal-Royo,et al.  The roles of complex mechanical stratigraphy and syn-kinematic sedimentation in fold development: insights from discrete-element modelling and application to the Pico del Águila anticline (External Sierras, Southern Pyrenees) , 2011 .

[39]  Alan H. Barr,et al.  Global and local deformations of solid primitives , 1984, SIGGRAPH.

[40]  Guillaume Caumon,et al.  Method for Stochastic Inverse Modeling of Fault Geometry and Connectivity Using Flow Data , 2012, Mathematical Geosciences.

[41]  J. Mallet,et al.  Editing faults within tetrahedral volume models in real time , 2007 .

[42]  J. Mallet Discrete modeling for natural objects , 1997 .