Abstract
The Sequential Spectral Turning Bands Method (S-STBM) builds Gaussian random fields (GRF) calibrated to desired response functions. An interesting application of S-STBM concerns the simulation of GRF subject to inequality constraints. S-STBM works by choosing the phase of each cosine function of the STBM algorithm instead of perturbating nodes of the GRF many thousand times using conditional distributions as in Gibbs sampler. Each chosen phase increasingly constrains the nodes to the desired inequalities. A method based on the sequential Gaussian simulation is introduced to accelerate convergence at the end of the process. Examples shown compare S-STBM approach to Gibbs sampler. Orders of magnitude reduction in computation time is achieved with our spectral method. Furthermore, examples show that the phase selection has no significant influence on the spatial correlation. Our approach is easily generalized to pluriGaussian simulations. Compared to Gibbs sampler, S-STBM is not limited to small systems (no memory limitation) and its complexity of O(n) makes it an efficient tool to simulate large GRF subject to inequality constraints.
Similar content being viewed by others
Availability of data and material
Not applicable
References
Allard D (1994) Simulating a geological lithofacies with respect to connectivity information using the truncated Gaussian model. In: Geostatistical Simulations, Springer Netherlands, pp 197–211. https://doi.org/10.1007/978-94-015-8267-4_16
Allard D, D’Or D, Biver P, Froidevaux R (2012) Non-parametric diagrams for pluriGaussian simulations of lithologies. In: Abrahamsen P, Hauge R, Kolbjornsen O (eds) vol Oslo 2012. Springer
Armstrong M, Galli A, Beucher H, Loc’h G, Renard D, Doligez B, Eschard R, Geffroy F (2011) PluriGaussian simulations in geosciences. Springer Berlin Heidelberg. https://doi.org/10.1007/978-3-642-19607-2
Arroyo D, Emery X (2017) Spectral simulation of vector random fields with stationary Gaussian increments in d-dimensional Euclidean spaces. Stoch Environ Res Risk Assess 31(7):1583–1592. https://doi.org/10.1007/s00477-016-1225-7
Beucher H, Renard D (2016) Truncated Gaussian and derived methods. Comptes Rendus Geosci 348(7):510–519. https://doi.org/10.1016/j.crte.2015.10.004
Blévec TL, Dubrule O, John CM, Hampson GJ (2018) Geostatistical modelling of cyclic and rhythmic facies architectures. Math Geosci 50(6):609–637. https://doi.org/10.1007/s11004-018-9737-y
Botev ZI (2016) The normal law under linear restrictions: simulation and estimation via minimax tilting. J R Stat Soc Ser B (Stat Methodol) 79(Part 1):1–24. https://doi.org/10.1111/rssb.12162
Chilès JP, Delfiner P (2012) Geostatistics: modeling spatial uncertainty. John Wiley & Sons, Inc., Hoboken https://doi.org/10.1002/9781118136188
Desassis N, Renard D, Pereira M, Freulon X (2019) PluriGaussian simulations with the stochastic partial differential equation (SPDE) approach. In: Petroleum Geostatistics 2019, EAGE Publications BV. https://doi.org/10.3997/2214-4609.201902174
Deutsch JL, Deutsch CV (2014) A multidimensional scaling approach to enforce reproduction of transition probabilities in truncated pluriGaussian simulation. Stoch Environ Res Risk Assess 28(3):707–716. https://doi.org/10.1007/s00477-013-0783-1
Emery X (2004) Testing the correctness of the sequential algorithm for simulating Gaussian random fields. Stoch Environ Res Risk Assess 18(6):401–413. https://doi.org/10.1007/s00477-004-0211-7
Emery X, Lantuéjoul C (2006) TBSIM: A computer program for conditional simulation of three-dimensional Gaussian random fields via the turning bands method. Comput Geosci 32(10):1615–1628. https://doi.org/10.1016/j.cageo.2006.03.001
Emery X, Peláez M (2011) Assessing the accuracy of sequential Gaussian simulation and cosimulation. Comput Geosci 15(4):673–689. https://doi.org/10.1007/s10596-011-9235-5
Emery X, Arroyo D, Peláez M (2014) Simulating large Gaussian random vectors subject to inequality constraints by Gibbs sampling. Math Geosci 46(3):265–283. https://doi.org/10.1007/s11004-013-9495-9
Emery X, Arroyo D, Porcu E (2015) An improved spectral turning-bands algorithm for simulating stationary vector Gaussian random fields. Stoch Environ Res Risk Assess 30(7):1863–1873. https://doi.org/10.1007/s00477-015-1151-0
Freulon X, de Fouquet C (1993) Conditioning a Gaussian model with inequalities. In: Quantitative Geology and Geostatistics, Springer Netherlands, pp 201–212. https://doi.org/10.1007/978-94-011-1739-5_17
Galli A, Beucher H, Le Loc’h G, Doligez B, Group H (1994) The pros and cons of the truncated Gaussian method. In: Geostatistical Simulations, Springer Netherlands, pp 217–233. https://doi.org/10.1007/978-94-015-8267-4_18
Galli A, Le Loc’h G, Geffroy F, Eschard R (2006) An application of the truncated pluri-Gaussian method for modeling geology. In: Stochastic Modeling and Geostatistics, American Association of Petroleum Geologists, pp 109–122. https://doi.org/10.1306/1063810ca53226
Hu L, Zhao Y, Liu Y, Scheepens C, Bouchard A (2013) Updating multipoint simulations using the ensemble Kalman filter. Comput Geosci 51:7–15. https://doi.org/10.1016/j.cageo.2012.08.020
Jakeman AJ, Barreteau O, Hunt RJ, Rinaudo JD, Ross A (eds) (2016) Integrated Groundwater Management. Springer International Publishing. https://doi.org/10.1007/978-3-319-23576-9
Lantuéjoul C (2002) Geostatistical Simulation. Springer Berlin Heidelberg, https://doi.org/10.1007/978-3-662-04808-5
Lantuéjoul C, Desassis N (2012) Simulation of a Gaussian random vector: a propagative version of the Gibbs sampler. In: Ninth international geostatistics Congress, Oslo, http://geostats2012.nr.no/pdfs/1747181.pdf
Lauzon D, Marcotte D (2019) Calibration of random fields by FFTMA-SA. Comput Geosci 127:99–110. https://doi.org/10.1016/j.cageo.2019.03.003
Lauzon D, Marcotte D (2020) Calibration of random fields by a sequential spectral turning bands method. Comput Geosci 135:104390. https://doi.org/10.1016/j.cageo.2019.104390
Lindgren F, Rue H, Lindström J (2011) An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach. J R Stat Soc Ser B (Stat Methodol) 73(4):423–498. https://doi.org/10.1111/j.1467-9868.2011.00777.x
Luenberger DG, Ye Y (2008) Linear and nonlinear programming (international series in operations research & management science). Springer, URL https://www.ebook.de/de/product/19292757/david_g_luenberger_yinyu_ye_linear_and_nonlinear_programming.html
Madani N, Emery X (2014) Simulation of geo-domains accounting for chronology and contact relationships: application to the Río Blanco copper deposit. Stoch Environ Res Risk Assess 29(8):2173–2191. https://doi.org/10.1007/s00477-014-0997-x
Marcotte D (2015) TASC3D: A program to test the admissibility in 3D of non-linear models of coregionalization. Comput Geosci 83:168–175. https://doi.org/10.1016/j.cageo.2015.07.012
Marcotte D (2016) Spatial turning bands simulation of anisotropic non-linear models of coregionalization with symmetric cross-covariances. Comput Geosci 89:232–238. https://doi.org/10.1016/j.cageo.2016.01.004
Marcotte D, Allard D (2017) Half-tapering strategy for conditional simulation with large datasets. Stoch Environ Res Risk Assess 32(1):279–294. https://doi.org/10.1007/s00477-017-1386-z
Marcotte D, Allard D (2018) Gibbs sampling on large lattice with GMRF. Comput Geosci 111:190–199. https://doi.org/10.1016/j.cageo.2017.11.012
Nussbaumer R, Mariethoz G, Gloaguen E, Holliger K (2018) Which path to choose in sequential Gaussian simulation. Math Geosci 50(1):97–120. https://doi.org/10.1007/s11004-017-9699-5
Ravalec-Dupin ML, Hu LY (2005) Gradual deformation of boolean simulations. In: Geostatistics Banff 2004, Springer Netherlands, pp 939–948. https://doi.org/10.1007/978-1-4020-3610-1_98
Refsgaard JC, Christensen S, Sonnenborg TO, Seifert D, Højberg AL, Troldborg L (2012) Review of strategies for handling geological uncertainty in groundwater flow and transport modeling. Adv Water Resour 36:36–50. https://doi.org/10.1016/j.advwatres.2011.04.006
Rezaee H, Marcotte D (2018) Calibration of categorical simulations by evolutionary gradual deformation method. Comput Geosci 22(2):587–605. https://doi.org/10.1007/s10596-017-9711-7
Rue H, Tjelmeland H (2002) Fitting Gaussian Markov random fields to Gaussian fields. Scand J Stat 29(1):31–49. https://doi.org/10.1111/1467-9469.00058
Safikhani M, Asghari O, Emery X (2017) Assessing the accuracy of sequential Gaussian simulation through statistical testing. Stoch Environ Res Risk Assess 31(2):523–533. https://doi.org/10.1007/s00477-016-1255-1
Shinozuka M (1971) Simulation of multivariate and multidimensional random processes. J Acoust Soc Am 49(1B):357–368. https://doi.org/10.1121/1.1912338
Silva DSF, Deutsch CV (2017) Multiple imputation framework for data assignment in truncated pluri-Gaussian simulation. Stoch Environ Res Risk Assess 31(9):2251–2263. https://doi.org/10.1007/s00477-016-1309-4
Talebi H, Asghari O, Emery X (2013) Application of pluriGaussian simulation to delineate the layout of alteration domains in Sungun copper deposit. Open Geosci. https://doi.org/10.2478/s13533-012-0146-3
Tran TT (1994) Improving variogram reproduction on dense simulation grids. Comput Geosci 20(7–8):1161–1168. https://doi.org/10.1016/0098-3004(94)90069-8
Zhang Y (2015) Ensemble methods of data assimilation in porous media flow for non-Gaussian prior probability density. PhD thesis, The University of Bergen
Funding
This research was made possible by a Canada Graduate Scholarships-Master’s (CGS M) from the Natural Sciences and Engineering Research Council of Canada (NSERC) and a B2X Scholarship from the Fonds de recherche du Québec – Nature et technologies (FRQNT) of D. Lauzon and a NSERC grant (RGPIN-2015-06653) of D. Marcotte.
Author information
Authors and Affiliations
Contributions
Not applicable
Corresponding author
Ethics declarations
Conflicts of interest
The authors declare that they have no conflict of interest.
Code availability
The Matlab computer codes are available at https://github.com/Danlauz/SSTBM-as-an-alternative-to-Gibbs-Sampler.git
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Lauzon, D., Marcotte, D. The sequential spectral turning band simulator as an alternative to Gibbs sampler in large truncated- or pluri- Gaussian simulations. Stoch Environ Res Risk Assess 34, 1939–1951 (2020). https://doi.org/10.1007/s00477-020-01850-9
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00477-020-01850-9