Effects of shear-thinning fluids on residual oil formation in microfluidic pore networks.

Two-phase immiscible displacement in porous media is controlled by capillary and viscous forces when gravitational effects are negligible. The relative importance of these forces is quantified through the dimensionless capillary number Ca and the viscosity ratio M between fluid phases. When the displacing fluid is Newtonian, the effects of Ca and M on the displacement patterns can be evaluated independently. However, when the injecting fluids exhibit shear-thinning viscosity behaviour the values of M and Ca are interdependent. Under these conditions, the effects on phase entrapment and the general displacement dynamics cannot be dissociated. In the particular case of shear-thinning aqueous polymer solutions, the degree of interdependence between M and Ca is determined by the polymer concentration. In this work, two-phase immiscible displacement experiments were performed in micromodels, using shear-thinning aqueous polymer solutions as displacing fluids, to investigate the effect of polymer concentration on the relationship between Ca and M, the recovery efficiency, and the size distribution of the trapped non-wetting fluid. Our results show that the differences in terms of magnitude and distribution of the trapped phase are related to the polymer concentration which influences the values of Ca and M.

Experimental study of crossover from capillary to viscous fingering for supercritical CO2-water displacement in a homogeneous pore network.

Carbon sequestration in saline aquifers involves displacing brine from the pore space by supercritical CO(2) (scCO(2)). The displacement process is considered unstable due to the unfavorable viscosity ratio between the invading scCO(2) and the resident brine. The mechanisms that affect scCO(2)-water displacement under reservoir conditions (41 °C, 9 MPa) were investigated in a homogeneous micromodel. A large range of injection rates, expressed as the dimensionless capillary number (Ca), was studied in two sets of experiments: discontinuous-rate injection, where the micromodel was saturated with water before each injection rate was imposed, and continuous-rate injection, where the rate was increased after quasi-steady conditions were reached for a certain rate. For the discontinuous-rate experiments, capillary fingering and viscous fingering are the dominant mechanisms for low (logCa ≤ -6.61) and high injection rates (logCa ≥ -5.21), respectively. Crossover from capillary to viscous fingering was observed for logCa = -5.91 to -5.21, resulting in a large decrease in scCO(2) saturation. The discontinuous-rate experimental results confirmed the decrease in nonwetting fluid saturation during crossover from capillary to viscous fingering predicted by numerical simulations by Lenormand et al. (J. Fluid Mech.1988, 189, 165-187). Capillary fingering was the dominant mechanism for all injection rates in the continuous-rate experiment, resulting in monotonic increase in scCO(2) saturation.

Correlation of oil-water and air-water contact angles of diverse silanized surfaces and relationship to fluid interfacial tensions.

The use of air-water, θ(wa), or air-liquid contact angles is customary in surface science, while oil-water contact angles, θ(ow), are of paramount importance in subsurface multiphase flow phenomena including petroleum recovery, nonaqueous phase liquid fate and transport, and geological carbon sequestration. In this paper we determine both the air-water and oil-water contact angles of silica surfaces modified with a diverse selection of silanes, using hexadecane as the oil. The silanes included alkylsilanes, alkylarylsilanes, and silanes with alkyl or aryl groups that are functionalized with heteroatoms such as N, O, and S. These silanes yielded surfaces with wettabilities from water wet to oil wet, including specific silanized surfaces functionalized with heteroatoms that yield intermediate wet surfaces. The oil-water contact angles for clean and silanized surfaces, excluding one partially fluorinated surface, correlate linearly with air-water contact angles with a slope of 1.41 (R = 0.981, n = 13). These data were used to examine a previously untested theoretical treatment relating air-water and oil-water contact angles in terms of fluid interfacial energies. Plotting the cosines of these contact angles against one another, we obtain the relationship cos θ(wa) = 0.667 cos θ(ow) + 0.384 (R = 0.981, n = 13), intercepting cos θ(ow) = -1 at -0.284, which is in excellent agreement with the linear assumption of the theory. The theoretical slope, based on the fluid interfacial tensions σ(wa), σ(ow), and σ(oa), is 0.67. We also demonstrate how silanes can be used to alter the wettability of the interior of a pore network micromodel device constructed in silicon/silica with a glass cover plate. Such micromodels are used to study multiphase flow phenomena. The contact angle of the resulting interior was determined in situ. An intermediate wet micromodel gave a contact angle in excellent agreement with that obtained on an open planar silica surface using the same silane.

Assessing performance and closure for soil vapor extraction: integrating vapor discharge and impact to groundwater quality.

