Geometric-mean reverse-time migration for elastic media

Author(s):  
Tong Bai ◽  
Bin Lyu ◽  
Fuchun Gao ◽  
Paul Williamson ◽  
Nori Nakata
Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. S355-S364 ◽  
Author(s):  
Jianhang Yin ◽  
Nori Nakata

Diffracted waves contain a great deal of valuable information about small-scale subsurface structure such as faults, pinch-outs, karsts, and fractures, which are closely related to hydrocarbon accumulation and production. Therefore, diffraction separation and imaging with high spatial resolution play an increasingly critical role in seismic exploration. We have applied the geometric-mean reverse time migration (GmRTM) method to diffracted waves for imaging only subsurface diffractors based on the difference of the wave phenomena between diffracted and reflected waves. Numerical tests prove the advantages of this method on diffraction imaging with higher resolution as well as fewer artifacts compared to conventional RTM even when we only have a small number of receivers. Then, we developed a workflow to extract diffraction information using a fully data-driven method, called common-reflection surface (CRS), before we applied GmRTM. Application of this workflow indicates that GmRTM further improves the quality of the image by combining with the diffraction-separation technique CRS in the data domain.


2016 ◽  
Vol 208 (2) ◽  
pp. 1103-1125 ◽  
Author(s):  
Zhiming Ren ◽  
Yang Liu ◽  
Mrinal K. Sen

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. KS51-KS60 ◽  
Author(s):  
Nori Nakata ◽  
Gregory C. Beroza

Time reversal is a powerful tool used to image directly the location and mechanism of passive seismic sources. This technique assumes seismic velocities in the medium and propagates time-reversed observations of ground motion at each receiver location. Assuming an accurate velocity model and adequate array aperture, the waves will focus at the source location. Because we do not know the location and the origin time a priori, we need to scan the entire 4D image (3D in space and 1D in time) to localize the source, which makes time-reversal imaging computationally demanding. We have developed a new approach of time-reversal imaging that reduces the computational cost and the scanning dimensions from 4D to 3D (no time) and increases the spatial resolution of the source image. We first individually extrapolate wavefields at each receiver, and then we crosscorrelate these wavefields (the product in the frequency domain: geometric mean). This crosscorrelation creates another imaging condition, and focusing of the seismic wavefields occurs at the zero time lag of the correlation provided the velocity model is sufficiently accurate. Due to the analogy to the active-shot reverse time migration (RTM), we refer to this technique as the geometric-mean RTM or GmRTM. In addition to reducing the dimension from 4D to 3D compared with conventional time-reversal imaging, the crosscorrelation effectively suppresses the side lobes and yields a spatially high-resolution image of seismic sources. The GmRTM is robust for random and coherent noise because crosscorrelation enhances signal and suppresses noise. An added benefit is that, in contrast to conventional time-reversal imaging, GmRTM has the potential to be used to retrieve velocity information by analyzing time and/or space lags of crosscorrelation, which is similar to what is done in active-source imaging.


Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S173-S183 ◽  
Author(s):  
Hejun Zhu

Divergence and curl operators used for the decomposition of P- and S-wave modes in elastic reverse time migration (RTM) change the amplitudes, units, and phases of extrapolated wavefields. I separate the P- and S-waves in elastic media based on the Helmholtz decomposition. The decomposed wavefields based on this approach have the same amplitudes, units, and phases as the extrapolated wavefields. To avoid expensive multidimensional integrals in the Helmholtz decomposition, I introduce a fast Poisson solver to efficiently solve the vector Poisson’s equation. This fast algorithm allows us to reduce computational complexity from [Formula: see text] to [Formula: see text], where [Formula: see text] is the total number of grid points. Because the decomposed P- and S-waves are vector fields, I use vector imaging conditions to construct PP-, PS-, SS-, and SP-images. Several 2D numerical examples demonstrate that this approach allows us to accurately and efficiently decompose P- and S-waves in elastic media. In addition, elastic RTM images based on the vector imaging conditions have better quality and avoid polarity reversal in comparison with images based on the divergence and curl separation or direct component-by-component crosscorrelation.


Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Tong Bai ◽  
Bin Lyu ◽  
Paul Williamson ◽  
Nori Nakata

Geometric-mean Reverse-time migration (GmRTM), a powerful cross-correlation-based imaging method, generates higher-resolution source images and is more robust to noise compared to conventional time-reversal imaging. The price to pay is the higher computational costs. Alternatively, we can adopt hybrid strategies by dividing the receivers into different groups. Conventional time reversal (i.e., wavefield summation) is performed inside each group, followed by the application of cross-correlation imaging condition among different groups. Such hybrid strategies can retain the advantages of both GmRTM and time-reversal, and are often more practical than pure GmRTM. Yet, designing appropriate grouping strategy is not trivial. Here, we propose two grouping strategies (adjacent and scattered) and use synthetic and field-data examples to evaluate their performance with various group numbers. In addition to the spatial resolution of the source image, robustness to random noise is another important assessment criterion, for which we consider two distribution patterns, such as concentrated and scattered, of traces contaminated with strong random noise. We also evaluated their effectiveness to visualize events (in the image domain) that are not completely recorded by all receivers. Our comprehensive tests illustrate the respective advantages of the two grouping strategies.


2020 ◽  
Vol 223 (3) ◽  
pp. 1935-1947
Author(s):  
Bin Lyu ◽  
Nori Nakata

SUMMARY Passive-seismic provides useful information for reservoir monitoring and structural imaging; for example, the locations of microseismic events are helpful to understand the extension of the hydraulic fracturing. However, passive-seismic imaging still faces some challenges. First, it is not easy to know where the passive-seismic events happened, which is known as passive-source locating. Additionally, the accuracy of the subsurface velocity model will influence the accuracy of the estimated passive-source locations and the quality of the structural imaging obtained from the passive-seismic data. Therefore the velocity inversion using the passive-seismic data is required to provide the velocity with higher accuracy. Focusing on these challenges, we develop an iterative passive-source location estimation and velocity inversion method using geometric-mean reverse-time migration (GmRTM) and full-waveform inversion (FWI). In each iteration, the source location is estimated using a high-resolution GmRTM method, which provides a better focusing of passive-source imaging compared to conventional wavefield scanning method. The passive-source FWI is then followed to optimize the velocity model using the estimated source location provided by GmRTM. The source location estimation and velocity inversion are implemented sequentially. We evaluate this iterative method using the Marmousi model data set. The experiment result and sensitivity analysis indicate that the proposed method is effective to locate the sources and optimize velocity model in the areas with complicated subsurface structures and noisy recordings.


Sign in / Sign up

Export Citation Format

Share Document