Geostatistical traveltime tomography in elliptically anisotropic media

In geological materials, anisotropy may arise due to different mechanisms and can be found at different scales. Neglecting anisotropy in traveltime tomographic reconstruction leads to artefacts that can obscure important subsurface features. In this paper, a geostatistical tomography algorithm to invert cross‐hole traveltime data in elliptically anisotropic media is presented. The advantages of geostatistical tomography are that the solution is regularized by the covariance of the model parameters, that known model parameters can be used as constraints and fitted exactly or within a prescribed variance and that stochastic simulations can be performed to appraise the variability of the solution space. The benefits of the algorithm to image anisotropic media are illustrated by two examples using synthetic georadar data and real seismic data. The first example confirms suspected electromagnetic anisotropy in the vadose zone caused by relatively rapid water content variations with respect to wavelength at georadar frequencies. The second presents how sonic log data can be used to constrain the inversion of cross‐well seismic data and how geostatistical simulations can be used to infer parameter uncertainty. Results of both examples show that considering anisotropy yields a better fit to the data at high ray angles and reduces reconstruction artefacts.     

Monitoring and volumetric estimation of injected CO2 using 4D seismic,petrophysical data,core measurements and well logging: a case study at Ketzin,Germany

More than 50 000 tons of CO₂ 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 CO₂ had been injected. We show here that rather complex time‐lapse signatures of this CO₂ 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 CO₂ saturation levels allow an estimate of the total amount of CO₂ visible in the seismic data to be made. These estimates are somewhat lower than the actual amount of CO₂ 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 CO₂ storage site using seismic methods.     

A genetic algorithm for filter design to enhance features in seismic images

We present a novel method to enhance seismic data for manual and automatic interpretation. We use a genetic algorithm to optimize a kernel that, when convolved with the seismic image, appears to enhance the internal characteristics of salt bodies and the sub‐salt stratigraphy. The performance of the genetic algorithm was validated by the use of test images prior to its application on the seismic data. We present the evolution of the resulting kernel and its convolved image. This image was analysed by a seismic interpreter, highlighting possible advantages over the original one. The effects of the kernel were also subject to an automatic interpretation technique based on principal component analysis. Statistical comparison of these results with those from the original image, by means of the Mann‐Whitney U‐test, proved the convolved image to be more appropriate for automatic interpretation.     

Non‐uniqueness with refraction inversion – the Mt Bulga shear zone

The tau‐p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non‐uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non‐uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable.