Estimation of anisotropy and anisotropic 3-D prestack depth migration, offshore Zaire

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1495-1513 ◽  
Author(s):  
Greg Ball

Chevron (operator), Teikoku Oil, and Unocal hold the offshore Zaire concession, West Africa. Recently, Unocal performed an anisotropic 3-D prestack depth migration in an attempt to optimally image tilted fault blocks in the presalt section. Depth migration was required because large‐displacement faults juxtaposed a high‐velocity massive carbonate and a low‐velocity marl. Large compressional‐wave anisotropy caused by thin layering was inferred for most of the depth section. Evidence on the nature and magnitude of the anisotropy came from comparing stacking velocities with well velocities, analyzing a multioffset vertical seismic profile (VSP), thin‐layer modeling in the massive carbonate using well log data, and rock physics laboratory velocity measurements on core from the marl. The subsurface was characterized as being locally transversely isotropic and the Thomsen parameters were used to describe the anisotropy. This required that five parameters be specified at each subsurface location to represent the compressional‐wave 3-D velocity field. The observed anisotropy deviated far from the elliptical condition, and modeling indicated that isotropic migration would exhibit misfocusing and mispositioning of events, and its quality would be dip‐dependent. The anisotropic migration was performed using Kirchhoff summation for which the traveltimes were determined using a finite‐difference scheme. Interpreted horizons in the migrated depth section accurately tied well depths.

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 1782-1791 ◽  
Author(s):  
M. Graziella Kirtland Grech ◽  
Don C. Lawton ◽  
Scott Cheadle

We have developed an anisotropic prestack depth migration code that can migrate either vertical seismic profile (VSP) or surface seismic data. We use this migration code in a new method for integrated VSP and surface seismic depth imaging. Instead of splicing the VSP image into the section derived from surface seismic data, we use the same migration algorithm and a single velocity model to migrate both data sets to a common output grid. We then scale and sum the two images to yield one integrated depth‐migrated section. After testing this method on synthetic surface seismic and VSP data, we applied it to field data from a 2D surface seismic line and a multioffset VSP from the Rocky Mountain Foothills of southern Alberta, Canada. Our results show that the resulting integrated image exhibits significant improvement over that obtained from (a) the migration of either data set alone or (b) the conventional splicing approach. The integrated image uses the broader frequency bandwidth of the VSP data to provide higher vertical resolution than the migration of the surface seismic data. The integrated image also shows enhanced structural detail, since no part of the surface seismic section is eliminated, and good event continuity through the use of a single migration–velocity model, obtained by an integrated interpretation of borehole and surface seismic data. This enhanced migrated image enabled us to perform a more robust interpretation with good well ties.


Geophysics ◽  
1990 ◽  
Vol 55 (9) ◽  
pp. 1235-1241 ◽  
Author(s):  
Jan Douma

Traveltime inversion of multioffset VSP data can be used to determine the depths of the interfaces in layered media. Many inversion schemes, however, assume isotropy and consequently may introduce erroneous structures for anisotropic media. Synthetic traveltime data are computed for layered anisotropic media and inverted assuming isotropic layers. Only the interfaces between these layers are inverted. For a medium consisting of a horizontal isotropic low‐velocity layer on top of a transversely isotropic layer with a horizontal axis of symmetry (e.g., anisotropy due to aligned vertical cracks), 2-D isotropic inversion results in an anticline. For a given axis of symmetry the form of this anticline depends on the azimuth of the source‐borehole direction. The inversion result is a syncline (in 3-D a “bowl” structure), regardless of the azimuth of the source‐borehole direction for a vertical axis of symmetry of the transversely isotropic layer (e.g., anisotropy due to horizontal bedding).


Geophysics ◽  
1999 ◽  
Vol 64 (5) ◽  
pp. 1630-1636 ◽  
Author(s):  
Ayon K. Dey ◽  
Larry R. Lines

In seismic exploration, statistical wavelet estimation and deconvolution are standard tools. Both of these processes assume randomness in the seismic reflectivity sequence. The validity of this assumption is examined by using well‐log synthetic seismograms and by using a procedure for evaluating the resulting deconvolutions. With real data, we compare our wavelet estimations with the in‐situ recording of the wavelet from a vertical seismic profile (VSP). As a result of our examination of the randomness assumption, we present a fairly simple test that can be used to evaluate the validity of a randomness assumption. From our test of seismic data in Alberta, we conclude that the assumption of reflectivity randomness is less of a problem in deconvolution than other assumptions such as phase and stationarity.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. C107-C123
Author(s):  
Ivan Pšenčík ◽  
Bohuslav Růžek ◽  
Petr Jílek

