1887

Abstract

Summary

We generate independent Gaussian random variables on a regular grid and use a spatial filter to smooth the independent random variables to obtain a spatially correlated Gaussian random field. The FFT is used to speed up the smoothing since convolution is a simple cell by-cell multiplication in the Fourier domain. A representation of the spatial convolution filter in the Fourier domain is efficiently obtained from the FFT of any stationary correlation function. Since FFT is cyclic, the grid must be padded to ensure that opposite sides are uncorrelated. The size of the padding is discussed in detail. Most standard covariance functions fail to be positive definite on finite cyclic domains. This causes striping artifacts in the final simulated realizations and failure to meet statistical properties such as variogram reproduction in the simulated realizations. These problems are addressed and solutions are provided to ensure near perfect statistical properties of the generated realizations. The method is fast and can generate a hundred million grid cell realization in approximately 1.5 minutes on a standard laptop PC. The method scales approximately linearly in the number of grid cells.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201802134
2018-09-03
2024-04-24
Loading full text...

Full text loading...

References

  1. Dietrich, C.R. and Newsam, G.N.
    [1993] A fast and exact method for multidimensional gaussian stochastic simulations. Water Resources Research, 29(8), 2861–2869.
    [Google Scholar]
  2. [1997] Fast and Exact Simulation of Stationary Gaussian Processes through Circulant Embedding of the Covariance Matrix. SIAM Journal on Scientific Computing, 18(4), 1088–1107.
    [Google Scholar]
  3. Frigo, M. and Johnson, S.G.
    [2005] The Design and Implementation of FFTW3. Proceedings of the IEEE, 93(2), 216–231.
    [Google Scholar]
  4. Gneiting, T, Ševčíková, H., Percival, D.B., Schlather, M. and Jiang, Y.
    [2006] Fast and Exact Simulation of Large Gaussian Lattice Systems in R2: Exploring the Limits. Journal of Computational and Graphical Statistics, 15(3), 483–501.
    [Google Scholar]
  5. Goff, J.H. and Jennings, J.WJr.,
    [1999] Improvement on Fourier-based unconditional and conditional simulations for band limited fractal (von Kármán) statistical models. Math. Geol., 31(6), 627–649.
    [Google Scholar]
  6. Jacquemyn, C., Melnikova, Y., Jackson, M., Hampson, G. and John, C.
    [2016] Geologic Modelling Using Parametric NURBS Surfaces. In: ECMOR XV, 15th European Conference on the Mathematics of Oil Recovery.Amsterdam, Netherlands.
    [Google Scholar]
  7. Lantuéjoul, C.
    [2002] Geostatistical Simulation: Models and Algorithms.Springer-Verlag Inc., Berlin.
    [Google Scholar]
  8. Minniakhmetov, I.R. and Pergament, A.H.
    [2012] A Spectal Approach to Conditional Simulation. In: 74th EAGE Conference and Exhibition incorporating SPE EUROPEC 2012. Society of Petroleum Engineers, Copenhagen, Denmark.
    [Google Scholar]
  9. Pardo-Igúzquiza, E. and Chica-Olmo, M.
    [1993] The Fourier Integral Method: An efficient spectral method for simulation of random fields. Math. Geol., 25(2), 177–217.
    [Google Scholar]
  10. Ravalec, M.L., Noetinger, B. and Hu, L.Y.
    [2000] The FFT Moving Average (FFT-MA) Generator: An Efficient Numerical Method for Generating and Conditioning Gaussian Simulations. Mathematical Geology, 32(6), 701–723.
    [Google Scholar]
  11. Ripley, B.D.
    [1987] Stochastic Simulation.John Wiley & Sons, New York.
    [Google Scholar]
  12. Vigsnes, M., Kolbjørnsen, O., Hauge, V.L., Dahle, P. and Abrahamsen, P.
    [2017] Fast and Accurate Approximation to Kriging Using Common Data Neighborhoods. Mathematical Geosciences, 49(5), 619–634.
    [Google Scholar]
  13. Wood, A.T.A. and Chan, G.
    [1994] Simulation of Stationary Gaussian Processes in [0, 1]d. Journal of Computational and Graphical Statistics, 3(4), 409–432.
    [Google Scholar]
  14. Yao, T., Calvert, C., Jones, T., Bishop, G., Ma, Y. and Foreman, L.
    [2006] Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling. Mathematical Geology, 38(1), 51–62.
    [Google Scholar]
  15. Zaytsev, V.N., Biver, P., Wackernagel, H. and Allard, D.
    [2016] Geostatistical Simulations on Irregular Reservoir Models. In: ECMOR XV, 15th European Conference on the Mathematics of Oil Recovery.Amsterdam, Netherlands.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201802134
Loading
/content/papers/10.3997/2214-4609.201802134
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error