Accelerating geostatistical simulations using graphics processing units (GPU)

Geostatistical simulations have become a widely used tool for modeling of oil and gas reservoirs and the assessment of uncertainty. One important current issue is the development of high-resolution models in a reasonable computational time. A possible solution is based on taking advantage of parallel computational strategies. In this paper we present a new methodology that exploits the benefits of graphics processing units (GPUs) along with the master-slave architecture for geostatistical simulations that are based on random paths. The methodology is a hybrid method in which different levels of master and slave processors are used to distribute the computational grid points and to maximize the use of multiple processors utilized in GPU. It avoids conflicts between concurrently simulated grid points, an important issue in high-resolution and efficient simulations. For the sake of comparison, two distinct parallelization methods are implemented, one of which is specific to pattern-based simulations. To illustrate the efficiency of the method, the algorithm for the simulation of pattern is adapted with the GPU. Performance tests are carried out with three large grid sizes. The results are compared with those obtained based on simulations with central processing units (CPU). The comparison indicates that the use of GPUs reduces the computation time by a factor of 26-85.

[1]  Grégoire Mariethoz,et al.  A general parallelization strategy for random path based geostatistical simulation methods , 2010, Comput. Geosci..

[2]  Xing Gao,et al.  Accelerating POCS interpolation of 3D irregular seismic data with Graphics Processing Units , 2010, Comput. Geosci..

[3]  Christian Obrecht,et al.  LBM based flow simulation using GPU computing processor , 2010, Comput. Math. Appl..

[4]  Jacek Strzelczyk,et al.  Analysis of ground deformations based on parallel geostatistical computations of PSInSAR data , 2009, 2009 17th International Conference on Geoinformatics.

[5]  Muhammad Sahimi,et al.  Development of optimal models of porous media by combining static and dynamic data: the permeability and porosity distributions. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  H. S. Vargas,et al.  Parallelization of Sequential Simulation Procedures , 2007 .

[7]  Pejman Tahmasebi,et al.  Multiple-point geostatistical modeling based on the cross-correlation functions , 2012, Computational Geosciences.

[8]  Sinopec Exploration.  New alternative to geophysical high performance computing: GPU Computing. , 2007 .

[9]  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.

[10]  E. I. Polyakova,et al.  The Nu Expression for Probabilistic Data Integration , 2007 .

[11]  Andre G. Journel,et al.  Conditional Indicator Simulation: Application to a Saskatchewan uranium deposit , 1984 .

[12]  Massimiliano Fatica,et al.  Medical Image Processing Using GPU-Accelerated ITK Image Filters , 2011 .

[13]  Muhammad Sahimi,et al.  Computer simulation of gas generation and transport in landfills. III: Development of lanfills' optimal model , 2007 .

[14]  J. Caers,et al.  Conditional Simulation with Patterns , 2007 .

[15]  Rúben Nunes,et al.  Parallelization of sequential Gaussian, indicator and direct simulation algorithms , 2010, Comput. Geosci..

[16]  J. Kärger,et al.  Flow and Transport in Porous Media and Fractured Rock , 1996 .

[17]  Tomas Svensson,et al.  Parallel computing with graphics processing units for high-speed Monte Carlo simulation of photon migration. , 2008, Journal of biomedical optics.

[18]  Muhammad Sahimi,et al.  Computer simulation of gas generation and transport in landfills. V: Use of artificial neural network and the genetic algorithm for short- and long-term forecasting and planning , 2011 .

[19]  G. Mariéthoz,et al.  An Improved Parallel Multiple-point Algorithm Using a List Approach , 2011 .

[20]  A. Journel.  Combining Knowledge from Diverse Sources: An Alternative to Traditional Data Independence Hypotheses , 2002 .

[21]  J. Xu.  OpenCL – The Open Standard for Parallel Programming of Heterogeneous Systems , 2009 .

[22]  Kenneth A. Hawick,et al.  Kriging Interpolation on High-Performance Computers , 1998, HPCN Europe.

[23]  Muhammad Sahimi,et al.  Development of optimal models of porous media by combining static and dynamic data: the porosity distribution. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  Roussos Dimitrakopoulos,et al.  Generalized Sequential Gaussian Simulation on Group Size ν and Screen-Effect Approximations for Large Field Simulations , 2004 .

[25]  J. Chilès,et al.  Geostatistics: Modeling Spatial Uncertainty , 1999 .

[26]  Sebastien Strebelle,et al.  Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics , 2002 .

[27]  R. Olea.  Geostatistics for Natural Resources Evaluation By Pierre Goovaerts, Oxford University Press, Applied Geostatistics Series, 1997, 483 p., hardcover, $65 (U.S.), ISBN 0-19-511538-4 , 1999 .

[28]  Albrecht Gebhardt,et al.  PVM kriging with R , 2003 .

[29]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[30]  Dan Cornford,et al.  Parallel Geostatistics for Sparse and Dense Datasets , 2010 .

[31]  Julián M. Ortiz,et al.  Multiple Point Geostatistical Simulation with Simulated Annealing: Implementation Using Speculative Parallel Computing , 2010 .

[32]  Jeffrey A. Pedelty,et al.  High Performance Geostatistical Modeling of Biospheric Resources , 2004 .

[33]  M. Sahimi,et al.  Generation of long-range correlations in large systems as an optimization problem. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  J. Caers,et al.  Stochastic Simulation of Patterns Using Distance-Based Pattern Modeling , 2010 .

[35]  Paul Switzer,et al.  Filter-Based Classification of Training Image Patterns for Spatial Simulation , 2006 .

[36]  Muhammad Sahimi,et al.  Efficient Computational Strategies for Solving Global Optimization Problems , 2010, Computing in Science & Engineering.

[37]  Kenneth Moreland,et al.  The FFT on a GPU , 2003, HWWS '03.

[38]  Gregoire Mariethoz,et al.  The Direct Sampling method to perform multiple‐point geostatistical simulations , 2010 .

[39]  Steve Mann,et al.  Computer vision signal processing on graphics processing units , 2004, 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[40]  Jef Caers,et al.  Sequential simulation with patterns , 2005 .

[41]  Alexandre Boucher,et al.  Applied Geostatistics with SGeMS: A User's Guide , 2009 .

[42]  R. Bordley.  A Multiplicative Formula for Aggregating Probability Assessments , 1982 .