Use of refraction, reflection, and wave-equation-based tomography for imaging beneath shallow gas: A Trinidad field data example

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE281-VE289 ◽  
Author(s):  
Nurul Kabir ◽  
Uwe Albertin ◽  
Min Zhou ◽  
Vishal Nagassar ◽  
Einar Kjos ◽  
...  

Shallow localized gas pockets cause challenging problems in seismic imaging because of sags and wipe-out zones they produce on imaged reflectors deep in the section. In addition, the presence of shallow gas generates strong surface-related and interbed multiples, making velocity updating very difficult. When localized gas pockets are very shallow, we have limited information to build a near-surface velocity model using ray-based reflection tomography alone. Diving-wave refraction tomography successfully builds a starting model for the very shallow part. Usual ray-based reflection tomography using single-parameter hyperbolic moveout might need many iterations to update the deeper part of the velocity model. In addition, the method generates a low-velocity anomaly in the deeper part of the model. We present an innovative method for building the final velocity model by combining refraction, reflection, and wave-equation-based tomography. Wave-equation-based tomography effectively generates a detailed update of a shallow velocity field, resolving the gas-sag problem. When applied as the last step, following the refraction and reflection tomography, it resolves the gas-sag problem but fails to remove the low-velocity anomaly generated by the reflection tomography in the deeper part of the model. To improve the methodology, we update the shallow velocity field using refraction tomography followed by wave-equation tomography before updating the deeper part of the model. This step avoids generating the low-velocity anomaly. Refraction and wave-equation-based tomography followed by reflection tomography generates a simpler velocity model, giving better focusing in the deeper part of the image. We illustrate how the methodology successfully improves resolution of gas anomalies and significantly reduces gas sag from an imaged section in the Greater Cassia area, Trinidad.

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE235-VE241 ◽  
Author(s):  
Juergen Fruehn ◽  
Ian F. Jones ◽  
Victoria Valler ◽  
Pranaya Sangvai ◽  
Ajoy Biswal ◽  
...  

Imaging in deep-water environments poses a specific set of challenges, both in data preconditioning and velocity model building. These challenges include scattered, complex 3D multiples, aliased noise, and low-velocity shallow anomalies associated with channel fills and gas hydrates. We describe an approach to tackling such problems for data from deep water off the east coast of India, concentrating our attention on iterative velocity model building, and more specifically the resolution of near-surface and other velocity anomalies. In the region under investigation, the velocity field is complicated by narrow buried canyons ([Formula: see text] wide) filled with low-velocity sediments, which give rise to severe pull-down effects; possible free-gas accumulation below an extensive gas-hydrate cap, causing dimming of the image below (perhaps as a result of absorption); and thin-channel bodies with low-velocity fill. Hybrid gridded tomography using a conjugate gradient solver (with [Formula: see text] vertical cell size) was applied to resolve small-scale velocity anomalies (with thicknesses of about [Formula: see text]). Manual picking of narrow-channel features was used to define bodies too small for the tomography to resolve. Prestack depth migration, using a velocity model built with a combination of these techniques, could resolve pull-down and other image distortion effects in the final image. The resulting velocity field shows high-resolution detail useful in identifying anomalous geobodies of potential exploration interest.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE93-VE100 ◽  
Author(s):  
James L. Simmons

A linear least-squares inversion is applied to the turning-ray first-arrival times of a shallow-marine seismic reflection data set to estimate the slowly varying (laterally and vertically) components of the near-surface velocity field. The velocity model is represented with a low-spatial-frequency parameterization (2D cubic B-splines) designed specifically for the predicted components of the data. This model parameterization effectively decouples the slowly varying background from the higher spatial-frequency component of the velocity field produced by shallow, low-velocity, gas-charged sands and allows the solution to be obtained in a single iteration. The observed first-arrival times (background and shallow anomaly-induced perturbations) and the slowly varying first-arrival times related to the background velocity are inverted separately. Similar velocity-model estimates result, demonstrating the decoupling imposed by the B-spline model parameterization. The background velocity and the low-velocity anomalies are best treated as separate inverse problems using very different model parameterizations. Ray tracing a synthetic model containing local low-velocity anomalies embedded in a smooth background does not accurately predict the anomaly-induced first-arrival time perturbations seen in the field data. Acoustic finite-difference waveform modeling shows that reflections and diffractions from the anomalies interfere with the diving-wave first arrivals. First-arrival times picked from the full-waveform synthetics more accurately predict the field data first-arrival times.


2021 ◽  
Vol 47 (2) ◽  
Author(s):  
Andrzej Michał Dalętka

Despite the increasing technological level of the reflection seismic method, the imaging of fold and thrust belts remains a demanding task, and usually leaves some questions regarding the dips, the shape of the subthrust structures or the most correct approach to velocity model building. There is no straightforward method that can provide structural representation of the near-surface geological boundaries and their velocities. The in-terpretation of refracted waves frequently remains the only available technique that may be used for this purpose, although one must be aware of its limitations which appear in the complex geological settings. In the presented study, the analysis of velocity values obtained in the shallow part of Carpathian orogenic wedge by means of various geophysical methods was carried out. It revealed the lack of consistency between the results of 3D refraction tomography and both the sonic log and uphole velocities. For that reason, instead of the indus-try-standard utilization of tomography, a novel, geologically-consistent method of velocity model building is pro-posed. In the near-surface part, the uphole velocities are assigned to the formations, documented by the surface geologic map. Interpreted time-domain horizons, supplemented by main thrusts, are used to make the velocity field fully-compatible with the litho-stratigraphic units of the Carpathians. The author demonstrates a retrospective overview of seismic data imaging in the area of the Polish Carpathian orogenic wedge and discusses the most recent global innovations in seismic methodology which are the key to successful hydrocarbon exploration in fold and thrust regions.


Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. B43-B53 ◽  
Author(s):  
Xianhuai Zhu ◽  
Kirk Wallace ◽  
Qingrong Zhu ◽  
Robert Hofer

