Estimation of thickness and velocity changes of injected carbon dioxide layers from prestack time-lapse seismic data

Abstract

In this study, we bring together the two main categories of time-lapse seismic analysis — amplitude analysis and time-shift analysis — to estimate simultaneously the changes in thickness and velocity of a 4D seismic anomaly. The methodology is applied to time-lapse seismic monitoring of carbon dioxide [Formula: see text] storage at Sleipner field, Norway, that shows significant 4D effects. The 4D anomalies resulting from [Formula: see text] injection appear as a multilayer reflection pattern within the relatively shallow Utsira Sand. This multireflective appearance within the sand layer is interpreted as [Formula: see text] layers trapped below thin shale layers. Because most of the [Formula: see text] layers are believed to be thin [Formula: see text], the interference between top and base of these layers needs to be taken into account in 4D seismic analysis. By studying the reflected event from a horizon below the Utsira Sand, we estimate 4D traveltime shifts caused by the presence of the [Formula: see text] layer above thishorizon. We then combine these traveltime shifts with measured amplitude changes for the top and base of the [Formula: see text] layer to estimate velocity and thickness changes for the thin [Formula: see text] layer. In 1999, after three years of injection, the most likely velocity change was around [Formula: see text] and the thickness of the [Formula: see text] layer was around [Formula: see text]. In 2001, the corresponding velocity change and thickness estimates were [Formula: see text] and [Formula: see text], respectively. Finally, in 2002, the most likely velocity change was [Formula: see text] and the thickness of the [Formula: see text] layer was [Formula: see text]. It is not straightforward to apply this method to a stack of [Formula: see text] layers because 4D time shifts below the Utsira Sand only provide information about the average time shift for all layers. The amplitude information for each individual [Formula: see text] layer cannot be resolved without knowing the velocity change within each layer. However, our result from a single [Formula: see text] layer may be used to constrain the velocity changes for the multilayer [Formula: see text] case.