A modeling study of preprocessing considerations for reverse-time migration

Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. T99-T106 ◽  
Author(s):  
Ian F. Jones

Much of the thinking behind conventional geophysical processing assumes that we wanted to image energy that propagates down from the surface of the earth, scatters from a reflector or diffractor, and then propagates back up to the recording surface without being reflected by any other feature. Such travel paths conform to the assumptions of one-way wave propagation, and most contemporary migration schemes are designed to image such data. In addition, the moveout behavior of these primary reflection events in the various prestack domains is well understood, and many of our standard data-preprocessing techniques relied on the assumption that this behavior adequately describes the events we wanted to preserve for imaging. As a corollary, events that do not conform to this prescribed behavior are classified as noise, and many of our standard preprocessing techniqueswere designed to remove them. We assessed the kinematics of moveout behavior of events that arise from two-way wave propagation and the effect of certain preprocessing techniques on those events. This was of interest to us because the recent rapid increase in available cost-effective computing power has enabled industrial implementation of migration algorithms—particularly reverse-time migration—that in principle can image events that reflect more than once on their way from source to receiver. We used 2D synthetic data to show that some conventional data-processing steps—particularly those used in suppression of complex reverberations (“multiples”)—remove nonreverberatory primary events from seismic reflection data. Specifically, they remove events that have repeated or turning reflections in the subsurface (such as double-bounce arrivals) but that otherwise are imageable using reverse-time migration.

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.


Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. S143-S149 ◽  
Author(s):  
Francisco A. da Silva Neto ◽  
Jessé C. Costa ◽  
Jörg Schleicher ◽  
Amélia Novais

Reverse-time migration (RTM) in 2.5D offers an alternative to improve resolution and amplitude when imaging 2D seismic data. Wave propagation in 2.5D assumes translational invariance of the velocity model. Under this assumption, we implement a finite-difference (FD) modeling algorithm in the mixed time-space/wavenumber domain to simulate the velocity and pressure fields for acoustic wave propagation and apply it in RTM. The 2.5D FD algorithm is truly parallel, allowing an efficient implementation in clusters. Storage and computing time requirements are strongly reduced compared to a full 3D FD simulation of the wave propagation. This feature makes 2.5D RTM much more efficient than 3D RTM, while achieving improved modeling of 3D geometrical spreading and phase properties of the seismic waveform in comparison to 2D RTM. Together with an imaging condition that compensates for uneven illumination and/or the obliquity factor, this allows recover of amplitudes proportional to the earth’s reflectivity. Numerical experiments using synthetic data demonstrate the better resolution and improved amplitude recovery of 2.5D RTM relative to 2D RTM.


Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. S37-S46 ◽  
Author(s):  
Bao D. Nguyen ◽  
George A. McMechan

An implicitly stable ratio imaging condition for prestack reverse-time migration (RTM) was defined using excitation criteria. Amplitude maxima and their corresponding occurrence times were saved at each grid point during forward source wavefield extrapolation. Application of the imaging condition involves dividing the amplitudes of the back-propagated receiver wavefield by the precomputed maximum source wavefield amplitude only at the grid points that satisfy the image time at each time step. The division normalizes by the source amplitude, so only the highest signal-to-noise ratio portion of the data is used. Provided that the source and receiver wavefield amplitudes are accurate at the reflection points, the peak wavelet amplitudes in the migrated image are the angle-dependent reflection coefficients and low wavenumber artifacts are significantly reduced compared to those in images calculated by crosscorrelation. Using excitation information and time-binning for the imaging condition improves computational and storage efficiency by three or more orders of magnitude when compared to crosscorrelation with the full source wavefield. Numerical tests with synthetic data for the Marmousi2 model have shown this method to be a cost-effective and practical imaging condition for use in prestack RTM.


2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo

