Prestack and poststack migration of crooked-line seismic reflection data: A case study from the South Portuguese Zone fold belt, southwestern Iberia

Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. B9-B18 ◽  
Author(s):  
C. Schmelzbach ◽  
C. Juhlin ◽  
R. Carbonell ◽  
J. F. Simancas
Keyword(s):  
Seismic Reflection ◽  
Fold Belt ◽  
The South ◽  
Prestack Migration ◽  
Time Migration ◽  
Straight Line ◽  
Prestack Time Migration ◽  
2D Imaging ◽  
The Common ◽  
Zero Offset

Crooked-line 2D seismic reflection survey geometries violate underlying assumptions of 2D imaging routines, affecting our ability to resolve the subsurface reliably. We compare three crooked-line imaging schemes involving prestack and poststack time migration using the 2D IBERSEIS deep seismic reflection profile running over the South Portuguese Zone thrust-and-fold belt to obtain crisp high-resolution images of the shallow crust. The crust is characterized by a complex subsurface geometry with conflicting dips of up to [Formula: see text]. In summary, the three schemes are (1) normal-moveout (NMO) corrections, dip-moveout (DMO) corrections, common-midpoint (CMP) stacking, CMP projection, and poststack time migration; (2) NMO corrections, DMO corrections, CMP projection, zero-offset time migration of the common-offset gathers, and CMP stacking; (3) CMP projection, prestack time migration in the common-offset domain, and CMP stacking. An essential element of all three schemes is a CMP projection routine, projecting the CMPs first binned along individual segments for preprocessing onto one straight line, which is parallel to the general dip direction of the subsurface structures. After CMP projection, the data satisfy the straight-line assumption of 2D imaging routines more closely. We observe that the prestack time-migration scheme yields comparable or more coherent synthetic and field-data images than the other two DMO-based schemes along the parts of the profile where the acquisition overall follows a straight line. However, the schemes involving DMO corrections are less plagued by migration artifacts than the prestack time-migration scheme along profile parts where the acquisition line is crooked. In particular, prominent migration artifacts on the prestack migrated synthetic data can be related to significant variations in source-receiver azimuths for which 2D prestack migration cannot account. Thus, the processing scheme including DMO corrections, CMP projection, and zero-offset migration of common-offset gathers offers a reliable and effective alternative to prestack migration for crooked-line 2D seismic reflection processing.

Prestack time migration by common-migrated-reflector-element stacking

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. S73-S82 ◽  
Author(s):  
Sergius Dell ◽  
Dirk Gajewski ◽  
Claudia Vanelle
Keyword(s):  
Initial Time ◽  
Depth Migration ◽  
Time Migration ◽  
Migration Method ◽  
Prestack Time Migration ◽  
The Common ◽  
A New Technique ◽  
New Formulation ◽  
Time Image ◽  
Migration Velocities

Time migration is an attractive tool to produce a subsurface image because it is faster and less sensitive to velocities errors than depth migration. However, a highly focused time image is only achievable with well-determined time-migration velocities. Therefore, a refinement of the initial time-migration velocities often is required. We introduced a new technique for prestack time migration, based on the common-migrated-reflector-element stack of common scatterpoint gathers, including an automatic update of time-migration velocities. The common scatterpoint gathers are generated using a new formulation of the double-square-root equation that is parametrized with the common-offset apex time. The common-migrated-reflector-element stack is a multiparameter stacking technique based on the Taylor expansion of traveltimes of time-migrated reflections in the paraxial vicinity of the image ray. Our 2D synthetic and field data examples demonstrated that the proposed method provides updated time-migration velocities that are more robust and have higher resolution compared with the initial time-migration velocities. The prestack time migration method also showed a clear improvement of the focusing of reflections for such geologic features as faults and salt structures.

scholarly journals Shot-gather time migration of planar reflectors without velocity model

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. S93-S101 ◽  
Author(s):  
Andrej Bóna
Keyword(s):  
Synthetic Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Time Migration ◽  
Second Derivatives ◽  
Migration Method ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
2D And 3D ◽  
Homogeneous Media

