Flow regimes during immiscible displacement

DOE Office of Scientific and Technical Information (OSTI.GOV)

Armstrong, Ryan T.; Mcclure, James; Berrill, Mark A.

Fractional ow of immiscible phases occurs at the pore scale where grain surfaces and phases interfaces obstruct phase mobility. However, the larger scale behavior is described by a saturation-dependent phenomenological relationship called relative permeability. As a consequence, pore-scale parameters, such as phase topology and/ or geometry, and details of the flow regime cannot be directly related to Darcy-scale flow parameters. It is well understood that relative permeability is not a unique relationship of wetting-phase saturation and rather depends on the experimental conditions at which it is measured. Herein we use fast X-ray microcomputed tomography to image pore-scale phase arrangements duringmore » fractional flow and then forward simulate the flow regimes using the lattice-Boltzmann method to better understand the underlying pore-scale flow regimes and their influence on Darcy-scale parameters. We find that relative permeability is highly dependent on capillary number and that the Corey model fits the observed trends. At the pore scale, while phase topologies are continuously changing on the scale of individual pores, the Euler characteristic of the nonwetting phase (NWP) averaged over a sufficiently large field of view can describe the bulk topological characteristics; the Euler characteristic decreases with increasing capillary number resulting in an increase in relative permeability. Lastly, we quantify the fraction of NWP that flows through disconnected ganglion dynamics and demonstrate that this can be a significant fraction of the NWP flux for intermediate wetting-phase saturation. Furthermore, rate dependencies occur in our homogenous sample (without capillary end effect) and the underlying cause is attributed to ganglion flow that can significantly influence phase topology during the fractional flow of immiscible phases.« less

Flow regimes during immiscible displacement

DOE PAGES

Armstrong, Ryan T.; Mcclure, James; Berrill, Mark A.; ...

2017-02-01

Fractional ow of immiscible phases occurs at the pore scale where grain surfaces and phases interfaces obstruct phase mobility. However, the larger scale behavior is described by a saturation-dependent phenomenological relationship called relative permeability. As a consequence, pore-scale parameters, such as phase topology and/ or geometry, and details of the flow regime cannot be directly related to Darcy-scale flow parameters. It is well understood that relative permeability is not a unique relationship of wetting-phase saturation and rather depends on the experimental conditions at which it is measured. Herein we use fast X-ray microcomputed tomography to image pore-scale phase arrangements duringmore » fractional flow and then forward simulate the flow regimes using the lattice-Boltzmann method to better understand the underlying pore-scale flow regimes and their influence on Darcy-scale parameters. We find that relative permeability is highly dependent on capillary number and that the Corey model fits the observed trends. At the pore scale, while phase topologies are continuously changing on the scale of individual pores, the Euler characteristic of the nonwetting phase (NWP) averaged over a sufficiently large field of view can describe the bulk topological characteristics; the Euler characteristic decreases with increasing capillary number resulting in an increase in relative permeability. Lastly, we quantify the fraction of NWP that flows through disconnected ganglion dynamics and demonstrate that this can be a significant fraction of the NWP flux for intermediate wetting-phase saturation. Furthermore, rate dependencies occur in our homogenous sample (without capillary end effect) and the underlying cause is attributed to ganglion flow that can significantly influence phase topology during the fractional flow of immiscible phases.« less

Disorder-induced stiffness degradation of highly disordered porous materials

NASA Astrophysics Data System (ADS)

Laubie, Hadrien; Monfared, Siavash; Radjaï, Farhang; Pellenq, Roland; Ulm, Franz-Josef

2017-09-01

The effective mechanical behavior of multiphase solid materials is generally modeled by means of homogenization techniques that account for phase volume fractions and elastic moduli without considering the spatial distribution of the different phases. By means of extensive numerical simulations of randomly generated porous materials using the lattice element method, the role of local textural properties on the effective elastic properties of disordered porous materials is investigated and compared with different continuum micromechanics-based models. It is found that the pronounced disorder-induced stiffness degradation originates from stress concentrations around pore clusters in highly disordered porous materials. We identify a single disorder parameter, φsa, which combines a measure of the spatial disorder of pores (the clustering index, sa) with the pore volume fraction (the porosity, φ) to scale the disorder-induced stiffness degradation. Thus, we conclude that the classical continuum micromechanics models with one spherical pore phase, due to their underlying homogeneity assumption fall short of addressing the clustering effect, unless additional texture information is introduced, e.g. in form of the shift of the percolation threshold with disorder, or other functional relations between volume fractions and spatial disorder; as illustrated herein for a differential scheme model representative of a two-phase (solid-pore) composite model material.

Two-Phase Solid/Fluid Simulation of Dense Granular Flows With Dilatancy Effects

NASA Astrophysics Data System (ADS)

Mangeney, Anne; Bouchut, Francois; Fernandez-Nieto, Enrique; Narbona-Reina, Gladys; Kone, El Hadj

2017-04-01

