Search

The undesired drill-bit deviation is a source of drilling risks and requires monitoring. Seismic-while-drilling is one potential method to achieve this and has been tested in a number of previous studies. In August 2020 in Örebro, Sweden, we conducted an experiment to test the feasibility of seismic-while-drilling drill-bit positioning and other applications of the method. We used the hammer drill-bit signal generated while drilling a 200 m deep well in hardrock conditions and implemented vertical stacking of the subsequent impulsive signals from the hammer source, generating enhanced direct arrivals from the drill-bit. Then, we used the relative arrival times to estimate the drill-bit position for selected bit depths, confirming our methodology with two numerical studies. We successfully estimated the position of the drill-bit for the numerical examples and for several of the real data examples, with the accuracy dependent on the receiver array geometry and the quality of the data. We conclude that this drill-bit positioning method shows potential for near-real-time monitoring in drilling operations that could be applicable for both impulsive and noise-retrieved drill-bit signals.

The injection of CO2 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 CO2 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 CO2 injection has a rather minor impact on the seismic amplitudes. Therefore, the observed amplitude anomaly is interpreted to be mainly due to CO2 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 CO2 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 CO2 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 CO2 saturation and used for its quantitative assessment. The estimated CO2 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.

The injection of CO2 at the Ketzin pilot CO2 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 CO2 was injected into a ...

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.

Seismic monitoring of an injected carbon dioxide (CO2) distribution at depth is an important issue in the geological storage of CO2. To help monitor changes in the subsurface during CO2 injection a series of 2D seismic surveys were acquired within the framework of the CO2SINK and CO2MAN projects at Ketzin, Germany at different stages of the injection process. Here we investigate using seismic full waveform inversion as a qualitative tool for time-lapse seismic monitoring given the constraints of the limited maximum offsets of the 2D seismic data. Prior to applying the inversion to the real data we first made a number of benchmark tests on synthetic data using a similar geometry as in the real data. Results from the synthetic benchmark tests show that it is difficult to recover the true value of the velocity anomaly due to the injection but that it is possible to qualitatively locate the distribution of the injected CO2. After the synthetic studies, we applied seismic full waveform inversion on the real time-lapse data from the Ketzin site along with conventional time-lapse processing. Both methods show a similar qualitative distribution of the injected CO2 and agree well with expectations based upon more extensive 3D time-lapse monitoring in the area.

A combination of seismic and geoelectric processing was studied by means of a structurally constrained inversion approach. Structural constraints were interpreted from the seismic data and integrated into the geoelectric inversion through a local regularization, which allowed inverted resistivities to behave discontinuously across defined boundaries. This arranged seismic processing and constrained resistivity inversion in a sequential workflow, making the generic assumption that the petrophysical parameters of both methods change across common lithostructural boundaries. We evaluated the approach using a numerical example and a real data example from the Ketzin [Formula: see text] pilot storage site, Germany. The latter demonstrated the efficiency of this approach for combining 4D seismic and surface-downhole geoelectric data. In consistence with the synthetic example, the constrained resistivity inversions produced clearer delineated images along the boundary between caprock and reservoir formation. Near the [Formula: see text]-flooded reservoir, the seismic and geoelectric time-lapse anomalies correlated well. At some distance to the downhole electrodes, however, the geoelectric images conveyed a notably lower resolution in comparison to the corresponding seismic images. Both methods confirm a northwesterly trend for the [Formula: see text] migration at the Ketzin site, although a rather northerly direction was initially expected. The results demonstrate the relevance of the presented approach for the combination of both methods for integrated geophysical [Formula: see text] storage monitoring.

A method for static correction of time-lapse differences in reflection arrival times of time-lapse prestack seismic data is presented. These arrival-time differences are typically caused by changes in the near-surface velocities between the acquisitions and had a detrimental impact on time-lapse seismic imaging. Trace-to-trace time shifts of the data sets from different vintages are determined by crosscorrelations. The time shifts are decomposed in a surface-consistent manner, which yields static corrections that tie the repeat data to the baseline data. Hence, this approach implies that new refraction static corrections for the repeat data sets are unnecessary. The approach is demonstrated on a 4D seismic data set from the Ketzin [Formula: see text] pilot storage site, Germany, and is compared with the result of an initial processing that was based on separate refraction static corrections. It is shown that the time-lapse difference static correction approach reduces 4D noise more effectively than separate refraction static corrections and is significantly less labor intensive.

