A prefetching technique for prediction of porous media flows

In many applications in flows through porous media, one needs to determine the properties of subsurface to detect, monitor, or predict the actions of natural or induced forces. Here, we focus on two important subsurface properties: rock permeability and porosity. A Bayesian approach using a Markov Chain Monte Carlo (MCMC) algorithm is well suited for reconstructing the spatial distribution of permeability and porosity, and quantifying associated uncertainty in these properties. A crucial step in this approach is the computation of a likelihood function, which involves solving a possibly nonlinear system of partial differential equations. The computation time for the likelihood function limits the number of MCMC iterations that can be performed in a practical period of time. This affects the consistency of the posterior distribution of permeability and porosity obtained by MCMC exploration. To speed-up the posterior exploration, we can use a prefetching technique, which relies on the fact that multiple likelihoods of possible states into the future in an MCMC chain can be computed ahead of time. In this paper, we show that the prefetching technique implemented on multiple processors can make the Bayesian approach computationally tractable for subsurface characterization and prediction of porous media flows.

[1]  Mike West,et al.  Markov Random Field Models for High-Dimensional Parameters in Simulations of Fluid Flow in Porous Media , 2002, Technometrics.

[2]  Anthony Brockwell.  Parallel Markov chain Monte Carlo Simulation by Pre-Fetching , 2006 .

[3]  S. E. Ahmed,et al.  Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference , 2008, Technometrics.

[4]  Christian B Allen,et al.  48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition , 2010 .

[5]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation. , 2008, Journal of chemical theory and computation.

[6]  Felipe Pereira,et al.  Rapid quantification of uncertainty in permeability and porosity of oil reservoirs for enabling predictive simulation , 2014, Math. Comput. Simul..

[7]  Robert A. van de Geijn,et al.  Making Programming Synonymous with Programming for Linear Algebra Libraries FLAME Working Note # 31 , 2008 .

[8]  Inanc Senocak,et al.  An MPI-CUDA Implementation for Massively Parallel Incompressible Flow Computations on Multi-GPU Clusters , 2010 .

[9]  M. Loève,et al.  Elementary Probability Theory , 1977 .

[10]  F. Pereira,et al.  A semi-discrete central scheme for the approximation of two-phase flows in three space dimensions , 2011, Math. Comput. Simul..

[11]  E. Tadmor,et al.  New High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection—Diffusion Equations , 2000 .

[12]  Yalchin Efendiev,et al.  An Efficient Two-Stage Sampling Method for Uncertainty Quantification in History Matching Geological Models , 2008 .

[13]  Eugene Wong,et al.  Stochastic processes in information and dynamical systems , 1979 .

[14]  L. Durlofsky.  Numerical calculation of equivalent grid block permeability tensors for heterogeneous porous media , 1991 .

[15]  Guillermo Rein,et al.  44th AIAA Aerospace Sciences Meeting and Exhibit , 2006 .

[16]  Amitabh Varshney,et al.  High-throughput sequence alignment using Graphics Processing Units , 2007, BMC Bioinformatics.

[17]  Yalchin Efendiev,et al.  Dynamic Data Driven Simulations in Stochastic Environments , 2006, Computing.

[18]  Jim Douglas,et al.  Operator splitting for three-phase flow in heterogeneous porous media , 2009 .

[19]  Felipe Pereira,et al.  Multiple Markov Chains Monte Carlo Approach for Flow Forecasting in Porous Media* , 2012, ICCS.

[20]  Yalchin Efendiev,et al.  An efficient two‐stage Markov chain Monte Carlo method for dynamic data integration , 2005 .

[21]  Felipe Pereira,et al.  On the numerical simulation of waterflooding of heterogeneous petroleum reservoirs , 1997 .

[22]  Michael Mascagni.  SPRNG: A Scalable Library for Pseudorandom Number Generation , 1999, PPSC.

[23]  Ingvar Strid.  Efficient parallelisation of Metropolis-Hastings algorithms using a prefetching approach , 2010, Comput. Stat. Data Anal..

[24]  Avi Bleiweiss,et al.  GPU accelerated pathfinding , 2008, GH '08.

[25]  Graham Pullan,et al.  Acceleration of a 3D Euler solver using commodity graphics hardware , 2008 .

[26]  T. LaForce,et al.  Bayesian Reservoir History Matching Considering Model and Parameter Uncertainties , 2012, Mathematical Geosciences.

[27]  Paulius Micikevicius,et al.  3D finite difference computation on GPUs using CUDA , 2009, GPGPU-2.

[28]  Yuanle Ma,et al.  Computational methods for multiphase flows in porous media , 2007, Math. Comput..

[29]  Inanc Senocak,et al.  CUDA Implementation of a Navier-Stokes Solver on Multi-GPU Desktop Platforms for Incompressible Flows , 2009 .

[30]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[31]  Manfred Krafczyk,et al.  TeraFLOP computing on a desktop PC with GPUs for 3D CFD , 2008 .

[32]  Geoff K. Nicholls,et al.  Sampling Conductivity Images via MCMC , 2007 .

[33]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[34]  Michael Presho,et al.  Application of the two-stage Markov chain Monte Carlo method for characterization of fractured reservoirs using a surrogate flow model , 2011 .

[35]  Weiguo Liu,et al.  Molecular Dynamics Simulations on Commodity GPUs with CUDA , 2007, HiPC.

[36]  G. Dagan.  Flow and transport in porous formations , 1989 .

[37]  Michael Mascagni,et al.  SPRNG: A Scalable Library for Pseudorandom Number Generation , 1999, PP.

[38]  C. Fox,et al.  Markov chain Monte Carlo Using an Approximation , 2005 .

[39]  Felipe Pereira,et al.  A MULTI-STAGE BAYESIAN PREDICTION FRAMEWORK FOR SUBSURFACE FLOWS , 2013 .

[40]  Mike Christie,et al.  Uncertainty quantification for porous media flows , 2006, J. Comput. Phys..

[41]  Charles J. Geyer,et al.  Practical Markov Chain Monte Carlo , 1992 .

[42]  Yalchin Efendiev,et al.  Preconditioning Markov Chain Monte Carlo Simulations Using Coarse-Scale Models , 2006, SIAM J. Sci. Comput..

[43]  Manish Vachharajani,et al.  GPU acceleration of numerical weather prediction , 2008, 2008 IEEE International Symposium on Parallel and Distributed Processing.