Analysis of permeability change in dissolving rough fractures using depth-averaged flow and reactive transport models
Many geophysical problems and engineering applications involve fluid flow through evolving faults and fractures. These situations are typical in CO 2 injection scenarios where supercritical CO 2 or CO 2 –saturated fluids are able to dissolve the solid matrix. Dissolution or precipitation within a fracture could lead to self-reinforcing or self-limiting behavior due to feedbacks from changes in permeability (k). There is limited predictive understanding of this phenomenon. Here, we investigated the case of fracture dissolution using a two-dimensional depth-averaged reactive transport model. The model simulated an acidic fluid dissolving a three-dimensional rough fracture. We observed different dissolution patterns across a wide range of geometric properties and different Peclet and Damkohler numbers (960 cases total). Despite this variation, a persistent linear increase in k with time (t) roughly manifested once the dissolution front advanced through the fracture. The time rate of change of k (k/t) is roughly constant for all simulations. The modeled rate is consistent with a theoretical value based on the parallel plates model and with experimentally derived values. Finally, the results confirm the common presence of linear time-dependence of fracture k due to self-reinforcement.