Abstract:This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO2 (scCO(2))-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 x 800 x 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO(2)). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO(2) and showed that the variations in the relative permeability curves for both water and scCO(2) increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. The impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

Theoretical and numerical investigation into roles of geofluid flow in ore forming systems: Integrated mass conservation and generic model approach

Theoretical Modeling of Stable Isotope Systems and Its Applications to Geochemistry of Hydrothermal ore Deposits

Theoretical Modelling of Water-Rock Δd-Δ 18 O Isotopic Exchange System and Source of Ore-Forming Fluid: A Case Study on Jinduicheng Superlarge-Scale Molybdenum Deposit, Central China

Investigating Effects of Structural Deformation Regimes on Mineralization Distributions in Fluid-Saturated Rocks: Computational Simulation Approach through Generic Models

Numerical Modeling of Structurally Controlled Ore Formation in Magmatic‐Hydrothermal Systems

Computational simulation of the hydrothermal ore-forming system associated with the Luziyuan Pb-Zn Deposit, Yunnan, China: An example of pore-fluid convection controlled mineralization

The Role of the Spatial Heterogeneity and Correlation Length of Surface Wettability on Two-Phase Flow in a CO2-water-rock System

FEM‐based dual length‐scale simulation of hydrothermal ore‐forming systems involving convective flow in fluid‐saturated porous media

Contrasting regimes of Cu, Zn and Pb transport in ore-forming hydrothermal fluids

Computationally simulating the hydrothermal mineralizing system involved in the Laochang Pb-Zn deposit, Gejiu ore district, Yunnan, China: an example of pore-fluid convection controlled mineralization

Combination model-based numerical simulation of the mineralizing processes within iron oxide-apatite systems

The Coupled THMC finite-element modeling of hydrothermal systems: Insights into the Jiama porphyry metallogenic system

Geochemistry of metamorphic fluids: From static-qualitative analysis to dynamic-quantitative study

Multiphase Reactive Transport Modeling – the Impacts of Two-Phase Flow on Mineral Reaction Rate

Novel Modeling of Non-Isothermal Flow-Induced Fine Particle Migration in Porous Media Based on the Derjaguin-Landau-Verwey-Overbeek Theory

Influence of different tectonic settings on fracture formation and fluid flow around upper-crustal magmatic intrusion: insights from numerical modelling

Ore-system as a geodynamic probe

A Perspective on Applied Geochemistry in Porous Media: Reactive Transport Modeling of Geochemical Dynamics and the Interplay with Flow Phenomena and Physical Alteration

Geochemical Kinetics of Oxygen Isotope Exchange During Fluid-Rock Interactions: A Review

Modeling of Multiphysical–Chemical Coupling for Coordinated Mining of Coal and Uranium in a Complex Hydrogeological Environment

A review of interaction mechanisms and microscopic simulation methods for CO 2 -water-rock system