Stochastic generation of explicit pore structures by thresholding Gaussian random fields

Jeffrey D. Hyman, C Larrabee Winter

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone.

Original languageEnglish (US)
Pages (from-to)16-31
Number of pages16
JournalJournal of Computational Physics
Volume277
DOIs
StatePublished - Nov 15 2014

Fingerprint

Pore structure
Topography
porosity
topography
Parameterization
Sandstone
Random variables
winter
thresholds
Topology
Geometry
random variables
sandstones
parameterization
functionals
solid phases
voids
topology
physical properties
geometry

Keywords

  • Direct numerical simulation
  • Minkowski functionals
  • Porous media
  • Stochastic methods

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)
  • Computer Science Applications

Cite this

Stochastic generation of explicit pore structures by thresholding Gaussian random fields. / Hyman, Jeffrey D.; Winter, C Larrabee.

In: Journal of Computational Physics, Vol. 277, 15.11.2014, p. 16-31.

Research output: Contribution to journalArticle

@article{f6e325ecb47145239d28b1ccb20f5f9f,
title = "Stochastic generation of explicit pore structures by thresholding Gaussian random fields",
abstract = "We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone.",
keywords = "Direct numerical simulation, Minkowski functionals, Porous media, Stochastic methods",
author = "Hyman, {Jeffrey D.} and Winter, {C Larrabee}",
year = "2014",
month = "11",
day = "15",
doi = "10.1016/j.jcp.2014.07.046",
language = "English (US)",
volume = "277",
pages = "16--31",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Stochastic generation of explicit pore structures by thresholding Gaussian random fields

AU - Hyman, Jeffrey D.

AU - Winter, C Larrabee

PY - 2014/11/15

Y1 - 2014/11/15

N2 - We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone.

AB - We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone.

KW - Direct numerical simulation

KW - Minkowski functionals

KW - Porous media

KW - Stochastic methods

UR - http://www.scopus.com/inward/record.url?scp=84906495432&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84906495432&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2014.07.046

DO - 10.1016/j.jcp.2014.07.046

M3 - Article

VL - 277

SP - 16

EP - 31

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -