Search

Production of methane from hydrate-bearing sediments requires hydrate dissociation for releasing mobile methane gas in sediments prior to gas production operations. Fines may migrate through or clog the pore space of sandy sediments depending on the geometry and topology of the pore space. Multiphase flow experiments were conducted on brine saturated uniform silica sand mixed with different percentages of kaolinite. Fines concentrations were correlated to computed tomography (CT) numbers in their host brine phase. Representative element volume (REV) procedure was used to study local porosity and local fines concentration. A cubical REV with a side length of 850 µm was selected. For low fines content (less than 6%), gas percolated through the specimens with no major particle translation or porosity change. For 6% kaolinite specimen, fractures were initiated, and propagated by densifying surrounding sediments. Local fines concentration study showed that fines migrated through the pores for specimens with fines content less than 6%. For the 6% kaolinite specimen, fines concentrated near fractures due to clogging. Clogging induced an increase in the capillary pressure at the pore throats and caused the capillary pressure to overcome the effective stress between sand particle, and fractures to be initiated. Sand particles experienced up to 200 µm translation and 15◦ rotation due to gas driven fracture in the 6% kaolinite specimen This research was made possible by the National Priorities Research Program (NPRP) grant # NPRP8-594-2-244 from Qatar National Research Fund (a member of Qatar Foundation). The findings achieved herein are solely the responsibility of the authors. The SMT images were collected using the X-ray Operations and Research Beamline Station 13-BMD at Argonne Photon Source (APS), Argonne National Laboratory. The authors thank Dr. Mark Rivers of APS for help in performing the SMT scans. They also acknowledge the support of GeoSoilEnviroCARS (Sector 13), which is supported by the National Science Foundation , USA, Earth Sciences ( EAR-1128799 ), and the US Department of Energy (DOE) , Geosciences ( DE-FG02-94ER14466 ). Use of the Advanced Photon Source, an Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory, was supported by DOE , USA under contract no. DE-AC02-06CH11357 . Scopus

The concept of representative elementary volume (REV) is critical to understand and predict the behavior of effective properties of complex heterogeneous materials in a multiscale manner. Porosity is usually used to define the REV of a given material, if porosity is stationary over consecutive volumetric increments, the sample is considered as REV regardless of the behavior of other microstructure properties. To investigate this common assumption, the change of key microstructure parameters such as local void ratio, sphericity, angularity, and coordination number of particles was computed over different volumes of REV for porosity. X-ray computed tomography was used to obtain three-dimensional volumes of natural sand packed at four different particle fragments. Efficient and robust algorithms were used to compute the microstructure properties from three-dimensional images. Findings revealed that the REVmin for porosity may not be adequate to be considered as a REV for parameters such as particle size dist...

The concept of representative elementary volume (REV) is critical to understand and predict the behaviour of effective parameters of complex heterogeneous media (e.g., soils) in a multiscale manner. Porosity is commonly used to define the REV of a given sample. In this paper we investigated whether the use of a REV for porosity can be used as a REV for other parameters such as particle size distribution, local void ratio and coordination number. X-ray computed tomography was used to obtain 3D images (i.e., volumes) of natural sand systems with different particle size distributions. 3D volumes of four different systems were obtained and a REV analysis was performed for these parameters utilizing robust 3D algorithms.Findings revealed that the REVmin for porosity may not be adequate to be considered as a REV for parameters such as particle size distribution, local void ratio and coordination number. The REVmin for these parameters was observed to be larger than the REVmin for porosity. Heterogeneity of the systems was found to be an important factor to determine the REV for the parameters analyzed in this paper. The REV analysis revealed that as the uniformity coefficient increased, a larger volume was required to obtain the REVmin for the distribution of particle sizes and coordination number whereas a smaller volume was required to obtain the REVmin for local void ratio. Therefore, determination of the REV for parameters described in this paper or any microscale parameter of concern should not be derived based on REV for porosity and should be determined based on their distributions over different volumes. © 2010 Elsevier B.V.

The objective of this paper was to investigate the impact of wettability of porous media on pore-scale characteristics of residual nonaqueous phase liquids (NAPLs). Synchrotron X-ray microtomography was used to obtain high-resolution three-dimensional images of fractionally wet sand systems with mean grain size of 250 microm. Pore-scale characteristics of NAPL blobs such as volume, lengths, interfacial areas, and sphericity index were computed using three-dimensional image processing algorithms. Four systems comprised of 100, 50, 25, and 0% NAPL-wet mass fractions containing the residual NAPL were imaged and analyzed. Findings indicate that spatial variation in wettability of porous media surfaces has a significant impact on pore-scale characteristics of residual NAPL blobs in saturated porous media systems. As the porous media comprises more water-wet surfaces, residual NAPL blobs increase in size and length due to the entrapment at large pore bodies. NAPL-water interfacial areas tend to increase as the NAPL-wet surface fractions increase in the systems. Overall residual NAPL saturations are less in fractionally wet systems and increase as the systems become more NAPL-wet or water-wet.

We present in this paper a methodology to calculate the distribution of local void ratio in porous media systems from high-resolution three-dimensional images. We introduce an algorithm to calculate the distribution of local void ratio from 3D images based on distance and watershed transforms. The watershed transform is used to segment touched or overlapped particles in an efficient way and the distance transform is used to calculate the boundaries of local void regions. The algorithm is validated using computer simulated 3D images of regular packing, irregular (non-spherical particles) packing, and random uniform spherical packing. Results show that the algorithm is robust, accurate and can be used to calculate local void ratio distribution of 3D systems regardless of irregularity in shapes, sizes, or arrangement of particles. X-ray microtomography images of different glass bead systems are used to calculate distributions of local void ratio. Parameters of distributions are function of porosity and particle-size distribution. The maximum local void ratio in each system is less than 3.0 and the minimum is greater than 0.2.

