Acoustic prestack migration of cross‐hole data

Geophysics ◽  
1988 ◽  
Vol 53 (8) ◽  
pp. 1015-1023 ◽  
Author(s):  
Liang‐Zie Hu ◽  
George A. McMechan ◽  
Jerry M. Harris
Keyword(s):  
Synthetic Data ◽  
Field Application ◽  
Scale Model ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Prestack Migration ◽  
Time Migration ◽  
Time Space ◽  
Low Pass

Subsurface imaging with common‐source cross‐hole data can be achieved using prestack reverse‐time migration. The algorithm consists of extrapolation of the recorded wave field, application of the excitation‐time imaging condition, and postprocessing of the resulting image with a low‐pass wavenumber filter. The wavenumber filter removes the artifact associated with the direct arrival; this artifact is not separable from the scattered data before migration because, in the cross‐hole geometry, they significantly overlap in time, space, and wavenumber. Migration of synthetic data produces the best possible results, but images produced by migration of scale‐model data are not greatly inferior. Apparently, acceptable images can be obtained from a surprisingly few sources, if these sources are located sufficiently far apart to give independent information and the recording aperture is sufficiently wide.

3-D prestack migration in anisotropic media

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 79-90 ◽  
Author(s):  
Zhengxin Dong ◽  
George A. McMechan
Keyword(s):  
A Priori ◽  
Three Dimensional ◽  
Synthetic Data ◽  
Anisotropic Media ◽  
P Wave ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Prestack Migration ◽  
Time Migration

A three‐dimensional (3-D) prestack reverse‐time migration algorithm for common‐source P‐wave data from anisotropic media is developed and illustrated by application to synthetic data. Both extrapolation of the data and computation of the excitation‐time imaging condition are implemented using a second‐order finite‐ difference solution of the 3-D anisotropic scalar‐wave equation. Poorly focused, distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation; well focused, clear images are obtained using anisotropic extrapolation. A priori estimation of the 3-D anisotropic velocity distribution is required. Zones of anomalous, directionally dependent reflectivity associated with anisotropic fracture zones are detectable in both the 3-D common‐ source data and the corresponding migrated images.

2.5D reverse-time migration

Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. S143-S149 ◽  
Author(s):  
Francisco A. da Silva Neto ◽  
Jessé C. Costa ◽  
Jörg Schleicher ◽  
Amélia Novais
Keyword(s):  
Wave Propagation ◽  
Computing Time ◽  
Synthetic Data ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Translational Invariance ◽  
Time Migration ◽  
Seismic Waveform ◽  
Time Space

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.

True-amplitude prestack depth migration

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. S155-S166 ◽  
Author(s):  
Feng Deng ◽  
George A. McMechan
Keyword(s):  
Velocity Ratio ◽  
Synthetic Data ◽  
Transmission Losses ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Prestack Migration ◽  
Time Migration ◽  
Geometric Spreading

Most current true-amplitude migrations correct only for geometric spreading. We present a new prestack depth-migration method that uses the framework of reverse-time migration to compensate for geometric spreading, intrinsic [Formula: see text] losses, and transmission losses. Geometric spreading is implicitly compensated by full two-way wave propagation. Intrinsic [Formula: see text] losses are handled by including a [Formula: see text]-dependent term in the wave equation. Transmission losses are compensated based on an estimation of angle-dependent reflectivity using a two-pass recursive reverse-time prestack migration. The image condition used is the ratio of receiver/source wavefield amplitudes. Two-dimensional tests using synthetic data for a dipping-layer model and a salt model show that loss-compensating prestack depth migration can produce reliable angle-dependent reflection coefficients at the target. The reflection coefficient curves are fitted to give least-squares estimates of the velocity ratio at the target. The main new result is a procedure for transmission compensation when extrapolating the receiver wavefield. There are still a number of limitations (e.g., we use only scalar extrapolation for illustration), but these limitations are now better defined.

Implicit static corrections in prestack migration of common‐source data

Geophysics ◽  
1990 ◽  
Vol 55 (6) ◽  
pp. 757-760 ◽  
Author(s):  
G. A. McMechan ◽  
H. W. Chen
Keyword(s):  
Surface Velocity ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Near Surface ◽  
Prestack Migration ◽  
Time Migration ◽  
Excitation Time ◽  
Source Data ◽  
Static Corrections

Static effects due to surface topography and near‐surface velocity variations may be accurately compensated for, in an implicit way, during prestack reverse‐time migration of common‐source gathers, obviating the need for explicit static corrections. Receiver statics are incorporated by extrapolating the observed data from the actual recorder positions; source statics are incorporated by computing the excitation‐time imaging conditions from the actual source positions.

Reverse time migration in tilted transversely isotropic media

2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo
Keyword(s):  
Seismic Anisotropy ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Computer Clusters ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Tti Media

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.

Reverse time migration

Geophysics ◽  
1983 ◽  
Vol 48 (11) ◽  
pp. 1514-1524 ◽  
Author(s):  
Edip Baysal ◽  
Dan D. Kosloff ◽  
John W. C. Sherwood
Keyword(s):  
Wave Amplitude ◽  
Synthetic Data ◽  
Data Sets ◽  
Field Observations ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Method ◽  
Zero Offset ◽  
Time Extrapolation

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.

An exact adjoint operation pair in time extrapolation and its application in least-squares reverse-time migration

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. H27-H33 ◽  
Author(s):  
Jun Ji
Keyword(s):  
Least Squares ◽  
Adjoint Operator ◽  
Synthetic Data ◽  
Forward Modeling ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Product Test ◽  
Least Squares Migration ◽  
Time Extrapolation

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.

Imaging reflection-blind areas using transmitted PS-waves

Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. S241-S250 ◽  
Author(s):  
Yi Luo ◽  
Qinglin Liu ◽  
Yuchun E. Wang ◽  
Mohammed N. AlFaraj
Keyword(s):  
Mode Conversion ◽  
Imaging Techniques ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Reflected Waves ◽  
Prestack Migration ◽  
Reflection Data ◽  
Time Migration ◽  
Transmitted Waves ◽  
Vsp Data

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.

Least-squares reverse time migration via linearized waveform inversion using a Wasserstein metric

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S411-S423
Author(s):  
Peng Yong ◽  
Jianping Huang ◽  
Zhenchun Li ◽  
Wenyuan Liao ◽  
Luping Qu
Keyword(s):  
Least Squares ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Gradient Algorithm ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wasserstein Metric ◽  
Time Migration ◽  
Imaging Results ◽  
Primal Dual

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.

Elastic reverse‐time migration

Geophysics ◽  
1987 ◽  
Vol 52 (10) ◽  
pp. 1365-1375 ◽  
Author(s):  
Wen‐Fong Chang ◽  
George A. McMechan
Keyword(s):  
Finite Difference ◽  
Vertical Component ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration ◽  
Two Component ◽  
Synthetic Test ◽  
Data Extrapolation

Elastic, prestack, reverse‐time, finite‐difference migration of two‐component seismic surface data requires data extrapolation and application of an imaging condition. Data extrapolation involves synchronous driving of the vertical‐component and horizontal‐component finite‐difference meshes with the time reverse of the recorded vertical and horizontal traces, respectively. Extrapolation uses the coupled elastic wave equation for variable velocity solved with a second‐order, explicit finite‐difference scheme. The imaging condition at any point in the grid is the one‐way traveltime from the source to that point. Elastic migrations of both synthetic test data and real two‐component common‐source gathers produce simpler images than acoustic migrations because of the coalescing of double reflections (compressional waves and shear waves) into single loci.

Sign in / Sign up

Export Citation Format

Share Document