Time filter

Source Type

Palo Alto, CA, United States

Cheng P.,Yale University | Bestehorn M.,TU Brandenburg | Firoozabadi A.,Yale University | Firoozabadi A.,Reservoir Engineering Research Institute
Water Resources Research | Year: 2012

Solubility trapping of carbon dioxide (CO2) in deep saline aquifers is considered one of the most effective methods for carbon sequestration. Dissolution of CO2 into the brine may create gravitational instabilities that lead to the onset of convection, which greatly enhances the storage efficiency and reduces the possibilities of leakage. Convection appears in the form of downward traveling fingers of relatively dense, CO2-dissolved fluid. Many natural aquifer formations display considerable permeability anisotropy, where the horizontal permeability k h may be several times greater than the vertical permeability k z. It has been previously found that increasing kh for a fixed kz reduces the critical time tc at which onset occurs and the critical wavelength λc with which the fingers initially form. We extend earlier work by showing how and why this occurs. Our results reveal new insights about λc. We have studied the behavior for times greater than tc using high-resolution numerical simulations. We show that the enhanced dissolution from convection can become significant much earlier in anisotropic media. Furthermore, the effects of anisotropy may be sustained for a long period of time. Our results suggest that permeability anisotropy can allow a wider range of aquifer formations to be considered for effective sequestration. © 2012. American Geophysical Union.

Shahraeeni E.,ETH Zurich | Shahraeeni E.,Reservoir Engineering Research Institute | Lehmann P.,ETH Zurich | Or D.,ETH Zurich
Water Resources Research | Year: 2012

Prediction of drying rates from porous media remains a challenge due to complex interactions between ambient conditions and porous medium properties. Evaporation from a gradually drying porous surface across air boundary layer exhibits nonlinear behavior due to enhanced diffusive fluxes from increasingly isolated active pores. These nonlinear interactions were quantified by modeling evaporation from surfaces composed of individual pores considering surface water content dynamics and internal transport within the medium. Wind tunnel experiments show that in contrast with nearly constant evaporation rates obtained at low atmospheric demand (typically <5 mm/d), evaporation fluxes under high atmospheric demand (high air velocities) exhibit a continuous decrease with surface drying even in the absence of internal capillary flow limitations. The isolated pore evaporation model captures surface drying dynamics for a range of atmospheric demands associated with air velocity and boundary layer thickness. As a surface dries under low atmospheric demand (low air speed, thick boundary layer), the remaining active pores become gradually isolated with a conforming vapor concentration field becoming increasingly three-dimensional thereby enhancing evaporative flux per pore. Such enhancement may fully compensate for reduced evaporative surface area leading to observed constant evaporation rate under low demand. For high evaporative demand, limitations to vapor field configuration within thin boundary layer limit flux compensation efficiency and leads to decreasing evaporative flux with surface drying irrespective of internal supply capacity. The model provides new insights into the intrinsic links between surface properties and atmospheric conditions in determining a range of evaporative dynamics for similar surface wetness conditions. © 2012 American Geophysical Union.

Lukanov B.,Reservoir Engineering Research Institute | Firoozabadi A.,Yale University
Langmuir | Year: 2014

The self-assembly of amphiphilic molecules is a key process in numerous biological and chemical systems. When salts are present, the formation and properties of molecular aggregates can be altered dramatically by the specific types of ions in the electrolyte solution. We present a molecular thermodynamic model for the micellization of ionic surfactants that incorporates quantum dispersion forces to account for specific ion effects explicitly through ionic polarizabilities and sizes. We assume that counterions are distributed in the diffuse region according to a modified Poisson-Boltzmann equation and can reach all the way to the micelle surface of charge. Stern layers of steric exclusion or distances of closest approach are not imposed externally; these are accounted for through the counterion radial distribution profiles due to the incorporation of dispersion potentials, resulting in a simple and straightforward treatment. There are no adjustable or fitted parameters in the model, which allows for a priori quantitative prediction of surfactant aggregation behavior based only on the initial composition of the system and the surfactant molecular structure. The theory is validated by accurately predicting the critical micelle concentration (CMC) for the well-studied sodium dodecyl sulfate (SDS) surfactant and its alkaline-counterion derivatives in mono- and divalent salts, as well as the molecular structure parameters of SDS micelles such as aggregation numbers and micelle surface potential. © 2014 American Chemical Society.

Myint P.C.,Yale University | Myint P.C.,Lawrence Livermore National Laboratory | Firoozabadi A.,Yale University | Firoozabadi A.,Reservoir Engineering Research Institute
Physics of Fluids | Year: 2013

The density increase from carbon dioxide (CO2) dissolution in water or hydrocarbons creates buoyancy-driven instabilities that may lead to the onset of convection. The convection is important for both CO2 sequestration in deep saline aquifers and CO2 improved oil recovery from hydrocarbon reservoirs. We perform linear stability analyses to study the effect of fluid compressibility and interface movement on the onset of buoyancy-driven convection in porous media. Compressibility relates to a non-zero divergence of the velocity field. The interface between the CO2 phase and the aqueous or hydrocarbon phase moves with time as a result of the volume change that occurs upon CO2 dissolution. Previous stability analyses have neglected these two aspects by assuming that the aqueous or hydrocarbon phase is incompressible and that the interface remains fixed in position. The stability analyses are used to compute two key quantities: (1) the critical time and (2) the critical wavenumber. Our results indicate that compressibility has a negligible effect on the critical time and the critical wavenumber in CO2-water mixtures. We use thermodynamics to derive an expression which shows that the two opposing physical processes which contribute to the divergence are comparable in magnitude and largely cancel each other. This result explains why compressibility does not significantly affect the onset, and it also demonstrates the link between compressibility and the volume change that causes movement of the interface. Compared to when the interface is fixed in position, a moving interface in CO2-water mixtures may reduce the critical time by up to around 10%, which can be significant in low permeability formations. The decrease in the critical time due to interface movement may be much more pronounced in hydrocarbons than in water. This could have important implications for CO2 improved oil recovery. © 2013 AIP Publishing LLC.

Moortgat J.,Reservoir Engineering Research Institute | Firoozabadi A.,Reservoir Engineering Research Institute | Firoozabadi A.,Yale University
Advances in Water Resources | Year: 2010

We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE-dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO2 injection into a vertical domain saturated with a multi-component petroleum fluid. © 2010 Elsevier Ltd.

Discover hidden collaborations