The wetting properties of pore walls have a strong effect on multiphase flow through porous media. However, the fluid flow behaviour in porous materials with both complex pore structures and non-uniform wettability are still unclear. Here, we performed unsteady-state quasi-static oil- and waterflooding experiments to study multiphase flow in two sister heterogeneous sandstones with variable wettability conditions (i.e. one natively water-wet and one chemically treated to be mixed-wet). The pore-scale fluid distributions during this process were imaged by laboratory-based X-ray micro-computed tomography (micro-CT). In the mixed-wet case, we observed pore filling events where the fluid interface appeared to be at quasi-equilibrium at every position along the pore body (13% by volume), in contrast to capillary instabilities typically associated with slow drainage or imbibition. These events corresponded to slow displacements previously observed in unsteady-state experiments, explaining the wide range of displacement time scales in mixed-wet samples. Our new data allowed us to quantify the fluid saturations below the image resolution, indicating that slow events were caused by the presence of microporosity and the wetting heterogeneity. Finally, we investigated the sensitivity of the multi-phase flow properties to the slow filling events using a state-of-the-art multi-scale pore network model. This indicated that pores where such events took place contributed up to 19% of the sample’s total absolute permeability, but that the impact on the relative permeability may be smaller. Our study sheds new light on poorly understood multiphase fluid dynamics in complex rocks, of interest to e.g. groundwater remediation and subsurface CO2 storage.

Representative elementary volumes (REVs) are an important concept in studying subsurface multiphase flow at the continuum scale. However, fluctuations in multiphase flow are currently not represented in continuum scale models, and their impact at the REV-scale is unknown. Previous pore-scale imaging studies on these fluctuations were limited to small samples with mm-scale diameters and volumes on the order of ~ 0.5 cm3. Here, we image steady-state co-injection experiments on a one-inch diameter core plug sample, with nearly two orders of magnitude larger volume (21 cm3), while maintaining a pore-scale resolution with X-ray micro-computed tomography. This was done for three total flow rates in a series of drainage fractional flow steps. Our observations differ markedly from those reported for mm-scale samples in two ways: the macroscopic fluid distribution was less ramified at low capillary numbers (Ca) of 10-7; and the volume fraction of intermittency initially increased with increasing Ca (similar to mm-scale observations), but then decreased at Ca of 10-7. Our results suggest that viscous forces may play a role in the cm-scale fluid distribution, even at such low Ca, dampening intermittent pathway flow. A REV study of the fluid saturation showed that this may be missed in smaller-scale samples. Pressure drop measurements suggest that the observed pore-scale fluctuations resulted in non-Darcy like upscaled behavior. Overall, we show the importance of large field-of-view high-resolution imaging to bridge the gap between pore- and continuum-scale multiphase flow studies, in particular of pore-scale fluctuations.

Many subsurface fluid flows, including the storage of CO underground or the production of oil, are transient processes incorporating multiple fluid phases. The fluids are not in equilibrium meaning macroscopic properties such as fluid saturation and pressure vary in space and time. However, these flows are traditionally modelled with equilibrium (or steady-state) flow properties, under the assumption that the pore scale fluid dynamics are equivalent. In this work, we used fast synchrotron X-ray tomography with 1s time resolution to image the pore scale fluid dynamics as the macroscopic flow transitioned to steady-state. For nitrogen or decane, and brine injected simultaneously into a porous rock we observed distinct pore scale fluid dynamics during transient flow. Transient flow was found to be characterised by intermittent fluid occupancy, whereby flow pathways through the pore space were constantly rearranging. The intermittent fluid occupancy was largest and most frequent when a fluid initially invaded the rock. But as the fluids established an equilibrium the dynamics decreased to either static interfaces between the fluids or small-scale intermittent flow pathways, depending on the capillary number and viscosity ratio. If the fluids were perturbed after an equilibrium was established, by changing the flow rate, the transition to a new equilibrium was quicker than the initial transition. Our observations suggest that transient flows require separate modelling parameters. The timescales required to achieve equilibrium suggest that several metres of an invading plume front will have flow properties controlled by transient pore scale fluid dynamics.

Image-based pore-scale modeling is an important method to study multiphase flow in permeable rocks. However, in many rocks, the pore size distribution is so wide that it cannot be resolved in a single pore-space image, typically acquired using micro-computed tomography (micro-CT). Recent multi-scale models therefore incorporate sub-voxel porosity maps, created by differential micro-CT imaging of a contrast fluid in the pores. These maps delineate different microporous flow zones in the model, which must be assigned petrophysical properties as input. The uncertainty on the pore scale physics in these models is therefore heightened by uncertainties on the representation of unresolved pores, also called sub-rock typing. Here, we address this by validating a multi-scale pore network model using a drainage experiment imaged with differential micro-CT on an Estaillades limestone sample. We find that porosity map-based sub-rock typing was unable to match the micrometer-scale experimental fluid distributions. To investigate why, we introduce a novel baseline sub-rock typing method, based on a 3D map of the experimental capillary pressure function. By incorporating this data, we successfully remove most of the sub-rock typing uncertainty from the model, obtaining a close fit to the experimental fluid distributions. Comparison between the two methods shows that in this sample, the porosity map is poorly correlated to the multiphase flow behavior of the microporosity. The validation method introduced in this paper serves to separate and address the uncertainties in multi-scale models, facilitating simulations in complex geological reservoir rocks important for e.g. geological storage of CO2 and renewable energy.