Abstract
The reconstruction of porous media is of great importance in predicting fluid transport properties, which are widely used in various fields such as catalysis, oil recovery, medicine and aging of building materials. The real three-dimensional structural data of porous media are helpful to describe the irregular topologic structures of porous media. By using multiple-point statistics (MPS) to extract the characteristics of real porous media acquired from micro computed tomography (micro-CT) scanning, the probabilities of structural characteristics of pore spaces are obtained first, and then reproduced in the reconstructed regions. One solution to overcome the anisotropy of training images is to use real 3D volume data as a training image (TI). The CPU cost and memory burden brought up by 3D simulations can be reduced greatly by selecting the optimal multiple-grid template size that is determined by the entropy of a TI. Moreover, both soft data and hard data are integrated in MPS simulation to improve the accuracy of reconstructed images. The variograms and permeabilities, computed by lattice Boltzmann method, of the reconstructed images and the target image obtained from real volume data are compared, showing that the structural characteristics of reconstructed porous media using our method are similar to those of real volume data.
Similar content being viewed by others
References
Auzerais FM, Dunsmuir J, Ferreol BB, Martys N, Olson J, Ramakrishnan TS, Rothman DH, Schwartz LM (1996) Transport in sandstone: a study based on three dimensional microtomography. Geophys Res Lett 23:705–708. doi:10.1029/96GL00776
Bear J (1972) Dynamics of fluids in porous media. American Elsevier Publishing Company, Inc., New York
Blunt MJ, Bijeljic B, Dong H, Gharbi O, Iglauer S, Mostaghimi P, Paluszny A, Pentland C (2013) Pore-scale imaging and modelling. Adv Water Resour 51:197–216
Comunian A, Renard P, Straubhaar J (2012) 3D multiple-point statistics simulation using 2D training images. Comput Geosci 40:49–65
Dong H (2007) Micro-CT imaging and pore network extraction. Ph.D. Dissertation, Imperial College London, London, UK
Flannery BP, Deckman HW, Roberge WG, Damico KL (1987) Three-dimensional X-ray microtomography. Science 237:1439–1444
Guardiano F, Srivastava RM (1993) Multivariate geostatistics: beyond bivariate moments. In: Soares A (ed) Geostatistics-Troia, vol 1. Kluwer, Dordrecht, pp 133–144
Hajizadeh A, Safekordi A, Farhadpour FA (2011) A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images. Adv Water Resour 34:1256–1267
Honarkhah M (2011) Stochastic simulation of patterns using distance-based pattern modeling. Ph.D. Dissertation, Stanford University, Stanford, CA
Honarkhah M, Caers J (2010) Stochastic simulation of patterns using distance-based pattern modeling. Math Geosci 42:487–517
Hsu KC, Chen KC (2010) Multiscale flow and transport model in three-dimensional fractal porous media. Stoch Environ Res Risk Assess 24(7):1053–1065
Huang T, Li X, Zhang T, Lu DT (2013a) GPU-accelerated Direct Sampling method for multiple-point statistical simulation. Comput Geosci 57:13–23
Huang T, Lu DT, Li X, Wang L (2013b) GPU-based SNESIM implementation for multiple-point statistical simulation. Comput Geosci 54:75–87
Journel AG (2002) Combining knowledge from diverse data sources: an alternative to traditional data independence hypothesis. Math Geol 34(5):573–596
Mariethoz G (2010) A general parallelization strategy for random path based geostatistical simulation methods. Comput Geosci 36:953–958
Okabe H, Blunt MJ (2004) Prediction of permeability for porous media reconstructed using multiple-point statistics. Phys Rev E 70(066135):1–10
Okabe H, Blunt MJ (2005) Pore space reconstruction using multiple-point statistics. J Petrol Sci Eng 46:121–137
Remy N, Boucher A, Wu JB (2009) Applied geostatistics with SGeMS: a users’ guide. Cambridge University Press, New York
Riva M, Guadagnini L, Guadagnini A (2010) Effects of uncertainty of lithofacies, conductivity and porosity distributions on stochastic interpretations of a field scale tracer test. Stoch Environ Res Risk Assess 24(7):955–970
Robert SM, Robert SB, Dary WG (1996) Boundary conditions for the lattice boltzmann method. Phys Fluids 8(7):1788–1800
Shannon CE (1948) A mathematical theory of communication. Bell Syst Tech J 27:379–423
Straubhaar J, Renard P, Mariethoz G, Froidevaux R, Besson O (2011) An improved parallel multiple-point algorithm using a list approach. Math Geosci 43(3):305–328
Straubhaar J, Walgenwitz A, Renard P (2013) Parallel multiple-point statistics algorithm based on list and tree structures. Math Geosci 45(2):131–147
Strebelle SB (2000) Sequential Simulation drawing structures from training images. Ph.D. Dissertation, Stanford University, Stanford, CA
Strebelle SB (2002) Conditional simulation of complex geological structures using multiple-point statistics. Math Geol 34(1):1–21
Tartakovsky AM (2010) Lagrangian simulations of unstable gravity-driven flow of fluids with variable density in randomly heterogeneous porous media. Stoch Environ Res Risk Assess 24(7):993–1002
Tran TT (1994) Improving variogram reproduction on dense simulation grids. Comput Geosci 20(7–8):1161–1168
Wildenschild D, Sheppard AP (2013) X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems. Adv Water Resour 51:217–246
Wu JB (2007) 4D seismic and multiple-point pattern data integration using geostatistics. Ph.D. Dissertation, Stanford University, Stanford, CA
Zhang, TF (2006) Filter-based training pattern classification for spatial pattern simulation. Ph.D. Dissertation, Stanford University, Stanford, CA
Acknowledgments
This work is supported by National Program on Key Basic Research Project of China (973 Program, No. 2011CB707305), the National Science and Technology Major Project (No. 2011ZX05009-006), CAS Strategic Priority Research Program (XDB10030402), Shanghai Municipal Natural Science Foundation (No. 11ZR1413700, 12ZR1412000), Talented People Introduction Foundation of Shanghai University of Electric Power (No. K2012-004, K-2013-019, K2014-020), the Excellent University Young Teachers Training Program of Shanghai Municipal Education Commission (No. ZZsdl12002, ZZsd113015), and Software Service Engineering Key Discipline Construction (fourth stage) of Shanghai Second Polytechnic University (No. XXKZD1301).
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Zhang, T., Du, Y., Huang, T. et al. Reconstruction of porous media using multiple-point statistics with data conditioning. Stoch Environ Res Risk Assess 29, 727–738 (2015). https://doi.org/10.1007/s00477-014-0947-7
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00477-014-0947-7