Standard migration techniques require a velocity model. A new and fast prestack time migration method is presented that does not require a velocity model as an input. The only input is a shot gather, unlike other velocity-independent migrations that also require input of data in other gathers. The output of the presented migration is a time-migrated image and the migration velocity model. The method uses the first and second derivatives of the traveltimes with respect to the location of the receiver. These attributes are estimated by computing the gradient of the amplitude in a shot gather. The assumptions of the approach are a laterally slowly changing velocity and reflectors with small curvatures; the dip of the reflector can be arbitrary. The migration velocity corresponds to the root mean square (rms) velocity for laterally homogeneous media for near offsets. The migration expressions for 2D and 3D cases are derived from a simple geometrical construction considering the image of the source. The strengths and weaknesses of the methods are demonstrated on synthetic data. At last, the applicability of the method is discussed by interpreting the migration velocity in terms of the Taylor expansion of the traveltime around the zero offset.

scholarly journals Deabsorption prestack time migration from rugged topography

2021 ◽  
Vol 18 (2) ◽  
pp. 291-303
Author(s):  
Changshan Han ◽  
Linong Liu ◽  
Zelin Liu ◽  
Zhengwei Li
Keyword(s):  
Signal To Noise Ratio ◽  
Three Dimensional ◽  
Signal To Noise ◽  
Time Migration ◽  
Q Model ◽  
Rugged Topography ◽  
Migration Method ◽  
Prestack Time Migration ◽  
The Common ◽  
Quality Factor Q

Abstract We developed a modified topography prestack time migration (PSTM) scheme that can improve the imaging resolution by applying effective Q to topography migration. The computation of the traveltime at each imaging location in the migration is based on the floating datum smoothed by rugged topography. Unlike the common quality factor Q, the effective Q only determines the frequency-dependent amplitude and the traveltime at a single imaging location, which enables us to establish a Q model in an inhomogeneous medium. Hence, we can acquire the effective Q using a scanning technology according to the width of the frequency band and signal-to-noise ratio of the imaging gathers. The proposed migration method can be integrated into the conventional topography migration workflow. Synthetic and three-dimensional (3D) field datasets indicate that the proposed deabsorption PSTM from rugged topography is effective.

Velocity continuation and the anatomy of residual prestack time migration

Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1650-1661 ◽  
Author(s):  
Sergey Fomel
Keyword(s):  
Continuous Process ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Image Transformation ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Seismic Image ◽  
Spectral Algorithms ◽  
Prestack Time Migration ◽  
Zero Offset

Velocity continuation is an imaginary continuous process of seismic image transformation in the postmigration domain. It generalizes the concepts of residual and cascaded migrations. Understanding the laws of velocity continuation is crucially important for a successful application of time‐migration velocity analysis. These laws predict the changes in the geometry and intensity of reflection events on migrated images with the change of the migration velocity. In this paper, I derive kinematic and dynamic laws for the case of prestack residual migration from simple geometric principles. The main theoretical result is a decomposition of prestack velocity continuation into three different components corresponding to residual normal moveout, residual dip moveout, and residual zero‐offset migration. I analyze the contribution and properties of each of the three components separately. This theory forms the basis for constructing efficient finite‐difference and spectral algorithms for time‐migration velocity analysis.

Migration velocity analysis with the Kirchhoff integral

Geophysics ◽  
1991 ◽  
Vol 56 (3) ◽  
pp. 365-370 ◽  
Author(s):  
Y. C. Kim ◽  
R. Gonzalez
Keyword(s):  
Maximum Amplitude ◽  
Migration Velocity ◽  
Downward Continuation ◽  
Velocity Analysis ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
Kirchhoff Integral ◽  
Migration Velocities

To obtain accurate migration velocities, we must estimate the velocity at migrated depth points. Wavefront focusing analysis with downward continuation yields the rms velocity at migrated depth points; however, the large amount of computation required for downward continuation limits use of this approach for routine processing. The purpose of this paper is to present an implementation of the Kirchhoff integral which makes the wavefront focusing analysis practical for time‐migration velocity analysis. Downward continuation focuses the wavefront to the zero offset at the depth controlled by the velocity used for the continuation. The migration velocity is then determined from the depth where the focused wavefront attains the maximum amplitude. The flexibility of the Kirchhoff integral allows us to compute only the zero‐offset trace at each depth point and lets us avoid most of the computation for the downward continuation of unstacked data. Furthermore, since the velocity is obtained from the location where the focused wavefront shows the maximum amplitude, prestack time migration with the velocity from this technique produces the maximum amplitude for the subsurface reflector.

Oriented prestack time migration using local slopes and predictive painting in the common-source domain for planar reflectors

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. S409-S418 ◽  
Author(s):  
M. Javad Khoshnavaz ◽  
Andrej Bóna ◽  
Aleksander Dzunic ◽  
Kevin Ung ◽  
Milovan Urosevic
Keyword(s):  
Seismic Data ◽  
Imaging Techniques ◽  
Velocity Model ◽  
Higher Order ◽  
Common Source ◽  
Velocity Analysis ◽  
Time Migration ◽  
Prestack Time Migration ◽  
The Common ◽  
Migration Technique

