Least-squares migration of multisource data with a deblurring filter

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. R135-R146 ◽  
Author(s):  
Wei Dai ◽  
Xin Wang ◽  
Gerard T. Schuster
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Computational Efficiency ◽  
Computational Cost ◽  
Processing Technique ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Method ◽  
Least Squares Migration

Least-squares migration (LSM) has been shown to be able to produce high-quality migration images, but its computational cost is considered to be too high for practical imaging. We have developed a multisource least-squares migration algorithm (MLSM) to increase the computational efficiency by using the blended sources processing technique. To expedite convergence, a multisource deblurring filter is used as a preconditioner to reduce the data residual. This MLSM algorithm is applicable with Kirchhoff migration, wave-equation migration, or reverse time migration, and the gain in computational efficiency depends on the choice of migration method. Numerical results with Kirchhoff LSM on the 2D SEG/EAGE salt model show that an accurate image is obtained by migrating a supergather of 320 phase-encoded shots. When the encoding functions are the same for every iteration, the input/output cost of MLSM is reduced by 320 times. Empirical results show that the crosstalk noise introduced by blended sources is more effectively reduced when the encoding functions are changed at every iteration. The analysis of signal-to-noise ratio (S/N) suggests that not too many iterations are needed to enhance the S/N to an acceptable level. Therefore, when implemented with wave-equation migration or reverse time migration methods, the MLSM algorithm can be more efficient than the conventional migration method.

Fast least-squares reverse time migration of VSP free-surface multiples with dynamic phase-encoding schemes

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. S321-S332 ◽  
Author(s):  
Xuejian Liu ◽  
Yike Liu ◽  
Majid Khan
Keyword(s):  
Free Surface ◽  
Least Squares ◽  
Computational Cost ◽  
Seismic Profile ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Multiple Reflections ◽  
Phase Encoding ◽  
Time Migration ◽  
Migration Method

For vertical seismic profile (VSP), free-surface multiples can provide much wider subsurface illumination when compared with primaries. However, migration of multiple reflections generates not only the desired image of reflection interfaces but also many crosstalk artifacts. Therefore, the least-squares reverse time migration method is used to image the VSP downgoing free-surface multiples (receiver-side ghosts) and iteratively suppress crosstalks, in which full downgoing data (including direct waves) and downgoing multiples are used as sources and observed data, respectively. To reduce the computational cost, we have developed the simultaneous imaging of different common-receiver gathers that are dynamically blended together with iterations through the altered realizations of the phase-encoding function. Relative to the popular encoding function with a combination of random time delays and polarities, only the random polarities can be applied for further increasing the computational efficiency. Synthetic experiments on Sigsbee2B and Pluto1.5 models indicate that the proposed method can effectively eliminate crosstalk artifacts and improve imaging resolution while calculated even more efficiently than reverse time migration of VSP ghosts.

From acoustic to elastic inverse extended Born modeling: a first insight in the marine environment

Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Milad Farshad ◽  
Hervé Chauris
Keyword(s):  
Least Squares ◽  
Marine Environment ◽  
Computational Cost ◽  
Physical Parameters ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
Physical Density ◽  
Acoustic Approximation

Elastic least-squares reverse time migration is the state-of-the-art linear imaging technique to retrieve high-resolution quantitative subsurface images. A successful application requires many migration/modeling cycles. To accelerate the convergence rate, various pseudoinverse Born operators have been proposed, providing quantitative results within a single iteration, while having roughly the same computational cost as reverse time migration. However, these are based on the acoustic approximation, leading to possible inaccurate amplitude predictions as well as the ignorance of S-wave effects. To solve this problem, we extend the pseudoinverse Born operator from acoustic to elastic media to account for the elastic amplitudes of PP reflections and provide an estimate of physical density, P- and S-wave impedance models. We restrict the extension to marine environment, with the recording of pressure waves at the receiver positions. Firstly, we replace the acoustic Green's functions by their elastic version, without modifying the structure of the original pseudoinverse Born operator. We then apply a Radon transform to the results of the first step to calculate the angle-dependent response. Finally, we simultaneously invert for the physical parameters using a weighted least-squares method. Through numerical experiments, we first illustrate the consequences of acoustic approximation on elastic data, leading to inaccurate parameter inversion as well as to artificial reflector inclusion. Then we demonstrate that our method can simultaneously invert for elastic parameters in the presence of complex uncorrelated structures, inaccurate background models, and Gaussian noisy data.

