Stress computation on unstructured geological grids using Finite Elements Method

Antoine Mazuyer and Arnaud Botella and Marianne Conin and Richard Giot. ( 2014 )
in: Proc. 34th Gocad Meeting, Nancy, ASGA

Abstract

Subsurface stress computation is useful for understanding and predicting mechanical problems such as hydraulic fracturation, borehole stability or initiation and propagation of seismic rupture. The Finite Element Method is a numerical technique for solving these mechanical problems. We present how a unstructured conformal geological grid combined with a software using Finite Element Method can compute the stress. To run these numerical simulations, geomechanical properties such as fault friction are needed. The different faults block have to be clearly defined in the mesh for setting up the numerical problem and visualizing the possibly discontinuous displacement field across the faults. The Finite Element Method also needs boundary conditions so the boundary surfaces of the geological model must be expressed as an input of the software. Then, the stress is computed using a Finite Element software. The computation performance for a tetrahedral grid and a hex-dominant grid is compared: for the same resolution, a hex-dominant grid produces results less precise than a tetrahedral grid, but the computation time is faster because there are less equations to solve.

Download / Links

BibTeX Reference

@inproceedings{RUNKJRM29,
 abstract = { Subsurface stress computation is useful for understanding and predicting mechanical problems such as hydraulic fracturation, borehole stability or initiation and propagation of seismic rupture. The Finite Element Method is a numerical technique for solving these mechanical problems. We present how a unstructured conformal geological grid combined with a software using Finite Element Method can compute the stress. To run these numerical simulations, geomechanical properties such as fault friction are needed. The different faults block have to be clearly defined in the mesh for setting up the numerical problem and visualizing the possibly discontinuous displacement field across the faults. The Finite Element Method also needs boundary conditions so the boundary surfaces of the geological model must be expressed as an input of the software. Then, the stress is computed using a Finite Element software.
The computation performance for a tetrahedral grid and a hex-dominant grid is compared: for the same resolution, a hex-dominant grid produces results less precise than a tetrahedral grid, but the computation time is faster because there are less equations to solve. },
 author = { Mazuyer, Antoine AND Botella, Arnaud AND Conin, Marianne AND Giot, Richard },
 booktitle = { Proc. 34th Gocad Meeting, Nancy },
 publisher = { ASGA },
 title = { Stress computation on unstructured geological grids using Finite Elements Method },
 year = { 2014 }
}