English

A matrix-free approach to geostatistical filtering

Methodology 2020-04-07 v1 Statistics Theory Applications Statistics Theory

Abstract

In this paper, we present a novel approach to geostatistical filtering which tackles two challenges encountered when applying this method to complex spatial datasets: modeling the non-stationarity of the data while still being able to work with large datasets. The approach is based on a finite element approximation of Gaussian random fields expressed as an expansion of the eigenfunctions of a Laplace--Beltrami operator defined to account for local anisotropies. The numerical approximation of the resulting random fields using a finite element approach is then leveraged to solve the scalability issue through a matrix-free approach. Finally, two cases of application of this approach, on simulated and real seismic data are presented.

Keywords

Cite

@article{arxiv.2004.02799,
  title  = {A matrix-free approach to geostatistical filtering},
  author = {Mike Pereira and Nicolas Desassis and Cédric Magneron and Nathan Palmer},
  journal= {arXiv preprint arXiv:2004.02799},
  year   = {2020}
}

Comments

25 pages, 8 figures

R2 v1 2026-06-23T14:41:24.749Z