Depth-variable frequency reverse time migration: Imaging results from Perdido fold belt and AVO analysis over SEAM model

Accurate seismic imaging can ensure safe and efficient tunnel construction under complex geological conditions. As a high-precision migration method, reverse time migration (RTM) has been introduced into tunnel seismic forward-prospecting. However, the resolution of traditional RTM imaging results may not meet the requirements in a complex tunnel environment, which affects the interpretation of tunnel seismic forward-prospecting results. In this study, we propose a least-squares RTM method based on the decoupled elastic wave equation in tunnels. The Born forward modeling operator and its exact adjoint migration imaging operator are derived to ensure a stable convergence of the conjugate gradient method. Moreover, a pseudo-Hessian based preconditioning operator is adopted to accelerate the convergence. Numerical examples are provided to verify the efficiency of the proposed scheme. A field test in a traffic tunnel construction site is performed to show the good application effect of the decoupled elastic least-squares RTM in practical situations.

Download Full-text

Phaseless Imaging by Reverse Time Migration: Acoustic Waves

Numerical Mathematics Theory Methods and Applications ◽

10.4208/nmtma.2017.m1617 ◽

2017 ◽

Vol 10 (1) ◽

pp. 1-21 ◽

Cited By ~ 23

Author(s):

Zhiming Chen ◽

Guanghui Huang

Keyword(s):

Acoustic Waves ◽

Numerical Experiments ◽

Cross Correlation ◽

Scattering Data ◽

Imaging Method ◽

Total Field ◽

Reverse Time ◽

Reverse Time Migration ◽

Time Migration ◽

Imaging Results

AbstractWe propose a reliable direct imaging method based on the reverse time migration for finding extended obstacles with phaseless total field data. We prove that the imaging resolution of the method is essentially the same as the imaging results using the scattering data with full phase information when the measurement is far away from the obstacle. The imaginary part of the cross-correlation imaging functional always peaks on the boundary of the obstacle. Numerical experiments are included to illustrate the powerful imaging quality

Download Full-text

An efficient local operator-based Q-compensated reverse time migration algorithm with multistage optimization

Geophysics ◽

10.1190/geo2017-0026.1 ◽

2018 ◽

Vol 83 (3) ◽

pp. S249-S259 ◽

Cited By ~ 5

Author(s):

Tong Zhou ◽

Wenyi Hu ◽

Jieyuan Ning

Keyword(s):

Seismic Attenuation ◽

Depth Information ◽

Frequency Noise ◽

Processing Unit ◽

Reverse Time ◽

Reverse Time Migration ◽

Pass Filter ◽

Low Pass Filter ◽

Time Migration ◽

Imaging Results

Most existing [Formula: see text]-compensated reverse time migration ([Formula: see text]-RTM) algorithms are based on pseudospectral methods. Because of the global nature of pseudospectral operators, these methods are not ideal for efficient parallelization, implying that they may suffer from high computational cost and inefficient memory usage for large-scale industrial problems. In this work, we reported a novel [Formula: see text]-RTM algorithm — the multistage optimized [Formula: see text]-RTM method. This [Formula: see text]-RTM algorithm uses a finite-difference method to compensate the amplitude and the phase simultaneously by uniquely combining two techniques: (1) a negative [Formula: see text] method for amplitude compensation and (2) a multistage dispersion optimization technique for phase correction. To prevent high-frequency noise from growing exponentially and ruining the imaging results, we apply a finite impulse response low-pass filter using the Kaiser window. The theoretical analyses and numerical experiments demonstrate that this [Formula: see text]-RTM algorithm precisely recovers the decayed amplitude and corrects the distorted phase caused by seismic attenuation effects, and hence produces higher resolution subsurface images with the correct structural depth information. This new method performs best in the frequency range of 10–70 Hz. Compared with pseudospectral [Formula: see text]-RTM methods, this [Formula: see text]-RTM approach offers nearly identical imaging quality. Based on local numerical differential operators, this [Formula: see text]-RTM method is very suitable for parallel computing and graphic processing unit implementation, an important feature for large 3D seismic surveys.

Download Full-text

Near-Borehole Imaging Using Full-Waveform Sonic Data

10.2118/204765-ms ◽

2021 ◽

Author(s):

Hala Alqatari ◽

Thierry-Laurent Tonellot ◽

Mohammed Mubarak

Keyword(s):

Seismic Data ◽

Reverse Time ◽

Reverse Time Migration ◽

Full Waveform ◽

Time Migration ◽

Imaging Results ◽

High Resolution Images ◽

Recorded Data ◽

The Impact ◽

Processing Steps

Abstract This work presents a full waveform sonic (FWS) dataset processing to generate high-resolution images of the near-borehole area. The dataset was acquired in a nearly horizontal well over a distance of 5400 feet. Multiple formation boundaries can be identified on the final image and tracked at up to 200 feet deep, along the wellbore's trajectory. We first present a new preprocessing sequence to prepare the sonic data for imaging. This sequence leverages denoising algorithms used in conventional surface seismic data processing to remove unwanted components of the recorded data that could harm the imaging results. We then apply a reverse time migration algorithm to the data at different processing stages to assess the impact of the main processing steps on the final image.

Download Full-text

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

Geophysics ◽

10.1190/geo2018-0619.1 ◽

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.

Download Full-text

Multi-scale high order finite different reverse time migration imaging

EURASIP Journal on Image and Video Processing ◽

10.1186/s13640-018-0290-x ◽

2018 ◽

Vol 2018 (1) ◽

Author(s):

Xiaodan Zhang ◽

Lei Zhu ◽

Hua Guo

Keyword(s):

High Order ◽

Reverse Time ◽

Reverse Time Migration ◽

Multi Scale ◽

Time Migration ◽

Migration Imaging

Download Full-text

An efficient vector elastic reverse time migration method in the hybrid time and frequency domain for anisotropic media

Geophysics ◽

10.1190/geo2018-0644.1 ◽

2019 ◽

Vol 84 (6) ◽

pp. S511-S522 ◽

Cited By ~ 1

Author(s):

Kai Gao ◽

Lianjie Huang

Keyword(s):

Frequency Domain ◽

Anisotropic Media ◽

Reverse Time ◽

Computationally Efficient ◽

Reverse Time Migration ◽

Computational Costs ◽

Time Migration ◽

Migration Imaging ◽

Efficient Vector ◽

Wave Migration

Vector elastic reverse time migration (ERTM) produces subsurface elastic images with correct polarities using multicomponent seismic data. However, the decomposition of elastic wavefields into vector P- and S-wavefields is computationally expensive, particularly in heterogeneous and complex anisotropic media. We have developed a computationally efficient vector ERTM method in the hybrid time and frequency domain by combining three existing techniques. Rather than decomposing elastic wavefields into vector qP- and qS-wavefields during time-domain wavefield propagation, we conduct the wavefield decomposition in the frequency domain for several selected frequencies. In general, the number of selected frequencies needed for migration imaging is much smaller than the number of time steps during forward and backward wavefield propagation, leading to greatly reduced computational costs associated with the qP-/qS-wavefield vector separation in complex heterogeneous anisotropic media. We further combine an implicit directional wavefield separation into the vector ERTM to enhance the image quality. The numerical results demonstrate that our method produces high-quality elastic-wave migration images with notably reduced computational costs compared to the conventional vector ERTM method.

Download Full-text