Ultrashallow depth imaging of a channel stratigraphy with first-arrival traveltime inversion and prestack depth migration: A case history from Bull Creek, Oklahoma

Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. B87-B96 ◽  
Author(s):  
Ammanuel Fesseha Woldearegay ◽  
Priyank Jaiswal ◽  
Alexander R. Simms ◽  
Hanna Alexander ◽  
Leland C. Bement ◽  
...  

Depth imaging in ultrashallow ([Formula: see text]) environments presents twofold challenge: (1) coda available for depth migration is very limited; and (2) conventional time processing with limited coda generally fails to estimate reliable velocity models for depth migration. We studied the combining of first-arrival traveltime inversion and prestack depth migration (PSDM) for depth imaging of ultrashallow paleochannel stratigraphy associated with the Bull Creek drainage system, Oklahoma. Restricted by a limited number of geophones (24) we acquired data for inversion and migration through two coincident profiles. The first profile for inversion has a wider survey-aperture (115-m maximum shot-receiver spacing) and consequently sparse CMP spacing (2.5 m), whereas the second profile for PSDM has denser CMP spacing (1 m) and consequently a narrower survey aperture (46-m maximum shot-receiver spacing). We also found that the velocity model from traveltime inversion of the wider-aperture data set is more preferable for depth-migration than the velocity model from time processing of the denser data set. The preferred depth image showed three episodes of incision whose chronological order is resolved through radio-carbon dating of terrace sediments. Results suggested that even with limited geophones, depth imaging of ultrashallow targets can be achieved by combining first-arrival traveltime inversion and PSDM through coincident wide- and narrow-aperture acquisitions.

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE269-VE280 ◽  
Author(s):  
Priyank Jaiswal ◽  
Colin A. Zelt

Imaging 2D multichannel land seismic data can be accomplished effectively by a combination of traveltime inversion and prestack depth migration (PSDM), referred to as unified imaging. Unified imaging begins by inverting the direct-arrival times to estimate a velocity model that is used in static corrections and stacking velocity analysis. The interval velocity model (from stacking velocities) is used for PSDM. The stacked data and the PSDM image are interpreted for common horizons, and the corresponding wide-aperture reflections are identified in the shot gathers. Using the interval velocity model, the stack interpretations are inverted as zero-offset reflections to constrain the corresponding interfaces in depth; the interval velocity model remains stationary. We define a coefficient of congruence [Formula: see text] that measures the discrepancy between horizons from the PSDM image andtheir counterparts from the zero-offset inversion. A value of unity for [Formula: see text] implies that the interpreted and inverted horizons are consistent to within the interpretational uncertainties, and the unified imaging is said to have converged. For [Formula: see text] greater than unity, the interval velocity model and the horizon depths are updated by jointly inverting the direct arrivals with the zero-offset and wide-aperture reflections. The updated interval velocity model is used again for both PSDM and a zero-offset inversion. Interpretations of the new PSDM image are the updated horizon depths. The unified imaging is applied to seismic data from the Naga Thrust and Fold Belt in India. Wide-aperture and zero-offset data from three geologically significant horizons are used. Three runs of joint inversion and PSDM are required in a cyclic manner for [Formula: see text] to converge to unity. A joint interpretation of the final velocity model and depth image reveals the presence of a triangle zone that could be promising for exploration.


Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 568-576 ◽  
Author(s):  
Young C. Kim ◽  
Worth B. Hurt, ◽  
Louis J. Maher ◽  
Patrick J. Starich

The transformation of surface seismic data into a subsurface image can be separated into two components—focusing and positioning. Focusing is associated with ensuring the data from different offsets are contributing constructively to the same event. Positioning involves the transformation of the focused events into a depth image consistent with a given velocity model. In prestack depth migration, both of these operations are achieved simultaneously; however, for 3-D data, the cost is significant. Prestack time migration is much more economical and focuses events well even in the presence of moderate velocity variations, but suffers from mispositioning problems. Hybrid migration is a cost‐effective depth‐imaging approach that uses prestack time migration for focusing; inverse migration for the removal of positioning errors; and poststack depth migration for proper positioning. When lateral velocity changes are moderate, the hybrid technique can generate a depth image that is consistent with a velocity field. For very complex structures that require prestack depth migration, the results of the hybrid technique can be used to create a starting velocity model, thereby reducing the number of iterations for velocity model building.


