Burying receivers for improved time‐lapse seismic repeatability: CO2CRC Otway field experiment

4D seismic is widely used to remotely monitor fluid movement in subsurface reservoirs. This technique is especially effective offshore where high survey repeatability can be achieved. It comes as no surprise that the first 4D seismic that successfully monitored the CO₂ sequestration process was recorded offshore in the Sleipner field, North Sea. In the case of land projects, poor repeatability of the land seismic data due to low S/N ratio often obscures the time‐lapse seismic signal. Hence for a successful on shore monitoring program improving seismic repeatability is essential. Stage 2 of the CO2CRC Otway project involves an injection of a small amount (around 15,000 tonnes) of CO₂/CH₄ gas mixture into a saline aquifer at a depth of approximately 1.5 km. Previous studies at this site showed that seismic repeatability is relatively low due to variations in weather conditions, near surface geology and farming activities. In order to improve time‐lapse seismic monitoring capabilities, a permanent receiver array can be utilised to improve signal to noise ratio and hence repeatability. A small‐scale trial of such an array was conducted at the Otway site in June 2012. A set of 25 geophones was installed in 3 m deep boreholes in parallel to the same number of surface geophones. In addition, four geophones were placed into boreholes of 1–12 m depth. In order to assess the gain in the signal‐to‐noise ratio and repeatability, both active and passive seismic surveys were carried out. The surveys were conducted in relatively poor weather conditions, with rain, strong wind and thunderstorms. With such an amplified background noise level, we found that the noise level for buried geophones is on average 20 dB lower compared to the surface geophones. The levels of repeatability for borehole geophones estimated around direct wave, reflected wave and ground roll are twice as high as for the surface geophones. Both borehole and surface geophones produce the best repeatability in the 30–90 Hz frequency range. The influence of burying depth on S/N ratio and repeatability shows that significant improvement in repeatability can be reached at a depth of 3 m. The level of repeatability remains relatively constant between 3 and 12 m depths.     

Post‐stack stratigraphic inversion workflow applied to carbon dioxide storage: application to the saline aquifer of Sleipner field‡

In the Norwegian North Sea, the Sleipner field produces gas with a high CO₂ content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO₂) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO₂ plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO₂ over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO₂ plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO₂ injection and the 2006 vintage carried out after a CO₂ injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO₂ in a reservoir. With the post‐stack inversion workflow applied to CO₂ storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.     

Full waveform inversion in the time lapse mode applied to CO2 storage at Sleipner

Carbon capture and storage is a viable greenhouse gas mitigation technology and the Sleipner CO₂ sequestration site in the North Sea is an excellent example. Storage of CO₂ at the Sleipner site requires monitoring over large areas, which can successfully be accomplished with time lapse seismic imaging. One of the main goals of CO₂ storage monitoring is to be able to estimate the volume of the stored CO₂ in the reservoir. This requires a parametrization of the subsurface as exact as possible. Here we use elastic 2D time‐domain full waveform inversion in a time lapse manner to obtain a P‐wave velocity constrain directly in the depth domain for a base line survey in 1994 and two post‐injection surveys in 1999 and 2006. By relating velocity change to free CO₂ saturation, using a rock physics model, we find that at the considered location the aquifer may have been fully saturated in some places in 1999 and 2006.     

Dual sensor streamer technology used in Sleipner CO2 injection monitoring

CO₂ has been injected into the saline aquifer Utsira Fm at the Sleipner field since 1996. In order to monitor the movement of the CO₂ in the sub‐surface, the seventh seismic monitor survey was acquired in 2010, with dual sensor streamers which enabled optimal towing depths compared to previous surveys. We here report both on the time‐lapse observations and on the improved resolution compared to the conventional streamer surveys. This study shows that the CO₂ is still contained in the subsurface, with no indications of leakage. The time‐lapse repeatability of the dual sensor streamer data versus conventional data is sufficient for interpreting the time‐lapse effects of the CO₂ at Sleipner, and the higher resolution of the 2010 survey has enabled a refinement of the interpretation of nine CO₂ saturated layers with improved thickness estimates of the layers. In particular we have estimated the thickness of the uppermost CO₂ layer based on an analysis of amplitude strength together with time‐separation of top and base of this layer and found the maximum thickness to be 11 m. This refined interpretation gives a good base line for future time‐lapse surveys at the Sleipner CO₂ injection site.     

