Phase diagram and permeability evolution for dissolving vertical fractures in a gravity field

Rock fractures always provide the main flow pathways in geological systems. When exposed to reactive flow, mineral dissolution causes different dissolution regimes and significantly impacts permeability change. Previous studies have examined dissolution regimes and permeability evolution in horizontal geometry, but dissolution dynamics and permeability changes of vertical fractures in a gravity field remain unexplored. Here, we conduct flow-dissolution experiments combined with a pore-scale modeling approach to investigate regime transition and permeability evolution in vertical fractures. We validate the modeling approach using experimental results and illustrate the importance of gravitational effects in dissolving vertical fractures. Our 3D numerical simulations reveal the regime transition, represented by a critical Richardson number Ri, from the forced convection regime to the buoyancy-driven convection regime as the importance of gravity increases. In the buoyancy-driven regime, gravitational effects promote the development of localized dissolution/flow channels. However, gravitational effects become negligible in the forced convection regime where wormholes and uniform patterns can be separated by the critical Pe number. Using the critical Ri and Pe numbers, we establish a phase diagram of dissolution regimes for vertical fractures in a gravity field, exhibiting good agreement with experiments. By correlating this phase diagram with permeability-aperture curves, we directly estimate the power law exponent of permeability evolution for each dissolution regime. This approach extends the classic phase diagram of dissolution regimes to consider gravitational effects on vertical fractures and provides a guide to directly estimate the parameters for permeability evolution, which is essential in reactive-transport modeling at the field scale.