ABSTRACT. Reverse time migration (RTM) is one of the most powerful methods used to generate images of the subsurface. The RTM was proposed in the early 1980s, but only recently it has been routinely used in exploratory projects involving complex geology – Brazilian pre-salt, for example. Because the method uses the two-way wave equation, RTM is able to correctly image any kind of geological environment (simple or complex), including those with anisotropy. On the other hand, RTM is computationally expensive and requires the use of computer clusters. This paper proposes to investigate the influence of anisotropy on seismic imaging through the application of RTM for tilted transversely isotropic (TTI) media in pre-stack synthetic data. This work presents in detail how to implement RTM for TTI media, addressing the main issues and specific details, e.g., the computational resources required. A couple of simple models results are presented, including the application to a BP TTI 2007 benchmark model.Keywords: finite differences, wave numerical modeling, seismic anisotropy. Migração reversa no tempo em meios transversalmente isotrópicos inclinadosRESUMO. A migração reversa no tempo (RTM) é um dos mais poderosos métodos utilizados para gerar imagens da subsuperfície. A RTM foi proposta no início da década de 80, mas apenas recentemente tem sido rotineiramente utilizada em projetos exploratórios envolvendo geologia complexa, em especial no pré-sal brasileiro. Por ser um método que utiliza a equação completa da onda, qualquer configuração do meio geológico pode ser corretamente tratada, em especial na presença de anisotropia. Por outro lado, a RTM é dispendiosa computacionalmente e requer o uso de clusters de computadores por parte da indústria. Este artigo apresenta em detalhes uma implementação da RTM para meios transversalmente isotrópicos inclinados (TTI), abordando as principais dificuldades na sua implementação, além dos recursos computacionais exigidos. O algoritmo desenvolvido é aplicado a casos simples e a um benchmark padrão, conhecido como BP TTI 2007.Palavras-chave: diferenças finitas, modelagem numérica de ondas, anisotropia sísmica.


Geophysics ◽  
1983 ◽  
Vol 48 (11) ◽  
pp. 1514-1524 ◽  
Author(s):  
Edip Baysal ◽  
Dan D. Kosloff ◽  
John W. C. Sherwood

Migration of stacked or zero‐offset sections is based on deriving the wave amplitude in space from wave field observations at the surface. Conventionally this calculation has been carried out through a depth extrapolation. We examine the alternative of carrying out the migration through a reverse time extrapolation. This approach may offer improvements over existing migration methods, especially in cases of steeply dipping structures with strong velocity contrasts. This migration method is tested using appropriate synthetic data sets.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. H27-H33 ◽  
Author(s):  
Jun Ji

To reduce the migration artifacts arising from incomplete data or inaccurate operators instead of migrating data with the adjoint of the forward-modeling operator, a least-squares migration often is considered. Least-squares migration requires a forward-modeling operator and its adjoint. In a derivation of the mathematically correct adjoint operator to a given forward-time-extrapolation modeling operator, the exact adjoint of the derived operator is obtained by formulating an explicit matrix equation for the forward operation and transposing it. The programs that implement the exact adjoint operator pair are verified by the dot-product test. The derived exact adjoint operator turns out to differ from the conventional reverse-time-migration (RTM) operator, an implementation of wavefield extrapolation backward in time. Examples with synthetic data show that migration using the exact adjoint operator gives similar results for a conventional RTM operator and that least-squares RTM is quite successful in reducing most migration artifacts. The least-squares solution using the exact adjoint pair produces a model that fits the data better than one using a conventional RTM operator pair.


Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. S241-S250 ◽  
Author(s):  
Yi Luo ◽  
Qinglin Liu ◽  
Yuchun E. Wang ◽  
Mohammed N. AlFaraj

We illustrate the use of mode-converted transmitted (e.g., PS- or SP-) waves in vertical seismic profiling (VSP) data for imaging areas above receivers where reflected waves cannot illuminate. Three depth-domain imaging techniques — move-out correction, common-depth-point (CDP) mapping, and prestack migration — are described and used for imag-ing the transmitted waves. Moveout correction converts an offset VSP trace into a zero-offset trace. CDP mapping maps each sample on an input trace to the location where the mode conversion occurs. For complex media, prestack migration (e.g., reverse-time migration) is used. By using both synthetic and field VSP data, we demonstrate that images derived from transmissions complement those from reflections. As an important application, we show that transmitted waves can illuminate zones above highly de-viated or horizontal wells, a region not imaged by reflection data. Because all of these benefits are obtained without extra data acquisition cost, we believe transmission imag-ing techniques will become widely adopted by the oil in-dustry.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S411-S423
Author(s):  
Peng Yong ◽  
Jianping Huang ◽  
Zhenchun Li ◽  
Wenyuan Liao ◽  
Luping Qu

