Enhancements to prestack frequency‐wavenumber (f-k) migration

Geophysics ◽  
1991 ◽  
Vol 56 (1) ◽  
pp. 27-40 ◽  
Author(s):  
Z. Li ◽  
W. Lynn ◽  
R. Chambers ◽  
Ken Larner ◽  
Ray Abma

Prestack frequency‐wavenumber (f-k) migration is a particularly efficient method of doing both full prestack time migration and migration velocity analysis. Conventional implementations of the method, however, can encounter several drawbacks: (1) poor resolution and spatial aliasing noise caused by insufficient sampling in the offset dimension, (2) poor definition of steep events caused by insufficient sampling in the velocity dimension, and (3) inadequate handling of ray bending for steep events. All three of these problems can be mitigated with modifications to the prestack f-k algorithm. The application of linear moveout (LMO) in the offset dimension prior to migration reduces event moveout and hence increases the bandwidth of non‐spatially aliased signals. To reduce problems of interpolation for steep events, the number of constant‐velocity migrations can be economically increased by performing residual poststack migrations. Finally, migration with a dip‐dependent imaging velocity addresses the issue of ray bending and thereby improves the positioning of steep events. None of these enhancements substantially increases the computational effort of f-k migration. Prestack f-k migration possesses a limitation for which no solution is readily available. Where lateral velocity variation is modest, steep events (such as fault‐plane reflections in sediments) may not be imaged as well as by other migration approaches. This shortcoming results from the restriction that, in the prestack f-k approach, a single velocity field must serve to perform two different functions: imaging and stacking. Nevertheless, in areas of strong velocity variation and gentle to moderate dip, the detailed velocity control afforded by the prestack f-k method is an excellent source of geologic information.

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C151-C162 ◽  
Author(s):  
Qi Hao ◽  
Alexey Stovas ◽  
Tariq Alkhalifah

The offset-midpoint traveltime pyramid describes the diffraction traveltime of a point diffractor in homogeneous media. We have developed an analytic approximation for the P-wave offset-midpoint traveltime pyramid for homogeneous orthorhombic media. In this approximation, a perturbation method and the Shanks transform were implemented to derive the analytic expressions for the horizontal slowness components of P-waves in orthorhombic media. Numerical examples were shown to analyze the proposed traveltime pyramid formula and determined its accuracy and the application in calculating migration isochrones and reflection traveltime. The proposed offset-midpoint traveltime formula is useful for Kirchhoff prestack time migration and migration velocity analysis for orthorhombic media.


Geophysics ◽  
1989 ◽  
Vol 54 (6) ◽  
pp. 701-717 ◽  
Author(s):  
Ken Larner ◽  
Craig J. Beasley ◽  
Walt Lynn

Primarily through synthetic and field data examples, this paper reviews the benefits of recent developments in time migration of seismic data and reveals limitations, some of them fundamental, that keep elusive the goal of imaging steep events with full accuracy. Even where velocity varies only with depth and, hence, time migration should suffice, accurate imaging of very steep events requires that the velocity structure be known with considerable precision and be finely sampled in depth. This sensitivity of migration accuracy to detail in velocity structure is attributable to the sensitivity of ray bending for wide‐angle rays to detail in the velocity structure. Also, interestingly, the presence of ray bending at interfaces is seen to enhance the steep‐event accuracy of some algorithms (e.g., phase shift and cascaded finite‐difference) while it degrades the accuracy of others (for example, conventional Kirchhoff summation and the frequency‐wavenumber domain method of Stolt). Of the various time‐migration schemes, a cascading of the Stolt method is the most efficient while having steep‐event accuracy in the presence of significant vertical velocity variation. Its behavior in the presence of even mild lateral velocity variation, however, differs greatly from that of the other methods and must be taken into account. A case involving 3-D migration of 3-D survey data shows how different issues in imaging of the subsurface (two‐pass versus single‐pass 3-D migration and algorithm choice in the presence of mild lateral velocity variation) can become intertwined in practice and lead to confusion as to which of the issues is essential for accurate imaging of the subsurface.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S317-S331 ◽  
Author(s):  
Jianfeng Zhang ◽  
Zhengwei Li ◽  
Linong Liu ◽  
Jin Wang ◽  
Jincheng Xu

We have improved the so-called deabsorption prestack time migration (PSTM) by introducing a dip-angle domain stationary-phase implementation. Deabsorption PSTM compensates absorption and dispersion via an actual wave propagation path using effective [Formula: see text] parameters that are obtained during migration. However, noises induced by the compensation degrade the resolution gained and deabsorption PSTM requires more computational effort than conventional PSTM. Our stationary-phase implementation improves deabsorption PSTM through the determination of an optimal migration aperture based on an estimate of the Fresnel zone. This significantly attenuates the noises and reduces the computational cost of 3D deabsorption PSTM. We have estimated the 2D Fresnel zone in terms of two dip angles through building a pair of 1D migrated dip-angle gathers using PSTM. Our stationary-phase QPSTM (deabsorption PSTM) was implemented as a two-stage process. First, we used conventional PSTM to obtain the Fresnel zones. Then, we performed deabsorption PSTM with the Fresnel-zone-based optimized migration aperture. We applied stationary-phase QPSTM to a 3D field data. Comparison with synthetic seismogram generated from well log data validates the resolution enhancements.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. U87-U96 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin

