An accelerated process-based method for the accurate computation of relative permeability from direct simulations of two-phase flow on micro-computed tomography images of porous media

Pore-scale forces have significant effects on the macroscopic behavior of multi-phase flow through porous media. We develop a robust and accurate accelerated process-based method for the computation of relative permeability from direct simulations of pore-scale two-phase flow on micro-computed tomography images. In the pressure drop calculation, we take advantage of an existing analysis that establishes a relationship between pore-scale forces and Darcy-scale pressure drops using an energy-conservation approach. We establish a thermodynamically consistent approximation of Darcy-scale viscous pressure drops as the rate of energy dissipation per unit flow rate of each flowing phase for the first time within the context of a free-energy lattice Boltzmann method (LBM). In addition, we propose and test a new computationally efficient partial-mirror periodic boundary condition for a fully coupled visco-capillary pore-scale flow simulator based on a free-energy LBM. The new boundary condition is imposed only in the main flow direction and significantly reduces the computational cost of the process-based relative-permeability computation protocol at a small compromise on accuracy.

We first compute primary-drainage and subsequent imbibition relative permeability curves for a reservoir sandstone rock sample. We use this real-reservoir dataset to validate the pressure drop computation method and the partial-mirror periodic boundary condition. We then simulate the entire drainage and imbibition cycle using an extensively studied Berea sandstone dataset. We quantitatively demonstrate that pore-scale direct numerical simulation-based relative permeability curves computed with our novel process-based method agree well with experimental steady-state relative permeability measurements. We also quantitatively demonstrate that the new partial-mirror periodic boundary condition accelerates the relative permeability computations 4 to 13 times.