Describing grain/fluid interaction in debris flows models is still an open and challenging issue with key impact on hazard assessment [1]. We present here a two-phase two-thin-layer model for fluidized debris flows that takes into account dilatancy effects. It describes the velocity of both the solid and the fluid phases, the compression/ dilatation of the granular media and its interaction with the pore fluid pressure [2]. The model is derived from a 3D two-phase model proposed by Jackson [3] and the mixture equations are closed by a weak compressibility relation. This relation implies that the occurrence of dilation or contraction of the granular material in the model depends on whether the solid volume fraction is respectively higher or lower than a critical value. When dilation occurs, the fluid is sucked into the granular material, the pore pressure decreases and the friction force on the granular phase increases. On the contrary, in the case of contraction, the fluid is expelled from the mixture, the pore pressure increases and the friction force diminishes. To account for this transfer of fluid into and out of the mixture, a two-layer model is proposed with a fluid or a solid layer on top of the two-phase mixture layer. Mass and momentum conservation are satisfied for the two phases, and mass and momentum are transferred between the two layers. A thin-layer approximation is used to derive average equations. Special attention is paid to the drag friction terms that are responsible for the transfer of momentum between the two phases and for the appearance of an excess pore pressure with respect to the hydrostatic pressure. Interestingly, when removing the role of water, our model reduces to a dry granular flow model including dilatancy. We first compare experimental and numerical results of dilatant dry granular flows. Then, by quantitatively comparing the results of simulation and laboratory experiments on submerged granular flows, we show that our model contains the basic ingredients making it possible to reproduce the interaction between the granular and fluid phases through the change in pore fluid pressure. In particular, we analyse the different time scales in the model and their role in granular/fluid flow dynamics. References [1] R. Delannay, A. Valance, A. Mangeney, O. Roche, P. Richard, J. Phys. D: Appl. Phys., in press (2016). [2] F. Bouchut, E. D. Fernández-Nieto, A. Mangeney, G. Narbona-Reina, J. Fluid Mech., 801, 166-221 (2016). [3] R. Jackson, Cambridges Monographs on Mechanics (2000).

Macroscopic modeling of heat and water vapor transfer with phase change in dry snow based on an upscaling method: Influence of air convection

NASA Astrophysics Data System (ADS)

Calonne, N.; Geindreau, C.; Flin, F.

2015-12-01

