True‐amplitude imaging and dip moveout

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 47-66 ◽  
Author(s):  
James L. Black ◽  
Karl L. Schleicher ◽  
Lin Zhang

True‐amplitude seismic imaging produces a three dimensional (3-D) migrated section in which the peak amplitude of each migrated event is proportional to the reflectivity. For a constant‐velocity medium, the standard imaging sequence consisting of spherical‐divergence correction, normal moveout (NMO), dip moveout (DMO), and zero‐offset migration produces a true‐amplitude image if the DMO step is done correctly. There are two equivalent ways to derive the correct amplitude‐preserving DMO. The first is to improve upon Hale’s derivation of F-K DMO by taking the reflection‐point smear properly into account. This yields a new Jacobian that simply replaces the Jacobian in Hale’s method. The second way is to calibrate the filter that appears in integral DMO so as to preserve the amplitude of an arbitrary 3-D dipping reflector. This latter method is based upon the 3-D acoustic wave equation with constant velocity. The resulting filter amounts to a simple modification of existing integral algorithms. The new F-K and integral DMO algorithms resulting from these two approaches turn out to be equivalent, producing identical outputs when implemented in nonaliased fashion. As dip increases, their output become progressively larger than the outputs of either Hale’s F-K method or the integral method generally associated with Deregowski and Rocca. This trend can be observed both on model data and field data. There are two additional results of this analysis, both following from the wave‐equation calibration on an arbitrary 3-D dipping reflector. The first is a proof that the entire imaging sequence (not just the DMO part) is true‐amplitude when the DMO is done correctly. The second result is a handy formula showing exactly how the zero‐phase wavelet on the final migrated image is a stretched version of the zero‐phase deconvolved source wavelet. This result quantitatively expresses the loss of vertical resolution due to dip and offset.

Geophysics ◽  
1990 ◽  
Vol 55 (1) ◽  
pp. 10-19 ◽  
Author(s):  
Martin Karrenbach

Three‐dimensional migration of zero‐offset data using a velocity varying with depth can be performed in one pass using Fourier transforms of time slices. The migration process is carried out entirely in the two‐dimensional spatial Fourier domain. The algorithm consecutively filters and adds time slices of the 3-D data volume in a way that is equivalent to summing energy over the diffraction surface of a point scatterer. The partial energy being distributed along a circle in a time slice is properly added in each summation step. Time‐slice migration is based on an integral solution of the acoustic wave equation known as the “Kirchhoff integral.” The wavelet shape in a 3-D data volume is preserved throughout the entire migration process. The frequency characteristics are maintained by summing weighted differences between time slices instead of summing the time slices themselves. Automatic weighting is achieved by time slicing at equal increments of diffraction radius. Tapering the summation operator reduces effects introduced by limiting the summation window. Time‐slice migration preserves the frequency content of a 3-D data volume during summation in a natural way. Since the migration scheme assumes a constant velocity within the entire time slice, it is a local process in time which migrates a 3-D data volume with a constant velocity or with a velocity which varies with depth. The migration algorithm is applied to numerical and physical model data. This method is especially suitable for a migration of a targeted subset of the 3-D data volume.


Geophysics ◽  
1993 ◽  
Vol 58 (7) ◽  
pp. 1030-1041 ◽  
Author(s):  
Hans A. Meinardus ◽  
Karl L. Schleicher

The standard seismic imaging sequence consists of normal moveout (NMO), dip moveout (DMO), stack, and zero‐offset migration. Conventional NMO and DMO processes remove much of the effect of offset from prestack data, but the constant velocity assumption in most DMO algorithms can compromise the ultimate results. Time‐variant DMO avoids the constant velocity assumption to create better stacks, especially for steeply dipping events. Time‐variant DMO can be implemented as a 3-D, f-k domain process using the dip decomposition method. Prestack data are moved out with a set of NMO velocities corresponding to discrete values of in‐line and crossline dips. The dip‐dependent NMO velocity is computed to remove the trace offset and azimuth dependence of event times for an arbitrary velocity function of depth. After stacking the moved out CMP gathers, a three‐dimensional (3-D) dip filter is applied to select the particular in‐line and crossline dip. The final zero‐offset image is obtained by summing all the dip‐filtered sections. This process generates a saddle‐shaped 3-D impulse response for a constant velocity gradient. The impulse response is more complicated for a general depth‐variable velocity function, where the response exhibits secondary branches, or triplications, at steeper dips. These complicated impulse responses, including amplitude and phase effects, are implicitly produced by the f-k process. The dip‐decomposition method of 3-D time‐variant DMO is an efficient and accurate process to correct for the effect of offset in the presence of an arbitrary velocity variation with depth. The impulse response of this process implicitly contains complex features like a 3-D saddle shape, triplications, amplitude, and phase. Field data from the Gulf of Mexico shows significant improvement on a steep salt flank event.


