Analyzing operators of 3-D imaging software with impulse responses

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1079-1092 ◽  
Author(s):  
William A. Schneider

No processing step changes seismic data more than 3-D imaging. Imaging techniques such as 3-D migration and dip moveout (DMO) generally change the position, amplitude, and phase of reflections as they are converted into reflector images. Migration and DMO may be formulated in many different ways, and various algorithms are available for implementing each formulation. These algorithms all make physical approximations, causing imaging software to vary with algorithm choice. Imaging software also varies because of additional implementation approximations, such as those that trade accuracy for efficiency. Imaging fidelity, then, generally depends upon algorithm, implementation, specific software parameters (such as aperture, antialias filter settings, and downward‐continuation step size), specific acquisition parameters (such as nominal x- and y-direction trace spacings and wavelet frequency range), and, of course, the velocity model. Successfully imaging the target usually requires using appropriate imaging software, parameters, and velocities. Impulse responses provide an easy way to quantitatively understand the operators of imaging software and then predict how specific imaging software will perform with the chosen parameters. (An impulse response is the image computed from a data set containing only one nonzero trace and one arrival on that trace.) I have developed equations for true‐amplitude impulse responses of 3-D prestack time migration, 3-D zero‐offset time migration, 3-D exploding‐reflector time migration, and DMO. I use these theoretical impulse responses to analyze the operators of actual imaging software for a given choice of software parameters, acquisition parameters, and velocity model. The procedure is simple: compute impulse responses of some software; estimate position, amplitude, and phase of the impulse‐response events; and plot these against the theoretical values. The method is easy to use and has proven beneficial for analyzing general imaging software and for parameter evaluation with specific imaging software.

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 ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. S71-S82 ◽  
Author(s):  
Chris Willacy ◽  
Maksym Kryvohuz

The imaging of steep salt boundaries has received much attention with the advent of improved wider azimuth acquisition designs and advanced imaging techniques such as reverse time migration (RTM), for example. However, despite these advancements in capability, there are cases in which the salt boundary is either poorly illuminated or completely absent in the migrated image. To provide a solution to this problem, we have developed two RTM methods for imaging salt boundaries, which use transmitted wavefields. In the first technique, downgoing waves, typically recorded in walkaway vertical seismic profile surveys, are used to image the salt flank via the generation of aplanatic isochrones. This image can be generated in the absence of an explicit interpretation of the salt flank using dual migration velocity models, as demonstrated on a 3D walkaway field data set from the Gulf of Mexico. In the second technique, we extend the basic theory to include imaging of upgoing source wavefields, which are transmitted at the base salt from below, as acquired by a surface acquisition geometry. This technique has similarities to the prism-imaging method, yet it uses transmitted instead of reflected waves at the salt boundary. Downgoing and upgoing methods are shown to satisfactorily generate an image of the salt flank; however, transmission imaging can create artifacts if reflection arrivals are included in the migration or the acquisition geometry is limited in extent. Increased wavelet stretch is also observed due to the higher transmission coefficient. An important benefit of these methods is that transmission imaging produces an opposite depth shift to errors in the velocity model compared with imaging of reflections. When combined with conventional seismic reflection surveys, this behavior can be used to provide a constraint on the accuracy of the salt and/or subsalt velocities.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB175-WB182 ◽  
Author(s):  
Yan Huang ◽  
Bing Bai ◽  
Haiyong Quan ◽  
Tony Huang ◽  
Sheng Xu ◽  
...  

The availability of wide-azimuth data and the use of reverse time migration (RTM) have dramatically increased the capabilities of imaging complex subsalt geology. With these improvements, the current obstacle for creating accurate subsalt images now lies in the velocity model. One of the challenges is to generate common image gathers that take full advantage of the additional information provided by wide-azimuth data and the additional accuracy provided by RTM for velocity model updating. A solution is to generate 3D angle domain common image gathers from RTM, which are indexed by subsurface reflection angle and subsurface azimuth angle. We apply these 3D angle gathers to subsalt tomography with the result that there were improvements in velocity updating with a wide-azimuth data set in the Gulf of Mexico.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB27-WB39 ◽  
Author(s):  
Zheng-Zheng Zhou ◽  
Michael Howard ◽  
Cheryl Mifflin