Application of diffracted wave analysis to time‐lapse seismic data for CO2 leakage detection

Time‐lapse seismic analysis is utilized in CO₂ geosequestration to verify the CO₂ containment within a reservoir. A major risk associated with geosequestration is a possible leakage of CO₂ from the storage formation into overlaying formations. To mitigate this risk, the deployment of carbon capture and storage projects requires fast and reliable detection of relatively small volumes of CO₂ outside the storage formation. To do this, it is necessary to predict typical seepage scenarios and improve subsurface seepage detection methods. In this work we present a technique for CO₂ monitoring based on the detection of diffracted waves in time‐lapse seismic data. In the case of CO₂ seepage, the migrating plume might form small secondary accumulations that would produce diffracted, rather than reflected waves. From time‐lapse data analysis, we are able to separate the diffracted waves from the predominant reflections in order to image the small CO₂ plumes. To explore possibilities to detect relatively small amounts of CO₂, we performed synthetic time‐lapse seismic modelling based on the Cooperative Research Centre for Greenhouse Gas Technologies (CO2CRC) Otway project data. The detection method is based on defining the CO₂ location by measuring the coherency of the signal along diffraction offset‐traveltime curves. The technique is applied to a time‐lapse stacked section using a stacking velocity to construct offset‐traveltime curves. Given the amount of noise found in the surface seismic data, the predicted minimum detectable amount of CO₂ is 1000–2000 tonnes. This method was also applied to real data obtained from a time‐lapse seismic physical model. The use of diffractions rather than reflections for monitoring small amounts of CO₂ can enhance the capability of subsurface monitoring in CO₂ geosequestration projects.     

Rough seas and time-lapse seismic

Time‐lapse seismic surveying has become an accepted tool for reservoir monitoring applications, thus placing a high premium on data repeatability. One factor affecting data repeatability is the influence of the rough sea‐surface on the ghost reflection and the resulting seismic wavelets of the sources and receivers. During data analysis, the sea‐surface is normally assumed to be stationary and, indeed, to be flat. The non‐flatness of the sea‐surface introduces amplitude and phase perturbations to the source and receiver responses and these can affect the time‐lapse image. We simulated the influence of rough sea‐surfaces on seismic data acquisition. For a typical seismic line with a 48‐fold stack, a 2‐m significant‐wave‐height sea introduces RMS errors of about 5–10% into the stacked data. This level of error is probably not important for structural imaging but could be significant for time‐lapse surveying when the expected difference anomaly is small. The errors are distributed differently for sources and receivers because of the different ways they are towed. Furthermore, the source wavelet is determined by the sea shape at the moment the shot is fired, whereas the receiver wavelet is time‐varying because the sea moves significantly during the seismic record.     

Time‐lapse seismic modelling of CO2 fluid substitution in the Devonian Redwater Reef,Alberta, Canada

Common shot ray tracing and finite difference seismic modelling experiments were undertaken to evaluate variations in the seismic response of the Devonian Redwater reef in the Alberta Basin, Canada after replacement of native pore waters in the upper rim of the reef with CO₂. This part of the reef is being evaluated for a CO₂ storage project. The input geological model was based on well data and the interpretation of depth‐converted, reprocessed 2D seismic data in the area. Pre‐stack depth migration of the ray traced and finite difference synthetic data demonstrate similar seismic attributes for the Mannville, Nisku, Ireton, Cooking Lake, and Beaverhill Lake formations and clear terminations of the Upper Leduc and Middle Leduc events at the reef margin. Higher amplitudes at the base of Upper‐Leduc member are evident near the reef margin due to the higher porosity of the foreslope facies in the reef rim compared to the tidal flat lagoonal facies within the central region of the reef. Time‐lapse seismic analysis exhibits an amplitude difference of about 14% for Leduc reflections before and after CO₂ saturation and a travel‐time delay through the reservoir of 1.6 ms. Both the ray tracing and finite difference approaches yielded similar results but, for this particular model, the latter provided more precise imaging of the reef margin. From the numerical study we conclude that time‐lapse surface seismic surveys should be effective in monitoring the location of the CO₂ plume in the Upper Leduc Formation of the Redwater reef, although the differences in the results between the two modelling approaches are of similar order to the effects of the CO₂ fluid replacement itself.     