Seismic imaging techniques often require an input velocity model. Velocity analysis is one of the most critical stages in seismic data processing. Standard ways to find the velocity model from seismic data in the time domain are constant velocity stack and semblance velocity analysis that may be time consuming and labor intensive. Oriented/velocity-less imaging using local event slopes is an alternative to the conventional imaging techniques. In some previous oriented techniques, seismic data must be sorted in two different domains, whereas seismic data are not always available in both domains and the use of interpolation is inevitable in such cases. Other methods are developed in terms of the higher order derivatives of traveltime with respect to offset, whereas estimation of the higher order derivatives is difficult to achieve with the required accuracy. We addressed the limitations by developing an oriented local slope based prestack time migration technique in only one domain: the common-source domain. The migration technique is developed for reflectors with small curvature. In the proposed approach, the need for the estimation of higher order derivatives is replaced by a point-to-point mapping of seismic data using the predictive painting technique. The theoretical contents of the proposed technique are tested on a simple synthetic data example and applied to a field data set.

Common‐reflection‐surface stack: Image and attributes

Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 97-109 ◽  
Author(s):  
Rainer Jäger ◽  
Jürgen Mann ◽  
German Höcht ◽  
Peter Hubral
Keyword(s):  
Optimal Parameter ◽  
Normal Wave ◽  
Normal Incidence ◽  
Prestack Migration ◽  
Computational Costs ◽  
Reflection Data ◽  
Common Reflection Surface ◽  
The Common ◽  
Zero Offset ◽  
Independent Parameters

The common‐reflection‐surface stack provides a zero‐offset simulation from seismic multicoverage reflection data. Whereas conventional reflection imaging methods (e.g. the NMO/dip moveout/stack or prestack migration) require a sufficiently accurate macrovelocity model to yield appropriate results, the common‐reflection‐surface (CRS) stack does not depend on a macrovelocity model. We apply the CRS stack to a 2-D synthetic seismic multicoverage dataset. We show that it not only provides a high‐quality simulated zero‐offset section but also three important kinematic wavefield attribute sections, which can be used to derive the 2-D macrovelocity model. We compare the multicoverage‐data‐derived attributes with the model‐derived attributes computed by forward modeling. We thus confirm the validity of the theory and of the data‐derived attributes. For 2-D acquisition, the CRS stack leads to a stacking surface depending on three search parameters. The optimum stacking surface needs to be determined for each point of the simulated zero‐offset section. For a given primary reflection, these are the emergence angle α of the zero‐offset ray, as well as two radii of wavefront curvatures [Formula: see text] and [Formula: see text]. They all are associated with two hypothetical waves: the so‐called normal wave and the normal‐incidence‐point wave. We also address the problem of determining an optimal parameter triplet (α, [Formula: see text], [Formula: see text]) in order to construct the sample value (i.e., the CRS stack value) for each point in the desired simulated zero‐offset section. This optimal triplet is expected to determine for each point the best stacking surface that can be fitted to the multicoverage primary reflection events. To make the CRS stack attractive in terms of computational costs, a suitable strategy is described to determine the optimal parameter triplets for all points of the simulated zero‐offset section. For the implementation of the CRS stack, we make use of the hyperbolic second‐order Taylor expansion of the stacking surface. This representation is not only suitable to handle irregular multicoverage acquisition geometries but also enables us to introduce simple and efficient search strategies for the parameter triple. In specific subsets of the multicoverage data (e.g., in the common‐midpoint gathers or the zero‐offset section), the chosen representation only depends on one or two independent parameters, respectively.

scholarly journals Zero-offset sections with a deblurring filter in the time domain

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. S239-S249
Author(s):  
Shihang Feng ◽  
Oz Yilmaz ◽  
Yuqing Chen ◽  
Gerard T. Schuster
Keyword(s):  
Time Domain ◽  
Constant Velocity ◽  
Field Data ◽  
Velocity Model ◽  
Velocity Models ◽  
Time Migration ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
Image Volume ◽  
The Time Domain