Various reverse time migration (RTM) angle gather generation techniques have been developed to address poor subsalt data quality and multiarrival induced problems in gathers from Kirchhoff migration. But these techniques introduce new problems, such as inaccuracies in 2D subsurface angle gathers and edge diffraction artifacts in 3D subsurface angle gathers. The unique rich-azimuth data set acquired over the Shenzi field in the Gulf of Mexico enabled the generally artifact-free generation of 3D subsurface angle gathers. Using this data set, we carried out suprasalt tomography and salt model building steps and then produced 3D angle gathers to update the subsalt velocity. We used tilted transverse isotropy RTM with extended image condition to generate full 3D subsurface offset domain common image gathers, which were subsequently converted to 3D angle gathers. The angle gathers were substacked along the subsurface azimuth axis into azimuth sectors. Residual moveout analysis was carried out, and ray-based tomography was used to update velocities. The updated velocity model resulted in improved imaging of the subsalt section. We also applied residual moveout and selective stacking to 3D angle gathers from the final migration to produce an optimized stack image.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA13-WA21 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin

Nonhyperbolic moveout analysis plays an increasingly important role in velocity model building because it provides valuable information for anisotropic parameter estimation. However, lateral heterogeneity associated with stratigraphic lenses such as channels and reefs can significantly distort the moveout parameters, even when the structure is relatively simple. We analyze the influence of a low-velocity isotropic lens on nonhyperbolic moveout inversion for horizontally layered VTI (transversely isotropic with a vertical symmetry axis) models. Synthetic tests demonstrate that a lens can cause substantial, laterally varying errors in the normal-moveout velocity [Formula: see text] and the anellipticity parameter [Formula: see text]. The area influenced by the lens can be identified using the residual moveout after the nonhyperbolic moveout correction as well as the dependence of errors in [Formula: see text] and [Formula: see text] on spreadlength. To remove such errors in [Formula: see text] and [Formula: see text], we propose a correction algorithm designed for a lens embedded in a horizontally layered overburden. This algorithm involves estimation of the incidence angle of the ray passing through the lens for each recorded trace. With the assumption that lens-related perturbation of the raypath is negligible, the lens-induced traveltime shifts are computed from the corresponding zero-offset time distortion (i.e., from “pull-up” or “push-down” anomalies). Synthetic tests demonstrate that this algorithm substantially reduces the errors in the effective and interval parameters [Formula: see text] and [Formula: see text]. The corrected traces and reconstructed “background” values of [Formula: see text] and [Formula: see text] are suitable for anisotropic time imaging and producing a high-quality stack.


Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1944-1946
Author(s):  
M. Tygel ◽  
J. Schleicher ◽  
P. Hubral

We highly appreciate the useful remarks of Dr. Barnes relating our work to well‐known practical seismic processing effects. This is of particular interest as normal‐moveout (NMO) correction and post‐stack time migration are still two very important processing steps. Most exploration geophysicists know about the significance of pulse distortions known as “NM0 stretch” and “frequency shifting due to zero‐offset time migration.” As a result of the discussion of Dr. Barnes, it should now be possible to better appreciate the importance of our very general formulas (27) describing the pulse distortion of seismic reflections from an arbitrarily curved subsurface reflector when subjected to a prestack depth migration in 3‐D laterally inhomogeneous media. This discussion thus relates in particular to such important questions as how to correctly sample signals in the time or depth domain in order to avoid spatial aliasing, or how to stack seismic data without loss of information due to destructive interference of wavelets of different lengths.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 546-556 ◽  
Author(s):  
Herman Chang ◽  
John P. VanDyke ◽  
Marcelo Solano ◽  
George A. McMechan ◽  
Duryodhan Epili

Portable, production‐scale 3-D prestack Kirchhoff depth migration software capable of full‐volume imaging has been successfully implemented and applied to a six‐million trace (46.9 Gbyte) marine data set from a salt/subsalt play in the Gulf of Mexico. Velocity model building and updates use an image‐driven strategy and were performed in a Sun Sparc environment. Images obtained by 3-D prestack migration after three velocity iterations are substantially better focused and reveal drilling targets that were not visible in images obtained from conventional 3-D poststack time migration. Amplitudes are well preserved, so anomalies associated with known reservoirs conform to the petrophysical predictions. Prototype development was on an 8-node Intel iPSC860 computer; the production version was run on an 1824-node Intel Paragon computer. The code has been successfully ported to CRAY (T3D) and Unix workstation (PVM) environments.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. S105-S111 ◽  
Author(s):  
Sheng Xu ◽  
Feng Chen ◽  
Bing Tang ◽  
Gilles Lambare