One of the most serious problems in anisotropic velocity analysis is the trade-off between anisotropy and lateral heterogeneity, especially if velocity varies on a scale smaller than the maximum offset. We have developed a P-wave MVA (migration velocity analysis) algorithm for transversely isotropic (TI) models that include layers with small-scale lateral heterogeneity. Each layer is described by constant Thomsen parameters [Formula: see text] and [Formula: see text] and the symmetry-direction velocity [Formula: see text] that varies as a quadratic function of the distance along the layer boundaries. For tilted TI media (TTI), the symmetry axis is taken orthogonal to the reflectors. We analyzed the influence of lateral heterogeneity on image gathers obtained after prestack depth migration and found that quadratic lateral velocity variation in the overburden can significantly distort the moveout of the target reflection. Consequently, medium parameters beneath the heterogeneous layer(s) are estimated with substantial error, even when borehole information (e.g., check shots or sonic logs) is available. Because residual moveout in the image gathers is highly sensitive to lateral heterogeneity in the overburden, our algorithm simultaneously inverts for the interval parameters of all layers. Synthetic tests for models with a gently dipping overburden demonstrate that if the vertical profile of the symmetry-direction velocity [Formula: see text] is known at one location, the algorithm can reconstruct the other relevant parameters of TI models. The proposed approach helps increase the robustness of anisotropic velocity model-building and enhance image quality in the presence of small-scale lateral heterogeneity in the overburden.


Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1118-1127 ◽  
Author(s):  
Dimitri Bevc ◽  
James L. Black ◽  
Gopal Palacharla

We analyze how time migration mispositions events in the presence of lateral velocity variation by examining the impulse response of depth modeling followed by time migration. By examining this impulse response, we lay the groundwork for the development of a remedial migration operator that links time and depth migration. A simple theory by Black and Brzostowski predicted that the response of zero‐offset time migration to a point diffractor in a v(x, z) medium would be a distinctive, cusp‐shaped curve called a plume. We have constructed these plumes by migrating synthetic data using several time‐migration methods. We have also computed the shape of the plumes by two geometrical construction methods. These two geometrical methods compare well and explain the observed migration results. The plume response is strongly influenced by migration velocity. We have studied this dependency by migrating synthetic data with different velocities. The observed velocity dependence is confirmed by geometrical construction. A simple first‐order theory qualitatively explains the behavior of zero‐offset time migration, but a more complete understanding of migration velocity dependence in a v(x, z) medium requires a higher order finite‐offset theory.


Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. S141-S152 ◽  
Author(s):  
Jeffrey Shragge

Migration of seismic data from topography using methods based on finite-difference (FD) approximation to acoustic wave propagation commonly suffers from a number of imaging drawbacks due to the difficulty of applying FD stencils to irregular computational meshes. Altering the computational geometry from Cartesian to a topographic coordinate system conformal to the data acquisition surface can circumvent many of these issues. The coordinate transformation approach allows for acoustic wave propagation and the crosscorrelation and inverse-scattering imaging conditions to be posed and computed directly in topographic coordinates. Resulting reverse time migration (RTM) images may then be interpolated back to the Cartesian domain using the known inverse mapping. Orthogonal 2D topographic coordinates can be developed using known conformal mapping transforms and serve as the computational mesh for performing migration from topography. Impulse response tests demonstrate the accuracy of the 2D generalized acoustic wave propagation. RTM imaging examples show the efficacy of performing migration from topography directly from the data acquisition surface on topographic meshes and the ability to image complex near-surface structure even in the presence of strong lateral velocity variation.


2014 ◽  
Vol 556-562 ◽  
pp. 592-596
Author(s):  
Su Zhen Shi ◽  
Juan Li ◽  
Yi Chen Zhao ◽  
Li Biao Yang ◽  
Yao Tang ◽  
...  

In order to improve imaging precision of small structures and small fault blocks of coal seams, the prestack time migration method is used for imaging. Preserved amplitude processing (PAP) is applied to prestack gather firstly after geological data and original seismic data of the exploration area are fully understood. Initial root mean square velocity field is established through the method of picking up root mean square velocity on CRP gather. Then, a precise root mean square velocity model is created after continuous iteration and modification. Meanwhile, appropriate algorithm and migration parameters are selected during the migration process. Finally, the imaging of small fault blocks and small faults in the prestack time migration section is clear and migration is highly coinciding with the case disclosed by boreholes. It’s proved that prestack time migration is especially suitable for processing 3D seismic data of small faults and small fault blocks in coal seams with complicated geological conditions.


Sign in / Sign up

Export Citation Format

Share Document