Geophysics ◽  
1981 ◽  
Vol 46 (8) ◽  
pp. 1116-1120 ◽  
Author(s):  
A. B. Weglein ◽  
W. E. Boyse ◽  
J. E. Anderson

We present a formalism for obtaining the subsurface velocity configuration directly from reflection seismic data. Our approach is to apply the results obtained for inverse problems in quantum scattering theory to the reflection seismic problem. In particular, we extend the results of Moses (1956) for inverse quantum scattering and Razavy (1975) for the one‐dimensional (1-D) identification of the acoustic wave equation to the problem of identifying the velocity in the three‐dimensional (3-D) acoustic wave equation from boundary value measurements. No a priori knowledge of the subsurface velocity is assumed and all refraction, diffraction, and multiple reflection phenomena are taken into account. In addition, we explain how the idea of slant stack in processing seismic data is an important part of the proposed 3-D inverse scattering formalism.


2001 ◽  
Vol 12 (10) ◽  
pp. 1497-1507 ◽  
Author(s):  
ERICH STOLL ◽  
STEFAN DANGEL

The acoustic wave equation is solved numerically for two and three-dimensional systems at the limit between near and far field propagation. Our results show that for large sound velocities, corresponding to wavelengths larger than the system, near field properties are dominant. When the near field conditions are no longer satisfied, standing waves close to the sound emitters and interference patterns between the near field and far field solutions appear. Our procedure is applied to sound sources, which broadcast coherent and continuous waves as well as to sources emitting bursts of incoherent and uncorrelated waves. Both cases can be used to simulate the spreading of low frequency seismic waves observed close to volcanoes and hydrocarbon reservoirs.


2015 ◽  
Vol 12 (3) ◽  
pp. 409-420 ◽  
Author(s):  
Xiao-Hui Cai ◽  
Yang Liu ◽  
Zhi-Ming Ren ◽  
Jian-Min Wang ◽  
Zhi-De Chen ◽  
...  

Geophysics ◽  
1988 ◽  
Vol 53 (5) ◽  
pp. 604-610 ◽  
Author(s):  
David Forel ◽  
Gerald H. F. Gardner

Prestack migration in a constant‐velocity medium spreads an impulse on any trace over an ellipsoidal surface with foci at the source and receiver positions for that trace. The same ellipsoid can be obtained by migrating a family of zero‐offset traces placed along the line segment from the source to the receiver. The spheres generated by migrating the zero‐offset impulses are arranged to be tangent to the ellipsoid. The resulting nonstandard moveout equation is equivalent to two consecutive moveouts, the first requiring no knowledge of velocity and the second being standard normal moveout (NMO). The first of these is referred to as dip moveout (DMO). Because this DMO-NMO algorithm converts any trace to an equivalent set of zero‐offset traces, it can be applied to any ensemble of traces no matter what the variations in azimuth and offset may be. In particular, this three‐dimensional perspective on DMO can be used with multifold inline data. Then it becomes clear that velocity‐independent DMO operates on radial‐trace profiles and not on constant‐offset profiles. Inline data over a three‐dimensional subsurface will be properly stacked by using DMO followed by NMO.


Geophysics ◽  
1988 ◽  
Vol 53 (9) ◽  
pp. 1194-1201 ◽  
Author(s):  
Jing Wen ◽  
George A. McMechan ◽  
Michael W. Booth

Programs for 3-D modeling and migration, using 3-D Fourier transforms to solve the scalar wave equation in frequency‐wavenumber space, are developed, implemented, tested, and applied to synthetic and scale‐model data. With microtasking to fully use four CRAY processors in parallel, we can solve a complete [Formula: see text] modeling problem in about 2.5 minutes (elapsed time); of this time, the two 3-D Fourier transforms take 1 minute and the wave‐equation calculations take 1.5 minutes. The corresponding migration also takes 2.5 minutes. Thus, even iterative 3-D processing is now feasible. The two main assumptions in our algorithm are that the earth has a constant velocity and that the data are zero‐offset (or stacked). Tests with model data verify that the algorithm produces the correct results when these assumptions are satisfied. Tests with scale‐model data show that approximate images may still be obtained when the assumptions are not strictly met; but the images contain a variety of distortions, primarily related to undermigration and overmigration, so caution is required in interpretation.


Sign in / Sign up

Export Citation Format

Share Document