Characterization of immiscible fluid displacement processes with various capillary numbers and viscosity ratios in 3D natural sandstone

To characterize the influence of reservoir conditions upon multiphase flow, we calculated fluid displacements (drainage processes) in 3D pore spaces of Berea sandstone using two-phase lattice Boltzmann (LB) simulations. The results of simulations under various conditions were used to classify the resulting two-phase flow behavior into three typical fluid displacement patterns on the diagram of capillary number (Ca) and viscosity ratio of the two fluids (M). In addition, the saturation of the nonwetting phase was calculated and mapped on the Ca–M diagram. We then characterized dynamic pore-filling events (i.e., Haines jumps) from the pressure variation of the nonwetting phase, and linked this behavior to the occurrence of capillary fingering. The results revealed the onset of capillary fingering in 3D natural rock at a higher Ca than in 2D homogeneous granular models, with the crossover region between typical displacement patterns broader than in the homogeneous granular model. Furthermore, saturation of the nonwetting phase mapped on the Ca–M diagram significantly depends on the rock models. These important differences between two-phase flow in 3D natural rock and in 2D homogeneous models could be due to the heterogeneity of pore geometry in the natural rock and differences in pore connectivity. By quantifying two-phase fluid behavior in the target reservoir rock under various conditions (e.g., saturation mapping on the Ca–M diagram), our approach could provide useful information for investigating suitable reservoir conditions for geo-fluid management (e.g., high CO₂ saturation in CO₂ storage).