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.

在数字岩石图像上直接对两种流体进行孔隙尺度模拟，为了解地表润湿现象在地质储层中对流动和传输的作用提供了有前途的工具。我们提供了模拟常规特殊核心分析实验室（SCAL）实验的计算协议，这些协议在开源LBPM软件包中实现。协议被描述为模拟非稳态位移，固定饱和状态下的稳态流以及模拟离心机实验。这些方法可用于推断相对渗透率和毛细曲线，以及另外基于第一原理理解双流体流动行为。应用形态学工具评估图像分辨率，建立初始条件并根据流体的分布实例化表面润湿图。描述了内部分析工具，该工具可测量双流体流动的基本方面，包括流体连通性和表面测量，用于跟踪在模拟过程中发生的流动行为的瞬态变化。通过将这些组件与双流体格子Boltzmann模型相结合来开发计算有效的工作流程，以定义可通过使用形态学工具逐步演化孔隙尺度流体分布来加速计算的混合方法。我们表明，所描述的方法可以应用于基于Benntheimer砂岩流模拟的表面润湿特性来恢复预期趋势。它们用于跟踪模拟过程中发生的流动行为的瞬态变化。通过将这些组件与双流体格子Boltzmann模型相结合来开发计算有效的工作流程，以定义可通过使用形态学工具逐步演化孔隙尺度流体分布来加速计算的混合方法。我们表明，所描述的方法可以应用于基于Benntheimer砂岩流模拟的表面润湿特性来恢复预期趋势。它们用于跟踪模拟过程中发生的流动行为的瞬态变化。通过将这些组件与双流体格子Boltzmann模型相结合来开发计算有效的工作流程，以定义可通过使用形态学工具逐步演化孔隙尺度流体分布来加速计算的混合方法。我们表明，所描述的方法可以应用于基于Benntheimer砂岩流模拟的表面润湿特性来恢复预期趋势。