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.

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.

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.

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

scholarly journals Plane-wave least-squares reverse-time migration

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. S165-S177 ◽  
Author(s):  
Wei Dai ◽  
Gerard T. Schuster
Keyword(s):  
Phase Shift ◽  
Plane Wave ◽  
Least Squares ◽  
Plane Waves ◽  
Migration Velocity ◽  
Stable Convergence ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Time Migration

A plane-wave least-squares reverse-time migration (LSRTM) is formulated with a new parameterization, where the migration image of each shot gather is updated separately and an ensemble of prestack images is produced along with common image gathers. The merits of plane-wave prestack LSRTM are the following: (1) plane-wave prestack LSRTM can sometimes offer stable convergence even when the migration velocity has bulk errors of up to 5%; (2) to significantly reduce computation cost, linear phase-shift encoding is applied to hundreds of shot gathers to produce dozens of plane waves. Unlike phase-shift encoding with random time shifts applied to each shot gather, plane-wave encoding can be effectively applied to data with a marine streamer geometry. (3) Plane-wave prestack LSRTM can provide higher-quality images than standard reverse-time migration. Numerical tests on the Marmousi2 model and a marine field data set are performed to illustrate the benefits of plane-wave LSRTM. Empirical results show that LSRTM in the plane-wave domain, compared to standard reverse-time migration, produces images efficiently with fewer artifacts and better spatial resolution. Moreover, the prestack image ensemble accommodates more unknowns to makes it more robust than conventional least-squares migration in the presence of migration velocity errors.

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 3-D Multi-Component Reverse Time Migration Method for Tunnel Seismic Data

Sensors ◽  
2021 ◽  
Vol 21 (9) ◽  
pp. 3244
Author(s):  
Peng Guan ◽  
Cuifa Shao ◽  
Yuyong Jiao ◽  
Guohua Zhang ◽  
Bin Li ◽  
...  
Keyword(s):  
Seismic Data ◽  
Cross Correlation ◽  
Measured Data ◽  
Ray Theory ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Kirchhoff Migration ◽  
Time Migration ◽  
Migration Imaging ◽  
Migration Method

Migration imaging is a key step in tunnel seismic data processing. Due to the limitation of tunnel space, tunnel seismic data are small-quantity, multi-component, and have a small offset. Kirchhoff migration based on the ray theory is limited to the migration aperture and has low migration imaging accuracy. Kirchhoff migration can no longer meet the requirements of high-precision migration imaging. The reverse time migration (RTM) method is used to realize cross-correlation imaging by reverse-time recursion principle of the wave equation. The 3-D RTM method cannot only overcome the effect of small offset, but also realize multi-component data imaging, which is the most accurate migration method for tunnel seismic data. In this paper, we will study the 3-D RTM method for multi-component tunnel seismic data. Combined with the modeled data and the measured data, the imaging accuracy of the 3-D Kirchhoff migration and 3-D RTM is analyzed in detail. By comparing single-component and multi-component Kirchhoff migration and RTM profile, the advantages of the multi-component RTM method are summarized. Compared with the Kirchhoff migration method, the 3-D RTM method has the following advantages: (1) it can overcome the effect of small offset and expand the range of migration imaging; (2) multi-component data can be realized to improve the energy of anomalous interface; (3) it can make full use of multiple waves to realize migration imaging and improve the resolution of the anomalous interface. The modeled data and the measured data prove the advantages of the 3-D multi-component RTM method.

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.

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