The investigation of fluid flow channeling and viscous fingering during immiscible two-phase displacement in subsurface porous media is crucial for optimizing CO2 geological sequestration and improving hydrocarbon recovery. In this study, we develop a pore-scale numerical framework for unsteady state immiscible displacement based on a body-centered cubic percolation network, which explicitly captures the coupled effects of pore-scale heterogeneity, capillary number, and unfavorable viscosity ratio on flow channeling and viscous fingering. The simulations reveal that viscous fingering and flow channeling preferentially occur along overlapping high conductivity pathways that conform to the minimum energy dissipation principle. Along these preferential routes, the local balance between viscous and capillary forces governs the stability of the two-phase interface and gives rise to distinct patterns and intensities of viscous fingering in the invading phase. Building on these insights, we establish a theoretical framework that quantifies how the critical pore radius and capillary number control the onset and growth of interfacial instabilities during immiscible displacement. The model demonstrates that lowering the injection rate, and hence, the effective capillary number, suppresses viscous fingering, leading to more stable displacement fronts. These findings provide practical guidance for the design of injection schemes, helping to enhance oil and gas recovery and improve the storage efficiency and security of CO2 geological sequestration projects.
Yang et al. (Mon,) studied this question.