1887
Volume 63, Issue 2
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images.

In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12156
2014-11-03
2024-04-26
Loading full text...

Full text loading...

References

  1. AdlerP.M., JacquinC.G. and ThovertJ.F.1992. The formation factor of reconstructed porous media. Journal of Water Resources Research28(6), 1571–1576.
    [Google Scholar]
  2. AndräH., CombaretN., DvorkinJ., GlattE., HanJ., KabelM. Et al.2013a. Digital rock physics benchmarks—Part I: Imaging and segmentation. Computers and Geosciences50, 25–32.
    [Google Scholar]
  3. AndräH., CombaretN., DvorkinJ., GlattE., HanJ., KabelM. Et al.2013b. Digital rock physics benchmarks—part II: Computing effective properties. Computers and Geosciences50, 33–43.
    [Google Scholar]
  4. ArnsC.H., KnackstedtM.A., PinczewskiW.V. and GarbocziE.J.2002. Computation of linear elastic properties from microtomographic images: methodology and agreement between theory and experiment. Geophysics67, 1396–1405.
    [Google Scholar]
  5. ArnsC.H., SheppardA.P., SokR.M., KnackstedtM.A. and PinczewskiW.V.2003. Relative permeability from tomographic images; effect of correlated heterogeneity. Journal of Petroleum Science and Engineering39(3–4), 247–259.
    [Google Scholar]
  6. ArnsC.H., KnackstedtM.A., PinczewskiW.V. and MartysN.S.2004. Virtual permeametry on microtomographic images. Journal of Petroleum Science and Engineering45(1–2)41–46.
    [Google Scholar]
  7. ArnsC.H., BaugetF., LimayeA., SakellariouA., SendenT.J., SheppardA.P. and SokR.M.2005. Pore‐Scale Characterization of Carbonates Using X‐Ray Microtomography. Society of Petroleum Engineering Journal10(4), 475–484.
    [Google Scholar]
  8. BatenburgK.J. and SijbersJ.2008. Optimal threshold selection for tomogram segmentation by projection distance minimization. IEEE Transaction in Medical Imaging28(5), 676–86.
    [Google Scholar]
  9. BelheineN., PlassiardJ.P., DonzéF.V., DarveF. and SeridiA.2009. Numerical simulation of drained triaxial test using 3d discrete element modeling. Computers and Geotechnics36 (1–2), 320–331.
    [Google Scholar]
  10. BieniekA. and MogaA.2000. An efficient watershed algorithm based on connected components. Pattern Recognition33, 907–916.
    [Google Scholar]
  11. BiswalB., OrenP.‐E., HeldR.J., BakkeS. and HilferR.Stochastic multiscale model for carbonate rocks. Physical Review E – Statistical, Nonlinear and Soft Matter . American Physical Society75(6).
    [Google Scholar]
  12. BohnR.B. and GarbocziE.J.2003. User manual for finite element and finite difference programs: a parallel version of NIST IR 6269. NIST Internal Report 6997.
  13. ChaoW.H, ChenY.Y, LinS.H., ShihY.Y. and TsangS.2009. Automatic segmentation of magnetic resonance images using a decision tree with spatial information, Computerized Medical Imaging and Graphics33(2), 111–121.
  14. ChenH., ChenS. and MatthaeusW.H.1992. Recovery of the Navier‐Stokes equations using a lattice‐gas Boltzmann method. Physical Review A45(8), 5339–5342.
    [Google Scholar]
  15. ClarkA.B.S., Al‐DayyaniT. and NurA.2009. Computations of porosity and permeability of sparic carbonate using multi‐scale CT images. Proceedings of the international symposium of the society of core analysts, SCA2009‐31.
  16. CundallP.A. and StrackO.D.L.1979. A discrete numerical model for granular assemblies. Geotechnique29, 47–65.
    [Google Scholar]
  17. DvorkinJ., NurA. and YinH.1994. Effective properties of cemented granular materials. Mechanics of Materials18(4), 351–366.
    [Google Scholar]
  18. EissaE.A. and KaziA.1988. Relation between static and dynamic Young's moduli of rocks. International Journal of Rock Mechanics and Mining Sciences and Geomechanics25(6), 479–482.
    [Google Scholar]
  19. FrijaJ., KervilerE., GéryS. and ZagdanskiA.M.1998. Computed radiography. Biomedicine and Pharmacotherapy52(2), 59–63.
    [Google Scholar]
  20. GarbocziE.J. and DayA.1995. Algorithm for computing the effective linear elastic properties of heterogeneous materials: Three dimensional results for composites with equal phase poisson ratios. Journal of the Mechanics and Physics of Solids43, 1349–1362.
    [Google Scholar]
  21. GelbJ., RothS., DongH., LiD., GuA., YunS. and YunW.2012. Non‐destructive local X‐Ray tomography for multi‐length scale analysis of carbonates: validations and observations. Society of Core Analysts, 24(2), 59.
    [Google Scholar]
  22. HollisC., VahrenkampV., TullS., MookerjeeA., TabernerC. and HuangY.2010. Pore system characterization in heterogeneous carbonates: An alternative approach to widely‐used rock‐typing methodologies. Marine and Petroleum Geology27(4), 772–793.
    [Google Scholar]
  23. HoltR.M., KjolaasJ., LarsenI., LiL., PillitteriA, G. and SonsteboE.F.2005. Comparison between controlled laboratory experiments and discrete particle simulations of the mechanical behaviour of rock. International Journal of Rock Mechanics and Mining Sciences42(7–8), 985–995.
    [Google Scholar]
  24. JouiniM.S. and VegaS.D.2011. Simulation of elastic properties in carbonates. The Leading Edge30, 1400–1407, doi:10.1190/1.3672485.
    [Google Scholar]
  25. KaczmarczykJ., DohnalikM. and ZalewskaJ.2010. Three‐Dimensional Pore Scale Fluid Flow Simulation Based on Computed microtomography carbonate rocks images. In Proceedings of the European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 Lisbon, Portugal.
  26. KalyaniS. and SwarupK.S.2011. Particle swarm optimization based K‐means clustering approach for security assessment in power systems. Expert Systems with Applications38 (9), 10839–10846.
    [Google Scholar]
  27. KangC.C., WangW.J. and KangC.H.2012. Image segmentation with complicated background by using seeded region growing. AEU – International Journal of Electronics and Communications66(9), 767–771.
    [Google Scholar]
  28. KeehmY. and MukerjiT.2004. Permeability and Relative Permeability from Digital Rocks: Issues on Grid Resolution and Representative Elementary Volume. Proceedings of the Annual Meeting of Society Exploration Geophysicists, Denver, USA. 10–15.
  29. KeehmY., MukerjiT. and NurA.2001. Computational rock physics at the pore scale: Transport properties and diagenesis in realistic pore geometries. The Leading Edge20(2), 180–183.
    [Google Scholar]
  30. KingM.S.1983. Static and dynamic elastic properties of rocks from the Canadian shield. International Journal of Rock Mechanics and Mining Sciences and Geomechanics20(5), 237–241.
    [Google Scholar]
  31. KnackstedtM.A., LathamS., MadadiM., SheppardA., VarslotT. and ArnsC.H.2009. Digital rock physics: 3D imaging of core material and correlations to acoustic and flow properties. The Leading Edge28(1), 28–33, doi: 10.1190/1.
    [Google Scholar]
  32. KoivuV., DecainM., GeindreauC., MattilaK., BlochJ.F. and KatajaM.2009. Transport properties of heterogeneous materials: Combining computerised X‐ray microtomography and direct numerical simulations. International Journal of Computational Fluid Dynamics23(10).
    [Google Scholar]
  33. LiL. and FjærE.2008. Investigation of the stress‐dependence of static and dynamic moduli of sandstones using a discrete element method. In Proceedings 42nd U.S. Rock Mechanics Symposium (USRMS), San Francisco, USA June 29 – July 2.
  34. MacQueenJ.B.1967. Some Methods for classification and Analysis of Multivariate Observations. Proceedings of 5th Berkeley Symposium on Mathematical Statistics and Probability. University of California Press. 281–297.
  35. MadonnaC., AlmqvistB.S.G. and SaengerE.H.2012. Digital rock physics: Numerical prediction of pressure‐dependent ultrasonic velocities using micro‐CT imaging. Geophysical Journal International189, 1475–1482.
    [Google Scholar]
  36. MaierR.S., BernardR.S. and GrunauD.W.1996. Boundary conditions for the lattice Boltzmann method. Physics of Fluids8(7) 1788–1801.
    [Google Scholar]
  37. MavkoG., MukerjiT. and DvorkinJ.2003. The Rock Physics Handbook, Cambridge University Press.
    [Google Scholar]
  38. MeyerF. and BeucherS.1990. Morphological segmentation. Journal of Visual Communication and Image Representation1(1), 21–46.
    [Google Scholar]
  39. OrenP.E. and BakkeS.2003. Reconstruction of Berea sandstone and pore‐scale modelling of wettability effects. Journal of Petroleum Science and Engineering39(3‐4), 177–199.
    [Google Scholar]
  40. OtsuN.1979. A threshold selection method from gray‐level histograms. Systems, Man, and Cybernetics Society9, 62–66.
    [Google Scholar]
  41. PanX., YuL. and KaoC.M., 2005. Spatial‐Resolution Enhancement in Computed Tomography.
  42. QifengD., YuefengS. and SullivanC.2011. Rock‐physics‐based carbonate pore type characterization and reservoir permeability heterogeneity evaluation, Upper San Andres reservoir, Permian Basin, west Texas. Journal of Applied Geophysics74(1), 8–18.
  43. QuintalB., FrehnerM., MadonnaC., TisatoN., KuteynikovaM. and SaengerE.H.2011. Integrated numerical and laboratory rock physics applied to seismic characterization of reservoir rocks. The Leading Edge30, 1360, DOI:10.1190/1.3672480.
    [Google Scholar]
  44. RambabuC. and ChakrabartiI.2007. An efficient immersion‐based watershed transform method and its prototype architecture. Journal of Systems Architecture53, 210–226.
    [Google Scholar]
  45. RedmondS.J. and HeneghanC.2007. A method for initialising the K‐means clustering algorithm using kd‐trees. Pattern Recognition Letters28(8), 965–973.
    [Google Scholar]
  46. RezaeiN. and ChatzisI.2011. Characterization of heterogeneities in porous media using constant rate air injection porosimetry. Journal of Petroleum Science and Engineering79(3–4), 113–124.
    [Google Scholar]
  47. RidlerT.W. and CalvardS.1978. Picture thresholding using an iterative selection method. IEEE Transaction on System, Man and Cybernetics, SMC‐8, 630–632.
  48. RosinP.L.2001. Unimodal thresholding. Pattern Recognition34(11), 2083–2096.
    [Google Scholar]
  49. SaengerE. H.2008. Numerical methods to determine effective elastic properties. International Journal of Engineering Science46(6), 598–605.
    [Google Scholar]
  50. SaengerE.H., EnzmannF., KeehmY. and SteebH.2011. Digital rock physics: Effect of fluid viscosity on effective elastic properties. Journal of Applied Geophysics74(4), 236–241.
    [Google Scholar]
  51. SchluterS., WellerU. and VogelH.J.2010. Segmentation of X‐ray microtomography images of soil using gradient masks. Computer and Geosciences36(10), 1246–1251.
    [Google Scholar]
  52. ShafeiB. and SteidlG.2012. Segmentation of images with separating layers by fuzzy c‐means and convex optimization. Journal of Visual Communication and Image Representation23 (4), 611–62.
    [Google Scholar]
  53. SokR., VarslotT. and GhousA.2010. Pore scale characterization of carbonates at multiple scales: integration of micro‐CT BSEM and FIBSEM. Petrophysics51(6), 379–387.
    [Google Scholar]
  54. SoneS., KasugaT., SakaiF.IzunoI. and OguchiM.1991. Digital Image Processing to Remove Blur from Linear Tomography of the Lung. Acta Radiologica32(5), 421.
    [Google Scholar]
  55. TremeauA. and BorelN.1997. A region growing and merging algorithm to color segmentation. Pattern Recognition30(7), 1191–1203.
    [Google Scholar]
  56. TurnerM., KnuefingL., ArnsC.H., SakellariouA., SendenT.J., SheppardA.P. et al 2004. Three‐dimensional imaging of multiphase flow in porous media. Physica A: Statistical Mechanics and its Applications339(1–2), 166–172.
    [Google Scholar]
  57. VillainN., GoussardY., IdierJ. and AllainM.2003. Three‐dimensional edge‐preserving image enhancement for computed tomography. Medical Imaging Journal22(10), 1275–1287.
    [Google Scholar]
  58. VincentL. and SoilleP.1991. Watersheds in Digital Spaces: An efficient algorithm based on immersion simulations. Pattern Analysis and Machine Intelligence13(6), 583–598.
    [Google Scholar]
  59. VogelH.J. and KretzschmarA.1996. Topological characterization of pore space in soil. Sample preparation and digital image processing. Geoderma73, 29–38.
    [Google Scholar]
  60. WhiteJ.A., BorjaR.I. and FredrichJ.T.2006. Calculating the effective permeability of sandstone with multiscale lattice Boltzmann/finite element simulations. Acta Geotechnica1(4), 195–209.
    [Google Scholar]
  61. WildenschildD. and SheppardA.P.2012. X‐ray imaging and analysis techniques for quantifying pore‐scale structure and processes in subsurface porous medium systems. Invited review paper for 35th issue of Advances in Water resources.
  62. WildenschildD., CulliganK.A. and ChristensenB.S.B.2004. Application of x‐ray microtomography to environmental fluid flow problems. Proceeding of SPIE 5535, Developments in X‐Ray Tomography IV, 432‐441. doi 10.1117/12.559056.
  63. YeomT., CalocaI.Z. and ChambersF.W.2008. Lattice Boltzmann Computations of Micro Channel and Micro Orifice Flows. In proceedings American Society of Mechanical Engineering Conference, Boston, USA, 907‐919. DOI:10.1115/IMECE2008-67449
  64. ZouQ. and HeX.1997. On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Physics of Fluids9(6), 1591.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12156
Loading
/content/journals/10.1111/1365-2478.12156
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): carbonate; elasticity; Microtomography; multiscale analysis; permeability; porosity

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error