Full text loading...
We introduce a novel workflow that integrates non-Euclidean shortest path distances (SPD) using locally varying geometric anisotropy (LVA) to define spatial covariance matrices. The resulting covariance can be used in both geophysical and geostatistical algorithms where spatial continuity is utilized. The use of LVA in proposed methodology bypasses the limitation of constant geometric anisotropy traditionally assumed in geostatistical methods; it produces a more accurate covariance and generates more geologically realistic models. The Hamiltonian Fast Marching (HFM) algorithm is more precise than other approaches to evaluate SPD, such as the Dijkstra algorithm. While robust spatial covariance estimation needs redundant measurements in any direction, in petroleum exploration only the vertical direction is sampled sufficiently. Lateral direction data scarcity is mitigated by extracting the local dips and anisotropy ratios along the major and minor continuity directions from the seismic data itself. We demonstrate the HFM method on Bayesian post-stack inversion and LVA ordinary kriging (OK) for estimating acoustic impedance. The proposed workflow is also valid for other inversion algorithms, varieties of kriging, and sequential Gaussian simulation.