Enhancing residual trapping of supercritical CO2 via cyclic injections

We utilize synchrotron X-ray tomographic imaging to investigate the pore-scale characteristics and residual trapping of supercritical CO2 (scCO2) over the course of multiple drainage-imbibition (D-I) cycles in Bentheimer sandstone cores. Capillary pressure measurements are paired with X-ray image-derived saturation and connectivity metrics which describe the extent of drainage and subsequent residual (end of imbibition) scCO2 trapping. For the first D-I cycle, residual scCO2 trapping is suppressed due to high imbibition capillary number (Ca ≈ 10−6); however, residual scCO2 trapping dramatically increases for subsequent D-I cycles carried out at the same Ca value. This behavior is not predicted by conventional multiphase trapping theory. The magnitude of scCO2 trapping increase is hysteretic and depends on the relative extent of the sequential drainage processes. The hysteretic pore-scale behavior of the scCO2-brine-sandstone system observed in this study suggests that cyclic multiphase flow could potentially be used to increase scCO2 trapping for sequestration applications

Three-Dimensional Imaging and Quantification of Biomass and Biofilms in Porous Media

A new method to resolve biofilms in three dimensions in porous media using high-resolution synchrotron-based x-ray computed microtomography (CMT) has been developed. Imaging biofilms in porous media without disturbing the natural spatial arrangement of the porous media and associated biofilm has been a challenging task, primarily because porous media generally precludes conventional imaging via optical microscopy; x-ray tomography offers a potential alternative. One challenge for using this method is that most conventional x-ray contrast agents are water-soluble and easily diffuse into biofilms. To overcome this problem, silver-coated microspheres were added to the fluid phase to create an x-ray contrast that does not diffuse into the biofilm mass. Using this approach, biofilm imaging in porous media was accomplished with sufficient contrast to differentiate between the biomass- and fluid-filled pore spaces. The method was validated by using a two-dimensional micro-model flow cell where both light microscopy and CMT imaging were used to im age the biofilm. The results of this work has been published in Water Resources Research (Iltis et al., 2010). Additional work needs to be done to optimize this imaging approach, specifically, we find that the quality of the images are highly dependent on the coverage of the biofilm with Ag particles, - which means that we may have issues in dead-end pore space and for very low density (fluffy) biofilms. What we can image for certain with this technique is the biofilm surface that is well-connected to flow paths and thus well-supplied with nutrients etc

The effect of non-uniform microscale distribution of sorption sites on solute diffusion in soil

Conventional models of solute transport in soil consider only soil volumes large enough to average over microscale heterogeneities, and it is assumed that microscale variations are unimportant at the macroscale. In this research we test this assumption for cases in which the microscale distribution of solute-sorbing sites is patchy. We obtain a set of equations at the macroscale that allow for the effect of the microscale distribution with the mathematical technique of homogenization. We combine these equations with an image-based model that describes the true microscale pore geometry in a real, structured soil measured with X-ray computed tomography. The resulting models are used to test the microscale averaging assumptions inherent in conventional models. We show that, in general, macroscale diffusion is little affected by microscale variation in the distribution of sorption sites. Therefore, for most purposes the assumption of microscale averaging used in conventional models is justified. The effects of microscale heterogeneity are noticeable only when (i) the rate of sorption is slow compared with diffusion, but still fast enough to affect macroscale transport and (ii) the defined macroscale volume approaches the microscale. We discuss the effects when these conditions are me

Incorporating bubble evolution and transport in constitutive relationships for quasi- and non-equilibrium two-phase flows in porous media

There is a need to better understand the presence and transport of bubbles in multi-phase subsurface porous media so that these processes can be accurately described, and more efficient engineered solutions can be developed. To this end, constitutive relationships between geometric state variables (fluid-fluid curvature, Jnw; non-wetting phase volume, Vn; fluid-fluid interfacial area, anw; and Euler characteristic, χn) have become increasingly more common in efforts to uniquely predict the state of a two-fluid flow system. Both lattice Boltzmann simulations and fast X-ray microtomography (μCT) imaging experiments have shown that a geometric state function using the non-dimensionalized invariant properties of saturation, specific interfacial area, and Euler characteristic can uniquely predict the mean curvature of the system for both quasi- and non-equilibrium conditions, however, the presence of bubble evolution and the subsequent transport phenomena have not been explored. This study investigates whether the geometric state function remains unique with the inclusion of bubble generation and transport under quasi- and non-equilibrium two-fluid flow. The data presented here suggests that bubble formation and entrapment occur in a manner that cannot be predicted by the more traditional capillary pressure-saturation-interfacial area, Pc(Sw, anw), relationship, and further extensions to the constitutive relationship are needed to fully capture these mechanisms

Investigating the pore-scale mechanisms of microbial enhanced oil recovery

Microbial Enhanced Oil Recovery (MEOR) is a process where microorganisms are used for tertiary oil recovery. Numerous mechanisms have been proposed in the literature through which microorganisms facilitate the mobilization of residual oil. Herein, we focus on the MEOR mechanisms of interfacial tension reduction (via biosurfactant) and bioclogging in water-wet micromodels, using Shewanella oneidensis (MR-1) that causes bioclogging and Bacillus mojavensis (JF-2) that produces biosurfactant and causes bioclogging. Micromodels were flooded with an assortment of flooding solutions ranging from metabolically active bacteria to nutrient limited bacteria to dead inactive biomass to asses the effectiveness of the proposed MEOR mechanisms of bioclogging and biosurfactant production. Results indicate tertiary flooding of the micromodel system with biomass and biosurfactant was optimal for oil recovery due to the combined effects of bioclogging of the pore-space and interfacial tension reduction. However, biosurfactant was able to recover oil in some cases dependent on wettability. Biomass without biosurfactant that clogged the pore-space also successfully produced additional oil recovery. When analyzing residual oil blob morphology, MEOR resulted in oil blob size and radius of curvature distributions similar to those obtained by an abiotic capillary desaturation test, where flooding rate was increased post secondary recovery. Furthermore, for the capillary number calculated during MEOR flooding with bioclogging and biosurfactant, lower residual oil saturation was measured than for the corresponding capillary number under abiotic conditions. These results suggest that bioclogging and biosurfactant MEOR is a potentially effective approach for pore morphology modification and thus flow alteration in porous media that can have a significant effect on oil recovery beyond that predicted by capillary number.Keywords: bioclogging, micromodel, multiphase flow, microbial enhanced oil recovery, water flooding, biosurfactant, interfacial curvatur

