Premium
An inverse method for estimating thickness and volume with time of a thin CO 2 ‐filled layer at the Sleipner Field, North Sea
Author(s) -
Cowton L. R.,
Neufeld J. A.,
White N. J.,
Bickle M. J.,
White J. C.,
Chadwick R. A.
Publication year - 2016
Publication title -
journal of geophysical research: solid earth
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.983
H-Index - 232
eISSN - 2169-9356
pISSN - 2169-9313
DOI - 10.1002/2016jb012895
Subject(s) - geology , inverse , reflection (computer programming) , layer (electronics) , mineralogy , amplitude , broadband , optics , materials science , composite material , geometry , computer science , physics , mathematics , programming language
Migration of CO 2 through storage reservoirs can be monitored using time lapse seismic reflection surveys. At the Sleipner Field, injected CO 2 is distributed throughout nine layers within the reservoir. These layers are too thin to be seismically resolvable by direct measurement of the separation between reflections from the top and bottom of each layer. Here we develop and apply an inverse method for measuring thickness changes of the shallowest layer. Our approach combines differences in traveltime down to a specific reflection together with amplitude measurements to determine layer thicknesses from time lapse surveys. A series of synthetic forward models were used to test the robustness of our inverse approach and to quantify uncertainties. In the absence of ambient noise, this approach can unambiguously resolve layer thickness. If a realistic ambient noise distribution is included, layer thicknesses of 1–6 m are accurately retrieved with an uncertainty of ±0.5 m. We used this approach to generate a thickness map of the shallowest layer. The fidelity of this result was tested using measurements of layer thickness determined from the 2010 broadband seismic survey. The calculated volume of CO 2 within the shallowest layer increases at a rate that is quadratic in time, despite an approximately constant injection rate into the base of the reservoir. This result is consistent with a diminished growth rate of the areal extent of underlying layers. Finally, the relationship between caprock topography and layer thickness is explored and potential migration pathways that charge this layer are identified.