We have developed a practical concept of compressional wave (P-wave) traveltime inversion in weakly to moderately anisotropic media of arbitrary symmetry and orientation. The concept provides sufficient freedom to explain and reproduce observed anisotropic seismic signatures to a high degree of accuracy. The key to this concept is the proposed P-wave anisotropy parameterization (A-parameters) that, together with the use of the weak-anisotropy approximation, leads to a significantly simplified theory. Here, as an example, we use a simple and transparent formula relating P-wave traveltimes to 15 P-wave A-parameters describing anisotropy of arbitrary symmetry. The formula is used in the inversion scheme, which does not require any a priori information about anisotropy symmetry and its orientation, and it is applicable to weak and moderate anisotropy. As the first step, we test applicability of the proposed scheme on a blind inversion of synthetic P-wave traveltimes generated in vertical seismic profile experiments in homogeneous models. Three models of varying anisotropy are used: tilted orthorhombic and triclinic models of moderate anisotropy (approximately 10%) and an orthorhombic model of strong anisotropy (>25%) with a horizontal plane of symmetry. In all cases, the inversion yields the complete set of 15 P-wave A-parameters, which make reconstruction of corresponding phase-velocity surfaces possible with high accuracy. The inversion scheme is robust with respect to noise and the source distribution pattern. Its quality depends on the angular illumination of the medium; we determine how the absence of nearly horizontal propagation directions affects inversion accuracy. The results of the inversion are applicable, for example, in migration or as a starting model for inversion methods, such as full-waveform inversion, if a model refinement is desired. A similar procedure could be designed for the inversion of S-wave traveltimes in anisotropic media of arbitrary symmetry.


Geophysics ◽  
2002 ◽  
Vol 67 (1) ◽  
pp. 98-109 ◽  
Author(s):  
John H. Bradford ◽  
D. S. Sawyer

It is common in shallow seismic studies for the compressional‐wave velocity in unconsolidated sediments to increase by a factor of four or more at the transition from dry or partial water saturation to full saturation. Under these conditions, conventional NMO velocity analysis fails and leads to large depth and layer thickness estimates if the Dix equation is assumed valid. Prestack depth migration (PSDM) is a means of improving image accuracy. A comparison of PSDM with conventional NMO processing for three field examples from differing hydrogeologic environments illustrates that PSDM can significantly improve image quality and accuracy.


2015 ◽  
Vol 3 (3) ◽  
pp. SW51-SW56
Author(s):  
Xiaomin Zhao ◽  
Shengwen Jin

Prestack Kirchhoff depth migration is commonly used in borehole seismic imaging, where there is uneven illumination due to the limitations of the source-receiver geometry. A new vertical seismic profile (VSP) migration/imaging workflow has been established that incorporates the structure-dip information derived from a newly developed structure tensor analysis into the existing VSP Kirchhoff migration/imaging technique. This allows us to better image the structures in the vicinity of a borehole and the far-field dipping events away from the borehole. We tested the workflow with the HESS salt model. The results were compared with those from reverse time migration, which found that Kirchhoff migration combined with structure-dip information not only reduced ambiguities of the imaging result but also allowed for imaging dip structures (e.g., fault) in the far region from the borehole. This allows for imaging dip structures and provides a useful extension of existing VSP imaging capabilities using Kirchhoff migration.


Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. S25-S36 ◽  
Author(s):  
Ernesto V. Oropeza ◽  
George A. McMechan

An efficient Kirchhoff-style prestack depth migration, called “parsimonious” migration, was developed a decade ago for isotropic 2D and 3D media by using measured slownesses to reduce the amount of ray tracing by orders of magnitude. It is conceptually similar to “map” migration, but its implementation has some differences. We have extended this approach to 2D tilted transversely isotropic (TTI) media and illustrated it with synthetic P-wave data. Although the framework of isotropic parsimonious may be retained, the extension to TTI media requires redevelopment of each of the numerical components, calculation of the phase and group velocity for TTI media, development of a new two-point anisotropic ray tracer, and substitution of an initial-angle isotropic shooting ray-trace algorithm for an anisotropic one. The model parameterization consists of Thomsen’s parameters ([Formula: see text], [Formula: see text], [Formula: see text]) and the tilt angle of the symmetry axis of the TI medium. The parsimonious anisotropic migration algorithm is successfully applied to synthetic data from a TTI version of the Marmousi2 model. The quality of the image improves by weighting the impulse response by the calculation of the anisotropic Fresnel radius. The accuracy and speed of this migration makes it useful for anisotropic velocity model building. The elapsed computing time for 101 shots for the Marmousi2 TTI model is 35 s per shot (each with 501 traces) in 32 Opteron cores.


Sign in / Sign up

Export Citation Format

Share Document