Non-stationary Geostatistical Modeling: A Case Study Comparing LVA Estimation Frameworks

Incorporating locally varying anisotropy (LVA) in geostatistical modeling improves estimates for structurally complex domains where a single set of anisotropic parameters modeled globally do not account for all geological features. In this work, the properties of two LVA-geostatistical modeling frameworks are explored through application to a complexly folded gold deposit in Ghana. The inference of necessary parameters is a significant requirement of geostatistical modeling with LVA; this work focuses on the case where LVA orientations, derived from expert geological interpretation, are used to improve the grade estimates. The different methodologies for inferring the required parameters in this context are explored. The results of considering different estimation frameworks and alternate methods of parameterization are evaluated with a cross-validation study, as well as visual inspection of grade continuity along select cross sections. Results show that stationary methodologies are outperformed by all LVA techniques, even when the LVA framework has minimal guidance on parameterization. Findings also show that additional improvements are gained by considering parameter inference where the LVA orientations and point data are used to infer the local range of anisotropy. Considering LVA for geostatistical modeling of the deposit considered in this work results in better reproduction of curvilinear geological features.

[1]  Jeff B. Boisvert,et al.  Iterative refinement of implicit boundary models for improved geological feature reproduction , 2017, Comput. Geosci..

[2]  T. J. O'Neil Application of Computers and Operations Research in the Mineral Industry , 1979 .

[3]  Wenlong Xu,et al.  Conditional curvilinear stochastic simulation using pixel-based algorithms , 1996 .

[4]  Grade Modelling with Local Anisotropy Angles: a Practical Point of View , .

[5]  F. Fouedjio Space Deformation Non-stationary Geostatistical Approach for Prediction of Geological Objects: Case Study at El Teniente Mine (Chile) , 2016, Natural Resources Research.

[6]  Y. Oshman,et al.  Averaging Quaternions , 2007 .

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

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

[9]  Jacques Rivoirard,et al.  A Generalized Convolution Model and Estimation for Non-stationary Random Functions , 2014, 1412.1373.

[10]  Francky Fouedjio,et al.  Second-order non-stationary modeling approaches for univariate geostatistical data , 2017, Stochastic Environmental Research and Risk Assessment.

[11]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[12]  Francky Fouedjio,et al.  Estimation of Space Deformation Model for Non-stationary Random Functions , 2014, 1412.1344.

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

[14]  J. Snepvangers,et al.  Mapping Curvilinear Structures with Local Anisotropy Kriging , 2005 .

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

[16]  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.

[17]  Clayton V. Deutsch,et al.  Geostatistics with Locally Varying Anisotropy , 2009 .