Adjoint-state reverse time migration of 4C data: Finite-frequency map migration for marine seismic imaging

Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. WA159-WA172 ◽  
Author(s):  
Clement Fleury ◽  
Ivan Vasconcelos

Recent advances in marine seismic acquisition allow for the recording of vector-acoustic ([VA] pressure and particle velocity) seismic data from dual-source configurations, i.e., using monopole as well as dipole sources. VA reverse time migration (RTM) can be custom designed to accurately handle amplitude and directivity information from 4C seismic data. We present a method for multicomponent RTM that is based on an adjoint-state formulation using the full VA wave equations for pressure and corresponding displacement fields. This method takes advantage of the directional finite-frequency information contained in the 4C acoustic fields by using source and receiver weighting operators in the adjoint-state imaging scheme. With this adjoint-state method, the source and receiver radiation properties are tailored by choosing specific weighting operators. Weighting operators were chosen so that source- and receiver-side ghost arrivals are jointly migrated with primary energy. Because the dipole field components (e.g., components of particle displacement or acceleration) are proportional to the spatial gradient components of the pressure field, our method is in fact a formulation for reverse-time map migration that images pressure fields while jointly using the directional information contained in its full 3C gradients. As a result, our reverse time 4C map migration method yields less aperture- and sampling-related artifacts when compared to imaging of the pressure-only or 2C seismic data. In addition, our method sets a framework for full-waveform inversion using dual-source 4C seismic data. We demonstrated our findings with synthetic data, including a subsalt imaging example.

Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. WA123-WA145 ◽  
Author(s):  
Ivan Vasconcelos

Novel technologies in seismic data acquisition allow for recording full vector-acoustic (VA) data: pointwise recordings of pressure and its multicomponent gradient, excited by pressure only as well as dipole/gradient sources. Building on recent connections between imaging and seismic interferometry, we present a wave-equation-based, nonlinear, reverse-time imaging approach that takes full advantage of dual-source multicomponent data. The method’s formulation relies on source-receiver scattering reciprocity, thus making proper use of VA fields in the wavefield extrapolation and imaging condition steps in a self-consistent manner. The VA imaging method is capable of simultaneously focusing energy from all in- and outgoing waves: The receiver-side up- and downgoing (receiver ghosts) fields are handled by the VA receiver extrapolation, whereas source-side in- and outgoing (source ghosts) arrivals are accounted for when combining dual-source data at the imaging condition. Additionally, VA imaging handles image amplitudes better than conventional reverse-time migration because it properly handles finite-aperture directivity directly from dual-source, 4C data. For nonlinear imaging, we provide a complete source-receiver framework that relies only on surface integrals, thus being computationally applicable to practical problems. The nonlinear image can be implicitly interpreted as a superposition of several nonlinear interactions between scattering components of data with those corresponding to the extrapolators (i.e., to the model). We demonstrate various features of the method using synthetic examples with complex subsurface features. The numerical results show, e.g., that the dual-source, VA image retrieves subsurface features with “super-resolution”, i.e., with resolution higher than the limits of Born imaging, but at the cost of introducing image artifacts not present in the linear image. Although the method does not require any deghosting as a preprocessing step, it can use separated up- and downgoing fields to generate independent subsurface images.


Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. T287-T300 ◽  
Author(s):  
Lasse Amundsen ◽  
Johan O. A. Robertsson

Methods for wavefield injection are used in, for instance, reverse time extrapolation of shot gathers in reverse time migration. For correct injection of recorded data without any ambiguity of the propagation direction, the wavefield-injection methodology requires pressure and particle velocity data such as multicomponent towed marine or seabed seismic recordings. We discovered that by carefully considering the models (medium parameters and boundary conditions) for injection, wavefield injection of multicomponent data can also be used to solve several long-standing challenges in marine seismic data processing by means of conventional time-space-domain finite-difference propagators. We outlined and demonstrated several of these important applications including up-down separation of wavefields (deghosting), direct-wave removal, source-signature estimation, multiple removal, and imaging using primaries and multiples. Only acoustic models are considered, but the concepts are straightforward to generalize to elastodynamic and electromagnetic models.


