Effects of cavern depth on surface subsidence and storage loss of oil-filled caverns

Abstract

Finite element analyses of oil-filled caverns were performed to investigate the effects of cavern depth on surface subsidence and storage loss, a primary performance criteria of SPR caverns. The finite element model used for this study was axisymmetric, approximating an infinite array of caverns spaced at 750 ft. The stratigraphy and cavern size were held constant while the cavern depth was varied between 1500 ft and 3000 ft in 500 ft increments. Thirty year simulations, the design life of the typical SPR cavern, were performed with boundary conditions modeling the oil pressure head applied to the cavern lining. A depth dependent temperature gradient of 0.012{degrees}F/ft was also applied to the model. The calculations were performed using ABAQUS, a general purpose of finite element analysis code. The user-defined subroutine option in ABAQUS was used to enter an elastic secondary creep model which includes temperature dependence. The calculations demonstrated that surface subsidence and storage loss rates increase with increasing depth. At lower depths the difference between the lithostatic stress and the oil pressure is greater. Thus, the effective stresses are greater, resulting in higher creep rates. Furthermore, at greater depths the cavern temperatures are higher which also produce higher creep rates. Together, these factors result in faster closure of the cavern. At the end of the 30 year simulations, a 1500 ft-deep cavern exhibited 4 percent storage loss and 4 ft of subsidence while a 3000 ft-deep cavern exhibited 33 percent storage loss and 44 ft of subsidence. The calculations also demonstrated that surface subsidence is directly related to the amount of storage loss. Deeper caverns exhibit more subsidence because the caverns exhibit more storage loss. However, for a given amount of storage loss, nearly the same magnitude of surface subsidence was exhibited, independent of cavern depth.