Wavelet‐based cepstrum decomposition of seismic data and its application in hydrocarbon detection

How to use cepstrum analysis for reservoir characterization and hydrocarbon detection is an initial question of great interest to exploration seismologists. In this paper, wavelet‐based cepstrum decomposition is proposed as a valid technology for enhancing geophysical responses in specific frequency bands, in the same way as traditional spectrum decomposition methods do. The calculation of wavelet‐based cepstrum decomposition, which decomposes the original seismic volume into a series of common quefrency volumes, employs a sliding window to move over each seismic trace sample by sample. The key factor in wavelet‐based cepstrum decomposition is the selection of the sliding‐window length as it limits the frequency ranges of the common quefrency section. Comparison of the wavelet‐based cepstrum decomposition with traditional spectrum decomposition methods, such as short‐time Fourier transform and wavelet transform, is conducted to demonstrate the effectiveness of the wavelet‐based cepstrum decomposition and the relation between these two technologies. In hydrocarbon detection, seismic amplitude anomalies are detected using wavelet‐based cepstrum decomposition by utilizing the first and second common quefrency sections. This reduces the burden of needing dozens of seismic volumes to represent the response to different mono‐frequency sections in the interpretation of spectrum decomposition in conventional spectrum decomposition methods. The model test and the application of real data acquired from the Sulige gas field in the Ordos Basin, China, confirm the effectiveness of the seismic amplitude anomaly section using wavelet‐based cepstrum decomposition for discerning the strong amplitude anomalies at a particular quefrency buried in the broadband seismic response. Wavelet‐based cepstrum decomposition provides a new method for measuring the instantaneous cepstrum properties of a reservoir and offers a new field of processing and interpretation of seismic reflection data.     

Feasibility of utilizing wavelet phase to map the CO2 plume at the Ketzin pilot site,Germany

Spectral decomposition is a powerful tool that can provide geological details dependent upon discrete frequencies. Complex spectral decomposition using inversion strategies differs from conventional spectral decomposition methods in that it produces not only frequency information but also wavelet phase information. This method was applied to a time‐lapse three‐dimensional seismic dataset in order to test the feasibility of using wavelet phase changes to detect and map injected carbon dioxide within the reservoir at the Ketzin carbon dioxide storage site, Germany. Simplified zero‐offset forward modelling was used to help verify the effectiveness of this technique and to better understand the wavelet phase response from the highly heterogeneous storage reservoir and carbon dioxide plume. Ambient noise and signal‐to‐noise ratios were calculated from the raw data to determine the extracted wavelet phase. Strong noise caused by rainfall and the assumed spatial distribution of sandstone channels in the reservoir could be correlated with phase anomalies. Qualitative and quantitative results indicate that the wavelet phase extracted by the complex spectral decomposition technique has great potential as a practical and feasible tool for carbon dioxide detection at the Ketzin pilot site.     

Possible thermal brightness temperature anomalies associated with the Lushan (China) MS7.0 earthquake on 20 April 2013

Lushan M_S7.0 earthquake occurred in Lushan county, Ya'an city, Sichuan province of China, on 20 April 2013, and caused 196 deaths, 23 people of missing and more than 12 thousand of people injured. In order to analyze the possible seismic brightness temperature anomalies which may be associated with Lushan earthquake, daily brightness temperature data for the period from 1 June 2011 to 31 May 2013 and the geographical extent of 25°E-35°N latitude and 98°E-108°E longitude are collected from Chinese geostationary meteorological satellite FY-2E. Continuous wavelet transform method which has good resolution both in time and frequency domains is used to analyze power spectrum of brightness temperature data. The results show that the relative wavelet power spectrum (RWPS) anomalies appeared since 24 January 2013 and still lasted on 19 April. Anomalies firstly appeared at the middle part of Longmenshan fault zone and gradually spread toward the southwestern part of Longmenshan fault. Anomalies also appeared along the Xianshuihe fault since about 1 March. Eventually, anomalies gathered at the intersection zone of Longmenshan and Xianshuihe faults. The anomalous areas and RWPS amplitude increased since the appearance of anomalies and reached maximum in late March. Anomalies attenuated with the earthquake approaching. And eventually the earthquake occurred at the southeastern edge of anomalous areas. Lushan earthquake was the only obvious geological event within the anomalous area during the time period, so the anomalous changes of RWPS are possibly associated with the earthquake.     

Seismic dynamic monitoring in CO<Subscript>2</Subscript> flooding based on characterization of frequency-dependent velocity factor

