Kriging in the Presence of Locally Varying Anisotropy Using Non-Euclidean Distances

A stationary specification of anisotropy does not always capture the complexities of a geologic site. In this situation, the anisotropy can be varied locally. Directions of continuity and the range of the variogram can change depending on location within the domain being modeled. Kriging equations have been developed to use a local anisotropy specification within kriging neighborhoods; however, this approach does not account for variation in anisotropy within the kriging neighborhood. This paper presents an algorithm to determine the optimum path between points that results in the highest covariance in the presence of locally varying anisotropy. Using optimum paths increases covariance, results in lower estimation variance and leads to results that reflect important curvilinear structures. Although CPU intensive, the complex curvilinear structures of the kriged maps are important for process evaluation. Examples highlight the ability of this methodology to reproduce complex features that could not be generated with traditional kriging.

[1]  David D. Pollard,et al.  Fundamentals of Structural Geology , 2005 .

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

[3]  Anders Løland,et al.  Spatial covariance modelling in a complex coastal domain by multidimensional scaling , 2003 .

[4]  H. Reading,et al.  Sedimentary Environments: Processes, Facies And Stratigraphy , 2011 .

[5]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[6]  S T Roweis,et al.  Nonlinear dimensionality reduction by locally linear embedding. , 2000, Science.

[7]  P. Guttorp,et al.  Nonparametric Estimation of Nonstationary Spatial Covariance Structure , 1992 .

[8]  J. Guilbert,et al.  The geology of ore deposits , 1985 .

[9]  Montserrat Fuentes,et al.  Modeling and predicting complex space–time structures and patterns of coastal wind fields , 2005 .

[10]  S. Mitra,et al.  Structural geology of fold and thrust belts , 1992 .

[11]  W. E. Galloway,et al.  Terrigenous Clastic Depositional Systems: Applications to Petroleum, Coal, and Uranium Exploration , 1996 .

[12]  T. Jones,et al.  Conditioning Geologic Models to Local Continuity Azimuth in Spectral Simulation , 2007 .

[13]  P. Ruffo,et al.  Depth, Dip and Gradient , 1993 .

[14]  Frank C. Curriero,et al.  On the Use of Non-Euclidean Distance Measures in Geostatistics , 2007 .

[15]  Michael Edward Hohn,et al.  An Introduction to Applied Geostatistics: by Edward H. Isaaks and R. Mohan Srivastava, 1989, Oxford University Press, New York, 561 p., ISBN 0-19-505012-6, ISBN 0-19-505013-4 (paperback), $55.00 cloth, $35.00 paper (US) , 1991 .

[16]  Clayton V. Deutsch,et al.  Geostatistical Reservoir Modeling , 2002 .

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

[18]  J. Tenenbaum,et al.  A global geometric framework for nonlinear dimensionality reduction. , 2000, Science.

[19]  Ya-Xiang Yuan,et al.  Optimization Theory and Methods: Nonlinear Programming , 2010 .

[20]  G. Bennison.  An introduction to geological structures and maps , 1975 .