3-D mechanical analysis of complex reservoirs: a novel mesh-free approach

Building geomechanical models for induced seismicity in complex reservoirs poses a major challenge, in particular if many faults need to be included. We developed a novel way of calculating induced stress changes and associated seismic moment response for structurally complex reservoirs with tens to hundreds of faults. Our specific target was to improve the predictive capability of stress evolution along multiple faults, and to use the calculations to enhance physics-based understanding of the reservoir seismicity. Our methodology deploys a mesh-free numerical and analytical approach for both the stress calculation and the seismic moment calculation. We introduce a high-performance computational method for high-resolution induced Coulomb stress changes along faults, based on a Green's function for the stress response to a nucleus of strain. One key ingredient is the deployment of an octree representation and calculation scheme for the nuclei of strain, based on the topology and spatial variability of the mesh of the reservoir flow model. Once the induced stress changes are evaluated along multiple faults, we calculate potential seismic moment release in a fault system supposing an initial stress field. The capability of the approach, dubbed as MACRIS (Mechanical Analysis of Complex Reservoirs for Induced Seismicity) is proven through comparisons with finite element models. Computational performance and suitability for probabilistic assessment of seismic hazards are demonstrated though the use of the complex, heavily faulted Gullfaks field.