Mesh-free stress solution for complex 3D reservoir simulation grids: a mixed nuclei of strain and analytical element approach

Summary An innovative method is presented for full 3D pressure and temperature dominated stress evaluations in the subsurface using tetrahedra-based analytical elements combined with the inflation point source solution. A mesh-free approach suitable to industry standard 3D flow simulation models (based on hexahedral cells) is obtained by representing the grid cells in tetrahedral elements, effectively preserving the 3D geometrical complexities of the reservoir. Contributions from neighboring grid cells are added via a Tartan grid representation of point sources that enables the spatial resolution to increase and the stress evaluations to be carried out in parallel, both of which significantly improve computational efficiency. The novel approach is demonstrated for synthetic low enthalpy geothermal models with clastic reservoir characteristics and varying degrees of geometrical and structural complexity. The method is shown able to accurately capture the effects of stress arching on complex faults causing reservoir throw and flow compartmentalization, and along the rim of the cold-water volume. Results of a synthetic geothermal development model of the heavily faulted Gullfaks field show the novel method to provide an accurate and computationally highly efficient approach for evaluating pressure and temperature dominated stress changes in structurally complex sedimentary reservoirs.