Modeling and simulating depositional sequences using latent Gaussian random fields

DOI10.1007/S11004-020-09875-0MaRDI QIDQ2040696zbMATH OpenOpenAlexFDO

Authors Paolo Fabbri, Carlo Gaetan, D. Allard

Publication date 14 July 2021

Published in Mathematical Geosciences (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2003.11383

zbMATH Keywords

conditional simulationdata augmentationPC priorstochastic 3D modelstratigraphic sequencesubsoil modeling

Mathematics Subject Classification ID

Random fields; image analysis (62M40) Geostatistics (86A32)

Abstract: Simulating a depositional (or stratigraphic) sequence conditionally on borehole data is a long-standing problem in hydrogeology and in petroleum geostatistics. This paper presents a new rule-based approach for simulating depositional sequences of surfaces conditionally on lithofacies thickness data. The thickness of each layer is modeled by a transformed latent Gaussian random field allowing for null thickness thanks to a truncation process. Layers are sequentially stacked above each other following the regional stratigraphic sequence. By choosing adequately the variograms of these random fields, the simulated surfaces separating two layers can be continuous and smooth. Borehole information is often incomplete in the sense that it does not provide direct information as to the exact layer some observed thickness belongs to. The latent Gaussian model proposed in this paper offers a natural solution to this problem by means of a Bayesian setting with a Markov Chain Monte Carlo (MCMC) algorithm that can explore all possible configurations compatible with the data. The model and the associated MCMC algorithm are validated on synthetic data and then applied to a subsoil in the Venetian Plain with a moderately dense network of cored boreholes.

Recommendations

• Geostatistical modelling of cyclic and rhythmic facies architectures
• scientific article; zbMATH DE number 1934834
• Conditional simulation of non-Gaussian random functions
• Conditional simulation of complex geological structures using multiple-point statistics
• Gradual deformation and iterative calibration of sequential stochastic simulations

Cites work

• A Latent Gaussian Markov Random-Field Model for Spatiotemporal Rainfall Disaggregation
• Computation of multivariate normal and t probabilities
• Conditional simulation of complex geological structures using multiple-point statistics
• Conditional simulation of multi-type non stationary Markov object models respecting specified proportions
• Conditioning of Marked Point Processes within a Bayesian Framework
• Constructing Priors that Penalize the Complexity of Gaussian Random Fields
• Geostatistical modelling of cyclic and rhythmic facies architectures
• Geostatistics. Modeling spatial uncertainty.
• Gibbs Sampling
• Inconsistent Estimation and Asymptotically Equal Interpolations in Model-Based Geostatistics
• Penalising model component complexity: a principled, practical approach to constructing priors
• Statistical analysis of zero-inflated nonnegative continuous data: a review
• Statistics for spatial data
• Tools for statistical inference. Methods for the exploration of posterior distributions and likelihood functions.
• Transition probability-based indicator geostatistics

Cited in

(3)

• Geostatistical modelling of cyclic and rhythmic facies architectures
• Stochastic models for the completeness of stratigraphic sections
• Geostatistical rock physics inversion for predicting the spatial distribution of porosity and saturation in the critical zone

Describes a project that uses

Uses Software

• mvtnorm
• spMC

This page was built for publication: Modeling and simulating depositional sequences using latent Gaussian random fields

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2040696)