The LBPM software package for simulating multiphase flow on digital images of porous rocks

TL;DR

The LBPM software package enables detailed pore-scale simulations of multiphase flow in digital rock images, replicating laboratory experiments to analyze flow behavior, wettability effects, and transport properties in porous media.

Contribution

This work introduces open-source protocols and tools within LBPM for simulating and analyzing multiphase flow on digital rocks, integrating morphological analysis and efficient workflows.

Findings

Successfully simulated flow in sandstone matching expected wettability trends.

Developed workflows that combine morphological tools with lattice Boltzmann models for efficiency.

Provided methods to infer relative permeability and capillary curves from simulations.

Abstract

Direct pore scale simulations of two-fluid flow on digital rock images provide a promising tool to understand the role of surface wetting phenomena on flow and transport in geologic reservoirs. We present computational protocols that mimic conventional special core analysis laboratory (SCAL) experiments, which are implemented within the open source LBPM software package. Protocols are described to simulate unsteady displacement, steady-state flow at fixed saturation, and to mimic centrifuge experiments. These methods can be used to infer relative permeability and capillary curves, and otherwise understand two-fluid flow behavior based on first principles. Morphological tools are applied to assess image resolution, establish initial conditions, and instantiate surface wetting maps based on the distribution of fluids. Internal analysis tools are described that measure essential aspects of two-fluid flow, including fluid connectivity and surface measures, which are used to track transient aspects of the flow behavior as they occur during simulation. Computationally efficient workflows are developed by combining these components with a two-fluid lattice Boltzmann model to define hybrid methods that can accelerate computations by using morphological tools to incrementally evolve the pore-scale fluid distribution. We show that the described methods can be applied to recover expected trends due to the surface wetting properties based on flow simulation in Benntheimer sandstone.