Using synchrotron-based X-Ray microtomography and functional contrast agents in environmental applications

Despite very rapid development in commercial X-ray tomography technology, synchrotron-based tomography facilities still have a number of advantages over conventional systems. The high photon flux inherent of synchrotron radiation sources allows for (i) high resolution to micro- or nanometer scales depending on the individual beamline, (ii) rapid acquisition times that allow for collection of sufficient data for statistically significant results in a short amount of time as well as prevention of temporal changes that would take place during longer scan times, and (iii) optimal implementation of contrast agents that allow us to resolve features that would not be decipherable in scans obtained with a polychromatic radiation source. This chapter highlights recent advances in capabilities at synchrotron sources, as well as implementation of synchrotron-based computed microtomography (CMT) to two topics of interest to researchers in the soil science, hydrology, and environmental engineering fields, namely multiphase flow in porous media and characterization of biofilm architecture in porous media. In both examples, we make use of contrast agents and photoelectric edge-specic scanning (single- or dual-energy type), in combination with advanced image processing techniques

Influence of microstructural properties on geophysical measurements in sand-clay mixtures

We have performed a series of laboratory experiments on saturated sand-clay mixtures. Measurements include frequency-dependent electrical properties using the four-electrode technique (10 niHz to 1 MHz), permeability, porosity, and acoustic velocities. We mixed clean Ottawa (quartz) sand with Na-montmorillonite (Wyoming bentonite) in a number of different configurations containing 0 to 10% clay: as a dispersed mixture, as discrete clay clusters, and arranged in distinct layers. Solutions of CaCl{sub 2} ranging from 0.0005 N to 0.75 N (0.05 to 64 mS/cm) and deionized water were used as saturating fluids. We found the electrical properties to be dependent on clay content, fluid conductivity, and microstructure in a complex fashion. Increasing fluid conductivity and increasing clay content generally resulted in higher electrical conductivity. For an individual sample, two main regions of conduction exist: a region dominated by surface conduction and a region where the ionic strength of the saturating fluid controlled conduction. The sample geometry (dispersed, nondispersed, or layered clay configuration) was found to greatly affect the magnitude of the surface conductance in the range of low fluid conductivity

Compressional and Shear Wave Velocities for Artificial Granular Media Under Simulated Near Surface Conditions

Laboratory ultrasonic experiments were made on artificial soil samples in order to observe the effects of slight overburden, sand/clay ratio and pore fluid saturation on compressional and shear wave velocities. Up to several meters of overburden were simulated by applying low uniaxial stress of about 0.1 MPa to a restrained sample. Samples were fabricated from Ottawa sand mixed with a swelling clay (Wyoming bentonite). The amount of clay added was 1 to 40 percent by mass. Most measurements were made under room-dry conditions, but some measurements were made for fully-saturated sand-clay mixtures and for partially-saturated sand samples. For the dry sand-clay samples, compressional (P) velocities were low, ranging from about 200 to 500 m/s for the mixtures at low stress. Shear (S) velocities were about half of the compressional velocity, about 70 to 250 m/s. Dramatic increases in all velocities occurred with small uniaxial loads, indicating strong nonlinearity. Composition and grain packing control the mechanical response at grain contacts and the resulting nonlinear response at low stresses. P and S velocities are sensitive to the amount of clay added, even at low concentrations. At these low equivalent overburden conditions, adhesion and capillarity at grain contacts affect wave amplitudes, velocities, and frequency content in the partial saturation case

Report on laboratory scale thermally-coupled processes experiments

Yucca Mountain Site Characterization Project (YMP) is studying Yucca Mountain, Nevada as a potential repository for high-level nuclear wastes. The studies include predictions of the quantity and composition of water in the repository near-field environment that will affect the release rate of radioactive nuclides from waste packages, and the transport of the nuclides through the rock mass adjacent to these packages. The radioactive decay heat from the high- level nuclear waste may increase the temperature in the rock mass to the extent that coupled thermal-mechanical-hydrological-chemical (TMHC) processes may exist in the originally -partially-saturated Topopah Spring tuff-the host rock for the potential repository in Yucca Mountain. Modeling the coupled TMHC processes is necessary to predict the quantity and quality of water in the near-field environment for the entire life span of a repository (tens of thousands of years). In situ thermal tests are required to build up the confidence level of the coupled TMHC models. The purposes of conducting the laboratory studies of the coupled TMHC processes are to enhance our understanding of those processes, and to assist the interpretation of the field test results. Laboratory experiments deal with controlled experimental and boundary conditions, smaller sample sizes, and simpler geometrical configurations (e.g., regular shape and single fracture). These characteristics make the laboratory results suitable for understanding the processes. This in turn will make incorporation of these processes in model calculations more manageable. However, it should be noted that small sample size and simple geometrical configuration make the results of the laboratory tests unsuitable for direct use in predicting behaviors of in situ rock mass. The laboratory tests included in this reporting period are summarized below, along with projection of future work. This report fulfills the level 4 Milestone ID: SPL7A5M4