More than 50 000 tons of CO2 have been injected at Ketzin into the Stuttgart Formation, a saline aquifer, at approximately 620 m depth, as of summer 2011. We present here results from the 1 st repeat 3D seismic survey that was performed at the site in autumn 2009, after about 22 000 tons of CO2 had been injected. We show here that rather complex time-lapse signatures of this CO2 can be clearly observed within a radius of about 300 m from the injection well. The highly irregular amplitude response within this radius is attributed to the heterogeneity of the injection reservoir. Time delays to a reflection below the injection level are also observed. Petrophysical measurements on core samples and geophysical logging of CO2 saturation levels allow an estimate of the total amount of CO2 visible in the seismic data to be made. These estimates are somewhat lower than the actual amount of CO2 injected at the time of the survey and they are dependent upon the choice of a number of parameters. In spite of some uncertainty, the close agreement between the amount injected and the amount observed is encouraging for quantitative monitoring of a CO2 storage site using seismic methods.

SUMMARY Recent seismic evidence suggested that most oceanic plate hydration is associated with trenchouter rise faulting prior to subduction. Hydration at trenches may have a significant impact on the subduction zone water cycle. Previous seismic experiments conducted to the northwest of Nicoya Peninsula, Northern Costa Rica, have shown that the subducting Cocos lithosphere is pervasively altered, which was interpreted to be due to both hydration (serpentinization) and fracturing of the crustal and upper-mantle rocks. New seismic wide-angle reflection and refraction data were collected along two profiles, running parallel to the Middle American trench axis offshore of central Nicaragua, revealing lateral changes of the seismic properties of the subducting lithosphere. Seismic structure along both profiles is characterized by low velocities both in the crust and upper mantle. Velocities in the uppermost mantle are found to be in the range 7.3–7.5 km s −1 ; thus are 8–10 per cent lower than velocities typical for unaltered peridotites and hence confirm the assumption that serpentinization is a common process at the trench-outer rise area offshore of Nicaragua. In addition, a prominent velocity anomaly occurred within the crust beneath two seamounts. Here, velocity reduction may indicate increased porosity and perhaps permeability, supporting the idea that seamounts serve as sites for water percolation and circulation.

A seismic wide‐angle and refraction experiment was conducted offshore of Nicaragua in the Middle American Trench to investigate the impact of bending‐related normal faulting on the seismic properties of the oceanic lithosphere prior to subduction. On the basis of the reflectivity pattern of multichannel seismic reflection (MCS) data it has been suggested that bending‐related faulting facilitates hydration and serpentinization of the incoming oceanic lithosphere. Seismic wide‐angle and refraction data were collected along a transect which extends from the outer rise region not yet affected by subduction into the trench northwest of the Nicoya Peninsula, where multibeam bathymetric data show prominent normal faults on the seaward trench slope. A tomographic joint inversion of seismic refraction and wide‐angle reflection data yield anomalously low seismic P wave velocities in the crust and uppermost mantle seaward of the trench axis. Crustal velocities are reduced by 0.2–0.5 km s−1 compared to normal mature oceanic crust. Seismic velocities of the uppermost mantle are 7.6–7.8 km s−1 and hence 5–7% lower than the typical velocity of mantle peridotite. These systematic changes in P wave velocity from the outer rise toward the trench axis indicate an evolutionary process in the subducting slab consistent with percolation of seawater through the faulted and fractured lithosphere and serpentinization of mantle peridotites. If hydration is indeed affecting the seismic properties of the mantle, serpentinization might be reaching 12–17% in the uppermost 3–4 km of the mantle, depending on the unknown degree of fracturing and its impact on the elastic properties of the subducting lithosphere.