scholarly journals Depth Migration Based on Two-Way Wave Equation to Image OBS Multiples: A Case Study in the South Shetland Margin (Antarctica)

Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-9
Author(s):  
Sha Song ◽  
Jiachun You ◽  
Qing Cao ◽  
Bin Chen ◽  
Xiaomeng Cao

With the development of marine seismic exploration, the ocean bottom seismometer (OBS) as a new seismic acquisition technology has been widely concerned. Although multiple waves are frequently viewed as noises, they may carry a wealth of subsurface information and produce a broader illumination than primary waves. To perform multiple wave imaging, we propose to utilize a two-way wave equation depth wavefield extrapolation method which is rarely used in this field. A simple dipping model is imaged by using primary and multiple waves, which proves the superiority of multiple waves in imaging over the primary waves and lays a foundation for practical application. Moreover, the comparison of multiple imaging results by reverse time migration and those by our proposed method demonstrates that our proposed method requires less storage space. In this study, we apply this migration method to actual OBS data collected in the South Shetland margin (Antarctica), where gas hydrates have been well documented. Firstly, the wavefield separation method is adopted to process the OBS data, so as to produce reliable primary and multiples waves; secondly, the ray-tracing method is used to derive the velocity field; and finally, the depth wavefield extrapolation method based on the two-way wave equation is applied to image primary and multiple waves. Migration results show that multiple waves provide a broader illumination and a clearer sediment structure than primary waves, especially for the highly shallow reflections.

2021 ◽  
Vol 11 (7) ◽  
pp. 3010
Author(s):  
Hao Liu ◽  
Xuewei Liu

The lack of an initial condition is one of the major challenges in full-wave-equation depth extrapolation. This initial condition is the vertical partial derivative of the surface wavefield and cannot be provided by the conventional seismic acquisition system. The traditional solution is to use the wavefield value of the surface to calculate the vertical partial derivative by assuming that the surface velocity is constant. However, for seismic exploration on land, the surface velocity is often not uniform. To solve this problem, we propose a new method for calculating the vertical partial derivative from the surface wavefield without making any assumptions about the surface conditions. Based on the calculated derivative, we implemented a depth-extrapolation-based full-wave-equation migration from topography using the direct downward continuation. We tested the imaging performance of our proposed method with several experiments. The results of the Marmousi model experiment show that our proposed method is superior to the conventional reverse time migration (RTM) algorithm in terms of imaging accuracy and amplitude-preserving performance at medium and deep depths. In the Canadian Foothills model experiment, we proved that our method can still accurately image complex structures and maintain amplitude under topographic scenario.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. S137-S155 ◽  
Author(s):  
Matteo Ravasi ◽  
Andrew Curtis

With the more widespread introduction of multicomponent recording devices in land and marine ocean-bottom seismic acquisition, elastic imaging may become mainstream in coming years. We have derived new, nonlinear, elastic imaging conditions. A correlation-type representation theorem for perturbed elastic media, commonly used in seismic interferometry to explain how a scattered wave response between two receivers/sources may be predicted given a boundary of sources/receivers, can be considered as a starting point for the derivation. Here, we use this theorem to derive and interpret imaging conditions for elastic migration by wavefield extrapolation (e.g., elastic reverse-time migration). Some approximations lead to a known, heuristically derived imaging condition that crosscorrelates P- and S-wave potentials that are separated in the subsurface after full-wavefield extrapolation. This formal connection reveals that the nonapproximated correlation-type representation theorem can be interpreted as a nonlinear imaging condition, that accounts also for multiply scattered and multiply converted waves, properly focusing such energy at each image point. We present a synthetic data example using either an ideal (acquisition on a full, closed boundary) or a real (partial boundary) seismic exploration survey, and we demonstrate the importance of nonlinearities in pure- and converted-mode imaging. In PP imaging, they result in better illumination and artifact reduction, whereas in PS imaging they show how zero time-lag and zero space-lag crosscorrelation imaging conditions are not ideal for imaging of converted-mode waves because no conversion arises from zero-offset experiments.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S557-S567 ◽  
Author(s):  
Yan Zhao ◽  
Ningbo Mao ◽  
Zhiming Ren