Q least-squares reverse time migration based on the first-order viscoacoustic quasidifferential equations

Geophysics ◽  
2021 ◽  
pp. 1-65
Author(s):  
Yingming Qu ◽  
Yixin Wang ◽  
Zhenchun Li ◽  
Chang Liu
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Surface Topography ◽  
Density Perturbation ◽  
Second Order ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
First Order ◽  
Time Migration ◽  
Grid Scheme

Seismic wave attenuation caused by subsurface viscoelasticity reduces the quality of migration and the reliability of interpretation. A variety of Q-compensated migration methods have been developed based on the second-order viscoacoustic quasidifferential equations. However, these second-order wave-equation-based methods are difficult to handle with density perturbation and surface topography. In addition, the staggered grid scheme, which has an advantage over the collocated grid scheme because of its reduced numerical dispersion and enhanced stability, works in first-order wave-equation-based methods. We have developed a Q least-squares reverse time migration method based on the first-order viscoacoustic quasidifferential equations by deriving Q-compensated forward-propagated operators, Q-compensated adjoint operators, and Q-attenuated Born modeling operators. Besides, our method using curvilinear grids is available even when the attenuating medium has surface topography and can conduct Q-compensated migration with density perturbation. The results of numerical tests on two synthetic and a field data sets indicate that our method improves the imaging quality with iterations and produces better imaging results with clearer structures, higher signal-to-noise ratio, higher resolution, and more balanced amplitude by correcting the energy loss and phase distortion caused by Q attenuation. It also suppresses the scattering and diffracted noise caused by the surface topography.

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.

Mini-batch Least-Squares Reverse Time Migration In A Deep Learning Framework

Geophysics ◽  
2020 ◽  
pp. 1-61
Author(s):  
Janaki Vamaraju ◽  
Jeremy Vila ◽  
Mauricio Araya-Polo ◽  
Debanjan Datta ◽  
Mohamed Sidahmed ◽  
...  
Keyword(s):  
Deep Learning ◽  
Least Squares ◽  
Spatial Resolution ◽  
Data Science ◽  
Computational Cost ◽  
Computation Time ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Learning Framework ◽  
Time Migration

Migration techniques are an integral part of seismic imaging workflows. Least-squares reverse time migration (LSRTM) overcomes some of the shortcomings of conventional migration algorithms by compensating for illumination and removing sampling artifacts to increase spatial resolution. However, the computational cost associated with iterative LSRTM is high and convergence can be slow in complex media. We implement pre-stack LSRTM in a deep learning framework and adopt strategies from the data science domain to accelerate convergence. The proposed hybrid framework leverages the existing physics-based models and machine learning optimizers to achieve better and cheaper solutions. Using a time-domain formulation, we show that mini-batch gradients can reduce the computation cost by using a subset of total shots for each iteration. Mini-batch approach does not only reduce source cross-talk but also is less memory intensive. Combining mini-batch gradients with deep learning optimizers and loss functions can improve the efficiency of LSRTM. Deep learning optimizers such as the adaptive moment estimation are generally well suited for noisy and sparse data. We compare different optimizers and demonstrate their efficacy in mitigating migration artifacts. To accelerate the inversion, we adopt the regularised Huber loss function in conjunction. We apply these techniques to 2D Marmousi and 3D SEG/EAGE salt models and show improvements over conventional LSRTM baselines. The proposed approach achieves higher spatial resolution in less computation time measured by various qualitative and quantitative evaluation metrics.

Least-squares reverse time migration in TTI media using a pure qP-wave equation

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. S199-S216
Author(s):  
Xinru Mu ◽  
Jianping Huang ◽  
Jidong Yang ◽  
Xu Guo ◽  
Yundong Guo
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Seismic Imaging ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
The Matrix ◽  
The Stability

