Restauration de bassins sédimentaires avec gocad Kine3D J-L. Rudkiewicz, M. Guiton, J-F. Lecomte, I. Moretti, S. Brocheray, M-C. Cacas, (IFP) Février 2009 23 février 2009
Summary Principles of geological restoration forward and backward reconstruction Kine3D plugins within Gocad Examples of restorations Limits when using Code-Aster 2
Backward restoration Present day at -65 Ma 3
Section in the Central Appenines 5 km 4
Forward modelling 5
Summary Principles of geological restoration Kine3D plugins within Gocad Examples of restorations Limits when using Code-Aster 6
Gocad Suite Base Module Seismic Interpretation Velocity Modeling Geologic Interpretation Reservoir Modeling Reservoir Production Drilling Planner http://www.pdgm.com 7
KINE3D Volumic restoration Kine3D 3: link Code Aster and Gocad Tetrahedra, Hexahedras Final geometry Gocad is used as a mesh generator, visualisation tool and post processing tool 8
Principle of 3D restoration The deformation energy of a rock volume is minimized, assuming displacement constraints on some parts Rocks follow linear elasticity Displacement constraints are used to input geological constraints Finite element meshes are used to discretize the equations Code Aster is used to perform computations 9
Restoration constraints 10
Flattening = in AFFE_CHAR_MECA(... DDL_IMPO=( _F( NOEUD='N1024',DZ=2220.3, ),... 11
Fault Lip Closure { } α α + + = + LIAISON_DDL=( _F ( NOEUD=('N777', 'N779' ), DDL=( 'DX', 'DX' ), COEF_MULT=( 1.0, -1 ), COEF_IMPO= -3.55859, ), 12
Fault Gliding Contact AFFE_CHAR_MECA( MODELE= mode_045, CONTACT=( _F( METHODE='CONTINUE', INTEGRATION='GAUSS', ITER_GEOM_MAXI=5, GROUP_MA_MAIT=('R2_5_MA'), GROUP_MA_ESCL=('R2_6_MA'), GLISSIERE='OUI', TOLE_PROJ_EXT=1000, CONTACT_INIT='OUI', ), = 13
Preventing Rigid Body Movements with Point Constraints (DDL_IMPO...) DDL_IMPO=( _F( NOEUD='N102',DX=0., ),... 14
Pin Faces = = DDL_IMPO=( _F( NOEUD='N102',DX=0., ),... 15
Summary Principles of geological restoration Kine3D plugins within Gocad Examples of restorations Limits when using Code-Aster 16
Données originales Agha Jari 20 x 16 x 7 km 17
Horizons 18
Maillages avec 6 grilles 19
Restauration des 6 grilles séparément Pin faces in SW and NE sides 20
Construction of missing parts 21
Putting the missing part in the right position Node contact on visible points Contact on pink fault 22
Summary Principles of geological restoration forward and backward reconstruction Kine3D plugins within Gocad Examples of restorations Limits when using Code-Aster the limited size of grids usable the behaviour with various types of COMP_INCR=(_F(RELATION='ELAS', DEFORMATION='GREEN',),), 23
Limits on grid sizes on a Linux 32 bits/2mo nx ny nz Noeuds SG Noeuds aster dont Lagrange Equations NOMBRE DE COEFFICIENT S NON NULS DANS LA MATRICE 20x22x6_C0plat 20 22 6 2640 5146 1660 12118 367038 40x50x6_C0plat 40 50 6 12000 20576 5680 50368 1621410 60x70x6_C0plat 60 70 6 25200 41491 10810 102853 3375720 40x50x12_C0plat 40 50 12 24000 34448 6784 89776 3134958 60x70x15_C0plat 60 70 15 63000 83314 13186 223570 8105472 24
Différence de comportement Green élastique Simo_Miehe/VMises_isot_line en élastique Contraintes : jonctions sur la partie jaune contact sur la faille rouge 25
Exemple "réel" visé à terme Il reste du chemin à parcourir! 26
Restoration with different rheology red grid: homogeneous stiff rheology of sandstone type green grid: soft rheology of shale type in the graben 27
Merci 28
Output data: dilatancy Local variation of area/volume induced by the restoration Dilatation is zero if one may unfold and close the fault lip without changing the area. Negative dilatation means global extension, positive gobal compression Dilatation is computed on area for surfaces and on volumes for 3D. θ = 29
Dilatancy depends upon boundary conditions 30
Dilatancy depends upon boundary conditions 31
Dilatancy depends upon boundary conditions 32
Dilatancy depends upon boundary conditions 33
Conclusions Fault sealing: contact to quality control the fault geometry, glide to get access to past past geometry or to define local deformation related to fault geometry Restoration geometry sensible to fault contact conditions not very sensible to rigid rotation preventing conditions not very sensitive to rheology Restoration derived parameters sensitive to rigid rotation preventing conditions sensitive to rheology contrast 34