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.

Upscaling of air–water system and capillary entrapment of methane in porous media

Some Challenges Posed by Coal Bed Methane Regional Assessment Modeling

Numerical Investigation of Wellbore Methane Leakage From a Dual‐Porosity Reservoir and Subsequent Transport in Groundwater

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

A Petroleum Vapor Intrusion Model Involving Upward Advective Soil Gas Flow Due to Methane Generation

In Situ Characterization of Three-Phase Flow in Mixed-Wet Porous Media Using Synchrotron Imaging

MATHEMATICAL SIMULATION FOR SUBMARINE GAS HYDRATE FORMATION: UPON THE ASSUMPTION OF UPWARD ADVECTION OF METHANE-BEARING POREWATER

Evaluation of the sealing capacity of overlying fine-grained sediments in the process of hydrate exploitation by thermal stimulation with warm water

Upscaling multiphase flow through heterogeneous porous media

Effect of Hydrate on Gas/Water Relative Permeability of Hydrate-Bearing Sediments: Pore-Scale Microsimulation by the Lattice Boltzmann Method

Analytical Model for Methane Migration Through Fractured Unsaturated Landfill Cover Soil

Non-isothermal, multi-phase, multi-component flows through deformable methane hydrate reservoirs

A Novel Theoretical Method for Upscaling Permeability in Hydrate‐Bearing Sediments

Role of capillary pressure on carbon dioxide sequestration: mathematical analysis

Permeability of Laboratory-Formed Porous Media Containing Methane Hydrate: Observations Using X-ray Computed Tomography and Simulations with Pore Network Models

Molecular simulation studies on the water/methane two-phase flow in a cylindrical silica nanopore: Formation mechanisms of water lock and implications for gas hydrate exploitation

A Coupled Model for Upscaling Water Flow in a Shale Matrix System from Pore Scale to Representative Elementary Area Scale

Simulation of bench-scale CO2 injection using a coupled continuum-discrete approach

Mechanisms for Upward Migration of Methane in Marine Sediments

Measurements of the Capillary Pressure-Saturation Relationship of Methane Hydrate Bearing Sediments

Role of heterogeneous surface wettability on dynamic immiscible displacement, capillary pressure, and relative permeability in a CO2-water-rock system