Stochastic generation of explicit pore structures by thresholding Gaussian random fields

Journal Article · Sat Nov 15 04:00:00 EST 2014 · Journal of Computational Physics

DOI:https://doi.org/10.1016/J.JCP.2014.07.046· OSTI ID:22382142

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. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.

• OSTI ID: 22382142
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 277; ISSN JCTPAH; ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
COMPUTERIZED SIMULATION
CONTROL
CORRELATIONS
FUNCTIONALS
KERNELS
LENGTH
MINKOWSKI SPACE
PORE STRUCTURE
POROUS MATERIALS
RANDOMNESS
SANDSTONES
SOLIDS
STOCHASTIC PROCESSES
TOPOGRAPHY
TOPOLOGY