Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spontaneous imbibition”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Spatial and Temporal Quantification of Spontaneous Imbibition

Spontaneous imbibition–the process of a wetting fluid displacing a non-wetting fluid purely by capillary forces–is a ubiquitous phenomenon in porous and fibrous materials. Here we present a new experimental method for quantification of spontaneous imbibition in geologic materials. This method makes it possible to perform spontaneous imbibition under elevated pressure conditions relevant to environmental and energy resource applications. In this work, computed tomography imaging reveals a new time-independent scaling relationship that describes local imbibition rates as a function of water saturation. Imbibition capillary pressure curves are calculated with this wetting-phase pressure and flowrate characterization, with no assumptions about the functional form, end point behavior, or scaling factors. Calculated end-point capillary pressure is nonzero, in agreement with recent pore-scale measurements of capillary pressure of trapped nonwetting phase. This work provides a new approach and insights into trapping and remobilization of non-wetting fluids in CO 2 storage reservoirs and contaminated groundwater aquifers.

58 GEOSCIENCES↗

Spontaneous imbibition in a square tube with corner films: theoretical model and numerical simulation

Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.

42 ENGINEERING↗

Preferential Solute Transport in Low Permeability Zones During Spontaneous Imbibition in Heterogeneous Porous Media

Multiphase flows in porous media|and the associated solute transport processes|are controlled by a combination of gravity, capillary, and viscous forces. Geologic heterogeneity influences flow and transport processes, yet gaps exist in our understanding of how heterogeneity impacts capillary-driven transport, such as during spontaneous imbibition. Here we use positron emission tomography, combined with a newly developed method for conducting spontaneous imbibition experiments, to observe solute transport into low permeability regions during imbibition. Using an experimentally parameterized 2D numerical model, we demonstrate how capillary-driven ow near the imbibition front can carry solutes into low permeability regions. This process displaces solutes from high permeability zones while the cumulative amount of solute in low permeability zones increases- opposite of what is observed under fully saturated solute transport conditions. These results highlight the complex flow and solute transport behavior that arise during multiphase displacements in heterogeneous geologic porous media.

58 GEOSCIENCES↗

Analytical solution and parameter estimation for heat of wetting and vapor adsorption during spontaneous imbibition in tuff

An analytical expression is derived for the thermal response observed during spontaneous imbibition of water into a dry core of zeolitic tuff. Sample tortuosity, thermal conductivity, and thermal source strength are estimated from fitting an analytical solution to temperature observations during a single laboratory test. Here, the closed-form analytical solution is derived using Green's functions for heat conduction in the limit of “slow” water movement; that is, when advection of thermal energy with the wetting front is negligible. The solution has four free fitting parameters and is efficient for parameter estimation. Laboratory imbibition data used to constrain the model include a time series of the mass of water imbibed, visual location of the wetting front through time, and temperature time series at six locations. The thermal front reached the end of the core hours before the visible wetting front. Thus, the predominant form of heating during imbibition in this zeolitic tuff is due to vapor adsorption in dry zeolitic rock ahead of the wetting front. The separation of the wetting front and thermal front in this zeolitic tuff is significant, compared to wetting front behavior of most materials reported in the literature. This work is the first interpretation of a thermal imbibition response to estimate transport (tortuosity) and thermal properties (including thermal conductivity) from a single laboratory test.

36 MATERIALS SCIENCE↗

Spontaneous Imbibition Tests and Parameter Estimation in Volcanic Tuff

We present a dynamic laboratory spontaneous imbibition test and interpretation method, demonstrated on volcanic tuff samples from the Nevada National Security Site. The method includes numerical inverse modeling to quantify uncertainty of estimated two-phase fluid flow properties. As opposed to other approaches requiring multiple different laboratory instruments, the dynamic imbibition method simultaneously estimates capillary pressure and relative permeability from one test apparatus.

36 MATERIALS SCIENCE↗

Parameter estimation from spontaneous imbibition into volcanic tuff

Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.