The phase velocity of seismic waves varies with the propagation frequency, and thus frequency-dependent phenomena appear when CO₂ gas is injected into a reservoir. By dynamically considering these phenomena with reservoir conditions it is thus feasible to extract the frequency-dependent velocity factor with the aim of monitoring changes in the reservoir both before and after CO₂ injection. In the paper, we derive a quantitative expression for the frequency-dependent factor based on the Robinson seismic convolution model. In addition, an inversion equation with a frequency-dependent velocity factor is constructed, and a procedure is implemented using the following four processing steps: decomposition of the spectrum by generalized S transform, wavelet extraction of cross-well seismic traces, spectrum equalization processing, and an extraction method for frequency-dependent velocity factor based on the damped least-square algorithm. An attenuation layered model is then established based on changes in the Q value of the viscoelastic medium, and spectra of migration profiles from forward modeling are obtained and analyzed. Frequency-dependent factors are extracted and compared, and the effectiveness of the method is then verified using a synthetic data. The frequency-dependent velocity factor is finally applied to target processing and oil displacement monitoring based on real seismic data obtained before and after CO₂ injection in the G89 well block within Shengli oilfield. Profiles and slices of the frequency-dependent factor determine its ability to indicate differences in CO₂ flooding, and the predicting results are highly consistent with those of practical investigations within the well block.     

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.     

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.     

CO2气体在地震孕育与发生过程中的作用及其机理研究

通过甘东南地震危险区地震宏观异常观测实例分析,系统研究CO_2气体在地震构造活动过程中的地球化学演化过程,以及发生的一系列酸碱平衡和氧化还原化学反应。结果表明地下深部CO_2气体不仅是其他微量气体Rn等向地表方向运移的载体,而且参与酸碱平衡和氧化还原反应,是地下深部生物化学反应的重要影响条件,地震前兆及宏观异常现象的发生常伴随着CO_2的异常变化。因此CO_2气体可以作为良好的示踪气体,在地震前兆观测及重大异常落实中尤其要重视对其进行监测。     

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.     

Measurements of the elastic and anelastic properties of sandstone flooded with supercritical CO2

The aim of this study was to investigate the effects of supercritical CO₂ (scCO₂) injection on the elastic and anelastic properties of sandstone at seismic and ultrasonic frequencies. We present the results of the low‐frequency and ultrasonic experiments conducted on water‐saturated sandstone (Donnybrook, Western Australia) flooded with scCO₂. The sandstone was cut in the direction perpendicular to a formation bedding plane and tested in a Hoek triaxial pressure cell. During the experiments with scCO_2, the low‐frequency and ultrasonic systems and the pump dispensing scCO₂ were held at a temperature of 42°C. The elastic parameters obtained for the sandstone with scCO₂ at seismic (0.1 Hz–100 Hz) and ultrasonic (～0.5 MHz) frequencies are very close to those for the dry rock. The extensional attenuation was also measured at seismic frequencies for the dry, water‐saturated, and scCO₂‐injected sandstones. The applicability of Gassmann's fluid substitution theory to obtained results was also tested during the experiments.     

地震复谱分解技术及其在烃类检测中的应用

谱分解技术在地震解释领域已得到广泛应用,但常用的谱分解方法存在两方面的不足.一是时间分辨率低,难以对薄层进行刻画;二是在烃类检测中多解性强,难以区分流体类型.为了改善该问题,本文提出一种基于地震复谱分解技术的烃类检测方法.复谱分解是指用一个包含多个不同频率Ricker子波的复子波库对地震道进行分解,从而得到时变子波频率和相位信息的过程.借助稀疏反演技术复谱分解可以获得高分辨率的时频能量谱和时频相位谱.本文首先通过拟合算例验证了复谱分解方法刻画薄层的能力以及求取子波频率和相位的准确性.然后利用基于Kelvin-Voigt模型的黏弹波动方程数值模拟对衰减引起子波相位改变的原因进行了分析.最后通过实际资料应用展示了本文方法在储层预测中的高时间分辨率优势,验证了利用子波相位信息识别气藏的有效性.     

Comparison of surface seismic sources at the CO2SINK site, Ketzin, Germany

In 2004 three seismic surface sources (VIBSIST, accelerated weight drop and MiniVib) were tested in a pilot study at the Ketzin test site, Germany, a study site for geological storage of CO₂ (EU project CO₂SINK). The main objectives of this pilot study were to 1) evaluate the response of the Ketzin site to reflection seismics, especially at the planned injection depth, 2) test different acquisition parameters and 3) use the results to guide the planning of the 3D survey. As part of these objectives, we emphasize the source performance comparison in this study. The sources were tested along two perpendicular lines of 2.4 km length each. Data were acquired by shooting at all stations (source and receiver spacing of 20 m) on both lines, allowing common‐midpoint stacked sections to be produced. The sources' signal characteristics based on signal‐to‐noise ratio, signal penetration and frequency content of raw shot records were analysed and stacked sections were compared. The results show that all three surface sources are suitable for reflection seismic studies down to a depth of about 1 km and provide enough bandwidth for resolving the geological targets at the site, i.e., the Weser and Stuttgart Formations. Near surface conditions, especially a thick weathering layer present in this particular area, strongly influence the data quality, as indicated by the difference in reflectivity and signal‐to‐noise ratio of the two common‐midpoint lines. The stacked sections of the MiniVib source show the highest frequency signals down to about 500 ms traveltime (approximately 500 m depth) but also the shallowest signal penetration depth. The VIBSIST source generates signals with the highest signal‐to‐noise ratio and greatest signal penetration depth of the tested sources. In particular, reflections below 900 ms (approximately 1 km depth) are best imaged by the VIBSIST source. The weight drop performance lies in between these two sources and might be recommended as an appropriate source for a 3D survey at this site because of the shorter production time compared to the VIBSIST and MiniVib sources.     

