Capillary desaturation curve for residual nonwetting phase in natural fractures

The displacement of a nonwetting phase by a wetting phase is characterized by the capillary number. Different forms of capillary number have been used in the literature for flow in porous media. A capillary number for a single rock fracture has been defined in the literature, using the mean aperture to characterize the trapping and mobilization in a fracture. We propose a new capillary-number definition for fractures that incorporates geometrical characterization of the fracture, dependent on the force balance on a trapped ganglion. The new definition is validated with laboratory experiments using five distinctive model fractures. The model fractures are made of glass plates, with a wide variety of hydraulic apertures, degrees of roughness, and correlation lengths of the roughness. The fracture surfaces were characterized in detail and statistically analyzed. The aperture distribution of each model fracture was represented as a 2D network of pore bodies connected by throats. The hydraulic aperture of each model fracture was measured experimentally. Capillary desaturation curves (CDCs) were generated experimentally using water/air in forced imbibition. The transparent nature of the system permits us to determine the residual air saturation as a function of pressure gradient from the captured images. The residual nonwetting saturation/capillary-number relationship obtained from different fractures varying in aperture and roughness can be represented approximately by a single curve in terms of the new definition of the capillary number. They do not fit a single trend using the conventional definition of the capillary number.