CO2 sequestration projects will in the coming years include both aquifer and depleted oil and gas field sites, with different stress paths and history. Stress changes and stress concentration effects on faults will have to be readily assessed, potentially endangering shallower permeable formations. Usually, a fault is modeled as a singularity with shear strength or friction properties, and simulations are run to determine whether the fault is reactivated through shear failure. In this paper, we model a simple rectilinear fault as a finite surface with lowered mechanical properties compared to elsewhere in the domain, which represents a fractured zone alongside the fault core. SINTEF’s Modified Discrete Element code is used coupled to the flow simulator TOUGH2, to model the fracture initiation and propagation, monitoring the permeability increase along the fault. A simplified scenario is simulated, with a sandstone storage reservoir bounded by a fault, penetrating a shale caprock to a shallower sandstone layer. The storage site either undergoes depletion before CO2 injection or has its pore pressure increased to simulate the case of aquifer storage. Results show that during depletion, shear stresses may develop such that fractures propagate alongside the fault to the upper aquifer. However, for the mirror fault orientation with regards to verticality, no such fractures develop. These results are reversed for the aquifer storage case.
Read full abstract