Additional Info

Links

Current Research and Scholarly Interests

ResearchI am interested in the recovery of unconventional hydrocarbon resources and mitigating carbon emissions from fossil fuels via geological sequestration of greenhouse gases. My research group and I examine the physics of flow through porous media at length scales that vary from the pore to the laboratory to the reservoir. The organizing themes are flow imaging to delineate the mechanisms of multiphase flow (oil, water, and gas) in porous media and the synthesis of models from experimental, theoretical, and field data. In all of our work, physical observations, obtained mainly from laboratory and field measurements, are interwoven with theory.

TeachingMy teaching interests center broadly around education of students to meet the energy challenges that we will face this century. I teach undergraduate courses that examine the interplay of energy use and environmental issues including renewable energy resources and sustainability. At the graduate level, I offer classes on enhanced oil recovery and the thermodynamics of hydrocarbon mixtures.

All Publications

Abstract

Silicon-based microfluidic devices, so-called micromodels in this application, are particularly useful laboratory tools for the direct visualization of fluid flow revealing pore-scale mechanisms controlling flow and transport phenomena in natural porous media. Current microfluidic devices with uniform etched depths, however, are limited when representing complex geometries such as the multiple-scale pore sizes common in carbonate rocks. In this study, we successfully developed optimized sequential photolithography to etch micropores (1.5 to 21 μm width) less deeply than the depth of wider macropores (>21 μm width) to improve the structural realism of an existing single-depth micromodel with a carbonate-derived pore structure. Surface profilimetry illustrates the configuration of the dual-depth dual-porosity micromodel and is used to estimate the corresponding pore volume change for the dual-depth micromodel compared to the equivalent uniform- or single-depth model. The flow characteristics of the dual-depth dual-porosity micromodel were characterized using micro-particle image velocimetry (μ-PIV), relative permeability measurements, and pore-scale observations during imbibition and drainage processes. The μ-PIV technique provides insights into the fluid dynamics within microfluidic channels and relevant fluid velocities controlled predominantly by changes in etching depth. In addition, the reduction of end-point relative permeability for both oil and water in the new dual-depth dual-porosity micromodel compared to the equivalent single-depth micromodel implies more realistic capillary forces occurring in the new dual-depth micromodel. Throughout the imbibition and drainage experiments, the flow behaviors of single- and dual-depth micromodels are further differentiated using direct visualization of the trapped non-wetting phase and the preferential mobilization of the wetting phase in the dual-depth micromodel. The visual observations agree with the relative permeability results. These findings indicate that dual-porosity and dual-depth micromodels have enhanced physical realism that is pertinent to oil recovery processes in complex porous media.

Understanding the role of brine ionic composition on oil recovery by assessment of wettability from colloidal forcesADVANCES IN COLLOID AND INTERFACE SCIENCEAlshakhs, M. J., Kovscek, A. R.2016; 233: 126-138

Abstract

The impact of injection brine salinity and ionic composition on oil recovery has been an active area of research for the past 25years. Evidence from laboratory studies and field tests suggests that implementing certain modifications to the ionic composition of the injection brine leads to greater oil recovery. The role of salinity modification is attributed to its ability to shift wettability of a rock surface toward water wetness. The amount of trapped oil released depends on the nature of rock, oil, and brine surface interactions. Reservoir rocks exhibit different affinities to fluids. Carbonates show stronger adsorption of oil films as opposed to the strongly water-wet and mixed-wet sandstones. The concentration of divalent ions and total salinity of the injection brine are other important factors to consider. Accordingly, this paper provides a review of laboratory and field studies of the role of brine composition on oil recovery from carbonaceous rock as well as rationalization of results using DLVO (Derjaguin, Landau, Verwey and Overbeek) theory of surface forces. DLVO evaluates the contribution of each component of the oil/brine/rock system to the wettability. Measuring zeta potential of each pair of surfaces by a charged particle suspension method is used to estimate double layer forces, disjoining pressure, and contact-angle. We demonstrate the applicability of the DLVO approach by showing a comprehensive experimental study that investigates the effect of divalent ions in carbonates, and uses disjoining pressure results to rationalize observations from core flooding and direct contact-angle measurements.

Abstract