58 GEOSCIENCES↗

Spontaneous Imbibition Velocities in Porous Metal Layers on Ceramic Substrates Calculated Using Microstructural Analyses

In this paper, the imbibition velocity of molten silver into a wettable porous nickel interlayer on a non-wetting sapphire substrate is investigated. Several imbibition models for porous media that use bulk characteristics and material properties are employed to predict the imbibition front location as a function of time. The model predictions are compared to experimental observations. It is found that pore size distribution is a better predictor of imbibition velocity than correlations based on permeability and porosity. The tortuosity of the microstructure is also found to have a significant effect and should be considered. A new model accounting for the different liquid contact angles on the underlying substrate and the porous interlayer material is proposed and achieves better agreement with the experimental observations.

36 MATERIALS SCIENCE↗

Mixed imbibition controls the advance of wetting fluid in multiscale geological media

Imbibition plays a central role during multiphase flows in geological media; yet, our understanding of imbibition dynamics across the multiscale features of rocks, including fractures, microcracks, and heterogeneous matrices is limited. Here, we deployed a unique combination of micro- and clinical X-ray CT imaging to probe directly the spatiotemporal evolution of imbibition in a natural, multiscale shale sample. The method enabled unraveling of mixed imbibition behavior, with spontaneous imbibition rates in the matrix pores and microcracks driven by capillarity proceeding as $\propto t^{-1/2}$, whereas pressure-driven flows in the macroscopic fractures yield apparent imbibition rates with $\propto t^{-1}$ scaling. Regardless of the aperture, we showed that spontaneous imbibition proceeds uniformly away from the fractures. The effects of microcracks are further investigated by comparing the dynamics of the imbibition fronts with and without microcracks. Finally, these results highlight the paramount role of microcracks on anomalous front roughening and the need to account for microcracks appropriately in the development of predictive theoretical models.

58 GEOSCIENCES↗

Wettability Alteration with Weak Acid-Assisted Surfactant Flood

Oil-wetness and heterogeneity are two key reasons for low oil recovery by waterflooding in carbonate reservoirs. Surfactants have been effective in altering the oil-wet matrix to a more water-wet condition and initiating spontaneous imbibition. Because it takes time for the surfactant to alter wettability, oil recovery from the tight matrix is slow and sometimes not economically feasible. Acids have the potential of dissolving minerals, which may alter wettability. In this study, the enhanced oil recovery (EOR) performance of an acid-assisted surfactant solution, a novel technique, was evaluated for low-temperature applications. A set of acids and their acetates were tested. Bulk rock-acid reaction, wettability alteration (WA) tests, and spontaneous imbibition measurements were conducted at reservoir temperature (35°C) to identify effective candidates. Coreflood tests were then performed to evaluate the selected acid-surfactant formulations. Before and after a coreflood test, the core was scanned using micro-computed tomography (CT) to investigate pore structure alteration. Bulk reaction measurements showed that weak acids, especially acetic acid (AA), have the desired low reaction rates at 35°C. WA tests showed that AA can remove the crude oil off the rock surface and alter wettability through mineral dissolution. The surfactant can reduce contact angles from 160° to 58°; adding acid into the surfactant can further reduce it to 52°. Spontaneous imbibition experiments showed the synergy between the acid and the surfactant; the AA-surfactant solution had the highest oil recovery (62.6%) among acid-surfactant formulations. The acid improves the WA efficiency by the surfactant through surface mineral dissolution and lower ζ-potential. The imbibition transports the acid-surfactant solutions into the matrix, which minimizes face dissolution. Coreflood tests show that the AA-surfactant flood can increase the oil recovery rate and recover about 8% more oil compared to the surfactant flood. Micro-CT showed that a few mineral particles were transported along the core and partially plugged pore throats, which reduced permeability and diverted flow leading to improved oil recovery. Here, the transport of the acid in reservoirs and the potential plugging issues have to be carefully evaluated in future studies.

