A Gradient-Based Blocking Markov Chain Monte Carlo Method for Stochastic Inverse Modeling

Inverse modeling for subsurface flow and transport in porous media is expected to improve the reliability of predictions in that the realizations generated are consistent with the observations of states. A gradient-based blocking Markov chain Monte Carlo (McMC) method is presented for stochastic inverse modeling. The method proposed effectively takes advantage of gradient information for tuning each realization to create a new “candidate” proposal, and hence it is capable of improving the performance of McMC. The gradients are efficiently computed by an adjoint method. The proposal mechanism is based on the optimization of a random seed field (or probability field), and thus it is able to preserve the prior model statistics. The method proposed has better performances than the single-component McMC and also avoids directly solving a difficult large-scale ill-conditioned optimization problem simply by turning it into a sampling procedure plus a sequence of well-conditioned optimization subproblems. A synthetic example demonstrates the method proposed.

[1]  James Martin,et al.  A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion , 2012, SIAM J. Sci. Comput..

[2]  J. Gómez-Hernández,et al.  Uncertainty assessment and data worth in groundwater flow and mass transport modeling using a blocking Markov chain Monte Carlo method. , 2009 .

[3]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

[4]  Binayak P. Mohanty,et al.  Efficient uncertainty quantification techniques in inverse problems for Richards’ equation using coarse-scale simulation models , 2009 .

[5]  J. Gómez-Hernández,et al.  Joint Sequential Simulation of MultiGaussian Fields , 1993 .

[6]  Philippe Renard,et al.  Blocking Moving Window algorithm: Conditioning multiple‐point simulations to hydrogeological data , 2010 .

[7]  M. Girolami,et al.  Riemann manifold Langevin and Hamiltonian Monte Carlo methods , 2011, Journal of the Royal Statistical Society: Series B (Statistical Methodology).

[8]  K. Hanson Use of probability gradients in hybrid MCMC and a new convergence test , 2002 .

[9]  Variational method for estimating the rate of convergence of Markov-chain Monte Carlo algorithms. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Jianlin Fu,et al.  Preserving spatial structure for inverse stochastic simulation using blocking Markov chain Monte Carlo method , 2008 .

[11]  H. Tchelepi,et al.  A multiscale adjoint method to compute sensitivity coefficients for flow in heterogeneous porous media , 2010 .

[12]  Hamdi A. Tchelepi,et al.  A multiscale method for subsurface inverse modeling: Single-phase transient flow , 2011 .

[13]  S. Duane,et al.  Hybrid Monte Carlo , 1987 .

[14]  A Blocking Markov Chain Monte Carlo Method for Inverse Stochastic Hydrogeological Modeling , 2009 .

[15]  Yuan Qi,et al.  Hessian-based Markov Chain Monte-Carlo Algorithms , 2017 .

[16]  John Shawe-Taylor,et al.  Advances in Neural Information Processing Systems 24: 25th Annual Conference on Neural Information Processing Systems 2011. Proceedings of a meeting held 12-14 December 2011, Granada, Spain , 2011, NIPS.