At the microscopic scale, i.e., pore scale, dry snow metamorphism is mainly driven by the heat and water vapor transfer and the sublimation-deposition process at the ice-air interface. Up to now, the description of these phenomena at the macroscopic scale, i.e., snow layer scale, in the snowpack models has been proposed in a phenomenological way. Here we used an upscaling method, namely, the homogenization of multiple-scale expansions, to derive theoretically the macroscopic equivalent modeling of heat and vapor transfer through a snow layer from the physics at the pore scale. The physical phenomena under consideration are steady state air flow, heat transfer by conduction and convection, water vapor transfer by diffusion and convection, and phase change (sublimation and deposition). We derived three different macroscopic models depending on the intensity of the air flow considered at the pore scale, i.e., on the order of magnitude of the pore Reynolds number and the Péclet numbers: (A) pure diffusion, (B) diffusion and moderate convection (Darcy's law), and (C) strong convection (nonlinear flow). The formulation of the models includes the exact expression of the macroscopic properties (effective thermal conductivity, effective vapor diffusion coefficient, and intrinsic permeability) and of the macroscopic source terms of heat and vapor arising from the phase change at the pore scale. Such definitions can be used to compute macroscopic snow properties from 3-D descriptions of snow microstructures. Finally, we illustrated the precision and the robustness of the proposed macroscopic models through 2-D numerical simulations.

Beyond Darcy's law: The role of phase topology and ganglion dynamics for two-fluid flow

DOE PAGES

Armstrong, Ryan T.; McClure, James E.; Berrill, Mark A.; ...

2016-10-27

Relative permeability quantifies the ease at which immiscible phases flow through porous rock and is one of the most well known constitutive relationships for petroleum engineers. It however exhibits troubling dependencies on experimental conditions and is not a unique function of phase saturation as commonly accepted in industry practices. The problem lies in the multi-scale nature of the problem where underlying disequilibrium processes create anomalous macroscopic behavior. Here we show that relative permeability rate dependencies are explained by ganglion dynamic flow. We utilize fast X-ray micro-tomography and pore-scale simulations to identify unique flow regimes during the fractional flow of immisciblemore » phases and quantify the contribution of ganglion flux to the overall flux of non-wetting phase. We anticipate our approach to be the starting point for the development of sophisticated multi-scale flow models that directly link pore-scale parameters to macro-scale behavior. Such models will have a major impact on how we recover hydrocarbons from the subsurface, store sequestered CO 2 in geological formations, and remove non-aqueous environmental hazards from the vadose zone.« less

High Fidelity Computational Analysis of CO2 Trapping at Pore Scales

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kumar, Vinod

2013-07-13

With an alarming rise in carbon dioxide (CO2) emission from anthropogenic sources, CO2 sequestration has become an attractive choice to mitigate the emission. Some popular storage media for CO{sub 2} are oil reservoirs, deep coal-bed, and deep oceanic-beds. These have been used for the long term CO{sub 2} storage. Due to special lowering viscosity and surface tension property of CO{sub 2}, it has been widely used for enhanced oil recovery. The sites for CO{sub 2} sequestration or enhanced oil recovery mostly consist of porous rocks. Lack of knowledge of molecular mobility under confinement and molecule-surface interactions between CO2 and naturalmore » porous media results in generally governed by unpredictable absorption kinetics and total absorption capacity for injected fluids, and therefore, constitutes barriers to the deployment of this technology. Therefore, it is important to understand the flow dynamics of CO{sub 2} through the porous microstructures at the finest scale (pore-scale) to accurately predict the storage potential and long-term dynamics of the sequestered CO{sub 2}. This report discusses about pore-network flow modeling approach using variational method and analyzes simulated results this method simulations at pore-scales for idealized network and using Berea Sandstone CT scanned images. Variational method provides a promising way to study the kinetic behavior and storage potential at the pore scale in the presence of other phases. The current study validates variational solutions for single and two-phase Newtonian and single phase non-Newtonian flow through angular pores for special geometries whose analytical and/or empirical solutions are known. The hydraulic conductance for single phase flow through a triangular duct was also validated against empirical results derived from lubricant theory.« less

The Backscattering Phase Function for a Sphere with a Two-Scale Relief of Rough Surface

NASA Astrophysics Data System (ADS)

Klass, E. V.

2017-12-01

The backscattering of light from spherical surfaces characterized by one and two-scale roughness reliefs has been investigated. The analysis is performed using the three-dimensional Monte-Carlo program POKS-RG (geometrical-optics approximation), which makes it possible to take into account the roughness of objects under study by introducing local geometries of different levels. The geometric module of the program is aimed at describing objects by equations of second-order surfaces. One-scale roughness is set as an ensemble of geometric figures (convex or concave halves of ellipsoids or cones). The two-scale roughness is modeled by convex halves of ellipsoids, with surface containing ellipsoidal pores. It is shown that a spherical surface with one-scale convex inhomogeneities has a flatter backscattering phase function than a surface with concave inhomogeneities (pores). For a sphere with two-scale roughness, the dependence of the backscattering intensity is found to be determined mostly by the lower-level inhomogeneities. The influence of roughness on the dependence of the backscattering from different spatial regions of spherical surface is analyzed.

Unfitted Two-Phase Flow Simulations in Pore-Geometries with Accurate

NASA Astrophysics Data System (ADS)

Heimann, Felix; Engwer, Christian; Ippisch, Olaf; Bastian, Peter

2013-04-01

The development of better macro scale models for multi-phase flow in porous media is still impeded by the lack of suitable methods for the simulation of such flow regimes on the pore scale. The highly complicated geometry of natural porous media imposes requirements with regard to stability and computational efficiency which current numerical methods fail to meet. Therefore, current simulation environments are still unable to provide a thorough understanding of porous media in multi-phase regimes and still fail to reproduce well known effects like hysteresis or the more peculiar dynamics of the capillary fringe with satisfying accuracy. Although flow simulations in pore geometries were initially the domain of Lattice-Boltzmann and other particle methods, the development of Galerkin methods for such applications is important as they complement the range of feasible flow and parameter regimes. In the recent past, it has been shown that unfitted Galerkin methods can be applied efficiently to topologically demanding geometries. However, in the context of two-phase flows, the interface of the two immiscible fluids effectively separates the domain in two sub-domains. The exact representation of such setups with multiple independent and time depending geometries exceeds the functionality of common unfitted methods. We present a new approach to pore scale simulations with an unfitted discontinuous Galerkin (UDG) method. Utilizing a recursive sub-triangulation algorithm, we extent the UDG method to setups with multiple independent geometries. This approach allows an accurate representation of the moving contact line and the interface conditions, i.e. the pressure jump across the interface. Example simulations in two and three dimensions illustrate and verify the stability and accuracy of this approach.

Comparison of theory and experiment for NAPL dissolution in porous media

NASA Astrophysics Data System (ADS)

Bahar, T.; Golfier, F.; Oltéan, C.; Lefevre, E.; Lorgeoux, C.

2018-04-01

Contamination of groundwater resources by an immiscible organic phase commonly called NAPL (Non Aqueous Phase Liquid) represents a major scientific challenge considering the residence time of such a pollutant. This contamination leads to the formation of NAPL blobs trapped in the soil and impact of this residual saturation cannot be ignored for correct predictions of the contaminant fate. In this paper, we present results of micromodel experiments on the dissolution of pure hydrocarbon phase (toluene). They were conducted for two values of the Péclet number. These experiments provide data for comparison and validation of a two-phase non-equilibrium theoretical model developed by Quintard and Whitaker (1994) using the volume averaging method. The model was directly upscaled from the averaged pore-scale mass balance equations. The effective properties of the macroscopic model were calculated over periodic unit cells designed from images of the experimental flow cell. Comparison of experimental and numerical results shows that the transport model predicts correctly - with no fitting parameters - the main mechanisms of NAPL mass transfer. The study highlights the crucial need of having a fair recovery of pore-scale characteristic lengths to predict the mass transfer coefficient with accuracy.

Pore-scale supercritical CO 2 dissolution and mass transfer under drainage conditions

DOE PAGES

Chang, Chun; Zhou, Quanlin; Oostrom, Mart; ...

2016-12-05

Recently, both core- and pore-scale imbibition experiments have shown non-equilibrium dissolution of supercritical CO 2 (scCO 2) and a prolonged depletion of residual scCO 2. In this paper, pore-scale scCO 2 dissolution and mass transfer under drainage conditions were investigated using a two-dimensional heterogeneous micromodel and a novel fluorescent water dye with a sensitive pH range between 3.7 and 6.5. Drainage experiments were conducted at 9 MPa and 40 °C by injecting scCO 2 into the sandstone-analogue pore network initially saturated by water without dissolved CO 2 (dsCO 2). During the experiments, time-lapse images of dye intensity, reflecting water pH,more » were obtained. These images show non-uniform pH in individual pores and pore clusters, with average pH levels gradually decreasing with time. Further analysis on selected pores and pore clusters shows that (1) rate-limited mass transfer prevails with slowly decreasing pH over time when the scCO 2-water interface area is low with respect to the volume of water-filled pores and pore clusters, (2) fast scCO 2 dissolution and phase equilibrium occurs when scCO 2 bubbles invade into water-filled pores, significantly enhancing the area-to-volume ratio, and (3) a transition from rate-limited to diffusion-limited mass transfer occurs in a single pore when a medium area-to-volume ratio is prevalent. The analysis also shows that two fundamental processes – scCO 2 dissolution at phase interfaces and diffusion of dsCO 2 at the pore scale (10–100 µm) observed after scCO 2 bubble invasion into water-filled pores without pore throat constraints – are relatively fast. The overall slow dissolution of scCO 2 in the millimeter-scale micromodel can be attributed to the small area-to-volume ratios that represent pore-throat configurations and characteristics of phase interfaces. Finally, this finding is applicable for the behavior of dissolution at pore, core, and field scales when water-filled pores and pore clusters of varying size are surrounded by scCO 2 at narrow pore throats.« less

Pore-scale supercritical CO 2 dissolution and mass transfer under drainage conditions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chang, Chun; Zhou, Quanlin; Oostrom, Mart

Recently, both core- and pore-scale imbibition experiments have shown non-equilibrium dissolution of supercritical CO 2 (scCO 2) and a prolonged depletion of residual scCO 2. In this paper, pore-scale scCO 2 dissolution and mass transfer under drainage conditions were investigated using a two-dimensional heterogeneous micromodel and a novel fluorescent water dye with a sensitive pH range between 3.7 and 6.5. Drainage experiments were conducted at 9 MPa and 40 °C by injecting scCO 2 into the sandstone-analogue pore network initially saturated by water without dissolved CO 2 (dsCO 2). During the experiments, time-lapse images of dye intensity, reflecting water pH,more » were obtained. These images show non-uniform pH in individual pores and pore clusters, with average pH levels gradually decreasing with time. Further analysis on selected pores and pore clusters shows that (1) rate-limited mass transfer prevails with slowly decreasing pH over time when the scCO 2-water interface area is low with respect to the volume of water-filled pores and pore clusters, (2) fast scCO 2 dissolution and phase equilibrium occurs when scCO 2 bubbles invade into water-filled pores, significantly enhancing the area-to-volume ratio, and (3) a transition from rate-limited to diffusion-limited mass transfer occurs in a single pore when a medium area-to-volume ratio is prevalent. The analysis also shows that two fundamental processes – scCO 2 dissolution at phase interfaces and diffusion of dsCO 2 at the pore scale (10–100 µm) observed after scCO 2 bubble invasion into water-filled pores without pore throat constraints – are relatively fast. The overall slow dissolution of scCO 2 in the millimeter-scale micromodel can be attributed to the small area-to-volume ratios that represent pore-throat configurations and characteristics of phase interfaces. Finally, this finding is applicable for the behavior of dissolution at pore, core, and field scales when water-filled pores and pore clusters of varying size are surrounded by scCO 2 at narrow pore throats.« less

Pore-scale supercritical CO 2 dissolution and mass transfer under drainage conditions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chang, Chun; Zhou, Quanlin; Oostrom, Mart

Abstract: Recently, both core- and pore-scale imbibition experiments have shown non-equilibrium dissolution of supercritical CO 2 (scCO 2) and a prolonged depletion of residual scCO 2. In this study, pore-scale scCO 2 dissolution and mass transfer under drainage conditions were investigated using a two-dimensional heterogeneous micromodel and a novel fluorescent water dye with a sensitive pH range between 3.7 and 6.5. Drainage experiments were conducted at 9 MPa and 40 °C by injecting scCO 2 into the sandstone-analogue pore network initially saturated by water without dissolved CO 2 (dsCO 2). During the experiments, time-lapse images of dye intensity, reflecting watermore » pH, were obtained. These images show non-uniform pH in individual pores and pore clusters, with average pH levels gradually decreasing with time. Further analysis on selected pores and pore clusters shows that (1) rate-limited mass transfer prevails with slowly decreasing pH over time when the scCO 2-water interface area is low with respect to the volume of water-filled pores and pore clusters, (2) fast scCO 2 dissolution and phase equilibrium occurs when scCO 2 bubbles invade into water-filled pores, significantly enhancing the area-to-volume ratio, and (3) a transition from rate-limited to diffusion-limited mass transfer occurs in a single pore when a medium area-to-volume ratio is prevalent. The analysis also shows that two fundamental processes – scCO 2 dissolution at phase interfaces and diffusion of dsCO 2 at the pore scale (10-100 µm) observed after scCO 2 bubble invasion into water-filled pores without pore throat constraints – are relatively fast. The overall slow dissolution of scCO 2 in the millimeter-scale micromodel can be attributed to the small area-to-volume ratios that represent pore-throat configurations and characteristics of phase interfaces. This finding is applicable for the behavior of dissolution at pore, core, and field scales when water-filled pores and pore clusters of varying size are surrounded by scCO 2 at narrow pore throats.« less

Influence of capillary end effects on steady-state relative permeability estimates from direct pore-scale simulations

NASA Astrophysics Data System (ADS)

Guédon, Gaël Raymond; Hyman, Jeffrey De'Haven; Inzoli, Fabio; Riva, Monica; Guadagnini, Alberto

2017-12-01

We investigate and characterize the influence of capillary end effects on steady-state relative permeabilities obtained in pore-scale numerical simulations of two-phase flows. Our study is motivated by the observation that capillary end effects documented in two-phase laboratory-scale experiments can significantly influence permeability estimates. While numerical simulations of two-phase flows in reconstructed pore-spaces are increasingly employed to characterize relative permeabilities, a phenomenon which is akin to capillary end effects can also arise in such analyses due to the constraints applied at the boundaries of the computational domain. We profile the relative strength of these capillary end effects on the calculation of steady-state relative permeabilities obtained within randomly generated porous micro-structures using a finite volume-based two-phase flow solver. We suggest a procedure to estimate the extent of the regions influenced by these capillary end effects, which in turn allows for the alleviation of bias in the estimation of relative permeabilities.

Pore-scale modeling of saturated permeabilities in random sphere packings.

PubMed

Pan, C; Hilpert, M; Miller, C T

2001-12-01

We use two pore-scale approaches, lattice-Boltzmann (LB) and pore-network modeling, to simulate single-phase flow in simulated sphere packings that vary in porosity and sphere-size distribution. For both modeling approaches, we determine the size of the representative elementary volume with respect to the permeability. Permeabilities obtained by LB modeling agree well with Rumpf and Gupte's experiments in sphere packings for small Reynolds numbers. The LB simulations agree well with the empirical Ergun equation for intermediate but not for small Reynolds numbers. We suggest a modified form of Ergun's equation to describe both low and intermediate Reynolds number flows. The pore-network simulations agree well with predictions from the effective-medium approximation but underestimate the permeability due to the simplified representation of the porous media. Based on LB simulations in packings with log-normal sphere-size distributions, we suggest a permeability relation with respect to the porosity, as well as the mean and standard deviation of the sphere diameter.

New insights on the complex dynamics of two-phase flow in porous media under intermediate-wet conditions.

PubMed

Rabbani, Harris Sajjad; Joekar-Niasar, Vahid; Pak, Tannaz; Shokri, Nima

2017-07-04

Multiphase flow in porous media is important in a number of environmental and industrial applications such as soil remediation, CO 2 sequestration, and enhanced oil recovery. Wetting properties control flow of immiscible fluids in porous media and fluids distribution in the pore space. In contrast to the strong and weak wet conditions, pore-scale physics of immiscible displacement under intermediate-wet conditions is less understood. This study reports the results of a series of two-dimensional high-resolution direct numerical simulations with the aim of understanding the pore-scale dynamics of two-phase immiscible fluid flow under intermediate-wet conditions. Our results show that for intermediate-wet porous media, pore geometry has a strong influence on interface dynamics, leading to co-existence of concave and convex interfaces. Intermediate wettability leads to various interfacial movements which are not identified under imbibition or drainage conditions. These pore-scale events significantly influence macro-scale flow behaviour causing the counter-intuitive decline in recovery of the defending fluid from weak imbibition to intermediate-wet conditions.

Pore-scale modeling of capillary trapping in water-wet porous media: A new cooperative pore-body filling model

NASA Astrophysics Data System (ADS)

Ruspini, L. C.; Farokhpoor, R.; Øren, P. E.

2017-10-01

We present a pore-network model study of capillary trapping in water-wet porous media. The amount and distribution of trapped non-wetting phase is determined by the competition between two trapping mechanisms - snap-off and cooperative pore-body filling. We develop a new model to describe the pore-body filling mechanism in geologically realistic pore-networks. The model accounts for the geometrical characteristics of the pore, the spatial location of the connecting throats and the local fluid topology at the time of the displacement. We validate the model by comparing computed capillary trapping curves with published data for four different water-wet rocks. Computations are performed on pore-networks extracted from micro-CT images and process-based reconstructions of the actual rocks used in the experiments. Compared with commonly used stochastic models, the new model describes more accurately the experimental measurements, especially for well connected porous systems where trapping is controlled by subtleties of the pore structure. The new model successfully predicts relative permeabilities and residual saturation for Bentheimer sandstone using in-situ measured contact angles as input to the simulations. The simulated trapped cluster size distributions are compared with predictions from percolation theory.

Toward direct pore-scale modeling of three-phase displacements

NASA Astrophysics Data System (ADS)

Mohammadmoradi, Peyman; Kantzas, Apostolos

2017-12-01

A stable spreading film between water and gas can extract a significant amount of bypassed non-aqueous phase liquid (NAPL) through immiscible three-phase gas/water injection cycles. In this study, the pore-scale displacement mechanisms by which NAPL is mobilized are incorporated into a three-dimensional pore morphology-based model under water-wet and capillary equilibrium conditions. The approach is pixel-based and the sequence of invasions is determined by the fluids' connectivity and the threshold capillary pressure of the advancing interfaces. In addition to the determination of three-phase spatial saturation profiles, residuals, and capillary pressure curves, dynamic finite element simulations are utilized to predict the effective permeabilities of the rock microtomographic images as reasonable representations of the geological formations under study. All the influential features during immiscible fluid flow in pore-level domains including wetting and spreading films, saturation hysteresis, capillary trapping, connectivity, and interface development strategies are taken into account. The capabilities of the model are demonstrated by the successful prediction of saturation functions for Berea sandstone and the accurate reconstruction of three-phase fluid occupancies through a micromodel.

Pore-scale simulation of CO2-water-rock interactions

NASA Astrophysics Data System (ADS)

Deng, H.; Molins, S.; Steefel, C. I.; DePaolo, D. J.

2017-12-01

In Geologic Carbon Storage (GCS) systems, the migration of scCO2 versus CO2-acidifed brine ultimately determines the extent of mineral trapping and caprock integrity, i.e. the long-term storage efficiency and security. While continuum scale multiphase reactive transport models are valuable for large scale investigations, they typically (over-)simplify pore-scale dynamics and cannot capture local heterogeneities that may be important. Therefore, pore-scale models are needed in order to provide mechanistic understanding of how fine scale structural variations and heterogeneous processes influence the transport and geochemistry in the context of multiphase flow, and to inform parameterization of continuum scale modeling. In this study, we investigate the interplay of different processes at pore scale (e.g. diffusion, reactions, and multiphase flow) through the coupling of a well-developed multiphase flow simulator with a sophisticated reactive transport code. The objectives are to understand where brine displaced by scCO2 will reside in a rough pore/fracture, and how the CO2-water-rock interactions may affect the redistribution of different phases. In addition, the coupled code will provide a platform for model testing in pore-scale multiphase reactive transport problems.

Pore-scale dynamics of enzyme adsorption, swelling and reactive dissolution determine sugar yield in hemicellulose hydrolysis for biofuel production

NASA Astrophysics Data System (ADS)

Dutta, Sajal Kanti; Chakraborty, Saikat

2016-12-01

Hemicelluloses are the earth’s second most abundant structural polymers, found in lignocellulosic biomass. Efficient enzymatic depolymerization of xylans by cleaving their β-(1 → 4)-glycosidic bonds to produce soluble sugars is instrumental to the cost-effective production of liquid biofuels. Here we show that the multi-scale two-phase process of enzymatic hydrolysis of amorphous hemicelluloses is dominated by its smallest scale-the pores. In the crucial first five hours, two to fourfold swelling of the xylan particles allow the enzymes to enter the pores and undergo rapid non-equilibrium adsorption on the pore surface before they hydrolyze the solid polymers, albeit non-competitively inhibited by the products xylose and xylobiose. Rapid pore-scale reactive dissolution increases the solid carbohydrate’s porosity to 80-90%. This tightly coupled experimental and theoretical study quantifies the complex temporal dynamics of the transport and reaction processes coupled across scales and phases to show that this unique pore-scale phenomenon can be exploited to accelerate the depolymerization of hemicelluloses to monomeric sugars in the first 5-6 h. We find that an ‘optimal substrate loading’ of 5 mg/ml (above which substrate inhibition sets in) accelerates non-equilibrium enzyme adsorption and solid hemicellulose depolymerization at the pore-scale, which contributes three-quarters of the soluble sugars produced for bio-alcohol fermentation.

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Li; Wang, Mengyi; Kang, Qinjun

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2-water two-phase flow, multicomponent (CO 2(aq), H +, HCO 3 –, CO 3 2 – and OH –) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2(aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is requiredmore » by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.« less

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

DOE PAGES

Chen, Li; Wang, Mengyi; Kang, Qinjun; ...

2018-04-26

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2-water two-phase flow, multicomponent (CO 2(aq), H +, HCO 3 –, CO 3 2 – and OH –) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2(aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is requiredmore » by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.« less

Pore scale study of multiphase multicomponent reactive transport during CO2 dissolution trapping

NASA Astrophysics Data System (ADS)

Chen, Li; Wang, Mengyi; Kang, Qinjun; Tao, Wenquan

2018-06-01

Solubility trapping is crucial for permanent CO2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO2-water two-phase flow, multicomponent (CO2(aq), H+, HCO3-, CO32- and OH-) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO2(aq) concentration, scCO2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is required by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Finally, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.

The effects of wettability and trapping on relationships between interfacial area, capillary pressure and saturation in porous media: A pore-scale network modeling approach

NASA Astrophysics Data System (ADS)

Raeesi, Behrooz; Piri, Mohammad

2009-10-01

SummaryWe use a three-dimensional mixed-wet random pore-scale network model to investigate the impact of wettability and trapping on the relationship between interfacial area, capillary pressure and saturation in two-phase drainage and imbibition processes. The model is a three-dimensional network of interconnected pores and throats of various geometrical shapes. It allows multiple phases to be present in each capillary element in wetting and spreading layers, as well as occupying the center of the pore space. Two different random networks that represent the pore space in Berea and a Saudi Arabia reservoir sandstone are used in this study. We allow the wettability of the rock surfaces contacted by oil to alter after primary drainage. The model takes into account both contact angle and trapping hystereses. We model primary oil drainage and water flooding for mixed-wet conditions, and secondary oil injection for a water-wet system. The total interfacial area for pores and throats are calculated when the system is at capillary equilibrium. They include contributions from the arc menisci (AMs) between the bulk and corner fluids, and from the main terminal menisci (MTMs) between different bulk fluids. We investigate hysteresis in these relationships by performing water injection into systems of varying wettability and initial water saturation. We show that trapping and contact angle hystereses significantly affect the interfacial area. In a strongly water-wet system, a sharp increase is observed at the beginning of water flood, which shifts the area to a higher level than primary drainage. As we change the wettability of the system from strongly water-wet to strongly oil-wet, the trapped oil saturation decreases significantly. Starting water flood from intermediate water saturations, greater than the irreducible water saturation, can also affect the non-wetting phase entrapment, resulting in different interfacial area behaviors. This can increase the interfacial area significantly in oil-wet systems. A qualitative comparison of our results with the experimental data available in literature for glass beads shows, with some expected differences, an encouraging agreement. Also, our results agree well with those generated by the previously developed models.

Pore-scale modeling of moving contact line problems in immiscible two-phase flow

NASA Astrophysics Data System (ADS)

Kucala, Alec; Noble, David; Martinez, Mario

2016-11-01

Accurate modeling of moving contact line (MCL) problems is imperative in predicting capillary pressure vs. saturation curves, permeability, and preferential flow paths for a variety of applications, including geological carbon storage (GCS) and enhanced oil recovery (EOR). Here, we present a model for the moving contact line using pore-scale computational fluid dynamics (CFD) which solves the full, time-dependent Navier-Stokes equations using the Galerkin finite-element method. The MCL is modeled as a surface traction force proportional to the surface tension, dependent on the static properties of the immiscible fluid/solid system. We present a variety of verification test cases for simple two- and three-dimensional geometries to validate the current model, including threshold pressure predictions in flows through pore-throats for a variety of wetting angles. Simulations involving more complex geometries are also presented to be used in future simulations for GCS and EOR problems. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Pore-scale dynamics of enzyme adsorption, swelling and reactive dissolution determine sugar yield in hemicellulose hydrolysis for biofuel production

PubMed Central

Dutta, Sajal Kanti; Chakraborty, Saikat

2016-01-01

Hemicelluloses are the earth’s second most abundant structural polymers, found in lignocellulosic biomass. Efficient enzymatic depolymerization of xylans by cleaving their β-(1 → 4)-glycosidic bonds to produce soluble sugars is instrumental to the cost-effective production of liquid biofuels. Here we show that the multi-scale two-phase process of enzymatic hydrolysis of amorphous hemicelluloses is dominated by its smallest scale–the pores. In the crucial first five hours, two to fourfold swelling of the xylan particles allow the enzymes to enter the pores and undergo rapid non-equilibrium adsorption on the pore surface before they hydrolyze the solid polymers, albeit non-competitively inhibited by the products xylose and xylobiose. Rapid pore-scale reactive dissolution increases the solid carbohydrate’s porosity to 80–90%. This tightly coupled experimental and theoretical study quantifies the complex temporal dynamics of the transport and reaction processes coupled across scales and phases to show that this unique pore-scale phenomenon can be exploited to accelerate the depolymerization of hemicelluloses to monomeric sugars in the first 5–6 h. We find that an ‘optimal substrate loading’ of 5 mg/ml (above which substrate inhibition sets in) accelerates non-equilibrium enzyme adsorption and solid hemicellulose depolymerization at the pore-scale, which contributes three-quarters of the soluble sugars produced for bio-alcohol fermentation. PMID:27905534

Persistent Homology to describe Solid and Fluid Structures during Multiphase Flow

NASA Astrophysics Data System (ADS)

Herring, A. L.; Robins, V.; Liu, Z.; Armstrong, R. T.; Sheppard, A.

2017-12-01

The question of how to accurately and effectively characterize essential fluid and solid distributions and structures is a long-standing topic within the field of porous media and fluid transport. For multiphase flow applications, considerable research effort has been made to describe fluid distributions under a range of conditions; including quantification of saturation levels, fluid-fluid pressure differences and interfacial areas, and fluid connectivity. Recent research has effectively used topological metrics to describe pore space and fluid connectivity, with researchers demonstrating links between pore-scale nonwetting phase topology to fluid mobilization and displacement mechanisms, relative permeability, fluid flow regimes, and thermodynamic models of multiphase flow. While topology is clearly a powerful tool to describe fluid distribution, topological metrics by definition provide information only on the connectivity of a phase, not its geometry (shape or size). Physical flow characteristics, e.g. the permeability of a fluid phase within a porous medium, are dependent on the connectivity of the pore space or fluid phase as well as the size of connections. Persistent homology is a technique which provides a direct link between topology and geometry via measurement of topological features and their persistence from the signed Euclidean distance transform of a segmented digital image (Figure 1). We apply persistent homology analysis to measure the occurrence and size of pore-scale topological features in a variety of sandstones, for both the dry state and the nonwetting phase fluid during two-phase fluid flow (drainage and imbibition) experiments, visualized with 3D X-ray microtomography. The results provide key insights into the dominant topological features and length scales of a media which control relevant field-scale engineering properties such as fluid trapping, absolute permeability, and relative permeability.

Toward multiscale modelings of grain-fluid systems

NASA Astrophysics Data System (ADS)

Chareyre, Bruno; Yuan, Chao; Montella, Eduard P.; Salager, Simon

2017-06-01

Computationally efficient methods have been developed for simulating partially saturated granular materials in the pendular regime. In contrast, one hardly avoid expensive direct resolutions of 2-phase fluid dynamics problem for mixed pendular-funicular situations or even saturated regimes. Following previous developments for single-phase flow, a pore-network approach of the coupling problems is described. The geometry and movements of phases and interfaces are described on the basis of a tetrahedrization of the pore space, introducing elementary objects such as bridge, meniscus, pore body and pore throat, together with local rules of evolution. As firmly established local rules are still missing on some aspects (entry capillary pressure and pore-scale pressure-saturation relations, forces on the grains, or kinetics of transfers in mixed situations) a multi-scale numerical framework is introduced, enhancing the pore-network approach with the help of direct simulations. Small subsets of a granular system are extracted, in which multiphase scenario are solved using the Lattice-Boltzman method (LBM). In turns, a global problem is assembled and solved at the network scale, as illustrated by a simulated primary drainage.

Comparison of Pore-Network and Lattice Boltzmann Models for Pore-Scale Modeling of Geological Storage of CO2 in Natural Reservoir Rocks

NASA Astrophysics Data System (ADS)

Kohanpur, A. H.; Chen, Y.; Valocchi, A. J.; Tudek, J.; Crandall, D.

2016-12-01

CO2-brine flow in deep natural rocks is the focus of attention in geological storage of CO2. Understanding rock/flow properties at pore-scale is a vital component in field-scale modeling and prediction of fate of injected CO2. There are many challenges in working at the pore scale, such as size and selection of representative elementary volume (REV), particularly for material with complex geometry and heterogeneity, and the high computational costs. These issues factor into trade-offs that need to be made in choosing and applying pore-scale models. On one hand, pore-network modeling (PNM) simplifies the geometry and flow equations but can provide characteristic curves on fairly large samples. On the other hand, the lattice Boltzmann method (LBM) solves Navier-Stokes equations on the real geometry but is limited to small samples due to its high computational costs. Thus, both methods have some advantages but also face some challenges, which warrants a more detailed comparison and evaluation. In this study, we used industrial and micro-CT scans of actual reservoir rock samples to characterize pore structure at different resolutions. We ran LBM models directly on the characterized geometry and PNM on the equivalent 3D extracted network to determine single/two-phase flow properties during drainage and imbibition processes. Specifically, connectivity, absolute permeability, relative permeability curve, capillary pressure curve, and interface location are compared between models. We also did simulations on several subsamples from different locations including different domain sizes and orientations to encompass analysis of heterogeneity and isotropy. This work is primarily supported as part of the Center for Geologic Storage of CO2, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science and partially supported by the International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) based at Kyushu University, Japan.

Investigation of Biogrout processes by numerical analysis at pore scale

NASA Astrophysics Data System (ADS)

Bergwerff, Luke; van Paassen, Leon A.; Picioreanu, Cristian; van Loosdrecht, Mark C. M.

2013-04-01

Biogrout is a soil improving process that aims to improve the strength of sandy soils. The process is based on microbially induced calcite precipitation (MICP). In this study the main process is based on denitrification facilitated by bacteria indigenous to the soil using substrates, which can be derived from pretreated waste streams containing calcium salts of fatty acids and calcium nitrate, making it a cost effective and environmentally friendly process. The goal of this research is to improve the understanding of the process by numerical analysis so that it may be improved and applied properly for varying applications, such as borehole stabilization, liquefaction prevention, levee fortification and mitigation of beach erosion. During the denitrification process there are many phases present in the pore space including a liquid phase containing solutes, crystals, bacteria forming biofilms and gas bubbles. Due to the amount of phases and their dynamic changes (multiphase flow with (non-linear) reactive transport), there are many interactions making the process very complex. To understand this complexity in the system, the interactions between these phases are studied in a reductionist approach, increasing the complexity of the system by one phase at a time. The model will initially include flow, solute transport, crystal nucleation and growth in 2D at pore scale. The flow will be described by Navier-Stokes equations. Initial study and simulations has revealed that describing crystal growth for this application on a fixed grid can introduce significant fundamental errors. Therefore a level set method will be employed to better describe the interface of developing crystals in between sand grains. Afterwards the model will be expanded to 3D to provide more realistic flow, nucleation and clogging behaviour at pore scale. Next biofilms and lastly gas bubbles may be added to the model. From the results of these pore scale models the behaviour of the system may be studied and eventually observations may be extrapolated to a larger continuum scale.

The effect of pore-scale geometry and wettability on two-phase relative permeabilities within elementary cells

NASA Astrophysics Data System (ADS)

Bianchi Janetti, Emanuela; Riva, Monica; Guadagnini, Alberto

2017-04-01

We study the relative role of the complex pore space geometry and wettability of the solid matrix on the quantification of relative permeabilities characterizing steady state immiscible two-phase flow in porous media. We do so by considering elementary cells, which are typically employed in upscaling frameworks based on, e.g., homogenization or volume averaging. In this context one typically relies on the solution of pore-scale physics at a scale which is much smaller than that of an investigated porous system. Pressure-driven two-phase flow following simultaneous co-current injection of water and oil is numerically solved for a suite of regular and stochastically generated two-dimensional explicit elementary cells with fixed porosity and sharing main topological/morphological features. We show that relative permeabilities of the randomly generated elementary cells are significantly influenced by the formation of preferential percolation paths (principal pathways), giving rise to a strongly nonuniform distribution of fluid fluxes. These pathways are a result of the spatially variable resistance that the random pore structures exert on the fluid. The overall effect on relative permeabilities of the diverse organization of principal pathways, as driven by a given random realization at the scale of the unit cell, is significantly larger than that of the wettability of the host rock. In contrast to what can be observed for the random cells analyzed, relative permeabilities of regular cells display a clear trend with contact angle at the investigated scale. Our findings suggest the need to perform systematic upscaling studies in a stochastic context, to propagate the effects of uncertain pore space geometries to a probabilistic description of relative permeability curves at the continuum scale.

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Oostrom, Martinus; Mehmani, Yashar; Romero Gomez, Pedro DJ

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing.more » Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.« less

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.