Modeling Immiscible Two‐Phase Flow in Rough Fractures From Capillary to Viscous Fingering

Abstract : We develop an efficient computational model for simulating fluid invasion patterns emerging in variable aperture fractures. This two-dimensional model takes into account the effect of capillary force on the fluid-fluid interfaces and viscous pressure drop in both fluid phases. The pressure distribution is solved at each time step based on mass balance and local cubic law, considering an imposed pressure jump condition at the fluid-fluid interface. This pressure jump corresponds to the Laplace pressure which includes both terms related to the out-of-plane (aperture-spanning) curvature and to the in-plane curvature. Simulating a configuration that emulates viscous fingering in two-dimensional random porous media confirms that the model accounts properly for the role of viscous forces. Furthermore, direct comparison with previously obtained experimental results shows that the model reproduces the observed drainage patterns in a rough fracture reasonably well. The evolutions of tip location, the inlet pressures, and the invading phase fractal dimensions are analyzed to characterize the transition from capillary fingering to viscous fingering regimes. A radial injection scenario of immiscible invasion is also studied with varying modified capillary number and viscosity ratio, showing displacement patterns ranging from capillary fingering to viscous fingering to stable displacement. Such simulations using two contact angles show that the invading phase becomes more compact when the wetting condition changes from strong to weak drainage, as already observed in 2-D porous media. The model can be used to bridge the gap in spatial scales of two-phase flow between pore-scale modeling approaches and the continuum Darcy-scale models.