When using seismic data to image complex structures, the reverse time migration (RTM) algorithm generally provides the best results when the velocity model is accurate. With an inexact model, moveouts appear in common image gathers (CIGs), which are either in the surface offset domain or in subsurface angle domain; thus, the stacked image is not well focused. In extended image gathers, the strongest energy of a seismic event may occur at non-zero-lag in time-shift or offset-shift gathers. Based on the operation of RTM images produced by the time-shift imaging condition, the non-zero-lag time-shift images exhibit a spatial shift; we propose an approach to correct them by a second pass of migration similar to zero-offset depth migration; the proposed approach is based on the local poststack depth migration assumption. After the proposed second-pass migration, the time-shift CIGs appear to be flat and can be stacked. The stack enhances the energy of seismic events that are defocused at zero time lag due to the inaccuracy of the model, even though the new focused events stay at the previous positions, which might deviate from the true positions of seismic reflection. With the stack, our proposed approach is also able to attenuate the long-wavelength RTM artifacts. In the case of tilted transverse isotropic migration, we propose a scheme to defocus the coherent noise, such as migration artifacts from residual multiples, by applying the original migration velocity model along the symmetry axis but with different anisotropic parameters in the second pass of migration. We demonstrate that our approach is effective to attenuate the coherent noise at subsalt area with two synthetic data sets and one real data set from the Gulf of Mexico.


Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. Q15-Q26 ◽  
Author(s):  
Giovanni Angelo Meles ◽  
Kees Wapenaar ◽  
Andrew Curtis

State-of-the-art methods to image the earth’s subsurface using active-source seismic reflection data involve reverse time migration. This and other standard seismic processing methods such as velocity analysis provide best results only when all waves in the data set are primaries (waves reflected only once). A variety of methods are therefore deployed as processing to predict and remove multiples (waves reflected several times); however, accurate removal of those predicted multiples from the recorded data using adaptive subtraction techniques proves challenging, even in cases in which they can be predicted with reasonable accuracy. We present a new, alternative strategy to construct a parallel data set consisting only of primaries, which is calculated directly from recorded data. This obviates the need for multiple prediction and removal methods. Primaries are constructed by using convolutional interferometry to combine the first-arriving events of upgoing and direct-wave downgoing Green’s functions to virtual receivers in the subsurface. The required upgoing wavefields to virtual receivers are constructed by Marchenko redatuming. Crucially, this is possible without detailed models of the earth’s subsurface reflectivity structure: Similar to the most migration techniques, the method only requires surface reflection data and estimates of direct (nonreflected) arrivals between the virtual subsurface sources and the acquisition surface. We evaluate the method on a stratified synclinal model. It is shown to be particularly robust against errors in the reference velocity model used and to improve the migrated images substantially.


2019 ◽  
Vol 219 (3) ◽  
pp. 1876-1884 ◽  
Author(s):  
M Glöckner ◽  
J Walda ◽  
S Dell ◽  
D Gajewski ◽  
J Karstens ◽  
...  

SUMMARY Standard seismic acquisition and processing require appropriate source–receiver offsets. P-cable technology represents the opposite, namely, very short source–receiver offsets at the price of increased spatial and lateral resolution with a high-frequency source. To use this advantage, a processing flow excluding offset information is required. This aim can be achieved with a processing tuned to diffractions because point diffractions scatter the same information in the offset and midpoint direction. Usually, diffractions are small amplitude events and a careful diffraction separation is required as a first step. We suggest the strategy to use a multiparameter stacking operator, for example, common-reflection surface, and stack along the midpoint direction. The obtained kinematic wave-front attributes are used to calculate time-migration velocities. A diffractivity map serves as a filter to refine the velocities. This strategy is applied to a 3-D P-cable data set to obtain a time-migrated image.


Sign in / Sign up

Export Citation Format

Share Document