Imaging and image processing in porous media research

European Journal of Soil Science, 2013

Water Resources Research, 2013

1] Lattice Boltzmann flow simulations provide a physics-based means of estimating intrinsic permeability from pore structure and accounting for inertial flow that leads to departures from Darcy's law. Simulations were used to compute intrinsic permeability where standard measurement methods may fail and to provide better understanding of departures from Darcy's law under field conditions. Simulations also investigated resolution issues. Computed tomography (CT) images were acquired at 0.8 mm interscan spacing for seven samples characterized by centimeter-scale biogenic vuggy macroporosity from the extremely transmissive sole-source carbonate karst Biscayne aquifer in southeastern Florida. Samples were as large as 0.3 m in length; 7-9 cm-scale-length subsamples were used for lattice Boltzmann computations. Macroporosity of the subsamples was as high as 81%. Matrix porosity was ignored in the simulations. Non-Darcy behavior led to a twofold reduction in apparent hydraulic conductivity as an applied hydraulic gradient increased to levels observed at regional scale within the Biscayne aquifer; larger reductions are expected under higher gradients near wells and canals. Thus, inertial flows and departures from Darcy's law may occur under field conditions. Changes in apparent hydraulic conductivity with changes in head gradient computed with the lattice Boltzmann model closely fit the Darcy-Forchheimer equation allowing estimation of the Forchheimer parameter. CT-scan resolution appeared adequate to capture intrinsic permeability ; however, departures from Darcy behavior were less detectable as resolution coarsened. (2013), Evaluation of permeability and non-Darcy flow in vuggy macroporous limestone aquifer samples with lattice Boltzmann methods, Water Resour. Res., 49,

Optics Express

The demand for high resolution neutron imaging has been steadily increasing over the past years. The number of facilities offering cutting edge resolution is however limited, due to (i) the design complexity of an optimized device able to reach a resolution in the order of ≈ 10 µm and (ii) limitations in available neutron flux. Here we propose a simple addition, based on a Fibre Optics Taper (FOT), that can be easily attached to an already existing scintillator-camera imaging detector in order to efficiently increase its spatial resolution and hence boost the capability of an instrument into high resolution applications.

Water Resources Research

Water drainage from porous media is a highly dynamic process often marked by rapid pistonlike air invasion events at the front and other rapid interfacial reconfigurations. Liquid phase entrapped behind the moving front drains at significantly slower rates often via gravity driven flow through corners and crevices. This distribution of slowly draining residual water phase determines the plant available water and biological functioning of soils. The study aims to determine the conditions for the flow regime transition from piston-like invasion at a drainage front to slower corner dominated flow at the pore and sample scale. This transition was observed experimentally for sand and glass beads with fast X-ray tomography, revealing water fragmentation into clusters of full pores interconnected by water corner films. The observed liquid morphology at the transition from piston to corner flow was reproduced by a quasi-static pore network model and predicted by percolation theory. The amount of capillary-retained water at flow transition controlling the subsequent drainage dynamics could be reproduced by an idealized star shaped pore whose geometry is deduced from macroscopic properties of the porous medium. Predictions of water content thresholds at flow transitions were in agreement with other critical saturation values associated with cessation of solute diffusion and of internal drainage (at field capacity) highlighting the criticality of water phase continuity disruption for formation of relatively stable unsaturated conditions controlled by slow corner flow that support life in soil.

Transport in Porous Media