The conventional common-midpoint stack is not equivalent to the zero-offset section due to the existence of velocity uncertainty. To obtain a zero-offset reflection section that preserves most reflections and diffractions, we have developed a velocity-independent workflow for reconstructing a high-quality zero-offset reflection section from prestack data with a deblurring filter. This workflow constructs a migration image volume by prestack time migration using a series of constant-velocity models. A deblurring filter for each constant-velocity model is applied to each time-migration image to get a deblurred image volume. To preserve all events in the image volume, each deblurred image panel is demigrated and then summed over the velocity axis. Compared with the workflow without a deblurring filter, the composite zero-offset reflection section has higher resolution and fewer migration artifacts. We evaluate applications of our method to synthetic and field data to validate its effectiveness.

A two‐pass approximation to 3-D prestack migration

Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 409-421 ◽  
Author(s):  
Anat Canning ◽  
Gerald H. F. Gardner
Keyword(s):  
Computation Time ◽  
Three Dimensions ◽  
Velocity Analysis ◽  
Depth Migration ◽  
Prestack Migration ◽  
Time Migration ◽  
Data Volume ◽  
Land Survey ◽  
Zero Offset ◽  
The One

A two‐pass approximation to 3-D Kirchhoff migration simplifies the migration procedure by reducing it to a succession of 2-D operations. This approach has proven very successful in the zero‐offset case. A two‐pass approximation to 3-D migration is described here for the prestack case. Compared to the one‐pass approach, the scheme presented here provides significant reduction in computation time and a relatively simple data manipulation scheme. The two‐pass method was designed using velocity independent prestack time migration (DMO‐PSI) applied in the crossline direction, followed by conventional prestack depth migration in the inline direction. Velocity analysis, an important part of prestack migration, is also included in the two‐pass scheme. It is carried out as a 2-D procedure after 3-D effects are removed from the data volume. The procedure presented here is a practical full volume 3-D prestack migration. One of its main benefits is a realistic and efficient iterative velocity analysis procedure in three dimensions. The algorithm was designed in the frequency domain and the computational scheme was optimized by processing individual frequency slices independently. Irregular trace distribution, a feature that characterizes most 3-D seismic surveys, is implicitly accounted for within the two‐pass algorithm. A numerical example tests the performance of the two‐pass 3-D prestack migration program in the presence of a vertical velocity gradient. A 3-D land survey from a fold and thrust belt region was used to demonstrate the algorithm in a complex geological setting. The results were compared with images from other 2-D and 3-D migration schemes and show improved resolution and higher signal content.

The equivalent offset method of prestack time migration

Geophysics ◽  
1998 ◽  
Vol 63 (6) ◽  
pp. 2042-2053 ◽  
Author(s):  
John C. Bancroft ◽  
Hugh D. Geiger ◽  
Gary F. Margrave
Keyword(s):  
Velocity Analysis ◽  
Square Root ◽  
Prestack Migration ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Scattered Energy ◽  
Nmo Correction ◽  
Velocity Information ◽  
Prestack Time Migration ◽  
Receiver Input

A prestack time migration is presented that is simple, efficient, and provides detailed velocity information. It is based on Kirchhoff prestack time migration and can be applied to both 2-D and 3-D data. The method is divided into two steps: the first is a gathering process that forms common scatterpoint (CSP) gathers; the second is a focusing process that applies a simplified Kirchhoff migration on the CSP gathers, and consists of scaling, filtering, normal moveout (NMO) correction, and stacking. A key concept of the method is a reformulation of the double square‐root equation (of source‐scatterpoint‐receiver traveltimes) into a single square root. The single square root uses an equivalent offset that is the surface distance from the scatterpoint to a colocated source and receiver. Input samples are mapped into offset bins of a CSP gather, without time shifting, to an offset defined by the equivalent offset. The single square‐root reformulation gathers scattered energy to hyperbolic paths on the appropriate CSP gathers. A CSP gather is similar to a common midpoint (CMP) gather as both are focused by NMO and stacking. However, the CSP stack is a complete Kirchhoff prestack migrated section, whereas the CMP stack still requires poststack migration. In addition, the CSP gather has higher fold in the offset bins and a much larger offset range due to the gathering of all input traces within the migration aperture. The new method gains computational efficiency by delaying the Kirchhoff computations until after the CSP gather has been formed. The high fold and large offsets of the CSP gather enables precise focusing of the velocity semblance and accurate velocity analysis. Our algorithm is formulated in the space‐time domain, which enables prestack migration velocity analysis to be performed at selected locations and permits prestack migration of a 3-D volume into an arbitrarily located 2-D line.

Sign in / Sign up

Export Citation Format

Share Document