Pore-scale numerical simulation of low salinity water flooding using the lattice Boltzmann method

Bridging molecular-scale interfacial science with continuum-scale models

Solid-water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid-water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

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.

Pore-scale modeling of two-phase flow: A comparison of the generalized network model to direct numerical simulation

Despite recent advances in pore-scale modeling of two-phase flow through porous media, the relative strengths and limitations of various modeling approaches have been largely unexplored. In this work, two-phase flow simulations from the generalized network model (GNM) [Phys. Rev. E 96, 013312 (2017)2470-004510.1103/PhysRevE.96.013312; Phys. Rev. E 97, 023308 (2018)2470-004510.1103/PhysRevE.97.023308] are compared with a recently developed lattice-Boltzmann model (LBM) [Adv. Water Resour. 116, 56 (2018)0309-170810.1016/j.advwatres.2018.03.014; J. Colloid Interface Sci. 576, 486 (2020)0021-979710.1016/j.jcis.2020.03.074] for drainage and waterflooding in two samples-a synthetic beadpack and a micro-CT imaged Bentheimer sandstone-under water-wet, mixed-wet, and oil-wet conditions. Macroscopic capillary pressure analysis reveals good agreement between the two models, and with experiments, at intermediate saturations but shows large discrepancy at the end-points. At a resolution of 10 grid blocks per average throat, the LBM is unable to capture the effect of layer flow which manifests as abnormally large initial water and residual oil saturations. Critically, pore-by-pore analysis shows that the absence of layer flow limits displacement to invasion-percolation in mixed-wet systems. The GNM is able to capture the effect of layers, and exhibits predictions closer to experimental observations in water and mixed-wet Bentheimer sandstones. Overall, a workflow for the comparison of pore-network models with direct numerical simulation of multiphase flow is presented. The GNM is shown to be an attractive option for cost and time-effective predictions of two-phase flow, and the importance of small-scale flow features in the accurate representation of pore-scale physics is highlighted.

Discontinuous dewetting dynamics of highly viscous droplets on chemically heterogeneous substrates

HYPOTHESIS

Droplet spreading on heterogeneous (chemical/structural) surfaces has revealed local disturbances that affect the advancing contact line. With droplet dewetting being less studied, we hypothesize that a receding droplet can be perturbed by localized heterogeneity which leads to irregular and discontinuous dewetting of the substrate.

EXPERIMENTS

The sessile drop method was used to study droplet dewetting at a wettability boundary. One-half of a hydrophilic surface was hydrophobically modified with either i) methyloctyldichlorosilane or ii) clustered macromolecules. A Lattice Boltzmann method (LBM) simulation was also developed to determine the effect of contact angle hysteresis and boundary conditions on the droplet dynamics.

FINDINGS

The two surface treatments were optimized to produce comparable water wetting characteristics. With a negative Gibbs free energy on the hydrophilic-half, the oil droplet receded to the hydrophobic-half. On the silanized surface, the droplet was pinned and the resultant droplet shape was a distorted spherical cap, having receded uniformly on the unmodified surface. Modifying the surface with clustered macromolecules, the droplet receded slightly to form a spherical cap. However, droplet recession was non-uniform and daughter droplets formed near the wettability boundary. The LBM simulation revealed that daughter droplets formed when θR > 164°, with the final droplet shape accurately described by imposing a diffuse wettability boundary condition.

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.

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.

Using energy balance to determine pore-scale wettability

HYPOTHESIS

Based on energy balance during two-phase displacement in porous media, it has recently been shown that a thermodynamically consistent contact angle can be determined from micro-tomography images. However, the impact of viscous dissipation on the energy balance has not been fully understood. Furthermore, it is of great importance to determine the spatial distribution of wettability. We use direct numerical simulation to validate the determination of the thermodynamic contact angle both in an entire domain and on a pore-by-pore basis.

SIMULATIONS

Two-phase direct numerical simulations are performed on complex 3D porous media with three wettability states: uniformly water-wet, uniformly oil-wet, and non-uniform mixed-wet. Using the simulated fluid configurations, the thermodynamic contact angle is computed, then compared with the input contact angles.

FINDINGS

The impact of viscous dissipation on the energy balance is quantified; it is insignificant for water flooding in water-wet and mixed-wet media, resulting in an accurate estimation of a representative contact angle for the entire domain even if viscous effects are ignored. An increasing trend in the computed thermodynamic contact angle during water injection is shown to be a manifestation of the displacement sequence. Furthermore, the spatial distribution of wettability can be represented by the thermodynamic contact angle computed on a pore-by-pore basis.

Evaluation of methods using topology and integral geometry to assess wettability

HYPOTHESIS

The development of high-resolution in situ imaging has allowed contact angles to be measured directly inside porous materials. We evaluate the use of concepts in integral geometry to determine contact angle. Specifically, we test the hypothesis that it is possible to determine an average contact angle from measurements of the Gaussian curvature of the fluid/fluid meniscus using the Gauss-Bonnet theorem.

THEORY AND SIMULATION

We show that it is not possible to unambiguously determine an average contact angle from the Gauss-Bonnet theorem. We instead present an approximate relationship: 2πn(1-cosθ)=4π-∫κ_G12dS₁₂, where n is the number of closed loops of the three-phase contact line where phases 1 and 2 contact the surface, θ is the average contact angle, while κ_G12 is the Gaussian curvature of the fluid meniscus which is integrated over its surface S₁₂. We then use the results of pore-scale lattice Boltzmann simulations to assess the accuracy of this approach to determine a representative contact angle for two-phase flow in porous media.

FINDINGS

We show that in simple cases with a flat solid surface, the approximate expression works well. When applied to simulations on pore space images, the equation provides a robust estimate of contact angle, accurate to within 3°, when averaged over many fluid clusters, although individual values can have significant errors because of the approximations used in the calculation.

A Laboratory Approach on the Hybrid-Enhanced Oil Recovery Techniques with Different Saline Brines in Sandstone Reservoirs

As foams are not thermodynamically stable and might be collapsed, foam stability is defined by interfacial properties and bulk solution. In this paper, we investigated foam injection and different salinity brines such as NaCl, CaCl2, KCl, and MgCl2 to measure cumulative oil production. According to the results of this experiment, it is concluded that sequential low-salinity water injections with KCl and foam flooding have provided the highest cumulative oil production in sandstone reservoirs. This issue is related to high wettability changes that had been caused by the KCl. As K+ is a monovalent cation, KCl has the highest wettability changes compared to other saline brines and formation water at 1000 ppm, which is due to the higher wettability changes of potassium (K+) over other saline ions. The interfacial tension for KCl at the lowest value is 1000 ppm and, for MgCl2, has the highest value in this concentration. Moreover, the formation brine, regarding its high value of salty components, had provided lower cumulative oil production before and after foam injection as it had mobilized more in the high permeable zones and, therefore, large volumes of oil would be trapped in the small permeable zones. This was caused by the low wettability alteration of the formation brine. Thereby, formation water flowed in large pores and the oil phase remained in small pores and channels. On the other hand, as foams played a significant role in the mobility control and sweep efficiency, at 2 pore volume, foam increased the pressure drop dramatically after brine injection. Consequently, foam injection after KCl brine injection had the maximum oil recovery factor of 63.14%. MgCl2 and formation brine had 41.21% and 36.51% oil recovery factor.