Monitoring CO2 saturation using time‐lapse amplitude versus offset analysis of 3D seismic data from the Ketzin CO2 storage pilot site,Germany

The injection of CO₂ at the Ketzin pilot site commenced in June 2008 and was terminated in August 2013 after 67 kT had been injected into a saline formation at a depth of 630–650 m. As part of the site monitoring program, four 3D surface seismic surveys have been acquired to date, one baseline and three repeats, of which two were conducted during the injection period, and one during the post‐injection phase. The surveys have provided the most comprehensive images of the spreading CO₂ plume within the reservoir layer. Both petrophysical experiments on core samples from the Ketzin reservoir and spectral decomposition of the 3D time‐lapse seismic data show that the reservoir pore pressure change due to CO₂ injection has a rather minor impact on the seismic amplitudes. Therefore, the observed amplitude anomaly is interpreted to be mainly due to CO₂ saturation. In this study, amplitude versus offset analysis has been applied to investigate the amplitude versus offset response from the top of the sandstone reservoir during the injection and post‐injection phases, and utilize it to obtain a more quantitative assessment of the CO₂ gaseous saturation changes. Based on the amplitude versus offset modelling, a prominent decrease in the intercept values imaged at the top of the reservoir around the injection well is indeed associated solely with the CO₂ saturation increase. Any change in the gradient values, which would, in case it was positive, be the only signature induced by the reservoir pressure variations, has not been observed. The amplitude versus offset intercept change is, therefore, entirely ascribed to CO₂ saturation and used for its quantitative assessment. The estimated CO₂ saturation values around the injection area in the range of 40%–60% are similar to those obtained earlier from pulsed neutron‐gamma logging. The highest values of 80% are found in the second seismic repeat in close vicinity to the injection and observation wells.     

Seismic time‐lapse effects of solution salt mining – a feasibility study

This article addresses the question whether time‐lapse seismic reflection techniques can be used to follow and quantify the effects of solution salt mining. Specifically, the production of magnesium salts as mined in the north of the Netherlands is considered. The use of seismic time‐lapse techniques to follow such a production has not previously been investigated. For hydrocarbon production and CO₂ storage, time‐lapse seismics are used to look at reservoir changes mainly caused by pressure and saturation changes in large reservoirs, while for solution mining salt is produced from caverns with a limited lateral extent, with much smaller production volumes and a fluid (brine) replacing a solid (magnesium salt). In our approach we start from the present situation of the mine and then study three different production scenarios, representing salt production both in vertical and lateral directions of the mine. The present situation and future scenarios have been transformed into subsurface models that were input to an elastic finite‐difference scheme to create synthetic seismic data. These data have been analysed and processed up to migrated seismic images, such that time‐lapse analyses of intermediate and final results could be done. From the analyses, it is found that both vertical and lateral production is visible well above the detection threshold in difference data, both at pre‐imaging and post‐imaging stages. In quantitative terms, an additional production of the mine of 6 m causes time‐shifts in the order of 2 ms (pre‐imaging) and 4 ms (post‐imaging) and amplitude changes of above 20% in the imaged sections. A laterally oriented production causes even larger amplitude changes at the edge of the cavern due to replacement of solid magnesium salt with brine introducing a large seismic contrast. Overall, our pre‐imaging and post‐imaging time‐lapse analysis indicates that the effects of solution salt mining can be observed and quantified on seismic data. The effects seem large enough to be observable in real seismic data containing noise.     