Synchrotron-based fast micro-tomography is the method of choice to observe in situ multiphase flow and displacement dynamics on the pore scale. However, the image processing workflow is sensitive to a suite of manually selected parameters which can lead to ambiguous results. In this work, the relationship between porosity and permeability in response to systematically varied gray-scale threshold values was studied for different segmentation approaches on a dataset of Berea sandstone at a voxel length of 3 µm. For validation of the image processing workflow, porosity, permeability, and capillary pressure were compared to laboratory measurements on a larger-sized core plug of the same material. It was found that for global thresholding, minor variations in the visually permissive range lead to large variations in porosity and even larger variations in permeability. The latter is caused by changes in the pore-scale flow paths. Pore throats were found to be open for flow at large thresholds but closed for smaller thresholds. Watershed-based segmentation was found to be significantly more robust to manually chosen input parameters. Permeability and capillary pressure closely match experimental values; for capillary pressure measurements, the plateau of calculated capillary pressure curves was similar to experimental curves. Modeling on structures segmented with hysteresis thresholding was found to overpredict experimental capillary pressure values, while calculated permeability showed reasonable agreement to experimental data. This demonstrates that a good representation of permeability or capillary pressure alone is not a sufficient quality criterion for appropriate segmentation, but the data should be validated with both parameters. However, porosity is the least reliable quality criterion. In Electronic supplementary material The online version of this article (

Journal of Intelligent Systems

This work intends to apprehend and emphasize the contribution of image-processing techniques and computer vision in the treatment of clay-based material known in Meknes region. One of the various characteristics used to describe clay in a qualitative manner is porosity, as it is considered one of the properties that with “kill or cure” effectiveness. For this purpose, we use scanning electron microscopy images, as they are considered the most powerful tool for characterising the quality of the microscopic pore structure of porous materials. We present various existing methods of segmentation, as we are interested only in pore regions. The results show good matching between physical estimation and Voronoi diagram-based porosity estimation.

Transport in Porous Media

Wet tensile testing is a common method to assess the stability of bentonite-bonded moulding sands. For wet tensile testing, a specimen is first heated from above in order to simulate heat-driven moisture transport induced by the casting process. Then, tensile stress is applied until rupture. In this study, neutron radiography imaging was applied to moulding sands in situ during heating and wet tensile testing in order to investigate the effects of water kinematics on the tensile strength. Neutron radiography allowed the localization of the rupture plane and the quantitative determination of the local water content with sub-mm resolution. Quantification of the temperature at the rupture plane and of the heat kinematics within the specimen was accomplished by temperature measurements both in situ and ex situ. In this way, experimental data correlating the wet tensile strength with the specific conditions of moulding sands at the rupture plane were obtained for the first time. Series of experiments with different initial sand moisture contents were conducted. The results show B Guntram Jordan

Physical review. E, Statistical, nonlinear, and soft matter physics, 2015

We explore a mean-field theory of fluid imbibition and drainage through permeable porous solids. In the limit of vanishing inertial and viscous forces, the theory predicts the hysteretic "retention curves" relating the capillary pressure applied across a connected domain to its degree of saturation in wetting fluid in terms of known surface energies and void space geometry. To avoid complicated calculations, we adopt the simplest statistical mechanics, in which a pore interacts with its neighbors through narrow openings called "necks," while being either full or empty of wetting fluid. We show how the main retention curves can be calculated from the statistical distribution of two dimensionless parameters λ and α measuring the specific areas of, respectively, neck cross section and wettable pore surface relative to pore volume. The theory attributes hysteresis of these curves to collective first-order phase transitions. We illustrate predictions with a porous dom...

Journal of Nondestructive Evaluation, 2018

X-ray computed tomography (CT) images of three-phase silica sand and glass bead specimens are analyzed and used to evaluate the segmentation performances of Otsu-, and recursion-based multilevel algorithms. A global image segmentation technique that combines iterative and recursive algorithms, namely a refined statistics-based global segmentation is proposed for segmenting multi-phase granular geomaterials. The performance of the proposed algorithm is tested by segmenting partially saturated silica sand and glass bead specimens. For the tested silica sand specimen, the refined statistics method estimated void ratio and degree of saturation were 0.67 and 39.35%. The estimates for the glass bead specimen yielded 0.64 and 43.49%, respectively. The true void ratio (0.66) and degree of saturation (37.71%) were determined with a user-controlled Image processing software package-Image-Pro. It was found that the proposed method estimated the void ratio and the degree of saturation with 1.52 and 4.35 percent errors for the silica sand and with 15.63 and 0.34 percent errors for the glass bead, respectively. The computational time of the proposed method was found to be shorter than other methods considered. Overall, it is concluded that the proposed technique performed better in segmenting three-phase granular geomaterials.

Journal of environmental management, 2021

If biochar is applied to soil or stormwater treatment media, the saturated hydraulic conductivity (K) may be altered, which is a critical property affecting media performance. While a significant number of studies document biochar's effect on a porous medium's K, predictive models are lacking. Herein models are advanced for predicting K for repacked natural soil and engineered media when amended with biochar of various particle sizes and application rates. Experiments were conducted using three repacked natural soils, two uniform sands, and a bioretention medium amended with a wood biochar sieved to seven different biochar particle size distributions and applied at three rates. Experimental measurements showed a strong positive correlation between the interporosity of each medium and K. Across all media, the classic Kozeny-Carman (K-C) model predicted K and the relative change in K because of biochar amendment for each medium best. For soils alone, a recently developed model...

Journal of Synchrotron Radiation, 2013

A new modular X-ray-transparent experimental cell enables tomographic investigations of fluid rock interaction under natural reservoir conditions (confining pressure up to 20 MPa, pore fluid pressure up to 15 MPa, temperature ranging from 296 to 473 K). The portable cell can be used at synchrotron radiation sources that deliver a minimum X-ray flux density of 109 photons mm−2s−1in the energy range 30–100 keV to acquire tomographic datasets in less than 60 s. It has been successfully used in three experiments at the bending-magnet beamline 2BM at the Advanced Photon Source. The cell can be easily machined and assembled from off-the-shelf components at relatively low costs, and its modular design allows it to be adapted to a wide range of experiments and lower-energy X-ray sources.

Journal of Geophysical Research: Solid Earth, 2018

 3D propagation of hydrofracturing cracks in heterogeneous rock is probed using computed tomography and image reconstruction methods.  Material heterogeneity and horizontal geostress ratios significantly affect the 3D hydrofracturing crack initiation and propagation.  Fractal models are proposed to assess the complexity of hydraulic fractures and the stimulation efficiency of heterogeneous reservoirs.

Journal of Radioanalytical and Nuclear Chemistry

We performed tomographic investigation of the most heavily perturbed (thus highly radioactive) rod and its pristine/unirradiated replicate from the target No. 12 of the Swiss neutron spallation source (SINQ). The tomographic dataset reveals the 3D re-distribution of the lead filling inside the irradiated Zircaloy tube. The change in the linear attenuation coefficient of both the lead filling and the Zircaloy tube of the irradiated rod (due to the presence of the entrapped spallation products) in comparison with the pristine/unirradiated material is quantified. The dataset provides valuable input for the enhancement of safety and efficiency of future spallation targets at SINQ.

PLOS ONE, 2019

We describe in this paper the experimental procedure, the data treatment and the quantification of the black body correction: an experimental approach to compensate for scattering and systematic biases in quantitative neutron imaging based on experimental data. The correction algorithm is based on two steps; estimation of the scattering component and correction using an enhanced normalization formula. The method incorporates correction terms into the image normalization procedure, which usually only includes open beam and dark current images (open beam correction). Our aim is to show its efficiency and reproducibility: we detail the data treatment procedures and quantitatively investigate the effect of the correction. Its implementation is included within the open source CT reconstruction software MuhRec. The performance of the proposed algorithm is demonstrated using simulated and experimental CT datasets acquired at the ICON and NEUTRA beamlines at the Paul Scherrer Institut.

Applied Sciences, 2021

A crucial criterion for the quality of the additively manufactured parts is the porosity content for achieving an acceptable final relative density. In addition, for jewelry applications, visible pores are unacceptable at or in the vicinity of the surface. In this study, non-destructive 3D neutron microtomography is applied to map the spatial distribution of pores in additively manufactured red-gold samples. The 3D imaging assessment underlines the high relative density of the printed red-gold sample and indicates residual pore sizes are predominantly below the limit of concern for jewelry applications. The 3D maps of pores within printed samples highlight the effect of the scanning strategy on the final quality and location of pores in the printed samples. These results confirm that neutron microtomography is a novel and precise tool to characterize residual porosity in additively manufactured gold alloys and other higher-Z materials where such investigation using other non-destruc...

Energies

Estimating porous media properties is a vital component of geosciences and the physics of porous media. Until now, imaging techniques have focused on methodologies to match image-derived flows or geomechanical parameters with experimentally identified values. Less emphasis has been placed on the compromise between image processing techniques and the consequences on topological and morphological characteristics and on computed properties such as permeability. The effects of some of the most popular image processing techniques (filtering and segmentation) available in open source on 3D X-ray Microscopy (micro-XRM) images are qualitatively and quantitatively discussed. We observe the impacts of various filters such as erosion-dilation and compare the efficiency of Otsu’s method of thresholding and the machine-learning-based software Ilastik for segmentation.

Hydrogeology Journal

Paleokarst and paleokarst aquifers are not as well-documented as shallow karst, and studies of pores and pore networks are rare. Paleokarst lacks the open conduits typical of karst due to compaction, infilling, and diagenetic recrystallization by deep burial, so groundwater movement is through matrix porosity. In this case study, porosity networks in two saline aquifers in carbonate paleokarst of the Silurian Guelph Formation and Salina A-1 Carbonate Unit in southern Ontario (Canada) have been studied in drill cores at microscopic to macroscopic scales, utilizing medical computed tomography (CT) scans, optical petrography, macroscopic core examination, and scanning electron microscopy. The CT scans provided nondestructive three-dimensional visualization and quantification of pore distribution, size and volume, pore connections, and estimates of total porosity, similar to gas porosimetry. In the A-1 Carbonate, 50% of the pore volume comprises layers of macropores associated with alga...