Iterative refinement of implicit boundary models for improved geological feature reproduction

Geological domains contain non-stationary features that cannot be described by a single direction of continuity. Non-stationary estimation frameworks generate more realistic curvilinear interpretations of subsurface geometries. A radial basis function (RBF) based implicit modeling framework using domain decomposition is developed that permits introduction of locally varying orientations and magnitudes of anisotropy for boundary models to better account for the local variability of complex geological deposits. The interpolation framework is paired with a method to automatically infer the locally predominant orientations, which results in a rapid and robust iterative non-stationary boundary modeling technique that can refine locally anisotropic geological shapes automatically from the sample data. The method also permits quantification of the volumetric uncertainty associated with the boundary modeling. The methodology is demonstrated on a porphyry dataset and shows improved local geological features. Partition of Unity is used to model non-stationary implicit geological shapes.Iterative refinements automatically extract non-stationary geological features.Complex geological modeling is improved considering local anisotropy.

[1]  Richard K. Beatson,et al.  Reconstruction and representation of 3D objects with radial basis functions , 2001, SIGGRAPH.

[2]  Clayton V. Deutsch,et al.  Non-stationary Geostatistical Modeling Based on Distance Weighted Statistics and Distributions , 2012, Mathematical Geosciences.

[3]  Determination of volumetric uncertainty for geological bodies detected by boreholes , 2015 .

[4]  Clayton V. Deutsch,et al.  Programs for kriging and sequential Gaussian simulation with locally varying anisotropy using non-Euclidean distances , 2011, Comput. Geosci..

[5]  C. Lajaunie,et al.  Foliation fields and 3D cartography in geology: Principles of a method based on potential interpolation , 1997 .

[6]  Salvatore Cuomo,et al.  Surface reconstruction from scattered point via RBF interpolation on GPU , 2013, 2013 Federated Conference on Computer Science and Information Systems.

[7]  Silva Maureira,et al.  Enhanced Geologic Modeling with Data-Driven Training Images for Improved Resources and Recoverable Reserves , 2015 .

[8]  Clayton V. Deutsch,et al.  Mineral Resource Estimation , 2013 .

[9]  Qi Xia,et al.  Implicit fitting and smoothing using radial basis functions with partition of unity , 2005, Ninth International Conference on Computer Aided Design and Computer Graphics (CAD-CG'05).

[10]  F. Fouedjio,et al.  Predictive Geological Mapping Using Closed-Form Non-stationary Covariance Functions with Locally Varying Anisotropy: Case Study at El Teniente Mine (Chile) , 2016, Natural Resources Research.

[11]  Giulio Casciola,et al.  Shape preserving surface reconstruction using locally anisotropic radial basis function interpolants , 2006, Comput. Math. Appl..

[12]  M. Jessell,et al.  Locating and quantifying geological uncertainty in three-dimensional models: Analysis of the Gippsland Basin, southeastern Australia , 2012 .

[13]  Richard K. Beatson,et al.  Fast fitting of radial basis functions: Methods based on preconditioned GMRES iteration , 1999, Adv. Comput. Math..

[14]  Richard K. Beatson,et al.  Fast Solution of the Radial Basis Function Interpolation Equations: Domain Decomposition Methods , 2000, SIAM J. Sci. Comput..

[15]  Hans-Peter Seidel,et al.  Sparse surface reconstruction with adaptive partition of unity and radial basis functions , 2006, Graph. Model..

[16]  Hans-Peter Seidel,et al.  A multi-scale approach to 3D scattered data interpolation with compactly supported basis functions , 2003, 2003 Shape Modeling International..

[17]  M. Shephard,et al.  Automatic generation of octree-based three-dimensional discretizations for Partition of Unity methods , 2000 .

[18]  E. Schetselaar,et al.  Three-Dimensional Modelling of Geological Surfaces Using Generalized Interpolation with Radial Basis Functions , 2014, Mathematical Geosciences.

[19]  Mauro Roisenberg,et al.  Spatial interpolation: an analytical comparison between kriging and RBF networks , 2013, ACM Symposium on Applied Computing.

[20]  C. Deutsch,et al.  Facies proportion uncertainty in presence of a trend , 2017 .

[21]  Joachim Pouderoux,et al.  Adaptive hierarchical RBF interpolation for creating smooth digital elevation models , 2004, GIS '04.

[22]  Jeff B. Boisvert,et al.  Stochastic Distance Based Geological Boundary Modeling with Curvilinear Features , 2013, Mathematical Geosciences.

[23]  P. Milanfar,et al.  Multiscale principal components analysis for image local orientation estimation , 2002, Conference Record of the Thirty-Sixth Asilomar Conference on Signals, Systems and Computers, 2002..

[24]  Jeff B. Boisvert,et al.  Inference of locally varying anisotropy fields from diverse data sources , 2015, Comput. Geosci..

[25]  A. Cruden,et al.  Regional dome evolution and its control on ore-grade distribution: Insights from 3D implicit modelling of the Navachab gold deposit, Namibia , 2015 .

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