Soil vapor extraction (SVE) is typically effective for removal of volatile contaminants from higher-permeability portions of the vadose zone. However, contamination in lower-permeability zones can persist due to mass transfer processes that limit the removal effectiveness. After SVE has been operated for a period of time and the remaining contamination is primarily located in lower-permeability zones, the remedy performance needs to be evaluated to determine whether the SVE system should be optimized, terminated, or transitioned to another technology to replace or augment SVE. Numerical modeling of vapor-phase contaminant transport was used to investigate the correlation between measured vapor-phase mass discharge, MF(r), from a persistent, vadose-zone contaminant source and the resulting groundwater contaminant concentrations. This relationship was shown to be linear, and was used to directly assess SVE remediation progress over time and to determine the level of remediation in the vadose zone necessary to protect groundwater. Although site properties and source characteristics must be specified to establish a unique relation between MF(r) and the groundwater contaminant concentration, this correlation provides insight into SVE performance and support for decisions to optimize or terminate the SVE operation or to transition to another type of treatment.

Liquid CO2 displacement of water in a dual-permeability pore network micromodel.

Permeability contrasts exist in multilayer geological formations under consideration for carbon sequestration. To improve our understanding of heterogeneous pore-scale displacements, liquid CO(2) (LCO(2))-water displacement was evaluated in a pore network micromodel with two distinct permeability zones. Due to the low viscosity ratio (logM = -1.1), unstable displacement occurred at all injection rates over 2 orders of magnitude. LCO(2) displaced water only in the high permeability zone at low injection rates with the mechanism shifting from capillary fingering to viscous fingering with increasing flow rate. At high injection rates, LCO(2) displaced water in the low permeability zone with capillary fingering as the dominant mechanism. LCO(2) saturation (S(LCO2)) as a function of injection rate was quantified using fluorescent microscopy. In all experiments, more than 50% of LCO(2) resided in the active flowpaths, and this fraction increased as displacement transitioned from capillary to viscous fingering. A continuum-scale two-phase flow model with independently determined fluid and hydraulic parameters was used to predict S(LCO2) in the dual-permeability field. Agreement with the micromodel experiments was obtained for low injection rates. However, the numerical model does not account for the unstable viscous fingering processes observed experimentally at higher rates and hence overestimated S(LCO2).

Enhanced remedial amendment delivery to subsurface using shear thinning fluid and aqueous foam.

A major issue with in situ subsurface remediation is the ability to achieve an even spatial distribution of remedial amendments to the contamination zones in an aquifer or vadose zone. Amendment delivery to the aquifer using shear thinning fluid and to the vadose zone using aqueous foam has the potential to enhance the distribution. 2-D saturated flow cell experiments were conducted to evaluate the enhanced fluid sweeping over heterogeneous system, improved contaminant removal, and extended amendment presence in low-permeability zones achieved by shear thinning fluid delivery. Unsaturated column and flow cell experiments were conducted to investigate the improvement on contaminant mobilization mitigation, amendment distribution, and lateral delivery implemented by foam delivery. It was demonstrated that the shear thinning fluid injection enhanced the fluid sweeping and increased the delivery of remedial amendment into low-perm zones. The presence of amendment distributed by the shear thinning fluid in the low-permeability zones was increased. Foam delivery was shown to mitigate the mobilization of highly mobile contaminant from sediments. It also achieved more uniform amendment distribution in a heterogeneous unsaturated system, and demonstrated remarkable increasing in lateral distribution of the injected liquid compared to direct liquid injection.

Pore-scale study of transverse mixing induced CaCO3 precipitation and permeability reduction in a model subsurface sedimentary system.

A microfluidic pore structure etched into a silicon wafer was used as a two-dimensional model subsurface sedimentary system (i.e., micromodel) to study mineral precipitation and permeability reduction relevant to groundwater remediation and geological carbon sequestration. Solutions containing CaCl(2) and Na(2)CO(3) at four different saturation states (Ω = [Ca(2+)][CO(3)(2-)]/K(spCaCO(3))) were introduced through two separate inlets, and they mixed by diffusion transverse to the main flow direction along the center of the micromodel resulting in CaCO(3) precipitation. Precipitation rates increased and the total amount of precipitates decreased with increasing saturation state, and only vaterite and calcite crystals were formed (no aragonite). The relative amount of vaterite increased from 80% at the lowest saturation state (Ω(v) = 2.8 for vaterite) to 95% at the highest saturation state (Ω(v) = 4.5). Fluorescent tracer tests conducted before and after CaCO(3) precipitation indicate that pore spaces were occluded by CaCO(3) precipitates along the transverse mixing zone, thus substantially reducing porosity and permeability, and potentially limiting transformation from vaterite to the more stable calcite. The results suggest that mineral precipitation along plume margins can decrease both reactant mixing during groundwater remediation, and injection and storage efficiency during CO(2) sequestration.