2014年10月7日云南景谷MS6.6地震热红外异常

以中国静止气象卫星FY-2C/2E红外遥感资料为依据,采用小波变换与功率谱相对变化法对2014年10月7日云南景谷MS6.6地震进行震例数据处理和分析。结果显示:震前热异常具有明显时空特征,2014年8月上旬震中及其周围区域出现热红外异常,异常区范围随时间逐渐扩大、幅度陡增;8月30日左右异常面积达到最大,异常持续时间近一个月,异常幅度最大值达到平均值的14倍多,最大值45天后发生地震;该次地震热红外异常特征显著,容易识别,可作为一种识别地震热异常信息的判据,进一步验证卫星热红外异常在地震预判方面的作用。     

Continuous wavelet transform,theoretical aspects and application to aeromagnetic data at the Huanghua Depression,Dagang Oilfield,China

We use the continuous wavelet transform based on complex Morlet wavelets, which has been developed to estimate the source distribution of potential fields. For magnetic anomalies of adjacent sources, they always superimpose upon each other in space and wavenumber, making the identification of magnetic sources problematic. Therefore, a scale normalization factor, a^?n, is introduced on the wavelet coefficients to improve resolution in the scalogram. By theoretical modelling, we set up an approximate linear relationship between the pseudo‐wavenumber and source depth. The influences of background field, random noise and magnetization inclination on the continuous wavelet transform of magnetic anomalies are also discussed and compared with the short‐time Fourier transform results. Synthetic examples indicate that the regional trend has little effect on our method, while the influence of random noise is mainly imposed on shallower sources with higher wavenumbers. The source horizontal position will be affected by the change of magnetization direction, whereas the source depth remains unchanged. After discussing the performance of our method by showing the results of various synthetic tests, we use this method on the aeromagnetic data of the Huanghua depression in central China to define the distribution of volcanic rocks. The spectrum slices in different scales are used to determine horizontal positions of volcanic rocks and their source depths are estimated from the modulus maxima of complex coefficients, which is in good accordance with drilling results.     

Anomalous CO2 content in the Gallicano thermo-mineral spring (Serchio Valley,Italy) before the 21 June 2013, Alpi Apuane earthquake (M = 5.2)

Since late 2002, a continuous automatic monitoring network is operating in Tuscany, Central Italy, to investigate the geochemical response of selected aquifers to local seismic activity. The monitoring is aimed at identifying possible earthquake geochemical precursors. The network is currently constituted by six stations, all equipped with sensors for the measurement of temperature, pH, redox potential, electrical conductivity, CO₂ and CH₄ dissolved concentration, that have been installed in the areas of highest seismic risk of the region. By combining geochemical data gathered from the automatic station of Gallicano (Garfagnana, Northern Tuscany), and obtained via chemical analyses of spring water samples collected during periodic field surveys in the area surrounding this station, the most significant aspects of the deep fluid circulation paths feeding the Gallicano thermo-mineral system have been investigated, and the geochemical baseline of the Gallicano spring defined. The CO₂ continuous signal recorded by the Gallicano automatic station has been then processed over the period 2003–2013 in the search for anomalies possibly related to local seismic activity. A substantial anomaly in CO₂ content has been observed at Gallicano in conjunction with the Alpi Apuane earthquake (M = 5.2) of 21 June 2013.     

基于布格重力异常研究运城盆地的构造特征

利用运城及周边地区的布格重力资料,对布格重力数据进行小波分解得到不同阶次的重力异常信息。依据各个阶次的小波变换结果并结合该区域的地质构造环境,对小波变换细节图中存在的特征较为明显的布格重力异常带进行深入分析,对异常所揭示的构造现象在地壳不同深度的分布特性等进行分析研究。结果表明:小波变换细节与区内地壳内部的隆起、凹陷及断裂构造具有很好的对应关系,可以对运城盆地的构造特征及地壳结构做更加深入的研究。     

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.     

时频域油气储层低频阴影检测

为了准确刻画地震信号的局部层次结构,实现高效率的三维地震资料瞬时谱分解,检测油气储层的低频阴影,构造了广义S变换.广义S变换通过引入两个参数,改造S变换的小波函数,使其小波能根据信号处理的具体应用需要而调整.模型信号的仿真发现,广义S变换具有更加优越和灵活可调的时频聚集性能.文中分析了低频阴影的机理,并利用广义S变换对三维实际地震资料进行了瞬时谱分解,它不仅能检测油气储层的低频阴影,而且可以刻画油气储层的岩性边界和空间展布,减小油气储层检测的多解性.     

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.     

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.     

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.