58 GEOSCIENCES↗

Surfactant-Aided Low-Salinity Waterflooding for Low-Temperature Carbonate Reservoirs

Carbonate reservoirs tend to be oil-wet/mixed-wet and heterogeneous because of mineralogy and diagenesis. The objective of this study is to improve oil recovery in low-temperature dolomite reservoirs using low-salinity and surfactant-aided spontaneous imbibition. The low-salinity brine composition was optimized using ζ-potential measurements, contact-angle (CA) experiments, and a novel wettability-alteration measure. Significant wettability alteration was observed on dolomite rocks at a salinity of 2,500 ppm. We evaluated 37 surfactants by performing CA, interfacial-tension (IFT), and spontaneous-imbibition experiments. Here, three (quaternary ammonium) cationic and one (sulfonate) anionic surfactants showed significant wettability alteration and produced 43–63% of original oil in place (OOIP) by spontaneous imbibition. At a low temperature (35°C), oil recovery by low-salinity effect is small compared with that by wettability-altering surfactants. Coreflood tests were performed with a selected low-salinity cationic surfactant solution. A novel coreflood was proposed that modeled heterogeneity and dynamic imbibition into low-permeability regions. The results of the “heterogeneous” coreflood were consistent with that of spontaneous-imbibition tests. These experiments demonstrated that a combination of low-salinity brine and surfactants can make originally oil-wet dolomite rocks more water-wet and improve oil recovery from regions bypassed by waterflood at a low temperature of 35°C.

02 PETROLEUM↗

Study on the Impacts of Capillary Number and Initial Water Saturation on the Residual Gas Distribution by NMR

The determination of microscopic residual gas distribution is beneficial for exploiting reservoirs to their maximum potential. In this work, both forced and spontaneous imbibition (waterflooding) experiments were performed on a high-pressure displacement experimental setup, which was integrated with nuclear magnetic resonance (NMR) to reveal the impacts of capillary number ( Ca ) and initial water saturation ( S wi ) on the residual gas distribution over four magnitudes of injection rates (Q = 0.001, 0.01, 0.1 and 1 mL/min), expressed as Ca (log Ca = –8.68, –7.68, –6.68 and –5.68), and three different S wi ( S wi = 0%, 39.34% and 62.98%). The NMR amplitude is dependent on pore volumes while the NMR transverse relaxation time ( T 2 ) spectrum reflects the characteristics of pore size distribution, which is determined based on a mercury injection (MI) experiment. Using this method, the residual gas distribution was quantified by comparing the T 2 spectrum of the sample measured after imbibition with the sample fully saturated by brine before imbibition. The results showed that capillary trapping efficiency increased with increasing S wi , and above 90% of residual gas existed in pores larger than 1 μm in the spontaneous imbibition experiments. The residual gas was trapped in pores by different capillary trapping mechanisms under different Ca , leading to the difference of residual gas distribution. The flow channels were mainly composed of micropores (pore radius, r < 1 μm) and mesopores ( r = 1–10 μm) at log Ca = –8.68 and –7.68, while of mesopores and macropores (r > 10 μm) at logCa = –5.68. At both S wi = 0% and 39.34%, residual gas distribution in macropores significantly decreased while that in micropores slightly increased with log Ca increasing to –6.68 and –5.68, respectively.

02 PETROLEUM↗

Influence of interfacial responses of Berea Sandstone in low-salinity waterflooding environments

