Structural barriers to complete homogenization and wormholing in dissolving porous and fractured rocks

TL;DR

This study investigates how structural heterogeneity in porous and fractured rocks influences dissolution patterns, revealing fundamental limits on flow homogenization and wormholing through network modeling and a unified flow focusing metric.

Contribution

It introduces a unified flow focusing metric to quantify dissolution patterns across different network models, highlighting the persistent influence of initial heterogeneity.

Findings

Uniform dissolution cannot fully homogenize flow due to inherent network heterogeneity.

Flow focusing propagates from the inlet in wormholing regimes.

Pre-existing flow paths enlarge during channeling regimes.

Abstract

Dissolution in porous media and fractured rocks alters both the chemical composition of the fluid and the physical properties of the solid. Depending on system conditions, reactive flow may enlarge pores uniformly, widen pre-existing channels, or trigger instabilities that form wormholes. The resulting pattern reflects feedbacks among advection, diffusion, surface reaction, and the initial heterogeneity of the medium. Porous and fractured media can exhibit distinct characteristics -- for example, the presence of large fractures can significantly alter the network topology and overall connectivity of the system. We quantify these differences with three network models -- a regular pore network, a disordered pore network, and a discrete fracture network -- evaluated with a unified metric: the flow focusing profile. This metric effectively captures evolution of flow paths across all systems: it reveals a focusing front that propagates from the inlet in the wormholing regime, a system-wide decrease in focusing during uniform dissolution, and the progressive enlargement of pre-existing flow paths in the channeling regime. The metric shows that uniform dissolution cannot eliminate heterogeneity resulting from the network topology. This structural heterogeneity -- rather than just pore-diameter or fracture-aperture variance -- sets a fundamental limit on flow homogenization and must be accounted for when upscaling dissolution kinetics from pore or fracture scale to the reservoir level.