Overburden complexity and repeatability of seismic data: Impacts of positioning errors at the Oseberg field, North Sea

Repeatability of seismic data plays a crucial role in time‐lapse seismic analysis. There are several factors that can decrease the repeatability, such as positioning errors, varying tide, source variations, velocity changes in the water layer (marine data) and undesired effects of various processing steps. In this work, the complexity of overburden structure, as an inherent parameter that can affect the repeatability, is studied. A multi‐azimuth three‐dimensional vertical‐seismic‐profiling data set with 10 000 shots is used to study the relationship between overburden structure and repeatability of seismic data. In most repeatability studies, two data sets are compared, but here a single data set has been used because a significant proportion of the 10 000 shots are so close to each other that a repeatability versus positioning error is possible. We find that the repeatability decreases by a factor of approximately 2 under an overburden lens. Furthermore, we find that the X‐ and Y‐components have approximately the same sensitivity to positioning errors as the Z‐component (for the same events) in this three‐dimensional vertical‐seismic‐profiling experiment. This indicates that in an area with complex overburden, positioning errors between monitor and base seismic surveys are significantly more critical than outside such an area. This study is based on a three‐dimensional three‐component vertical‐seismic‐profiling data set from a North Sea reservoir and care should be taken when extrapolating these observations into a general four‐dimensional framework.     

CO2 storage in the high Arctic: efficient modelling of pre‐stack depth‐migrated seismic sections for survey planning

The sequestration of CO₂ in subsurface reservoirs constitutes an immediate counter‐measure to reduce anthropogenic emissions of CO₂, now recognized by international scientific panels to be the single most critical factor driving the observed global climatic warming. To ensure and verify the safe geological containment of CO₂ underground, monitoring of the CO₂ site is critical. In the high Arctic, environmental considerations are paramount and human impact through, for instance, active seismic surveys, has to be minimized. Efficient seismic modelling is a powerful tool to test the detectability and imaging capability prior to acquisition and thus improve the characterization of CO₂ storage sites, taking both geological setting and seismic acquisition set‐up into account. The unique method presented here avoids the costly generation of large synthetic data sets by employing point spread functions to directly generate pre‐stack depth‐migrated seismic images. We test both a local‐target approach using an analytical filter assuming an average velocity and a full‐field approach accounting for the spatial variability of point spread functions. We assume a hypothetical CO₂ plume emplaced in a sloping aquifer inspired by the conditions found at the University of Svalbard CO₂ lab close to Longyearbyen, Svalbard, Norway, constituting an unconventional reservoir–cap rock system. Using the local‐target approach, we find that even the low‐to‐moderate values of porosity (5%–18%) measured in the reservoir should be sufficient to induce significant change in seismic response when CO₂ is injected. The sensitivity of the seismic response to changes in CO₂ saturation, however, is limited once a relatively low saturation threshold of 5% is exceeded. Depending on the illumination angle provided by the seismic survey, the quality of the images of five hypothetical CO₂ plumes of varying volume differs depending on the steepness of their flanks. When comparing the resolution of two orthogonal 2D surveys to a 3D survey, we discover that the images of the 2D surveys contain significant artefacts, the CO₂‐brine contact is misplaced and an additional reflector is introduced due to the projection of the point spread function of the unresolvable plane onto the imaging plane. All of these could easily lead to a misinterpretation of the behaviour of the injected CO₂. Our workflow allows for testing the influence of geological heterogeneities in the target aquifer (igneous intrusions, faults, pervasive fracture networks) by utilizing increasingly complex and more realistic geological models as input as more information on the subsurface becomes available.     

The first post‐injection seismic monitor survey at the Ketzin pilot CO2 storage site: results from time‐lapse analysis

