Application of waveform tomography to marine seismic reflection data from the Queen Charlotte Basin of western Canada

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. B55-B70 ◽  
Author(s):  
E. M. Takam Takougang ◽  
A. J. Calvert

To obtain a higher resolution quantitative P-wave velocity model, 2D waveform tomography was applied to seismic reflection data from the Queen Charlotte sedimentary basin off the west coast of Canada. The forward modeling and inversion were implemented in the frequency domain using the visco-acoustic wave equation. Field data preconditioning consisted of f-k filtering, 2D amplitude scaling, shot-to-shot amplitude balancing, and time windowing. The field data were inverted between 7 and 13.66 Hz, with attenuation introduced for frequencies ≥ 10.5 Hz to improve the final velocity model; two different approaches to sampling the frequencies were evaluated. The limited maximum offset of the marine data (3770 m) and the relatively high starting frequency (7 Hz) were the main challenges encountered during the inversion. An inversion strategy that successively recovered shallow-to-deep structures was designed to mitigate these issues. The inclusion of later arrivals in the waveform tomography resulted in a velocity model that extends to a depth of approximately 1200 m, twice the maximum depth of ray coverage in the ray-based tomography. Overall, there is a good agreement between the velocity model and a sonic log from a well on the seismic line, as well as between modeled shot gathers and field data. Anomalous zones of low velocity in the model correspond to previously identified faults or their upward continuation into the shallow Pliocene section where they are not readily identifiable in the conventional migration.

Geophysics ◽  
2011 ◽  
Vol 76 (1) ◽  
pp. SA19-SA34 ◽  
Author(s):  
Shohei Minato ◽  
Toshifumi Matsuoka ◽  
Takeshi Tsuji ◽  
Deyan Draganov ◽  
Jürg Hunziker ◽  
...  

Crosswell reflection method is a high-resolution seismic imaging method that uses recordings between boreholes. The need for downhole sources is a restrictive factor in its application, for example, to time-lapse surveys. An alternative is to use surface sources in combination with seismic interferometry. Seismic interferometry (SI) could retrieve the reflection response at one of the boreholes as if from a source inside the other borehole. We investigate the applicability of SI for the retrieval of the reflection response between two boreholes using numerically modeled field data. We compare two SI approaches — crosscorrelation (CC) and multidimensional deconvolution (MDD). SI by MDD is less sensitive to underillumination from the source distribution, but requires inversion of the recordings at one of the receiver arrays from all the available sources. We find that the inversion problem is ill-posed, and propose to stabilize it using singular-value decomposition. The results show that the reflections from deep boundaries are retrieved very well using both the CC and MDD methods. Furthermore, the MDD results exhibit more realistic amplitudes than those from the CC method for downgoing reflections from shallow boundaries. We find that the results retrieved from the application of both methods to field data agree well with crosswell seismic-reflection data using borehole sources and with the logged P-wave velocity.


2016 ◽  
Vol 4 (3) ◽  
pp. SH1-SH9
Author(s):  
Steven D. Sloan ◽  
J. Tyler Schwenk ◽  
Robert H. Stevens

Variability of material properties in the shallow subsurface presents challenges for near-surface geophysical methods and exploration-scale applications. As the depth of investigation decreases, denser sampling is required, especially of the near offsets, to accurately characterize the shallow subsurface. We have developed a field data example using high-resolution shallow seismic reflection data to demonstrate how quickly near-surface properties can change over short distances and the effects on field data and processed sections. The addition of a relatively thin, 20 cm thick, low-velocity layer can lead to masked reflections and an inability to map shallow reflectors. Short receiver intervals, on the order of 10 cm, were necessary to identify the cause of the diminished data quality and would have gone unknown using larger, more conventional station spacing. Combined analysis of first arrivals, surface waves, and reflections aided in determining the effects and extent of a low-velocity layer that inhibited the identification and constructive stacking of the reflection from a shallow water table using normal-moveout-based processing methods. Our results also highlight the benefits of using unprocessed gathers to pragmatically guide processing and interpretation of seismic data.


2016 ◽  
Author(s):  
David K. Smythe

