Diffraction imaging for fault-karst structure by least-squares reverse time migration

2020 ◽  
pp. 1-40
Author(s):  
Xinru Mu ◽  
Jianping Huang ◽  
Liyun Fu ◽  
Shikai Jian ◽  
Bing Hu ◽  
...  
Keyword(s):  
Least Squares ◽  
Western China ◽  
Small Scale ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Imaging Results ◽  
Heterogeneous Structures ◽  
Karst Reservoir

The fault-karst reservoir, which evolved from the deformation and karstification of carbonate rock, is one of the most important reservoir types in western China. Along the deep-seated fault zones, there are a lot widely spread and densely distributed fractures and vugs. The energy of the diffractions generated by heterogeneous structures, such as faults, fractures and vugs, are much weaker than that of the reflections produced by continuous formation interface. When using conventional full wavefield imaging method, the imaging results of continuous layers usually cover small-scale heterogeneities. Given that, we use plane-wave destruction (PWD) filter to separate the diffractions from the full data and image the separated diffractions using least-squares reverse time migration (LSRTM) method. We use several numerical examples to demonstrate that the newly developed diffractions LSRTM (D-LSRTM) can improve the definition of the heterogeneous structures, characterize the configuration and internal structure of the fault-karst structure well and enhance the interpretation accuracy for fault-karst reservoir.

scholarly journals Phaseless Imaging by Reverse Time Migration: Acoustic Waves

2017 ◽  
Vol 10 (1) ◽  
pp. 1-21 ◽  
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

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.

scholarly journals Attenuating multiple-related imaging artifacts using combined imaging conditions

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. S469-S475 ◽  
Author(s):  
Carlos Alberto da Costa Filho ◽  
Andrew Curtis
Keyword(s):  
Small Scale ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Depth Migration ◽  
Time Migration ◽  
Migration Method ◽  
Combined Imaging ◽  
Spatial Locations ◽  
Imaging Conditions

The objective of prestack depth migration is to position reflectors at their correct subsurface locations. However, migration methods often also generate artifacts along with physical reflectors, which hamper interpretation. These spurious reflectors often appear at different spatial locations in the image depending on which migration method is used. Therefore, we have devised a postimaging filter that combines two imaging conditions to preserve their similarities and to attenuate their differences. The imaging filter is based on combining the two constituent images and their envelopes that were obtained from the complex vertical traces of the images. We have used the method to combine two images resulting from different migration schemes, which produce dissimilar artifacts: a conventional migration method (equivalent to reverse time migration) and a deconvolution-based imaging method. We show how this combination may be exploited to attenuate migration artifacts in a final image. A synthetic model containing a syncline and stochastically generated small-scale heterogeneities in the velocity and density distributions was used for the numerical example. We compared the images in detail at two locations where spurious events arose and also at a true reflector. We found that the combined imaging condition has significantly fewer artifacts than either constituent image individually.

Amplitude-preserving scalar PP and PS imaging condition for elastic reverse time migration based on a wavefield decoupling method

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. S113-S125 ◽  
Author(s):  
Xiyan Zhou ◽  
Xu Chang ◽  
Yibo Wang ◽  
Zhenxing Yao
Keyword(s):  
Wave Vector ◽  
Numerical Models ◽  
P Wave ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Decoupling Method ◽  
Time Migration ◽  
Imaging Results

To eliminate crosstalk within the imaging results of elastic reverse time migration (ERTM), we can separate the coupled P- and S-waves from the forward source wavefield and the backpropagated receiver wavefield. The P- and S-wave decoupling method retains the original phase, amplitude, and physical meaning in the separated wavefields. Thus, it is a vital wavefield separation method in ERTM. However, because these decomposed wavefields are vectors, we could consider how to retrieve scalar images that reveal the real reflectivity of the subsurface. For this purpose, we derive a scalar P-wave equation from the velocity-stress relationship for PP imaging. The phase and amplitude of this scalar P-wave are consistent with the scalarized P-wave. Therefore, this scalar P-wave can be exploited to perform PP imaging directly, with the imaging result retaining the amplitude characteristics. For PS imaging, it is difficult to calculate a dynamic preserved scalar S-wave. However, we have developed a scalar PS imaging method that divides the PS image into energy and sign components according to the geometric relationship between the wavefield vibration and propagation directions. The energy is calculated through the amplitude crosscorrelation of the forward P-wave and backpropagated S-wave from the receivers. The sign is obtained from the dot product of the forward P-wave vector and the backpropagated S-wave vector. These PP and PS imaging methods are suitable for 2D and 3D isotropic media and maintain the correct amplitude information while eliminating polarity-reversal phenomena. Several numerical models are used to verify the robustness and effectiveness of our method.

scholarly journals Research on the Improved Least Squares Reverse Time Migration Imaging Method

2021 ◽  
Vol 236 ◽  
pp. 04017
Author(s):  
Xiaodan Zhang ◽  
Dongxiao Liu ◽  
Guizhong Liu ◽  
Haiyang Gou ◽  
Rui Li ◽  
...  
Keyword(s):  
Least Squares ◽  
Iteration Process ◽  
Convergence Speed ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Imaging ◽  
And Performance ◽  
Key Techniques