The injection of CO₂ at the Ketzin pilot CO₂ storage site started in June 2008 and ended in August 2013. During the 62 months of injection, a total amount of about 67 kt of CO₂ was injected into a saline aquifer. A third repeat three‐dimensional seismic survey, serving as the first post‐injection survey, was acquired in 2015, aiming to investigate the recent movement of the injected CO₂. Consistent with the previous two time‐lapse surveys, a predominantly west–northwest migration of the gaseous CO₂ plume in the up‐dip direction within the reservoir is inferred in this first post‐injection survey. No systematic anomalies are detected through the reservoir overburden. The extent of the CO₂ plume west of the injection site is almost identical to that found in the 2012 second repeat survey (after injection of 61 kt); however, there is a significant decrease in its size east of the injection site. Assessment of the CO₂ plume distribution suggests that the decrease in the size of the anomaly may be due to multiple factors, such as limited vertical resolution, CO₂ dissolution, and CO₂ migration into thin layers, in addition to the effects of ambient noise. Four‐dimensional seismic modelling based on dynamic flow simulations indicates that a dynamic balance between the newly injected CO₂ after the second repeat survey and the CO₂ migrating into thin layers and being dissolved was reached by the time of the first post‐injection survey. In view of the significant uncertainties in CO₂ mass estimation, both patchy and non‐patchy saturation models for the Ketzin site were taken into consideration.     

Joint cross-well and single-well seismic studies of CO2 injection in an oil reservoir

A series of time‐lapse seismic cross‐well and single‐well experiments were conducted in a diatomite reservoir to monitor the injection of CO₂ into a hydrofracture zone, based on P‐ and S‐wave data. A high‐frequency piezo‐electric P‐wave source and an orbital‐vibrator S‐wave source were used to generate waves that were recorded by hydrophones as well as 3‐component geophones. During the first phase the set of seismic experiments was conducted after the injection of water into the hydrofractured zone. The set of seismic experiments was repeated after a time period of seven months during which CO₂ was injected into the hydrofractured zone. The questions to be answered ranged from the detectability of the geological structure in the diatomic reservoir to the detectability of CO₂ within the hydrofracture. Furthermore, it was intended to determine which experiment (cross‐well or single‐well) is best suited to resolve these features. During the pre‐injection experiment, the P‐wave velocities exhibited relatively low values between 1700 and 1900 m/s, which decreased to 1600–1800 m/s during the post‐injection phase (?5%). The analysis of the pre‐injection S‐wave data revealed slow S‐wave velocities between 600 and 800 m/s, while the post‐injection data revealed velocities between 500 and 700 m/s (?6%). These velocity estimates produced high Poisson's ratios between 0.36 and 0.46 for this highly porous (～50%) material. Differencing post‐ and pre‐injection data revealed an increase in Poisson's ratio of up to 5%. Both velocity and Poisson's ratio estimates indicate the dissolution of CO₂ in the liquid phase of the reservoir accompanied by an increase in pore pressure. The single‐well data supported the findings of the cross‐well experiments. P‐ and S‐wave velocities as well as Poisson's ratios were comparable to the estimates of the cross‐well data. The cross‐well experiment did not detect the presence of the hydrofracture but appeared to be sensitive to overall changes in the reservoir and possibly the presence of a fault. In contrast, the single‐well reflection data revealed an arrival that could indicate the presence of the hydrofracture between the source and receiver wells, while it did not detect the presence of the fault, possibly due to out‐of‐plane reflections.     

Seismic time‐lapse imaging using interferometric least‐squares migration: Case study

Seismic time‐lapse surveys are susceptible to repeatability errors due to varying environmental conditions. To mitigate this problem, we propose the use of interferometric least‐squares migration to estimate the migration images for the baseline and monitor surveys. Here, a known reflector is used as the reference reflector for interferometric least‐squares migration, and the data are approximately redatumed to this reference reflector before imaging. This virtual redatuming mitigates the repeatability errors in the time‐lapse migration image. Results with synthetic and field data show that interferometric least‐squares migration can sometimes reduce or eliminate artifacts caused by non‐repeatability in time‐lapse surveys and provide a high‐resolution estimate of the time‐lapse change in the reservoir.     