Wetting conditions and fluid-fluid interfacial properties are essential to understanding oil recovery in the context of low-salinity and smart waterflooding. Many have deemed wettability alteration to be the dominant interaction, while some highlight the importance of fluid-fluid interactions. In some rock-fluid systems, benefits from the use of additives can arise to favor hydrocarbon production. Understanding the balance between the fluid-fluid and rock-fluid interactions when utilizing additives is crucial to evaluate benefits in improved efficiency. Ultimately, this understanding will shed light on recovery mechanisms. Regarding oil recovery, previous studies indicate that the injection of a naphthenic acid (NA) blend dissolved in an injection brine improves oil recovery in Berea Sandstone, particularly under low-salinity conditions. Interfacial rheology results for the brine-oil interface correlated well with the recovery response. This observation was speculated to serve as a potential improved recovery mechanism through oil-trapping suppression. However, the said studies did not evaluate potentially concurrent rock wettability alteration during recovery or the initial wetting condition. The purpose of this work is to determine the production behavior in systems with NAs, namely NA-10 (3-Cyclohexanepropinoic acid), NA-14 (Cyclopentane carboxylic acid), NA-18 (Cyclohexane carboxylic acid), and NA-mix (Naphthenic acid mixture), at different initial wetting conditions under spontaneous imbibition. In addition, we tested the NAs ability to alter wettability in Berea Sandstone and determine each “additive’s” capacity to improve recovery. Here, time-domain nuclear magnetic resonance (TD-NMR), contact angle, and Amott cell tests, i.e. spontaneous imbibition, were used to analyze the influence of the selected naphthenic acids on oil displacement under different initial wetting conditions. Our results indicate that individual naphthenic acids indeed have potential as additives in low-salinity waterflooding for oil-wet systems while negatively impacting recovery in water-wet systems. However, results suggest that the evolution of the interfacial properties and wettability conditions are equally important on the recovery behavior.

54 ENVIRONMENTAL SCIENCES↗

The Impact of Naphthenic Acids on Dynamic Fluid–Fluid Interactions: Implication for Enhanced Oil Recovery

Previous coreflooding results and wettability analyses in our group show that injection of naphthenic-acid-enriched water can improve oil recovery over traditional waterflooding. This observation is still a subject of research efforts without a definitive explanation. Naphthenic acids (NA) have been reported to drive wettability alteration and increase the water–oil interface elasticity. These alterations depend on the NA carbon number and aqueous-phase salinity, among other conditions, as reported in the literature. Smart-water flooding (SWF) research often links recovery to the initial wettability condition, being higher for initially oil-wet rock. SWF refers to a technique in which the aqueous-phase ion composition or/and salinity are changed to maximize oil recovery. Given NAs’ complex solution behavior, selecting acid combinations that prompt oil recovery is a difficult objective. The aim of this research is to determine the effects of select naphthenic acids on the oil–water interfacial rheology and wettability alteration and how these interfacial effects are associated with oil recovery under spontaneous imbibition. NAs were selected based on their carbon number, molecular structure, and solubility in the saline solution used in this research. We aimed at exploring which NAs should be used to regulate interfacial properties so as to either increase oil recovery or accelerate production. Time-domain nuclear magnetic resonance, interfacial dilatational rheology, and liquid-bridge experiments, i.e., proxy of snap-off, were conducted. A baseline was established using results obtained with a previously tested sulfate-rich aqueous phase, shown to be effective in recovering oil. Results show that NA14 and N18 increase the water–oil interfacial viscoelasticity and induce interfacial healing but led to different recovery factors. N10, while effective at inducing water wetness in oil-wet rock, is ineffective at increasing the recovery factor. We concluded that wettability and oil–water interfacial rheology are not exclusive, and instead they can synergistically favor EOR benefits. Moreover, oil recovery benefits under spontaneous imbibition are shown to depend strongly on the initial wettability conditions.

