Induced seismicity constraints on subsurface geological structure, Paradox Valley, Colorado

• Published in: Geophysical journal international 2015-02, Vol.200 (2), p.1172-1195
• Main Authors: Block, Lisa V., Wood, Christopher K., Yeck, William L., King, Vanessa M.
• Format: Article
• Language: English
• ISSN: 0956-540X; 1365-246X
• DOI: 10.1093/gji/ggu459

Precise relative hypocentres of seismic events induced by long-term fluid injection at the Paradox Valley Unit (PVU) brine disposal well provide constraints on the subsurface geological structure and compliment information available from deep seismic reflection and well data. We use the 3-D spatial distribution of the hypocentres to refine the locations, strikes, and throws of subsurface faults interpre­ted previously from geophysical surveys and to infer the existence of previously unidentified subsurface faults. From distinct epicentre lineations and focal mechanism trends, we identify a set of conjugate fracture orientations consistent with shear-slip reactivation of late-Palaeozoic fractures over a widespread area, as well as an additional fracture orientation present only near the injection well. We propose simple Mohr–Coulomb fracture models to explain these observations. The observation that induced seismicity preferentially occurs along one of the identified conjugate fracture orientations can be explained by a rotation in the direction of the regional maximum compressive stress from the time when the fractures were formed to the present. Shear slip along the third fracture orientation observed near the injection well is inconsistent with the current regional stress field and suggests a local rotation of the horizontal stresses. The detailed subsurface model produced by this analysis provides important insights for anticipating spatial patterns of future induced seismicity and for evaluation of possible additional injection well sites that are likely to be seismically and hydrologically isolated from the current well. In addition, the interpreted fault patterns provide constraints for estimating the maximum magnitude earthquake that may be induced, and for building geomechanical models to simulate pore pressure diffusion, stress changes and earthquake triggering.