When 4D seismic is not applicable: Alternative monitoring scenarios for the Arab-D reservoir in the Ghawar Field

Ghawar, the largest oilfield in the world, produces oil from the Upper Jurassic Arab‐D carbonate reservoir. The high rigidity of the limestone–dolomite reservoir rock matrix and the small contrast between the elastic properties of the pore fluids, i.e. oil and water, are responsible for the weak 4D seismic effect due to oil production. A feasibility study was recently completed to quantify the 4D seismic response of reservoir saturation changes as brine replaced oil. The study consisted of analysing reservoir rock physics, petro‐acoustic data and seismic modelling. A seismic model of flow simulation using fluid substitution concluded that time‐lapse surface seismic or conventional 4D seismic is unlikely to detect the floodfront within the repeatability of surface seismic measurements. Thus, an alternative approach to 4D seismic for reservoir fluid monitoring is proposed. Permanent seismic sensors could be installed in a borehole and on the surface for passive monitoring of microseismic activity from reservoir pore‐pressure perturbations. Reservoir production and injection operations create these pressure or stress perturbations. Reservoir heterogeneities affecting the fluid flow could be mapped by recording the distribution of epicentre locations of these microseisms or small earthquakes. The permanent borehole sensors could also record repeated offset vertical seismic profiling surveys using a surface source at a fixed location to ensure repeatability. The repeated vertical seismic profiling could image the change in reservoir properties with production.     

Quantitative detection of fluid distribution using time-lapse seismic

Although previous seismic monitoring studies have revealed several relationships between seismic responses and changes in reservoir rock properties, the quantitative evaluation of time‐lapse seismic data remains a challenge. In most cases of time‐lapse seismic analysis, fluid and/or pressure changes are detected qualitatively by changes in amplitude strength, traveltime and/or Poisson's ratio. We present the steps for time‐lapse seismic analysis, considering the pressure effect and the saturation scale of fluids. We then demonstrate a deterministic workflow for computing the fluid saturation in a reservoir in order to evaluate time‐lapse seismic data. In this approach, we derive the physical properties of the water‐saturated sandstone reservoir, based on the following inputs: V_P, V_S, ρ and the shale volume from seismic analysis, the average properties of sand grains, and formation‐water properties. Next, by comparing the in‐situ fluid‐saturated properties with the 100% formation‐water‐saturated reservoir properties, we determine the bulk modulus and density of the in‐situ fluid. Solving three simultaneous equations (relating the saturations of water, oil and gas in terms of the bulk modulus, density and the total saturation), we compute the saturation of each fluid. We use a real time‐lapse seismic data set from an oilfield in the North Sea for a case study.     

Application of the continuous wavelet transform on seismic data for mapping of channel deposits and gas detection at the CO2SINK site, Ketzin, Germany

Conventional seismic data are band limited and therefore, provide limited geological information. Every method that can push the limits is desirable for seismic data analysis. Recently, time‐frequency decomposition methods are being used to quickly extract geological information from seismic data and, especially, for revealing frequency‐dependent amplitude anomalies. Higher frequency resolution at lower frequencies and higher temporal resolution at higher frequencies are the objectives for different time‐frequency decomposition methods. Continuous wavelet transform techniques, which are the same as narrow‐band spectral analysis methods, provide frequency spectra with high temporal resolution without the windowing process associated with other techniques. Therefore, this technique can be used for analysing geological information associated with low and high frequencies that normally cannot be observed in conventional seismic data. In particular, the continuous wavelet transform is being used to detect thin sand bodies and also as a direct hydrocarbon indicator. This paper presents an application of the continuous wavelet transform method for the mapping of potential channel deposits, as well as remnant natural gas detection by mapping low‐frequency anomalies associated with the gas. The study was carried out at the experimental CO₂ storage site at Ketzin, Germany (CO₂SINK). Given that reservoir heterogeneity and faulting will have significant impact on the movement and storage of the injected CO₂, our results are encouraging for monitoring the migration of CO₂ at the site. Our study confirms the efficiency of the continuous wavelet transform decomposition method for the detection of frequency‐dependent anomalies that may be due to gas migration during and after the injection phase and in this way, it can be used for real‐time monitoring of the injected CO₂ from both surface and borehole seismics.     

