Statistical fusion of two-scale images of porous media

The reconstruction of the architecture of void space in porous media is a challenging task, since porous media contain pore structures at multiple scales. Whereas past methods have been limited to producing samples with matching statistical behavior, the patterns of grey-level values in a measured sample actually say something about the unresolved details, thus we propose a statistical fusion framework for reconstructing high-resolution porous media images from low-resolution measurements. The proposed framework is based on a posterior sampling approach in which information obtained by low-resolution (MRI or X-ray) measurements is combined with prior models inferred from high-resolution microscopic data, typically 2D. In this paper, we focus on two-scale reconstruction tasks in which the measurements resolve only the large scale structures, leaving the small-scale to be inferred. The evaluation of the results generated by the proposed method shows the strong ability of the proposed method in reconstructing fine-scale structures positively correlated with the underlying ground truth. Comparing our method with the recent method of Okabe and Blunt [12], in which the measurements are also used in the reconstruction, we conclude that our method is more robust to the resolution of the measurement, and more closely matches the underlying fine-scale field.

[1]  Leonhard Held,et al.  Gaussian Markov Random Fields: Theory and Applications , 2005 .

[2]  Marios A. Ioannidis,et al.  Time-Dependent Diffusion and Surface-Enhanced Relaxation in Stochastic Replicas of Porous Rock , 2004 .

[3]  Paul W. Fieguth Hierarchical posterior sampling for images and random fields , 2003, Proceedings 2003 International Conference on Image Processing (Cat. No.03CH37429).

[4]  Christoph H. Arns,et al.  Developing a virtual materials laboratory , 2007 .

[5]  Paul W. Fieguth,et al.  Posterior Sampling of Scientific Images , 2006, ICIAR.

[6]  Edward R. Vrscay,et al.  Hierarchical annealing for synthesis of binary porous media images , 2005 .

[7]  J. G. Berryman Porous media—geometry and transports: by Pierre M. Adler. Butterworth-Heineman, Boston, MA (1992). 544 pp., US $95 , 1993 .

[8]  S. Torquato,et al.  Random Heterogeneous Materials: Microstructure and Macroscopic Properties , 2005 .

[9]  Anil K. Jain,et al.  Markov random fields : theory and application , 1993 .

[10]  Ioannidis,et al.  On the Geometry and Topology of 3D Stochastic Porous Media. , 2000, Journal of colloid and interface science.

[11]  Siegfried Stapf,et al.  NMR Imaging in Chemical Engineering , 2006 .

[12]  Marios A. Ioannidis,et al.  Pore size distribution in multiscale porous media as revealed by DDIF–NMR, mercury porosimetry and statistical image analysis , 2007 .

[13]  A. Radlinski,et al.  Angstrom-to-millimeter characterization of sedimentary rock microstructure. , 2004, Journal of colloid and interface science.

[14]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[15]  B. Inkson,et al.  Nanoscale tomography in materials science , 2007 .

[16]  Paul W. Fieguth,et al.  Hierarchical Sampling with Constraints , 2009, ICIAR.

[17]  S. Chakraborty,et al.  Induced pressure gradients due to entrance and exit effects in electroosmotically driven flows through nanopores within the continuum regime , 2008 .

[18]  C. Tsakiroglou,et al.  Dual‐porosity modelling of the pore structure and transport properties of a contaminated soil , 2008 .

[19]  S. Torquato,et al.  Reconstructing random media. II. Three-dimensional media from two-dimensional cuts , 1998 .

[20]  O. Torsæter,et al.  Stochastic reconstruction of particulate media from two-dimensional images. , 2002, Journal of colloid and interface science.

[21]  Gerhard Winkler,et al.  Image Analysis, Random Fields and Markov Chain Monte Carlo Methods: A Mathematical Introduction , 2002 .

[22]  Hilfer,et al.  Stochastic reconstruction of sandstones , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[23]  Yi-Qiao Song,et al.  Spatial Heterogeneity Length Scales in Carbonate Rocks , 2007 .

[24]  Anders Kaestner,et al.  Imaging and image processing in porous media research , 2008 .

[25]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[26]  D. Geman Random fields and inverse problems in imaging , 1990 .

[27]  M. Blunt,et al.  Pore space reconstruction of vuggy carbonates using microtomography and multiple‐point statistics , 2007 .

[28]  Paul W. Fieguth,et al.  Image Resolution Enhancement with Hierarchical Hidden Fields , 2009, ICIAR.

[29]  Simon K. Alexander,et al.  Hierarchical Annealing for Synthesis of Binary Images , 2009 .

[30]  Anthony Roberts Statistical reconstruction of three-dimensional porous media from two-dimensional images , 1997 .

[31]  Jun Yao,et al.  A new stochastic method of reconstructing porous media , 2007 .