Pore-scale processes govern fundamental behavior in multiphase porous media systems. A high-resolution, three-dimensional image of the interior of a multiphase porous media system was obtained using synchrotron X-ray tomography. The system was imaged at a resolution of 12.46 mum following entrapment of the nonwetting phase at residual saturation. First, the physically representative network structure of the porous media system is extracted from the void space. This provides a direct mapping of the pore bodies and throats and enables pore-level calculations of coordination numbers, aspect ratios, and pore body and throat correlations. Next, algorithms developed to calculate properties of the entrapped nonwetting phase, such as volume, sphericity, interfacial area, and orientation, are applied to the residual nonwetting phase blobs. Finally, correlations between the pore network structure and nonwetting phase characteristics are examined. As expected, it was found that the nonwetting phase was trapped primarily in the largest pore spaces, the pore bodies with the highest aspect ratios, and the pore bodies with the highest coordination numbers. This work shows that, while there may be limitations related to the ability to capture REV-sized domains for some of the multiphase flow properties and phenomena, high-resolution X-ray tomography is able to provide the high quality datasets needed to observe and quantify the pore-scale phenomena and processes that govern multiphase flow in unconsolidated porous media systems.

This paper presents application of a series of algorithms used to extract pore network structure from high-resolution three-dimensional synchrotron microtomography images of unconsolidated porous media systems. These algorithms are based on the three-dimensional skeletonization that simplifies the pore space to networks in the form of nodes connected to paths. Dilation algorithms were developed to generate inscribed spheres on the nodes and paths of the medial axis to represent pore-bodies and pore-throats of the network, respectively. The end result is a physically representative pore network structure, i.e. three-dimensional spatial distribution (i.e. x-, y-, and z-coordinates) of pore-bodies and pore-throats, pore-body size distribution, pore-throat size distribution, and the connectivity. Systems analyzed in this study include different glass bead systems and natural marine sand. The media ranged in size from 0.123 to 1.0 mm, while the image volumes ranged between 7.7 and 108.9 mm3. In addition to extracting the pore network structure, the porosity, specific surface area, and representative elementary volume analysis on the porosity were calculated. Spatial correlation between pore-body sizes in the network was investigated using semivariograms and integral scale concepts. The impact of resolution on the calculated property was also investigated. In this work, we show that microtomography is an effective tool to non-destructively extract the structure of many systems. The quality of the datasets depends on photon energy, photon flux, size of the sample, type of the sample, and size of the sample ‘features’. Results show that the developed method of extracting pore network structure is applicable to ideal and natural porous media systems. The impact of resolution on the quantification of the network structure properties varies in its significance based on feature size of the system and the properties being calculated. Therefore, a thorough resolution sensitivity analysis should be carried out to determine the degree of error associated with a system imaged at a given resolution.

While network models of porous materials have traditionally been constructed using regular or disordered lattices, recent developments allow the direct modeling of more realistic structures such as sphere packings, microtomographic images, or computer-simulated materials. One of the obstacles in these newer approaches is the generation of network structures that are physically representative of the real systems. In this paper, we present and compare two different algorithms to extract pore network parameters from three-dimensional images of unconsolidated porous media systems. The first approach, which utilizes a pixelized image of the pore space, is an extension to unconsolidated systems of a medial-axis based approach (MA). The second approach uses a modified Delaunay tessellation (MDT) of the grain locations. The two algorithms are validated using theoretical packings with known properties and then the networks generated from random packing are compared. For the regular packings, both methods are able to provide the correct pore network structure, including the number, size, and location of inscribed pore bodies, the number, size, and location of inscribed pore throats, and the connectivity. Despite the good agreement for the regular packings, there were differences in both the spatial mapping and statistical distributions in network properties for the random packings. The discrepancies are attributed to the pixelization at low resolution, non-uniqueness of the inscribed pore-body locations, and differences in merging processes used in the algorithms, and serve to highlight the difficulty in creating a unique network from a complex, continuum pore space.

The objective of this work was to investigate the impact of grain geometry (size and shape) of porous media on the morphology of residual NAPL. Synchrotron microtomography was used to obtain maps of residual NAPL in multiphase systems. High-resolution, three-dimensional images of natural sand systems, comprising a range of grain sizes and shapes were imaged and analyzed. Findings indicate that residual NAPL saturation is influenced by the shapes of grains of the porous medium more than their sizes. In systems composed of grains with similar sphericity and angularity, residual saturations are independent of median grain sizes at the same operating regime (capillary-controlled regime in this work). Residual saturations tend to increase as the system comprised more angular or non-spherical grains where relatively large NAPL blobs are entrapped in such systems. While volumes of individual blobs tend to decrease as grain size decreases, grain geometry has more profound effects on the morphology of the residual NAPL blobs. Within a system composed of grains with similar shape characteristics, total NAPL-water interfacial area increases as grain sizes decrease where a large number of small blobs are trapped. Total meniscus NAPL-water interfacial area exhibits a linear relation with total interfacial area where it tends to increase as grain sizes decrease. However, while meniscus interfacial areas of individual blobs are highly influenced by pore geometry; residual blobs trapped in pores with complex geometry tend to have higher meniscus interfacial areas due to their branched nature which increases contacts with the wetting phase. U.S. Department of Energy, the National Science Foundation (Contract No. HRD-0503362) and Louisiana Board of Regents (Contract No. LINK-45). Scopus