Least-squares reverse time migration (LSRTM), an effective tool for imaging the structures of the earth from seismograms, can be characterized as a linearized waveform inversion problem. We have investigated the performance of three minimization functionals as the [Formula: see text] norm, the hybrid [Formula: see text] norm, and the Wasserstein metric ([Formula: see text] metric) for LSRTM. The [Formula: see text] metric used in this study is based on the dynamic formulation of transport problems, and a primal-dual hybrid gradient algorithm is introduced to efficiently compute the [Formula: see text] metric between two seismograms. One-dimensional signal analysis has demonstrated that the [Formula: see text] metric behaves like the [Formula: see text] norm for two amplitude-varied signals. Unlike the [Formula: see text] norm, the [Formula: see text] metric does not suffer from the differentiability issue for null residuals. Numerical examples of the application of three misfit functions to LSRTM on synthetic data have demonstrated that, compared to the [Formula: see text] norm, the hybrid [Formula: see text] norm and [Formula: see text] metric can accelerate LSRTM and are less sensitive to non-Gaussian noise. For the field data application, the [Formula: see text] metric produces the most reliable imaging results. The hybrid [Formula: see text] norm requires tedious trial-and-error tests for the judicious threshold parameter selection. Hence, the more automatic [Formula: see text] metric is recommended as a robust alternative to the customary [Formula: see text] norm for time-domain LSRTM.


Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. S47-S64
Author(s):  
Yang Zhao ◽  
Tao Liu ◽  
Xueyi Jia ◽  
Hongwei Liu ◽  
Zhiguang Xue ◽  
...  

Angle-domain common-image gathers (ADCIGs) from elastic reverse time migration (ERTM) are valuable tools for seismic elastic velocity estimation. Traditional ADCIGs are based on the concept of common-offset domains, but common-shot domain implementations are often favored for computational cost considerations. Surface-offset gathers (SOGs) built from common-offset migration may serve as an alternative to the common-shot ADCIGs. We have developed a theoretical kinematic framework between these two domains, and we determined that the common SOG gives an alternative measurement of kinematic correctness in the presence of incorrect velocity. Specifically, we exploit analytical expressions for the image misposition between these two domains, with respect to the traveltime perturbation caused by velocity errors. Four formulations of the PP and PS residual moveout functions are derived and provide insightful information of the velocity error, angle, and PS velocity ratio contained in ERTM gathers. The analytical solutions are validated with homogeneous examples with a series of varied parameters. We found that the SOGs may perform in the way of simplicity and linearity as an alternative to the common-shot migration. To make a full comparison with ADCIGs, we have developed a cost-effective workflow of ERTM SOGs. A fast vector P- and S-wave decomposition can be obtained via spatial gradients at selected time steps. A selected ERTM imaging condition is then modified in which the migration is done by offset groups between each source and receiver pair for each P- and S-wave decomposition. Two synthetic (marine and land) examples are used to demonstrate the feasibility of our methods.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. T167-T174 ◽  
Author(s):  
Dan Kosloff ◽  
Reynam C. Pestana ◽  
Hillel Tal-Ezer

A new scheme for the calculation of spatial derivatives has been developed. The technique is based on recursive derivative operators that are generated by an [Formula: see text] fit in the spectral domain. The use of recursive operators enables us to extend acoustic and elastic wave simulations to shorter wavelengths. The method is applied to the numerical solution of the 2D acoustic wave equation and to the solution of the equations of 2D dynamic elasticity in an isotropic medium. An example of reverse-time migration of a synthetic data set shows that the numerical dispersion can be significantly reduced with respect to schemes that are based on finite differences. The method is tested for the solutions of the equations of dynamic elasticity by comparing numerical and analytic solutions to Lamb’s problem.


Sign in / Sign up

Export Citation Format

Share Document