A new compressive scheme to simulate species transfer across fluid interfaces using the Volume-Of-Fluid method

Pore-scale simulation of low-salinity waterflooding in mixed-wet systems: effect of corner flow, surface heterogeneity and kinetics of wettability alteration

The initial wettability state of the candidate oil reservoirs for low-salinity waterflooding (LSWF) is commonly characterized as mixed-wet. In mixed-wet systems, both the two-phase flow dynamics and the salt transport are significantly influenced by the corner flow of the wetting phase. Thus this study aims at comprehensive evaluation of LSWF efficiency by capturing the effect of corner flow and non-uniform wettability distribution. In this regard, direct numerical simulations under capillary-dominated flow regime were performed using the OpenFOAM Computational Fluid Dynamics toolbox. The results indicate that corner flow results in the transport of low-salinity water ahead of the primary fluid front and triggers a transition in the flow regime from a piston-like to multi-directional displacement. This then makes a substantial difference of 22% in the ultimate oil recovery factors between the 2D and quasi-3D models. Furthermore, the interplay of solute transport through corners and wettability alteration kinetics can lead to a new oil trapping mechanism, not reported in the literature, that diminishes LSWF efficiency. While the findings of this study elucidate that LSWF does exhibit improved oil recovery compared to high-salinity waterflooding, the complicating phenomena in mixed-wet systems can significantly affect the efficiency of this method and make it less successful.

Numerical Simulation on Insoluble Surfactant Mass Transfer on Deformable Bubble Interface in a Couette Flow by Phase-Field Lattice Boltzmann Method-Finite-Difference Method Hybrid Approach

Elaborate management on bubble shape and transportations depends on the balance between multiple physical behaviors for two-phase flow with Marangoni stress and the interface mass transfer. In this paper, a new model combining PFLBM (phase-field lattice Boltzmann method) and FDM ( finite-difference method) coupling with the ghost-cell method was built. The PFLBM-FDM was validated for the high accuracy, less computational cost, and low mass loss compared to other methods. Based on the PFLBM-FDM, a surfactant-laden bubble deformed and transported in a laminar Couette flow was investigated. The deformation ratio and transportation velocity were explored with different density ratios, surface tensions, shear velocities, and diffusion coefficients. The numerical results showed that the equilibrium state of the bubble deformation was decided only by the dimensionless numbers when the Sh number was higher than 100. Moreover, the transportation velocity of the bubble can be controlled by the balance between the Marangoni stress and shear velocity. When the Sh is lower than 100, the Marangoni stress from the surfactant is not a long-range force, which only works at the early flow. Otherwise, the Marangoni stress will be a long-range force that provides a persistent force to accelerate the bubble by ∼10%. Increasing ReH will further intensify the effect. Based on all the data, a correlation of the bubble deformation including with the densities of two fluids was obtained and the error range is less 5%.

Effect of Pore Space Stagnant Zones on Interphase Mass Transfer in Porous Media, for Two-Phase Flow Conditions

Interphase mass transfer is an important solute transport process in two-phase flow in porous media. During two-phase flow, hydrodynamically stagnant and flowing zones are formed, with the stagnant ones being adjacent to the interfaces through which the interphase mass transfer happens. Due to the existence of these stagnant zones in the vicinity of the interface, the mass transfer coefficient decreases to a certain extent. There seems to be a phenomenological correlation between the mass transfer coefficient and the extent of the stagnant zone which, however, is not yet fully understood. In this study, the phase-field method-based continuous species transfer model is applied to simulate the interphase mass transfer of a dissolved species from the immobile, residual, non-aqueous phase liquid (NAPL) to the flowing aqueous phase. Both scenarios, this of a simple cavity and this of a porous medium, are investigated. The effects of flow rates on the mass transfer coefficient are significantly reduced when the stagnant zone and the diffusion length are larger. It is found that the stagnant zone saturation can be a proxy of the overall diffusion length of the terminal menisci in the porous medium system. The early-stage mass transfer coefficient continuously decreases due to the depletion of the solute in the small NAPL clusters that are in direct contact with the flowing water. The long-term mass transfer mainly happens on the interfaces associated with large NAPL clusters with larger diffusion lengths, and the mass transfer coefficient is mainly determined by the stagnant zone saturation.

Comparative investigation of a lattice Boltzmann boundary treatment of multiphase mass transport with heterogeneous chemical reactions