An improved Least Squares Reverse Time Migration (LSRTM) method is proposed in the paper, which can effectively improve convergence speed and imaging accuracy. Firstly, the key techniques in the implementation of LSRTM are discussed. Secondly, a condition factor is introduced in the iteration process of conjugate gradient method. Finally, the imaging effect and performance of the algorithm are analyzed. The experiment results indicate that it can speed up the convergence speed and improve the convergence accuracy, so as to improve the imaging effect. Compared with the conventional LSRTM, the data residual of improved LSRTM can be reduced by about 5%.

scholarly journals Elastic least-squares reverse-time migration with density variation for flaw imaging in heterogeneous structures

2020 ◽  
Vol 29 (3) ◽  
pp. 035017 ◽  
Author(s):  
Jing Rao ◽  
Jizhong Yang ◽  
Jiaze He ◽  
Ming Huang ◽  
Ernst Rank
Keyword(s):  
Least Squares ◽  
Density Variation ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Heterogeneous Structures

Viscoacoustic least-squares reverse time migration using a time-domain complex-valued wave equation

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S479-S499 ◽  
Author(s):  
Jidong Yang ◽  
Hejun Zhu
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Time Domain ◽  
Inverse Operator ◽  
Total Variation Regularization ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Complex Valued

With limited recording apertures, finite-frequency source functions, and irregular subsurface illuminations, traditional imaging methods have been insufficient to produce satisfactory reflectivity images with high resolution and amplitude fidelity. This is because most traditional imaging approaches are commonly formulated as the adjoint instead of the inverse operator with respect to the forward-modeling operator. In addition, intrinsic attenuation introduces amplitude loss and phase dispersion during wave propagation. Without considering these effects, migrated images might be kinematically and dynamically incorrect. We have developed a viscoacoustic least-squares reverse time migration (LSRTM) method based on a time-domain complex-valued wave equation. According to the Born approximation, we first linearized the viscoacoustic wave equation and derived a demigration operator. Then, using the complex-valued Lagrange multiplier method, we derived the adjoint viscoacoustic wave equation and corresponding sensitivity kernel. With the forward and adjoint operators, a linear inverse problem is formulated to estimate the subsurface reflectivity model. A total-variation regularization scheme is introduced to enhance the robustness of our viscoacoustic LSRTM, and a diagonal Hessian is used as the preconditioner to accelerate the convergence. Three synthetic examples are used to demonstrate that our approach enables us to compensate attenuation effects, improve imaging resolution, and enhance amplitude fidelity in comparison with the adjoint imaging method.

Sparsity-promoting multi-parameter pseudoinverse Born inversion in acoustic media

Geophysics ◽  
2021 ◽  
pp. 1-62
Author(s):  
Milad Farshad ◽  
Hervé Chauris
Keyword(s):  
Least Squares ◽  
Cross Talk ◽  
Radon Transform ◽  
Weighted Least Squares ◽  
Variable Density ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Acoustic Media

Least-squares reverse-time migration has become the method of choice for quantitative seismic imaging. The main drawback of such scheme is that it requires many migration/modeling cycles. The convergence of least-squares reverse-time migration can be accelerated by using a suitable preconditioner. In the context of extended domain in a variable density acoustic media, the pseudoinverse Born operator is the recommended preconditioner, providing quantitative results within a single iteration. This method consists of two steps: application of the pseudoinverse Born operator, and inversion of two parameters using an efficient weighted least-squares approach based on the Radon transform. As expected, cross-talk artifacts are generated in the second step due to limited acquisition. We present a variable density pseudoinverse Born operator constrained with the ℓ1-norm for each model parameter to suppress the artifacts. The fast iterative shrinkage-thresholding algorithm is used to carry out the optimization problem. In classical iterative least-squares migration, the ℓ1-norm constraints would affect the whole imaging process. As the imaging method is split into two steps, only the Radon transform part is modified, where no wave-based operators are involved. Through numerical experiments, we verify the robustness of the proposed method against different migration artifacts including the parameter cross-talk, interfaces with abrupt truncations, sparse shot acquisition geometry, noisy data and high contrast complex structures.

Decoupled elastic least-squares reverse time migration and its application in tunnel geological forward-prospecting

Geophysics ◽  
2021 ◽  
pp. 1-136
Author(s):  
Bin Liu ◽  
Jiansen Wang ◽  
Yuxiao Ren ◽  
Xu Guo ◽  
Lei Chen ◽  
...  
Keyword(s):  
Least Squares ◽  
Tunnel Construction ◽  
Stable Convergence ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Imaging ◽  
Imaging Results ◽  
Geological Conditions ◽  
Migration Method

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.

scholarly journals A GPU implementation of least-squares reverse time migration

2021 ◽  
Vol 1719 (1) ◽  
pp. 012030
Author(s):  
Phudit Sombutsirinun ◽  
Chaiwoot Boonyasiriwat
Keyword(s):  
Least Squares ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Gpu Implementation
Sign in / Sign up

Export Citation Format

Share Document