The stability of foam in porous media is extremely important for realizing the advantages of foamed gas on gas mobility reduction. Foam texture (i.e., bubbles per volume of gas) achieved is dictated by foam generation and coalescence processes occurring at the pore-level. For foam injection to be widely applied during gas injection projects, we need to understand these pore-scale events that lead to foam stability/instability so that they are modeled accurately. Foam flow has been studied for decades, but most efforts focused on studying foam generation and coalescence in the absence of oil. Here, the extensive existing literature is reviewed and analyzed to identify open questions. Then, we use etched-silicon micromodels to observe foam generation and coalescence processes at the pore-level. Special emphasis is placed on foam coalescence in the presence of oil. For the first time, lamella pinch-off as described by Myers and Radke [40] is observed in porous media and documented. Additionally, a new mechanism coined "hindered generation" is found. Hindered generation refers to the role oil plays in preventing the successful formation of a lamella following snap-off near a pore throat.

Abstract

Sandstone formations are ubiquitous in both aquifers and petroleum reservoirs, of which clay is a major constituent. The release of clay particles from pore surfaces as a result of reduced injection fluid salinity can greatly modify the recovery of hydrocarbons from subsurface formations by shifting the wettability properties of the rock. In this paper we demonstrate a microfluidic approach whereby kaolinite is deposited into a two-dimensional microfluidic network (micromodel) to enable direct pore-scale, real-time visualization of fluid-solid interactions with representative pore-geometry and realistic surface interactions between the reservoir fluids and the formation rock. Structural characterization of deposited kaolinite particles agrees well with natural modes of occurrence in Berea sandstones; hence, the clay deposition method developed in this work is validated. Specifically, more than 90% of the deposited clay particles formed pore-lining structures and the remainder formed pore bridging structures. Further, regions of highly concentrated clay deposition likely leading to so-called Dalmatian wetting properties were found throughout the micromodel. Two post-deposition treatments are described whereby clay is adhered to the silicon surface reversibly and irreversibly resulting in microfluidic systems that are amenable to studies on (i) the fundamental mechanisms governing the increased oil recovery during low salinity waterfloods and (ii) the effect of a mixed-wet surface on oil recovery, respectively. The reversibly functionalized platform is used to determine the conditions at which stably adhered clay particles detach. Specifically, injection brine salinity below 6000 ppm of NaCl induced kaolinite particle release from the silicon surface. Furthermore, when applied to an aged system with crude oil, the low salinity waterflood recovered an additional 14% of the original oil in place compared to waterflooding with the formation brine.

Abstract

In a conventional ramped temperature oxidation kinetics cell experiment, an electrical furnace is used to ramp temperature at a prescribed rate. Thus, the heating rate of a kinetics cell experiment is limited by furnace performance to heating rates of about 0.5-3 °C/min. A new reactor has been designed to overcome this limit. It uses an induction heating method to ramp temperature. Induction heating is fast and easily controlled. The new reactor covers heating rates from 1 to 30 °C/min. This is the first time that the oxidation profiles of a crude oil are available over such a wide range of heating rate. The results from an induction reactor and a conventional kinetics cell at roughly 2 °C/min are compared to illustrate consistency between the two reactors. The results at low heating rate are the same as the conventional kinetics cell. As presented in the paper, the new reactor couples well with the isoconversional method for interpretation of reaction kinetics.

Abstract

The motivation for this work is a dramatically improved understanding of the fluid mechanics of drainage processes with applications such as CO_{2} storage in saline aquifers and water-alternating-gas injection as an enhanced oil recovery method. In this paper we present in situ distributions of wetting and nonwetting fluids obtained during core-scale two-phase immiscible drainage experiments. The ratio of the viscosity of the resident fluid to that of the invading fluid varies across a range of 0.43 to 150. Saturation distributions observed during dynamic displacement experiments are surprisingly smooth and do not display only one or a few dominant fingers, contrary to the indications of the current literature. The analysis of the saturation distribution using the fractal dimensions of the dynamic three-dimensional saturation distributions suggests that the constitutive relationships for porous media, namely, the relative permeability functions, are history dependent. Accordingly, it is suggested that the nonlinear, unstable flow regime is the regime where efforts to improve physical understanding must be focused.

Abstract

A well-characterized set of large-scale laboratory experiments is presented, illustrating forced imbibition displacements in the presence of irreducible wetting phase saturation in a cylindrical, homogeneous Berea sandstone rock. Experiments are designed to operate in the regime of compact microscopic flows and large-scale viscous instability. The distribution of fluid phases during the flow process is visualized by high-resolution computed tomography imaging. Linear stability analysis and high-accuracy numerical simulations are employed to analyze the ability of macroscopic continuum equations to provide a consistent approximation of the displacement process. The validity of the equilibrium relative permeability functions, which form the basis for the continuum model, is fundamentally related to the stability of the displacement process. It is shown that not only is the stable flow regime modeled accurately by existing continuum models, but the onset of instability as well as the initial unstable modes are also determined with reasonable accuracy for unstable flows. However, the continuum model is found to be deficient in the case of fully developed unstable flows.