Summary
Subsurface injection of fluids in a stress‐sensitive naturally fractured rock faces the problem of near‐wellbore fracture evolution and associated changes in rock properties. Numerical modeling of the changes in permeability and poroelastic properties of the near‐wellbore region is challenging due to the coupling between fracture dynamics and poromechanics across multiple length scales of fractures and the host rock. We present a numerical framework to model anisotropic and dynamic evolution in rock properties based on a coupled formulation of fluid flow, rock mechanics, and fracturing, where fracturing‐induced damage is used to update the rock properties. A generalized fixed‐stress method, which accounts for damage‐induced anisotropy in flow and deformation processes, is developed to sequentially solve the equations of flow, mechanics, and fracture evolution. We demonstrate the usefulness of our framework in quantifying the effects of injection rate variability, initial fracture distribution, and in‐situ stress state on the evolution in permeability, elastic stiffness, and the Biot parameters. Our framework does not require the computational mesh to conform to existing or future fractures, allows simultaneous growth of multiple randomly distributed fractures, and can be implemented relatively easily in existing coupled flow‐geomechanics simulators to extend them to model fracturing at the reservoir scale.