Geophysics ◽  
1993 ◽  
Vol 58 (6) ◽  
pp. 873-882 ◽  
Author(s):  
Roelof Jan Versteeg

To get a correct earth image from seismic data acquired over complex structures it is essential to use prestack depth migration. A necessary condition for obtaining a correct image is that the prestack depth migration is done with an accurate velocity model. In cases where we need to use prestack depth migration determination of such a model using conventional methods does not give satisfactory results. Thus, new iterative methods for velocity model determination have been developed. The convergence of these methods can be accelerated by defining constraints on the model in such a way that the method only looks for those components of the true earth velocity field that influence the migrated image. In order to determine these components, the sensitivity of the prestack depth migration result to the velocity model is examined using a complex synthetic data set (the Marmousi data set) for which the exact model is known. The images obtained with increasingly smoothed versions of the true model are compared, and it is shown that the minimal spatial wavelength that needs to be in the model to obtain an accurate depth image from the data set is of the order of 200 m. The model space that has to be examined to find an accurate velocity model from complex seismic data can thus be constrained. This will increase the speed and probability of convergence of iterative velocity model determination methods.


Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1903-1912 ◽  
Author(s):  
Igor B. Morozov ◽  
Alan Levander

Wide‐aperture, prestack depth migration requires application of challenging and time‐consuming velocity analysis and depth focusing, collectively referred to here as depth focusing. We present an approach to depth focusing using (1) a detailed starting velocity model obtained by a 1‐D transformation of the first‐arrival times, followed iteratively by (2) interactive analysis of the common‐image gathers, (3) computation of coherency attributes of the wavefield in the depth domain, and (4) 2‐D traveltime tomography to update the background velocity model. We employ two interactive method of migration velocities refinement. In the first method (similar to the common‐midpoint velocity spectrum approach), residual velocity updates are picked directly from the common‐image gathers. In another method (analogous to the common velocity stacks), we pick the velocity updates from the areas of maximum coherency in depth sections that are migrated using rescaled traveltime maps. Both types of migration velocity picks, optionally combined with the first arrivals, are inputs for a 2‐D traveltime inversion scheme that uses either the infinite‐frequency or a finite‐bandwidth approximation. This flexible and versatile depth focusing approach is implemented for several prestack depth migration algorithms and illustrated on an application to a real, ultrashallow seismic data set. The technique resolves overburden velocity variations and facilitates reliable high‐resolution reflection imaging of a paleochannel that was the target of the study.


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 ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE13-VE23 ◽  
Author(s):  
Frank Adler ◽  
Reda Baina ◽  
Mohamed Amine Soudani ◽  
Pierre Cardon ◽  
Jean-Baptiste Richard

Velocity-model estimation with seismic reflection tomography is a nonlinear inverse problem. We present a new method for solving the nonlinear tomographic inverse problem using 3D prestack-depth-migrated reflections as the input data, i.e., our method requires that prestack depth migration (PSDM) be performed before tomography. The method is applicable to any type of seismic data acquisition that permits seismic imaging with Kirchhoff PSDM. A fundamental concept of the method is that we dissociate the possibly incorrect initial migration velocity model from the tomographic velocity model. We take the initial migration velocity model and the residual moveout in the associated PSDM common-image gathers as the reference data. This allows us to consider the migrated depth of the initial PSDM as the invariant observation for the tomographic inverse problem. We can therefore formulate the inverse problem within the general framework of inverse theory as a nonlinear least-squares data fitting between observed and modeled migrated depth. The modeled migrated depth is calculated by ray tracing in the tomographic model, followed by a finite-offset map migration in the initial migration model. The inverse problem is solved iteratively with a Gauss-Newton algorithm. We applied the method to a North Sea data set to build an anisotropic layer velocity model.


Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 594-603 ◽  
Author(s):  
Robert J. Ferguson ◽  
Gary F. Margrave