Seismic imaging is challenging in the Bohai Bay PL19-3 area, offshore China. Bohai Bay Field is seismically obscured, making well penetrations the only reliable source of data for subsurface interpretation. With the help of turning-ray tomography, we are able to obtain a reliable near-surface velocity model approximately down to 700 m below sea level, using the first arrivals picked from streamer data. Depth migration using velocities estimated from turning-ray tomography has improved shallow structures and fault definition. However, reservoir level structures from 800 to 1500 m are still poorly imaged. A viscoelastic modeling study with assigned variable Q and shallow velocity profiles, with and without shallow gas-induced scatterers, demonstrates that scattering is the primary controlling phenomenon causing imaging difficulty within the obscured zone. Due to scattering, imaging tests at the target level were unsuccessful even with precise knowledge of velocity.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. U39-U47 ◽  
Author(s):  
Hui Liu ◽  
Hua-wei Zhou ◽  
Wenge Liu ◽  
Peiming Li ◽  
Zhihui Zou

First-arrival traveltime tomography is a popular approach to building the near-surface velocity models for oil and gas exploration, mining, geoengineering, and environmental studies. However, the presence of velocity-inversion interfaces (VIIs), across which the overlying velocity is higher than the underlying velocity, might corrupt the tomographic solutions. This is because most first-arrival raypaths will not traverse along any VII, such as the top of a low-velocity zone. We have examined the impact of VIIs on first-arrival tomographic velocity model building of the near surface using a synthetic near-surface velocity model. This examination confirms the severe impact of VIIs on first-arrival tomography. When the source-to-receiver offset is greater than the lateral extent of the VIIs, good near-surface velocity models can still be established using a multiscale deformable-layer tomography (DLT), which uses a layer-based model parameterization and a multiscale scheme as regularization. Compared with the results from a commercial grid-based tomography, the DLT delivers much better near-surface statics solutions and less error in the images of deep reflectors.


2018 ◽  
Vol 6 (4) ◽  
pp. SM63-SM70 ◽  
Author(s):  
Tian Jun ◽  
Peng Gengxin ◽  
Junru Jiao ◽  
Grace (Yan) Yan ◽  
Xianhuai Zhu

A special challenge for land seismic exploration is estimating velocities, in part due to complex near-surface structures, and in some instances because of rugose topography over foothills. We have developed an integrated turning-ray and reflection-tomographic method to face this challenge. First, turning-ray tomography is performed to derive a near-surface velocity-depth model. Then, we combine the near-surface model with the initial-subsurface model. Taking the combined model as starting model, we go through a reflection tomographic process to build the model for imaging. During reflection tomography, the near-surface model and subsurface models are jointly updated. Our method has been successfully applied to a 2D complex synthetic data example and a 3D field data example. The results demonstrate that our method derives a very decent model even when there is no reflection information available in a few hundred meters underneath the surface. Joint tomography can lead to geologic plausible models and produce subsurface images with high fidelity.


2018 ◽  
Vol 6 (4) ◽  
pp. SM27-SM37 ◽  
Author(s):  
Jing Li ◽  
Kai Lu ◽  
Sherif Hanafy ◽  
Gerard Schuster

Two robust imaging technologies are reviewed that provide subsurface geologic information in challenging environments. The first one is wave-equation dispersion (WD) inversion of surface waves and guided waves (GW) for the shear-velocity (S-wave) and compressional-velocity (P-wave) models, respectively. The other method is traveltime inversion for the velocity model, in which supervirtual refraction interferometry (SVI) is used to enhance the signal-to-noise ratio of far-offset refractions. We have determined the benefits and liabilities of both methods with synthetic seismograms and field data. The benefits of WD are that (1) there is no layered-medium assumption, as there is in conventional inversion of dispersion curves. This means that 2D or 3D velocity models can be accurately estimated from data recorded by seismic surveys over rugged topography, and (2) WD mostly avoids getting stuck in local minima. The liability is that WD for surface waves is almost as expensive as full-waveform inversion (FWI) and, for Rayleigh waves, only recovers the S-velocity distribution to a depth no deeper than approximately 1/2 to 1/3 wavelength of the lowest-frequency surface wave. The limitation for GW is that, for now, it can estimate the P-velocity model by inverting the dispersion curves from GW propagating in near-surface low-velocity zones. Also, WD often requires user intervention to pick reliable dispersion curves. For SVI, the offset of usable refractions can be more than doubled, so that traveltime tomography can be used to estimate a much deeper model of the P-velocity distribution. This can provide a more effective starting velocity model for FWI. The liability is that SVI assumes head-wave first arrivals, not those from strong diving waves.


Sign in / Sign up

Export Citation Format

Share Document