Anisotropy is a common phenomenon in subsurface strata and should be considered in seismic imaging and inversion. Seismic imaging in a vertical transversely isotropic (VTI) medium does not take into account the effects of the tilt angles, which can lead to degraded migrated images in areas with strong anisotropy. To correct such waveform distortion, reduce related image artifacts, and improve migration resolution, a tilted transversely isotropic (TTI) least-squares reverse time migration (LSRTM) method is presented. In the LSRTM, a pure qP-wave equation is used and solved with the finite-difference method. We have analyzed the stability condition for the pure qP-wave equation using the matrix method, which is used to ensure the stability of wave propagation in the TTI medium. Based on this wave equation, we derive a corresponding demigration (Born modeling) and adjoint migration operators to implement TTI LSRTM. Numerical tests on the synthetic data show the advantages of TTI LSRTM over VTI RTM and VTI LSRTM when the recorded data contain strong effects caused by large tilt angles. Our numerical experiments illustrate that the sensitivity of the adopted TTI LSRTM to the migration velocity errors is much higher than that to the anisotropic parameters (including epsilon, delta, and tilted angle parameters), and its sensitivity to the epsilon model and tilt angle is higher than that to the delta model.

Full Wavefield Least-Squares Reverse Time Migration

Geophysics ◽  
2021 ◽  
pp. 1-48
Author(s):  
Mikhail Davydenko ◽  
Eric Verschuur
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Input Data ◽  
Waveform Inversion ◽  
Secondary Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Internal Multiples ◽  
The One

Waveform inversion based on Least-Squares Reverse Time Migration (LSRTM) usually involves Born modelling, which models the primary-only data. As a result the inversion process handles only primaries and corresponding multiple elimination pre-processing of the input data is required prior to imaging and inversion. Otherwise, multiples left in the input data are mapped as false reflectors, also known as crosstalk, in the final image. At the same time the developed Full Wavefield Migration (FWM) methodology can handle internal multiples in an inversion-based imaging process. However, because it is based on the framework of the one-way wave equation, it cannot image dips close to and beyond 90 degrees. Therefore, we aim at upgrading LSRTM framework by bringing functionality of FWM to handle internal multiples. We have discovered that the secondary source term, used in the original formulation of FWM to define a wavefield relationship that allows to model multiple scattering via reflectivity, can be injected into a pressure component when simulating the two-way wave equation using finite-difference modelling. We use this modified forward model for estimating the reflectivity model with automatic crosstalk supression and validate the method on both synthetic and field data containing visible internal multiples.

Least-squares reverse time migration method using the factorization of the Hessian matrix

2021 ◽  
Vol 18 (1) ◽  
pp. 94-100
Author(s):  
Sun Xiao-Dong ◽  
Teng Hou-Hua ◽  
Ren Li-Juan ◽  
Wang Wei-Qi ◽  
Li Zhen-Chun
Keyword(s):  
Least Squares ◽  
Hessian Matrix ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Method

Preconditioned acoustic least-squares two-way wave-equation migration with exact adjoint operator

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. S1-S13 ◽  
Author(s):  
Linan Xu ◽  
Mauricio D. Sacchi
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Adjoint Operator ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Product Test ◽  
Dot Product ◽  
The Time Domain ◽  
Adjoint Operators

We have investigated the problem of designing the forward operator and its exact adjoint for two-way wave-equation least-squares migration. We study the problem in the time domain and pay particular attention to the individual operators that are required by the algorithm. We derive our algorithm using the language of linear algebra and establish a simple path to design forward and adjoint operators that pass the dot-product test. We also found that the exact adjoint operator is not equal to the classic reverse time migration algorithm. For instance, one must pay particular attention to boundary conditions to compute the exact adjoint that accurately passes the dot-product test. Forward and adjoint operators are adopted to solve the so-called least-squares reverse time migration problem via the method of conjugate gradients. We also examine a preconditioning strategy to invert extended images.

Sign in / Sign up

Export Citation Format

Share Document