A review of non-invasive imaging methods and applications in contaminant hydrogeology research.

Contaminant hydrogeological processes occurring in porous media are typically not amenable to direct observation. As a result, indirect measurements (e.g., contaminant breakthrough at a fixed location) are often used to infer processes occurring at different scales, locations, or times. To overcome this limitation, non-invasive imaging methods are increasingly being used in contaminant hydrogeology research. Four of the most common methods, and the subjects of this review, are optical imaging using UV or visible light, dual-energy gamma radiation, X-ray microtomography, and magnetic resonance imaging (MRI). Non-invasive imaging techniques have provided valuable insights into a variety of complex systems and processes, including porous media characterization, multiphase fluid distribution, fluid flow, solute transport and mixing, colloidal transport and deposition, and reactions. In this paper we review the theory underlying these methods, applications of these methods to contaminant hydrogeology research, and methods' advantages and disadvantages. As expected, there is no perfect method or tool for non-invasive imaging. However, optical methods generally present the least expensive and easiest options for imaging fluid distribution, solute and fluid flow, colloid transport, and reactions in artificial two-dimensional (2D) porous media. Gamma radiation methods present the best opportunity for characterization of fluid distributions in 2D at the Darcy scale. X-ray methods present the highest resolution and flexibility for three-dimensional (3D) natural porous media characterization, and 3D characterization of fluid distributions in natural porous media. And MRI presents the best option for 3D characterization of fluid distribution, fluid flow, colloid transport, and reaction in artificial porous media. Obvious deficiencies ripe for method development are the ability to image transient processes such as fluid flow and colloid transport in natural porous media in three dimensions, the ability to image many reactions of environmental interest in artificial and natural porous media, and the ability to image selected processes over a range of scales in artificial and natural porous media.

Estimation of interfacial tension between organic liquid mixtures and water.

Knowledge of IFT values for chemical mixtures helps guide the design and analysis of various processes, including NAPL remediation with surfactants or alcohol flushing, enhanced oil recovery, and chemical separation technologies, yet available literature values are sparse. A comprehensive comparison of thermodynamic and empirical models for estimating interfacial tension (IFT) of organic chemical mixtures with water is conducted, mainly focusing on chlorinated organic compounds for 14 ternary, three quaternary, and one quinary systems. Emphasis is placed on novel results for systems with three and four organic chemical compounds, and for systems with composite organic compounds like lard oil and mineral oil. Seven models are evaluated: the ideal and nonideal monolayer models (MLID and MLNID), the ideal and nonideal mutual solubility models (MSID and MSNID), an empirical model for ternary systems (EM), a linear mixing model based on mole fractions (LMMM), and a newly developed linear mixing model based on volume fractions of organic mixtures (LMMV) for higher order systems. The two ideal models (MLID and MSID) fit ternary systems of chlorinated organic compounds without surface active compounds relatively well. However, both ideal models did not perform well for the mixtures containing a surface active compound. However, for these systems, both the MLNID and MSNID models matched the IFT data well. It is shown that the MLNID model with a surface coverage value (0.00341 mmol/m2) obtained in this study can practically be used for chlorinated organic compounds. The LMMM results in poorer estimates of the IFT as the difference in IFT values of individual organic compounds in a mixture increases. The EM, with two fitting parameters, provided accurate results for all 14 ternarysystems including composite organic compounds. The new LMMV method for quaternary and higher component systems was successfully tested. This study shows that the LMMV may be able to be used for higher component systems and it can be easily incorporated into compositional multiphase flow models using only parameters from ternary systems.

Numerical and experimental investigation of DNAPL removal mechanisms in a layered porous medium by means of soil vapor extraction.

The purpose of this work is to identify the mechanisms that govern the removal of carbon tetrachloride (CT) during soil vapor extraction (SVE) by comparing numerical and analytical model simulations with a detailed data set from a well-defined intermediate-scale flow cell experiment. The flow cell was packed with a fine-grained sand layer embedded in a coarse-grained sand matrix. A total of 499 mL CT was injected at the top of the flow cell and allowed to redistribute in the variably saturated system. A dual-energy gamma radiation system was used to determine the initial NAPL saturation profile in the fine-grained sand layer. Gas concentrations at the outlet of the flow cell and 15 sampling ports inside the flow cell were measured during subsequent CT removal using SVE. Results show that CT mass was removed quickly in coarse-grained sand, followed by a slow removal from the fine-grained sand layer. Consequently, effluent gas concentrations decreased quickly at first, and then started to decrease gradually, resulting in long-term tailing. The long-term tailing was mainly due to diffusion from the fine-grained sand layer to the coarse-grained sand zone. An analytical solution for a one-dimensional advection and a first-order mass transfer model matched the tailing well with two fitting parameters. Given detailed knowledge of the permeability field and initial CT distribution, we were also able to predict the effluent concentration tailing and gas concentration profiles at sampling ports using a numerical simulator assuming equilibrium CT evaporation. The numerical model predictions were accurate within the uncertainty of independently measured or literature derived parameters. This study demonstrates that proper numerical modeling of CT removal through SVE can be achieved using equilibrium evaporation of NAPL if detailed fine-scale knowledge of the CT distribution and physical heterogeneity is incorporated into the model. However, CT removal could also be fit by a first-order mass transfer analytical model, potentially leading to an erroneous conclusion that the long-term tailing in the experiment was kinetically controlled due to rate-limited NAPL evaporation.

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction.

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.