A new depth migration method suitable for heterogeneous media is presented. The well‐known phase shift plus interpolation (PSPI) method and the recently introduced nonstationary phase‐shift (NSPS) method are combined into a single symmetric operator with improved accuracy and stability and with similar computational effort. For prestack depth migration, the symmetric operator is used in a recursive wavefield extrapolation to compute incident and reflected wavefields at any desired depth, and the ratio of the incident and reflected wavefields at a particular depth is used to estimate seismic reflectivity. When the velocity model is made piecewise constant laterally, the symmetric extrapolation operator can be computed efficiently using ordinary phase‐shift extrapolation for a series of reference velocities and appropriate spatial windowing. Migration of the Marmousi synthetic data set by symmetric nonstationary phase shift (SNPS) provides an image that compares favorably with an image of the zero‐offset reflectivity derived from the Marmousi velocity model.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. S185-S197 ◽  
Author(s):  
Bertrand Duquet ◽  
Patrick Lailly

Full-volume seismic imaging is essential for a sound interpretation of structurally complex geologies. Prestack depth imaging is the most appropriate tool for such imaging, but it requires a precise and often complex velocity model. In such situations, 3D Kirchhoff prestack depth migration can be quite expensive. On the other hand, a wavefield approach, although generally tremendously expensive, is not affected by the complexity of the velocity model. We propose an affordable 3-D wavefield prestack depth-migration technique. It is designed for marine surveys for which the source-receiver azimuth is approximately constant. The technique applies a plane-wave migration algorithm to time-shifted data — quite a surprising approach when we realize that marine surveys do not allow the synthesis of genuine plane-wave data. Additionally, the imaging principle has to be modified to give results consistent with shot-record migration. Our technique also produces image gathers that allow an update of the velocity model by means of migration velocity analysis. Results from synthetics and conventional marine data demonstrate the effectiveness of the method.


2021 ◽  
Vol 225 (2) ◽  
pp. 1020-1031
Author(s):  
Huachen Yang ◽  
Jianzhong Zhang ◽  
Kai Ren ◽  
Changbo Wang

SUMMARY A non-iterative first-arrival traveltime inversion method (NFTI) is proposed for building smooth velocity models using seismic diving waves observed on irregular surface. The new ray and traveltime equations of diving waves propagating in smooth media with undulant observation surface are deduced. According to the proposed ray and traveltime equations, an analytical formula for determining the location of the diving-wave turning points is then derived. Taking the influence of rough topography on first-arrival traveltimes into account, the new equations for calculating the velocities at turning points are established. Based on these equations, a method is proposed to construct subsurface velocity models from the observation surface downward to the bottom using the first-arrival traveltimes in common offset gathers. Tests on smooth velocity models with rugged topography verify the validity of the established equations, and the superiority of the proposed NFTI. The limitation of the proposed method is shown by an abruptly-varying velocity model example. Finally, the NFTI is applied to solve the static correction problem of the field seismic data acquired in a mountain area in the western China. The results confirm the effectivity of the proposed NFTI.


Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1226-1237 ◽  
Author(s):  
Irina Apostoiu‐Marin ◽  
Andreas Ehinger

Prestack depth migration can be used in the velocity model estimation process if one succeeds in interpreting depth events obtained with erroneous velocity models. The interpretational difficulty arises from the fact that migration with erroneous velocity does not yield the geologically correct reflector geometries and that individual migrated images suffer from poor signal‐to‐noise ratio. Moreover, migrated events may be of considerable complexity and thus hard to identify. In this paper, we examine the influence of wrong velocity models on the output of prestack depth migration in the case of straight reflector and point diffractor data in homogeneous media. To avoid obscuring migration results by artifacts (“smiles”), we use a geometrical technique for modeling and migration yielding a point‐to‐point map from time‐domain data to depth‐domain data. We discover that strong deformation of migrated events may occur even in situations of simple structures and small velocity errors. From a kinematical point of view, we compare the results of common‐shot and common‐offset migration. and we find that common‐offset migration with erroneous velocity models yields less severe image distortion than common‐shot migration. However, for any kind of migration, it is important to use the entire cube of migrated data to consistently interpret in the prestack depth‐migrated domain.


Sign in / Sign up

Export Citation Format

Share Document