Amplitude energy attenuation and phase distortion of seismic waves caused by formation viscoelasticity reduce the resolution of reverse time migration (RTM) images. Q-RTM often is used to compensate the attenuation effects and improve the resolution of seismic imaging. However, serious high-frequency noise and tremendous amplitude will be produced during the wavefield extrapolation of Q-RTM, resulting in its inability to be imaged. Many Q-RTM algorithms solve the problem of instability through low-pass filtering in the wavenumber domain, but the method is less efficient in computation and has a truncation effect in the wavefield. We have developed a stable and efficient Q-RTM method, in which a regularization term was introduced into the viscoacoustic wave equation to suppress the high-frequency noise, and the finite-difference method was used to solve the viscoacoustic wave equation with a regularization term. We used the model example to visually demonstrate the instability of wavefield extrapolation in Q-RTM and compared the effect and computational efficiency of the two stabilization processing methods, low-pass filtering and regularization. Meanwhile, our method is not involved in solving the fractional derivatives by using the pseudo-spectral method, the computational efficiency also can be improved. We tested the Q-RTM approach on a simple layered model, Marmousi model, and real seismic data. The results of numerical examples demonstrated that the Q-RTM method can solve the problem of instability effectively and obtain a higher resolution image with lower computational cost.


Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. A29-A33 ◽  
Author(s):  
Yu Zhang ◽  
Guanquan Zhang

We have proposed a new method, a one-step extrapolation algorithm, to solve the acoustic wave equation. By introducing a square-root operator, the two-way wave equation can be formulated as a first-order partial differential equation in time, which is similar to the one-way wave equation. To solve the new wave equation, we used a stable explicit extrapolation method in the time direction and handled lateral velocity variations in the space and wavenumber domains. Unlike the conventional explicit finite-difference schemes, the new method does not suffer from numerical instability or numerical dispersion problems. It can be used to design cost-effective and high-quality reverse time migration or modeling code.


Geophysics ◽  
2021 ◽  
pp. 1-65
Author(s):  
Yingming Qu ◽  
Yixin Wang ◽  
Zhenchun Li ◽  
Chang Liu

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.


2018 ◽  
Vol 6 (4) ◽  
pp. SM39-SM50
Author(s):  
Jingbo Wang ◽  
Zhongshan Qi ◽  
Penggui Jing ◽  
Tianfa Zheng ◽  
Yanqi Li ◽  
...  

Geologic studies indicate that the platform-margin reef-shallow facies in Permo-Triassic marine strata in the Micang-Dabashan foothill belt in the Sichuan Basin are favorable exploration targets for oil and gas exploration. However, the typical dual-complexity problem (complex surface condition and subsurface structure) brings a great challenge for seismic technology targeting of those potential oil and gas reservoirs. To overcome this problem, varieties of advanced seismic acquisition and processing methods have been used to improve the imaging quality of piedmont seismic data since 2000. Some improvements have been achieved: The reflection waves from the far offset and deep layer can be acquired in shot gathers from limestone outcropped areas, and the signal-to-noise ratio (S/N) of reflection and diffraction waves in the stack section has been enhanced significantly so as to reveal amounts of valuable geologic information. The resolution and the S/N of seismic migration imaging for the strong fold zone in marine strata have been improved partially, so that the structure of the step-fault zone and the enveloping of gypsum rock are clearer than those revealed by the old seismic section. Even so, actual drilling data demonstrate that the subsurface structures of the foothill belt are far more complex than those revealed by the current seismic imaging results. Therefore, postdrilling evaluation for the validity of seismic techniques implemented in the Nanjiang and Zhenba piedmont zone has been carried out. The results indicate that the current acquisition scheme and processing workflow cannot completely fulfill the requirements of high-precision velocity modeling and migration imaging of complex structures (such as footwalls of thrust fault and small-scale fault blocks) in the piedmont zone, especially when the rugged surface and the widespread limestone outcrop appear simultaneously. Finally, we have developed some potential needs of seismic theories and techniques in the foothill belt, including seismic wave propagation, acquisition, and processing technology.


2017 ◽  
Vol 10 (1) ◽  
pp. 1-21 ◽  
Author(s):  
Zhiming Chen ◽  
Guanghui Huang

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


Sign in / Sign up

Export Citation Format

Share Document