2007 ◽  
Author(s):  
Marco A. Barsottelli Botelho* ◽  
Antônio V. Moura Lima

2018 ◽  
Vol 8 (2) ◽  
Author(s):  
Juan Guillermo Paniagua Castrillón ◽  
Olga Lucía Quintero- Montoya

Low-frequency artifacts in reverse time migration result from unwanted cross-correlation of the source and receiver wavefields at non-reflecting points along ray-paths. These artifacts can hide important details in migrated models and increase poor interpretation risk. Some methods have been proposed to avoid or reduce the number of these artifacts, preserving reflections, and improving model quality, implementing other strategies such as modification of the wave equation, proposing other imaging conditions, and using image filtering techniques. One of these methods uses wavefield decomposition, correlating components of the wavefields that propagate in opposite directions. We propose a method for extracting directional information from the RTM imaging condition wavefields to obtain characteristics allowing for better, more refined imaging. The method works by separating directional information about the wavefields based on the continuous wavelet transform (CWT), and the analysis of the main changes on the frequency content revealed within the scalogram obtained by a Gaussian wavelet family. Through numerical applications, we demonstrate that this method can effectively remove undesired artifacts in migrated images. In addition, we use the Laguerre-Gauss filtering to improve the results obtained with the proposed method.


2021 ◽  
Author(s):  
Hala Alqatari ◽  
Thierry-Laurent Tonellot ◽  
Mohammed Mubarak

Abstract This work presents a full waveform sonic (FWS) dataset processing to generate high-resolution images of the near-borehole area. The dataset was acquired in a nearly horizontal well over a distance of 5400 feet. Multiple formation boundaries can be identified on the final image and tracked at up to 200 feet deep, along the wellbore's trajectory. We first present a new preprocessing sequence to prepare the sonic data for imaging. This sequence leverages denoising algorithms used in conventional surface seismic data processing to remove unwanted components of the recorded data that could harm the imaging results. We then apply a reverse time migration algorithm to the data at different processing stages to assess the impact of the main processing steps on the final image.


2021 ◽  
Author(s):  
Pavlo Kuzmenko ◽  
Viktor Buhrii ◽  
Carlo D'Aguanno ◽  
Viktor Maliar ◽  
Hrigorii Kashuba ◽  
...  

Abstract Processing of the seismic data acquired in areas of complex geology of the Dnieper-Donets basin, characterized by the salt tectonics, requires special attention to the salt dome interpretation. For this purpose, Kirchhoff Depth Imaging and Reverse Time Migration (RTM) were applied and compared. This is the first such experience in the Dnieper-Donets basin. According to international experience, RTM is the most accurate seismic imaging method for steep and vertical geological (acoustic contrast) boundaries. Application of the RTM on 3D WAZ land data is a great challenge in Dnieper-Donets Basin because of the poor quality of the data with a low signal-to-noise ratio and irregular spatial sampling due to seismic acquisition gaps and missing traces. The RTM algorithm requires data, organized to native positions of seismic shots. For KPSDM we used regularized data after 5D interpolation. This affects the result for near salt reflection. The analysis of KPSDM and RTM results for the two areas revealed the same features. RTM seismic data looked more smoothed, but for steeply dipping reflections, lateral continuity of reflections was much improved. The upper part (1000 m) of the RTM has shadow zones caused by low fold. Other differences between Kirchhoff data and RTM are in the spectral content, as the former is characterized by the full range of seismic frequency spectrum. Conversely, beneath the salt, the RTM has reflections with steep dips which are not observed on the KPSDM. It is possible to identify new prospects using the RTM seismic image. Reverse Time Migration of 3D seismic data has shown geologically consistent results and has the potential to identify undiscovered hydrocarbon traps and to improve salt flank delineation in the complex geology of the Dnieper-Donets Basin's salt domes.


Sign in / Sign up

Export Citation Format

Share Document