Speeding up the flash calculations in two-phase compositional flow simulations - The application of sparse grids

Flash calculations have become a performance bottleneck in the simulation of compositional flow in subsurface reservoirs. We apply a sparse grid surrogate model to substitute the flash calculation and thus try to remove the bottleneck from the reservoir simulation. So instead of doing a flash calculation in each time step of the simulation, we just generate a sparse grid approximation of all possible results of the flash calculation before the reservoir simulation. Then we evaluate the constructed surrogate model to approximate the values of the flash calculation results from this surrogate during the simulations. The execution of the true flash calculation has been shifted from the online phase during the simulation to the offline phase before the simulation. Sparse grids are known to require only few unknowns in order to obtain good approximation qualities. In conjunction with local adaptivity, sparse grids ensure that the accuracy of the surrogate is acceptable while keeping the memory usage small by only storing a minimal amount of values for the surrogate. The accuracy of the sparse grid surrogate during the reservoir simulation is compared to the accuracy of using a surrogate based on regular Cartesian grids and the original flash calculation. The surrogate model improves the speed of the flash calculations and the simulation of the whole reservoir. In an experiment, it is shown that the speed of the online flash calculations is increased by about 2000 times and as a result the speed of the reservoir simulations has been enhanced by 21 times in the best conditions.

[1]  Dirk Pflüger,et al.  Fastsg: A Fast Routines Library for Sparse Grids , 2012, ICCS.

[2]  Dongbin Xiu,et al.  Characterization of discontinuities in high-dimensional stochastic problems on adaptive sparse grids , 2011, J. Comput. Phys..

[3]  Vasile Gradinaru,et al.  Fourier transform on sparse grids: Code design and the time dependent Schrödinger equation , 2007, Computing.

[4]  Dirk Pflüger,et al.  Spatially Adaptive Sparse Grids for High-Dimensional Problems , 2010 .

[5]  Robert D. Falgout,et al.  hypre: A Library of High Performance Preconditioners , 2002, International Conference on Computational Science.

[6]  Wei Yan,et al.  Comparison of Two Methods for Speeding up Flash Calculations in Compositional Simulations , 2011, ANSS 2011.

[7]  Michael Griebel,et al.  Adaptive sparse grid multilevel methods for elliptic PDEs based on finite differences , 1998, Computing.

[8]  Abbas Firoozabadi,et al.  Thermodynamics of Hydrocarbon Reservoirs , 1999 .

[9]  Dirk Pflüger,et al.  A Parallel and Distributed Surrogate Model Implementation for Computational Steering , 2012, 2012 11th International Symposium on Parallel and Distributed Computing.

[10]  K. Hallatschek.  Fouriertransformation auf dünnen Gittern mit hierarchischen Basen , 1992 .

[11]  M. Michelsen.  The isothermal flash problem. Part I. Stability , 1982 .

[12]  Michael L. Michelsen,et al.  Increasing the Computational Speed of Flash Calculations With Applications for Compositional, Transient Simulations , 2006 .

[13]  Fabio Nobile,et al.  A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..

[14]  L. Dumas,et al.  Global Optimization Based on Sparse Grid Surrogate Models for Black-box Expensive Functions , 2014 .

[15]  Russell T. Johns,et al.  A Rapid and Robust Method To Replace Rachford-Rice in Flash Calculations , 2007 .

[16]  Michael L. Michelsen,et al.  Simplified flash calculations for Cubic Equations of State , 1986 .

[17]  Benjamin Peherstorfer,et al.  Fast Insight into High-Dimensional Parametrized Simulation Data , 2012, 2012 11th International Conference on Machine Learning and Applications.

[18]  M. Michelsen.  The isothermal flash problem. Part II. Phase-split calculation , 1982 .

[19]  James Nance,et al.  Reaction Path Following with Sparse Interpolation. , 2014, Journal of chemical theory and computation.

[20]  Michael Griebel,et al.  Fast Discrete Fourier Transform on Generalized Sparse Grids , 2014 .

[21]  Benjamin Peherstorfer,et al.  Density Estimation with Adaptive Sparse Grids for Large Data Sets , 2014, SDM.

[22]  Jochen Garcke,et al.  Sparse Grids in a Nutshell , 2012 .