Beyond Stationary Simulation; Modern Approaches to Stochastic Modelling

TL;DR

This paper explores modern non-stationary stochastic simulation methods, including GANs, for more realistic spatial modeling, demonstrating their effectiveness through synthetic experiments and a real geochemical case study.

Contribution

It introduces three advanced stochastic simulation techniques, incorporating non-stationary covariance, multipoint simulation, and GANs, surpassing classical stationary methods in realism and speed.

Findings

Non-stationary methods better preserve complex spatial structures.

GAN-based simulations achieve faster generation times.

Synthetic and real data experiments validate improved accuracy.

Abstract

Stochastic and conditional simulation methods have been effective towards producing realistic realizations and simulations of spatial numerical models that share equal probability of occurrence. Application of these methods are valuable throughout the domain of earth science for their ability to simulate sampled study data. Such stochastic methods have also been adopted into other fields outside of geostatistics domains, especially within the computing and data science community. Classical techniques for stochastic simulation have primarily consisted of stationary methods due to their brisk simulation speed and mathematical simplicity. However, advances in modern computing now allow for the implementation of more advanced non-stationary simulation methods, consisting of multiple varying structures, and allowing for much more accurate and realistic simulations. As some of these calculations may still be slow, the application of machine learning techniques, namely the Generative Adversarial Network (GAN) can be used to surpass previous simulation generation speeds and allows for greater parameterization. This work presents three stochastic simulation methods: stochastic simulation using non-stationary covariance, multipoint simulation, and conditional GANs. An SPDE method was used as a benchmark comparison. Experiments using synthetic data are used to showcase the effectiveness of each of these methods at maintaining non-stationary structures and conditioned data. A case study implementing non-stationary covariances is also presented on real geochemical samples coming from La Ceinture de roches vertes de la Haute-Eastmain, located in the Superior Province.