Utilizing spectral decomposition to determine the distribution of injected CO2 at the Snøhvit Field

Time‐lapse 3D seismic reflection data, covering the CO₂ storage operation at the Snøhvit gas field in the Barents Sea, show clear amplitude and time‐delay differences following injection. The nature and extent of these changes suggest that increased pore fluid pressure contributes to the observed seismic response, in addition to a saturation effect. Spectral decomposition using the smoothed pseudo‐Wigner–Ville distribution has been used to derive discrete‐frequency reflection amplitudes from around the base of the CO₂ storage reservoir. These are utilized to determine the lateral variation in peak tuning frequency across the seismic anomaly as this provides a direct proxy for the thickness of the causative feature. Under the assumption that the lateral and vertical extents of the respective saturation and pressure changes following CO₂ injection will be significantly different, discrete spectral amplitudes are used to distinguish between the two effects. A clear spatial separation is observed in the distribution of low‐ and high‐frequency tuning. This is used to discriminate between direct fluid substitution of CO₂, as a thin layer, and pressure changes that are distributed across a greater thickness of the storage reservoir. The results reveal a striking correlation with findings derived from pressure and saturation discrimination algorithms based on amplitude versus offset analysis.     

Analysis of time‐lapse travel‐time and amplitude changes to assess reservoir compartmentalization

Fluid depletion within a compacting reservoir can lead to significant stress and strain changes and potentially severe geomechanical issues, both inside and outside the reservoir. We extend previous research of time‐lapse seismic interpretation by incorporating synthetic near‐offset and full‐offset common‐midpoint reflection data using anisotropic ray tracing to investigate uncertainties in time‐lapse seismic observations. The time‐lapse seismic simulations use dynamic elasticity models built from hydro‐geomechanical simulation output and a stress‐dependent rock physics model. The reservoir model is a conceptual two‐fault graben reservoir, where we allow the fault fluid‐flow transmissibility to vary from high to low to simulate non‐compartmentalized and compartmentalized reservoirs, respectively. The results indicate time‐lapse seismic amplitude changes and travel‐time shifts can be used to qualitatively identify reservoir compartmentalization. Due to the high repeatability and good quality of the time‐lapse synthetic dataset, the estimated travel‐time shifts and amplitude changes for near‐offset data match the true model subsurface changes with minimal errors. A 1D velocity–strain relation was used to estimate the vertical velocity change for the reservoir bottom interface by applying zero‐offset time shifts from both the near‐offset and full‐offset measurements. For near‐offset data, the estimated P‐wave velocity changes were within 10% of the true value. However, for full‐offset data, time‐lapse attributes are quantitatively reliable using standard time‐lapse seismic methods when an updated velocity model is used rather than the baseline model.     

Calendar time interpolation of amplitude maps from 4D seismic data

A calendar time interpolation method for 2D seismic amplitude maps, done in two steps, is presented. The contour interpolation part is formulated as a quadratic programming problem, whereas the amplitude value interpolation is based on a conditional probability formulation. The method is applied on field data from the Sleipner CO₂ storage project. The output is a continuous image (movie) of the CO₂ plume. Besides visualization, the output can be used to better couple 4D seismic to other types of data acquired. The interpolation uncertainty increases with the time gap between consecutive seismic surveys and is estimated by leaving a survey out (blind test). Errors from such tests can be used to identify problems in understanding the flow and possibly improve the interpolation scheme for a given case. Field‐life cost of various acquisition systems and repeat frequencies are linked to the time‐lapse interpolation errors. The error in interpolated amplitudes increased by 3%‐4% per year of interpolation gap for the Sleipner case. Interpolation can never fully replace measurements.