Empirical or theoretical extensions of Darcy's law for immiscible two‐phase flow have shown significant limitations in properly modeling the flow at the continuum scale. We tackle this problem by proposing a set of upscaled equations based on pore‐scale flow regimes, that is, the topology of flowing phases. The incompressible Navier‐Stokes equation is upscaled by means of multiple‐scale expansions and its closures derived from the mechanical energy balance for different flow regimes at the pore scale. We also derive the applicability conditions of the upscaled equations based on the order of magnitude of relevant dimensionless numbers, that is, Eotvos, Reynolds, capillary, Froude numbers, and the viscosity and density ratio of the system, as well as a set of closures valid for the basic flow regimes of low Eotvos number systems, that is, core‐annular and plug and drop traffic flows. We provide analytical expressions for the relative permeability of the wetting and nonwetting phases in different flow regimes and demonstrate that the effect of the flowing‐phases topology on the relative permeabilities is significant. Finally, we show that the classical two‐phase Darcy law is recovered for a limited range of operative conditions, while specific terms accounting for interfacial and wall interactions should be incorporated to accurately model ganglia or drop traffic flow.

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). We injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

The employment of 2D models to investigate the properties of 3D flows in porous media is ubiquitous in the literature. The limitations of such approaches are often overlooked. Here, we assess to which extent 2D flows in porous media are suitable representations of 3D flows. To this purpose, we compare representative elementary volume (REV) scales obtained by 2D and 3D numerical simulations of flow in porous media. The stationarity of several quantities, namely porosity, permeability, mean and variance of velocity, is evaluated in terms of both classical and innovative statistics. The variance of velocity, strictly connected to the hydrodynamic dispersion, is included in the analysis in order to extend conclusions to transport phenomena. Pore scale flow is simulated by means of a Lattice Boltzmann model. The results from pore scale simulations point out that the 2D approach often leads to inconsistent results, due to the profound difference between 2D and 3D flows through porous media. We employ the error in the evaluation of REV as a quantitative measure for the reliability of a 2D approach. Moreover, we show that the acceptance threshold for a 2D representation to be valid strongly depends on which flow/transport quantity is sought.

Relative permeability measurements are commonly fitted with the Corey and Brooks‐Corey correlations. Despite such correlations fit experimental data generally well, they are mostly of an empirical nature. Here, we propose a semiempirical model to determine relative permeabilities of the wetting and the nonwetting phases in real 3‐D porous media that accounts for pore‐scale flow regimes. The starting point is the homogenization framework proposed by Picchi and Battiato (2018, https://doi.org/10.1029/2018WR023172), where the upscaling is conducted for different spatial distributions of the flowing phases in the capillary tube setting. First, we extend the approach to realistic media by allowing pore‐scale flow regimes to coexist in a complex geometry while accounting for capillary and viscous limits in the dynamics. Then, we discuss the scaling behavior of normalized relative permeabilities in terms of the phases viscosity ratio and identify three classes which govern their scaling. We also derive an analytical expression for the fractional flow. Finally, we provide a detailed validation of the proposed model for both relative permeabilities and fractional flow against data from numerical simulations and experiments available in the literature. The data set used for validation covers a wide range of systems, ranging from brine‐CO₂ to oil‐water flows. The equations derived capture well the trend of both numerical and experimental data.