Mesh Generation for 3D Geological Reservoirs with Arbitrary Stratigraphic Surface Constraints

With the advanced field observation, image and drilling technology, geological structure of reservoirs can be described in more details. A novel 3D mesh generation method for geological reservoir models is proposed and implemented with arbitrary stratigraphical surface constraints , which ensures the detailed geological structure and material properties are well described and analysed. Specifically, the stratigraphic interfaces are firstly extracted and meshed, and then the tetrahedron mesh is generated with the constraints of the meshed surfaces. The proposed approach includes the following five steps: (1) extracting stratum interfaces; (2) creating a background mesh with size field on the interfaces; (3) constructing geodesic isolines from the interface boundaries to the interior; (4) employing a geodesicbased approach to create triangles on the area between adjacent isolines and to combine them together; (5) generating tetrahedron mesh for 3D reservoirs with constraints of the above surface triangular mesh. Th is approach has been implemented and applied to the Lawn Hill reservoir model to demonstrate its effectiveness and usefulness.

[1]  Hua Xu,et al.  An approach to computer modeling and visualization of geological faults in 3D , 2003 .

[2]  J A Sethian,et al.  Computing geodesic paths on manifolds. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Chandrajit L Bajaj,et al.  An Automatic 3D Mesh Generation Method for Domains with Multiple Materials. , 2010, Computer methods in applied mechanics and engineering.

[4]  S. H. Lo,et al.  OPTIMIZATION OF TETRAHEDRAL MESHES BASED ON ELEMENT SHAPE MEASURES , 1997 .

[5]  Huilin Xing,et al.  Finite element modelling of interacting fault systems , 2007 .

[6]  Yehoshua Y. Zeevi,et al.  The farthest point strategy for progressive image sampling , 1997, IEEE Trans. Image Process..

[7]  Michela Spagnuolo,et al.  Shape Analysis and Structuring , 2008 .

[8]  Marco Attene,et al.  Recent Advances in Remeshing of Surfaces , 2008, Shape Analysis and Structuring.

[9]  Mariette Yvinec,et al.  Feature preserving Delaunay mesh generation from 3D multi‐material images , 2009, Comput. Graph. Forum.

[10]  A. Ord,et al.  Fluid Flow during Deformation Associated with Structural Closure of the Isa Superbasin at 1575 Ma in the Central and Northern Lawn Hill Platform, Northern Australia , 2006 .

[11]  Mark Meyer,et al.  Interactive geometry remeshing , 2002, SIGGRAPH.

[12]  M. Glikson,et al.  Timing and Chemistry of Fluid-Flow Events in the Lawn Hill Platform, Northern Australia , 2006 .

[13]  Andrea Tonini,et al.  Three-dimensional reconstruction of the Carrara Syncline (Apuane Alps, Italy): An approach to reconstruct and control a geological model using only field survey data , 2009, Comput. Geosci..

[14]  R E Guldberg,et al.  The accuracy of digital image-based finite element models. , 1998, Journal of biomechanical engineering.

[15]  John M. Schreiner,et al.  Inter-surface mapping , 2004, SIGGRAPH 2004.

[16]  L. Feltrin.  Predictive modelling of prospectivity for Pb-Zn deposits in the Lawn Hill Region, Queensland, Australia , 2008 .

[17]  Huilin Xing,et al.  An indirect approach for automatic generation of quadrilateral meshes with arbitrary line constraints , 2011 .

[18]  Qiang Du,et al.  Constrained boundary recovery for three dimensional Delaunay triangulations , 2004 .

[19]  Pierre Alliez,et al.  Isotropic Remeshing of Surfaces: A Local Parameterization Approach , 2003, IMR.

[20]  Nicolas S. Holliman,et al.  Integration of regional to outcrop digital data: 3D visualisation of multi-scale geological models , 2009, Comput. Geosci..

[21]  L. Cohen,et al.  Surface segmentation using geodesic centroidal tesselation , 2004 .

[22]  P. George,et al.  3D Delaunay mesh generation coupled with an advancing-front approach , 1998 .

[23]  Arthur W. Toga,et al.  The generation of tetrahedral mesh models for neuroanatomical MRI , 2011, NeuroImage.

[24]  Laurent D. Cohen,et al.  Geodesic Remeshing Using Front Propagation , 2003, International Journal of Computer Vision.

[25]  Jonathan Richard Shewchuk,et al.  What is a Good Linear Element? Interpolation, Conditioning, and Quality Measures , 2002, IMR.

[26]  Vladislav Kraevoy,et al.  Cross-parameterization and compatible remeshing of 3D models , 2004, SIGGRAPH 2004.

[27]  Peter Schröder,et al.  Consistent mesh parameterizations , 2001, SIGGRAPH.

[28]  Martin Isenburg,et al.  Isotropic surface remeshing , 2003, 2003 Shape Modeling International..

[29]  Alla Sheffer,et al.  Geodesic-based Surface Remeshing , 2003, IMR.

[30]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[31]  Christopher C. Pain,et al.  A systematic approach to unstructured mesh generation for ocean modelling using GMT and Terreno , 2008, Comput. Geosci..

[32]  Ziji Wu,et al.  Multiple material marching cubes algorithm , 2003 .

[33]  Chandrajit L. Bajaj,et al.  Surface Smoothing and Quality Improvement of Quadrilateral/Hexahedral Meshes with Geometric Flow. , 2005 .

[34]  Huilin Xing,et al.  A boundary focused quadrilateral mesh generation algorithm for multi-material structures , 2013, J. Comput. Phys..

[35]  Cheng-Chung Lin,et al.  Consistent parametrization by quinary subdivision for remeshing and mesh metamorphosis , 2003, GRAPHITE '03.