MicroCT Analysis

  • Science Problem
  • Approach
  • Results
  • Impact
  • Future Developments
  • Publications
  • Science Problem

    Demanding for tools to better understand processes in fluid-rock systems related to geologic sequestration of CO2 by using quantitative image analysis methods for automated measurements of experimental data as the calculation of material interfaces, quantification of CO2 during flow, recovery of material structures, etc. Other activities involve material interface reconstruction in simulation and deriving measurements (wetted surface area etc.) and methods for tracking boundaries to analyze pore region evolution.

    Permeability (K) is a fundamental measure for Earth scientists to determine the fluid flow through a porous material, and standard equations to estimate K lacks pore geometry. Previous work by Hazen included empirical formulation of K for saturated sands based on empirical coefficients and particle size. The main problem of this approach is that the particles are supposed to be quite homogeneous and to present fine grains. First Kozeny, then Carman proposed improvements to permeability calculation and postulated a new semiempirical, semitheoretical formulation known as Kozeny-Carman (KC) - this model considers geometry by including surface area per unit volume of particles in the formula. Further improvements are necessary to cope with limitations of KC formula, e.g. KC is only valid for laminar flow. Moreover, permeability is often associated to hydraulic conductivity because the water properties, e.g. viscosity, density, are well known and instruments use the flow of water through the material to calculate permeability. More about instruments . Information about permeability importance to oil industry.

    In a nutshell, our research consists in building algorithms that takes the INPUT: hundreds/thousands of slices, encoding tridimensional represention of solid objects and provides the OUTPUT: measurements of material porosity, permeability estimates, granularity of material, visualization of structures as pore network


    One of the most challenge questions posed to material sciences is to determine permeability of materials with complex structure. Can we invent new formulations to calculate permeability regarding geometrical information about the structure?
    Synthetic beads: microCT images of porous materials are seldom composed by homogeneous particles, and they also present artifacts that can interfere in the calculation of permeability. In order to isolate permeability models from issues associated to CT imaging and granularity, we generate synthetic datasets of randomly distributed sphere packing using numerical simulation of Stokes flow. The simulation algorithm creates numerical constructions of jammed backed bead beds, with identical and nonoverlapping spheres, and following algorithms proposed by Salvatore Torquato, later modified by Todd Weisgraber (LLNL). Link to further information and C++ code used in our experiments. The use of synthetic beads is only the first step toward describing porous media, using a new descriptor based on topological analysis: max-flow curves

    Pore network:

    Derived measures: a) effective pore radii: radius of a sphere of equivalent volume; b) effective throat radii = radius of a circle of equivalente area

    Fluid permeability curves: More about permeability curves


    The extraction of Reeb graphs through which we represent the pore network allowed the extraction two main descriptors: flow curve and capacity.
    Synthetic beads:

    Flow diagram

  • Analysis of 12 synthetic bead packed columns - download the source data used to generate columns HERE

  • Created ImageJ plugin to transform the output files from the C++ code into an image stack. Our code is HERE. Why JAVA?

  • Amazing visualizations for better scientific exploration - check Vis/Analytics website

    Experimental data:

  • Impact

  • Check this work in CRD News [HERE]

  • Software accessible to experts in experimental data (and not necessarily programmers) [HERE]
  • Future Developments


  • Ushizima, D.M., Ajo-Franklin, J., Macdowell, A., Morozov, D., Nico, P., Parkinson, Bethel E.W, Sethian J.A., "Statistical segmentation and porosity quantification of 3D x-ray microtomography", in SPIE Optics and Photonics: XXXIV Applications of Digital Image Processing, Vol.8135-1, pp.1-14, Aug 2011, San Diego, CA, USA. LBNL publication number pending.[pdf]

  • Ushizima, D.M., Weber, G.H., Ajo-Franklin, J., Kim, Y., Macdowell, A., Morozov, D., Nico, P., Parkinson, D., Trebotich, D., Wan, J., and Bethel E.W., "Analysis and visualization for multiscale control of geologic CO2", in: Journal of Physics: Conference Series, Proceedings of SciDAC 2011, July 2011, Denver, CO, USA. LBNL publication number pending.[pdf]
  • Use Quant-CT