Abstract:In Quaternary deposits, the characterization of subsurface heterogeneity and its associated uncertainty is critical when dealing with groundwater resource management. The combination of different data types through joint inversion has proven to be an effective way to reduce final model uncertainty. Moreover, it allows the final model to be in agreement with a wider spectrum of data available on site. However, integrating them stochastically through an inversion is very time-consuming and resource expensive, due to the important number of physical simulations needed. The use of multi-fidelity models, by combining low-fidelity inexpensive and less accurate models with high-fidelity expensive and accurate models, allows one to reduce the time needed for inversion to converge. This multiscale logic can be applied for the generation of Quaternary models. Most Quaternary sedimentological models can be considered as geological units (large scale), populated with facies (medium scale), and finally completed by physical parameters (small scale). In this paper, both approaches are combined. A simple and fast time-domain EM 1D geophysical direct problem is used to first constrain a simplified stochastic geologically consistent model, where each stratigraphic unit is considered homogeneous in terms of facies and parameters. The ensemble smoother with multiple data assimilation (ES-MDA) algorithm allows generating an ensemble of plausible subsurface realizations. Fast identification of the large-scale structures is the main point of this step. Once plausible unit models are generated, high-fidelity transient groundwater flow models are incorporated. The low-fidelity models are populated stochastically with heterogeneous facies and their associated parameter distribution. ES-MDA is also used for this task by directly inferring the property values (hydraulic conductivity and resistivity) from the generated model. To preserve consistency, geophysical and hydrogeological data are inverted jointly. This workflow ensures that the models are geologically consistent and are therefore less subject to artifacts due to localized poor-quality data. It is able to robustly estimate the associated uncertainty with the final model. Finally, due to the simplification of both the direct problem and the geology during the low-fidelity part of the inversion, it greatly reduces the time required to converge to an ensemble of complex models while preserving consistency.

Multifidelity adaptive sequential Monte Carlo for geophysical inversion

Adaptive sequential Monte Carlo for posterior inference and model selection among complex geological priors

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems

Surrogate-Based Inverse Modeling of the Hydrological System: An Adaptive Approach Considering Surrogate Structural Error.

Bayesian seismic inversion: a fast sampling Langevin dynamics Markov chain Monte Carlo method

Inverse modeling of hydrologic systems with adaptive multi-fidelity Markov chain Monte Carlo simulations

Efficiency and robustness in Monte Carlo sampling of 3-D geophysical inversions with Obsidian v0.1.2: Setting up for success

Adaptive Multifidelity Data Assimilation for Nonlinear Subsurface Flow Problems

An Efficient Bayesian Inversion of a Geothermal Prospect Using a Multivariate Adaptive Regression Spline Method

DE-MCMC Based Stochastic Seismic Inversion Incorporating Gaussian Mixture Model

Geological structure-guided hybrid MCMC and Bayesian linearized inversion methodology

Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes

Sequential Monte Carlo with Highly Informative Observations

Adaptivity in Bayesian Inverse Finite Element Problems: Learning and Simultaneous Control of Discretisation and Sampling Errors

Efficient History Matching Using the Markov-Chain Monte Carlo Method by Means of the Transformed Adaptive Stochastic Collocation Method

Sequential pCN-MCMC, an efficient MCMC method for Bayesian inversion of high-dimensional multi-Gaussian priors

Adaptive Multi-Fidelity Polynomial Chaos Approach to Bayesian Inference in Inverse Problems

Stochastic multi-fidelity joint hydrogeophysical inversion of consistent geological models

Integration of Active Learning and MCMC Sampling for Efficient Bayesian Calibration of Mechanical Properties

Multilevel Sequential Monte Carlo Samplers

Bayesian Inversion of Seismic and Electromagnetic Data for Marine Gas Reservoir Characterization Using Multi-Chain Markov Chain Monte Carlo Sampling