Abstract. North American shale basins differ from their European counterparts in that the latter are one to two orders of magnitude smaller in area, but correspondingly thicker, and are cut or bounded by normal faults penetrating from the shale to the surface. There is thus an inherent risk of groundwater resource contamination via these faults during or after unconventional resource appraisal and development. US shale exploration experience cannot simply be transferred to the UK. The Bowland Basin, with 1900 m of Lower Carboniferous shale, is in the vanguard of UK shale gas development. A vertical appraisal well to test the shale by hydraulic fracturing (fracking), the first such in the UK, triggered earthquakes. Re-interpretation of the 3D seismic reflection data, and independently the well casing deformation data, both show that the well was drilled through the earthquake fault, and did not avoid it, as concluded by the exploration operator. Faulting in this thick shale is evidently difficult to recognise. The Weald Basin is a shallower Upper Jurassic unconventional oil play with stratigraphic similarities to the Bakken play of the Williston Basin, USA. Two Weald licensees have drilled, or have applied to drill, horizontal appraisal wells based on inadequate 2D seismic reflection data coverage. I show, using the data from the one horizontal well drilled to date, that one operator failed identify two small but significant through-going normal faults. The other operator portrayed a seismic line as an example of fault-free structure, but faulting had been smeared out by reprocessing. The case histories presented show that: (1) UK shale exploration to date is characterised by a low degree of technical competence, and (2) regulation, which is divided between four separate authorities, is not up to the task. If UK shale is to be exploited safely: (1) more sophisticated seismic imaging methods need to be developed and applied to both basins, to identify faults in shale with throws as small as 4–5 m, and (2) the current lax and inadequate regulatory regime must be overhauled, unified, and tightened up.


2017 ◽  
Vol 90 (2) ◽  
pp. 187-195
Author(s):  
A. I. Opara ◽  
C. C. Agoha ◽  
C. N. Okereke ◽  
U. P. Adiela ◽  
C. N. Onwubuariri ◽  
...  

Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. R59-R80 ◽  
Author(s):  
Michael Warner ◽  
Andrew Ratcliffe ◽  
Tenice Nangoo ◽  
Joanna Morgan ◽  
Adrian Umpleby ◽  
...  

We have developed and implemented a robust and practical scheme for anisotropic 3D acoustic full-waveform inversion (FWI). We demonstrate this scheme on a field data set, applying it to a 4C ocean-bottom survey over the Tommeliten Alpha field in the North Sea. This shallow-water data set provides good azimuthal coverage to offsets of 7 km, with reduced coverage to a maximum offset of about 11 km. The reservoir lies at the crest of a high-velocity antiformal chalk section, overlain by about 3000 m of clastics within which a low-velocity gas cloud produces a seismic obscured area. We inverted only the hydrophone data, and we retained free-surface multiples and ghosts within the field data. We invert in six narrow frequency bands, in the range 3 to 6.5 Hz. At each iteration, we selected only a subset of sources, using a different subset at each iteration; this strategy is more efficient than inverting all the data every iteration. Our starting velocity model was obtained using standard PSDM model building including anisotropic reflection tomography, and contained epsilon values as high as 20%. The final FWI velocity model shows a network of shallow high-velocity channels that match similar features in the reflection data. Deeper in the section, the FWI velocity model reveals a sharper and more-intense low-velocity region associated with the gas cloud in which low-velocity fingers match the location of gas-filled faults visible in the reflection data. The resulting velocity model provides a better match to well logs, and better flattens common-image gathers, than does the starting model. Reverse-time migration, using the FWI velocity model, provides significant uplift to the migrated image, simplifying the planform of the reservoir section at depth. The workflows, inversion strategy, and algorithms that we have used have broad application to invert a wide-range of analogous data sets.


Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1339-1347 ◽  
Author(s):  
Kate C. Miller ◽  
Steven H. Harder ◽  
Donald C. Adams ◽  
Terry O’Donnell

Shallow seismic reflection surveys commonly suffer from poor data quality in the upper 100 to 150 ms of the stacked seismic record because of shot‐associated noise, surface waves, and direct arrivals that obscure the reflected energy. Nevertheless, insight into lateral changes in shallow structure and stratigraphy can still be obtained from these data by using first‐arrival picks in a refraction analysis to derive a near‐surface velocity model. We have used turning‐ray tomography to model near‐surface velocities from seismic reflection profiles recorded in the Hueco Bolson of West Texas and southern New Mexico. The results of this analysis are interval‐velocity models for the upper 150 to 300 m of the seismic profiles which delineate geologic features that were not interpretable from the stacked records alone. In addition, the interval‐velocity models lead to improved time‐to‐depth conversion; when converted to stacking velocities, they may provide a better estimate of stacking velocities at early traveltimes than other methods.


Sign in / Sign up

Export Citation Format

Share Document