Medina-Rodriguez, Bryan X. (ORCID:0000000262922831↗

Pore Structure and Wettability of Bossier Shale, East Texas, United States: Insights from Integrated Porosimetry, Scattering, and Imbibition Approaches

In this study, ascertaining the pore geometry and wettability characteristics of tight shales is of great significance for revealing the mechanisms of occurrence, migration, and production of shale resources. Taking the Bossier Shale collected from Well A, East Texas, USA, as an example, the pore geometry and wettability behavior in a broad nm-& μm scale pore spectrum were quantified based on integrated techniques, including water immersion porosimetry (WIP), mercury intrusion porosimetry(MIP), (ultra)small-angle X-ray scattering [(U)SAXS], contact angle, and liquid spontaneous imbibition (SI). Mainly owing to differences in the sample sizes used, data interpretation, as well as the detectable pore type (e.g., connected or not) and pore diameter ranges among different methods, sample porosities derived from WIP, MIP, and (U)SAXS were quite different, which range from 7.70 to 13.34%, 5.12 to 9.52%, and 2.48 to 6.44%, respectively. Therefore, a comparison of the porosities derived from different methods must be handled with utmost care. Additionally, although (U)SAXS is an effective technique for detecting both connected and non-connected pores, the fraction of non-connected pores could be underestimated by directly comparing the pore size distribution derived from (U)SAXS and MIP tests. This is because the sample sizes used for and pore information reflected by these two methods are different. Furthermore, the studied Bossier Shales exhibited various contact angle values and imbibition behaviors when using differently polarized liquids, revealing their mixed-wet characteristics. Comprehensively considering the differences in contact angle, imbibition slope, and imbibed liquid volume, three sub-categories of wettability behavior, with respect to more oil-wet, more water-wet, and intermediate mixed-wet, were qualitatively identified.

(U)SAXS↗

A Dynamic Pore Network Model for Imbibition Simulation Considering Corner Film Flow

Wetting films can develop in the corners of angular pores under strong wetting conditions. Modeling the dynamics of corner film remains elusive using direct numerical simulations because of the significant scale difference between main meniscus and corner film flow. In this paper, the modified interacting capillary bundle model (ICB), developed in our previous work to describe accurately corner film dynamics in a single square tube, is incorporated into a single-pressure dynamic pore network model (DPNM) to simulate imbibition in strongly wetting porous media with corner film flow. The traditional pore network is decomposed into several layers of interacting subpore networks where the 0th layer of subpore network simulates the main meniscus flow and higher layers the corner film flow. The fluid flow between different layers is captured by interlayer throats. In addition, the snap-off mechanism caused by the thickening of wetting corner film is considered. The accuracy of the developed model is validated for four cases: spontaneous imbibition in a single square tube, wetting fluid redistribution through corner films under a capillary pressure difference, snap off in a narrow throat connecting two large pores, and imbibition dynamics in a real microfluidic porous geometry. The validated model is then used to simulate both spontaneous and controlled imbibition in a pore network with random pore size distribution. Finally, the interaction between corner film and main meniscus flow in porous media is analyzed from a pore-scale perspective.

58 GEOSCIENCES↗

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of a nanoparticle on wettability alteration and wettability retainment of carbonate reservoirs

The oil recovery in many carbonate reservoirs is low due to oil-wetness and heterogeneity. In this study, nanoparticles are evaluated for altering the wettability of oil-wet limestones and their properties are compared with those of an anionic surfactant. A surface-modified silica nanoparticle (SiNP) with a negative zeta potential was found to be aqueous stable in brines. Wettability studies showed that the SiNP cannot change the wettability of an initially oil-wet calcite plate, but SiNP treated calcite surface can remain water-wet after aging in oil. Spontaneous imbibition tests confirmed the observations of wettability tests and showed that a SiNP treated carbonate core retained its water-wettability during oil injection and aging. The imbibition into a SiNP treated and oil aged core was comparable to that in a water-wet core. However, SiNP was not able to remove the oil layers in initially oil-wet cores and imbibe water into these cores. In contrast, the anionic surfactant could alter the oil-wet carbonate rocks to a more water-wet condition, but failed to prevent the water-wet surface from getting oil-wet during oil-aging. The results suggest that the SiNP can be injected into wettability-altered (water-wet) reservoirs/regions to help retain the water-wettability during the long-term oil production. The 0.5 wt% SiNP dispersions transported through 18 mD limestone cores without any plugging. The SiNP retention was measured to be 2.5 mg/m 2 , much higher than that of the anionic surfactant.

02 PETROLEUM↗