Multiphase reactive transport in porous media is an important component of many natural and engineering processes. In the present study, boundary schemes for the continuum species transport-lattice Boltzmann (CST-LB) mass transport model and the multicomponent pseudopotential model are proposed to simulate heterogeneous chemical reactions in a multiphase system. For the CST-LB model, a lattice-interface-tracking scheme for the heterogeneous chemical reaction boundary is provided. Meanwhile, a local-average virtual density boundary scheme for the multicomponent pseudopotential model is formulated based on the work of Li et al. [Li, Yu, and Luo, Phys. Rev. E 100, 053313 (2019)10.1103/PhysRevE.100.053313]. With these boundary treatments, a numerical implementation is put forward that couples the multiphase fluid flow, interfacial species transport, heterogeneous chemical reactions, and porous matrix structural evolution. A series of comparison benchmark cases are investigated to evaluate the numerical performance for different pseudopotential wetting boundary treatments, and an application case of multiphase dissolution in porous media is conducted to validate the present models' ability to solve complex problems. By applying the present LB models with reasonable boundary treatments, multiphase reactive transport in various natural or engineering scenarios can be simulated accurately.

New type of pore-snap-off and displacement correlations in imbibition

HYPOTHESIS

Imbibition of a fluid into a porous material involves the invasion of a wetting fluid in the pore space through piston-like displacement, film and corner flow, snap-off and pore bypassing. These processes have been studied extensively in two-dimensional (2D) porous systems; however, their relevance to three-dimensional (3D) natural porous media is poorly understood. Here, we investigate these pore-scale processes in a natural rock sample using time-resolved 3D (i.e., four-dimensional or 4D) X-ray imaging.

EXPERIMENTS

We performed a capillary-controlled drainage-imbibition experiment on an initially brine-saturated carbonate rock sample. The sample was imaged continuously during imbibition using 4D X-ray imaging to visualize and analyze fluid displacement and snap-off processes at the pore-scale.

FINDINGS

We discover a new type of snap-off that occurs in pores, resulting in the entrapment of a small portion of the non-wetting phase in pore corners. This contrasts with previously-observed snap-off in throats which traps the non-wetting phase in pore centers. We relate the new type of pore-snap-off to the pinning of fluid-fluid interfaces at rough surfaces, creating contact angles close to 90°. Subsequently, we provide correlations for displacement events as a function of pore-throat geometry. Our findings indicate that having a small throat does not necessarily favor snap-off: the key criterion is the throat radius in relation to the pore radius involved in a displacement event, captured by the aspect ratio.

Lattice Boltzmann modeling of interfacial mass transfer in a multiphase system

In the present study, a numerical model based on the lattice Boltzmann method (LBM) is proposed to simulate multiphase mass transfer, referred to as the CST-LB model. This model introduced continuum species transfer (CST) formulation by an additional collision term to model the mass transfer across the multiphase interface. The boundary condition treatment of this model is also discussed. In order to verify the applicability, the CST-LB model is combined with the pseudopotential multiphase model to simulate a series of benchmark cases, including concentration jump near the interface, gas dissolution in a closed system, species transport during drainage in a capillary tube, and multiphase species transport in the porous media. This CST-LB model can also be coupled with other multiphase LBMs since the model depends on the phase fraction field, which is not explicitly limited to specified multiphase models.

GeoChemFoam: Direct Modelling of Multiphase Reactive Transport in Real Pore Geometries with Equilibrium Reactions

GeoChemFoam is an open-source OpenFOAM-based toolbox that includes a range of additional packages that solve various flow processes from multiphase transport with interface transfer, to single-phase flow in multiscale porous media, to reactive transport with mineral dissolution. In this paper, we present a novel multiphase reactive transport solver for simulations on complex pore geometries, including microfluidic devices and micro-CT images, and its implementation in GeoChemFoam. The geochemical model includes bulk and surface equilibrium reactions. Multiphase flow is solved using the Volume-Of-Fluid method, and the transport of species is solved using the continuous species transfer method. The reactive transport equations are solved using a sequential operator splitting method, with the transport step solved using GeoChemFoam, and the reaction step solved using Phreeqc, the US geological survey’s geochemical software. The model and its implementation are validated by comparison with analytical solutions in 1D and 2D geometries. We then simulate multiphase reactive transport in two test pore geometries: a 3D pore cavity and a 3D micro-CT image of Bentheimer sandstone. In each case, we show the pore-scale simulation results can be used to develop upscaled models that are significantly more accurate than standard macro-scale equilibrium models.