Determination of NAPL-water interfacial areas in well-characterized porous media.

The nonaqueous-phase liquid (NAPL)-water interfacial area is an important parameter which influences the rate of NAPL dissolution in porous media. The aim of this study was to generate a set of baseline data for specific interfacial area for a two-phase-entrapped NAPL-water system in well-characterized porous media and subsequently use these data to evaluate two current theoretical models. The first model tested distributes entrapped NAPL over the pore classes based on Land's algorithm and assumes the resulting blobs to be spherical. The other model is thermodynamically based, assuming that reversible work done on the system results in an increase in interfacial area, such that the area between drainage and imbibition retention curves can be related to the interfacial area. Interfacial tracer tests (IFTT) were used to measure specific entrapped NAPL (hexadecane)-water interfacial areas in columns packed with four grades (12/20, 20/30, 30/40, 40/50) of silica sand. By use of the anionic surfactant dihexylsulfosuccinate (Aerosol MA80), IFTT gave specific interfacial areas between 58 cm(-1) for the finest sand and 16 cm(-1) for the coarsest, compared to values of between 33 and 7 cm(-1) for the first model and between 19 and 5 cm(-1) for the thermodynamic model. Results from the literature suggest that nonspherical blobs shapes occur relatively frequently; hence it is reasonable to suggest that the assumption of spherical NAPL blobs may explain the underprediction by the first model. The thermodynamic model underestimates the interfacial area because it assumes that entrapment occurs only within the largest pores. A modified version of the latter model, allowing entrapment across all pore classes, yielded values between 58 and 13 cm(-1). Of the models tested the modified thermodynamic model best predicts the interfacial area.

Dissolution of nonuniformly distributed immiscible liquid: intermediate-scale experiments and mathematical modeling.

The purpose of this work is to examine the effect of nonuniform distributions of immiscible organic liquid on dissolution behavior, with a specific focus on the condition dependency of dissolution (i.e., mass transfer) rate coefficients associated with applying mathematical models of differing complexities to measured data. Dissolution experiments were conducted using intermediate-scale flow cells packed with sand in which well-characterized zones of residual trichloroethene (TCE) and 1,2-dichloroethane (DCA) saturation were emplaced. A dual-energy gamma radiation system was used for in-situ measurement of NAPL saturation. Aqueous concentrations of TCE and DCA measured in the flow-cell effluent were significantly less than solubility, due primarily to dilution associated with the nonuniform immiscible-liquid distribution and bypass flow effects associated with physical heterogeneity. A quantitative analysis of flow and transport was conducted using a three-dimensional mathematical model wherein immiscible-liquid distribution, permeability variability, and sampling effects were explicitly considered. Independent values for the initial dissolution rate coefficients were obtained from dissolution experiments conducted using homogeneously packed columns. The independent predictions obtained from the model provided good representations of NAPL dissolution behavior and of total TCE/DCA mass removed, signifying model robustness. This indicates that for the complex three-dimensional model, explicit consideration of the larger scale factors that influenced immiscible-liquid dissolution in the flow cells allowed the use of a dissolution rate coefficient that represents only local-scale mass transfer processes. Conversely, the use of simpler models that did not explicitly consider the nonuniform immiscible-liquid distribution required the use of dissolution rate coefficients that are approximately 3 orders of magnitude smaller than the values obtained from the column experiments. The rate coefficients associated with the simpler models represent composite or lumped coefficients that incorporate the effects of the larger scale dissolution processes associated with the nonuniform immiscible-liquid distribution, which are not explicitly represented in the simpler models, as well as local-scale mass transfer. These results demonstrate that local-scale dissolution rate coefficients, such as those obtained from column experiments, can be used in models to successfully predict dissolution and transport of immiscible-liquid constituents at larger scales when the larger scale factors influencing dissolution behavior are explicitly accounted for in the model.