scholarly journals Data-driven wavefield focusing and imaging with multidimensional deconvolution: Numerical examples for reflection data with internal multiples

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. WA107-WA115 ◽  
Author(s):  
Filippo Broggini ◽  
Roel Snieder ◽  
Kees Wapenaar

Standard imaging techniques rely on the single scattering assumption. This requires that the recorded data do not include internal multiples, i.e., waves that have bounced multiple times between reflectors before reaching the receivers at the acquisition surface. When multiple reflections are present in the data, standard imaging algorithms incorrectly image them as ghost reflectors. These artifacts can mislead interpreters in locating potential hydrocarbon reservoirs. Recently, we introduced a new approach for retrieving the Green’s function recorded at the acquisition surface due to a virtual source located at depth. We refer to this approach as data-driven wavefield focusing. Additionally, after applying source-receiver reciprocity, this approach allowed us to decompose the Green’s function at a virtual receiver at depth in its downgoing and upgoing components. These wavefields were then used to create a ghost-free image of the medium with either crosscorrelation or multidimensional deconvolution, presenting an advantage over standard prestack migration. We tested the robustness of our approach when an erroneous background velocity model is used to estimate the first-arriving waves, which are a required input for the data-driven wavefield focusing process. We tested the new method with a numerical example based on a modification of the Amoco model.

Author(s):  
K. Wapenaar ◽  
E. Slob ◽  
F. Broggini ◽  
R. Snieder ◽  
J. Thorbecke ◽  
...  

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. Q23-Q37 ◽  
Author(s):  
Satyan Singh ◽  
Roel Snieder

Recent papers show that imaging with the retrieved Green’s function constructed by the Marchenko equations, called Marchenko imaging, reduces artifacts from internal and free-surface multiples compared with standard imaging techniques. Even though artifacts are reduced, they can still be present in the image, depending on the imaging condition used. We have found that when imaging with the up- and downgoing Green’s functions, the multidimensional deconvolution (MDD) imaging condition yields better images than correlation and deconvolution. “Better” in this case means improved resolution, fewer artifacts, and a closer match with the true reflection coefficient of the model. We have determined that the MDD imaging condition only uses primaries to construct the image, whereas multiples are implicitly subtracted in the imaging step. Consequently, combining the first arrival of the downgoing Green’s function with the complete upgoing Green’s function produces superior (or at least equivalent) images than using the one-way Green’s functions because the first arrival of the downgoing Green’s function excludes all the downgoing multiply reflected waves. We also find that standard imaging algorithms which use the redatumed reflection response, constructed with the one-way Green’s functions, produce images with reduced artifacts from multiples compared with standard imaging conditions, which use surface reflection data. All imaging methods that rely on the Marchenko equations require the same inputs as standard imaging techniques: the reflection response at the surface and a smooth estimate of the subsurface velocities.


Geophysics ◽  
2021 ◽  
pp. 1-59
Author(s):  
Evert Slob ◽  
Lele Zhang ◽  
Eric Verschuur

Marchenko multiple elimination schemes are able to attenuate all internal multiple reflections in acoustic reflection data. These can be implemented with and without compensation for two-way transmission effects in the resulting primary reflection dataset. The methods are fully automated and run without human intervention, but require the data to be properly sampled and pre-processed. Even when several primary reflections are invisible in the data because they are masked by overlapping primaries, such as in the resonant wedge model, all missing primary reflections are restored and recovered with the proper amplitudes. Investigating the amplitudes in the primary reflections after multiple elimination with and without compensation for transmission effects shows that transmission effects are properly accounted for in a constant velocity model. When the layer thickness is one quarter of the wavelength at the dominant frequency of the source wavelet, the methods cease to work properly. Full wavefield migration relies on a velocity model and runs a non-linear inversion to obtain a reflectivity model which results in the migration image. The primary reflections that are masked by interference with multiples in the resonant wedge model, are not recovered. In this case, minimizing the data misfit function leads to the incorrect reflector model even though the data fit is optimal. This method has much lower demands on data sampling than the multiple elimination schemes, but is prone to get stuck in a local minimum even when the correct velocity model is available. A hybrid method that exploits the strengths of each of these methods could be worth investigating.


2014 ◽  
Vol 199 (3) ◽  
pp. 1367-1371 ◽  
Author(s):  
Kees Wapenaar ◽  
Evert Slob

Abstract Recent work on the Marchenko equation has shown that the scalar 3-D Green's function for a virtual source in the subsurface can be retrieved from the single-sided reflection response at the surface and an estimate of the direct arrival. Here, we discuss the first steps towards extending this result to multicomponent data. After introducing a unified multicomponent 3-D Green's function representation, we analyse its 1-D version for elastodynamic waves in more detail. It follows that the main additional requirement is that the multicomponent direct arrival, needed to initiate the iterative solution of the Marchenko equation, includes the forward-scattered field. Under this and other conditions, the multicomponent Green's function can be retrieved from single-sided reflection data, and this is demonstrated with a 1-D numerical example.


Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. S241-S250 ◽  
Author(s):  
Yi Luo ◽  
Qinglin Liu ◽  
Yuchun E. Wang ◽  
Mohammed N. AlFaraj

We illustrate the use of mode-converted transmitted (e.g., PS- or SP-) waves in vertical seismic profiling (VSP) data for imaging areas above receivers where reflected waves cannot illuminate. Three depth-domain imaging techniques — move-out correction, common-depth-point (CDP) mapping, and prestack migration — are described and used for imag-ing the transmitted waves. Moveout correction converts an offset VSP trace into a zero-offset trace. CDP mapping maps each sample on an input trace to the location where the mode conversion occurs. For complex media, prestack migration (e.g., reverse-time migration) is used. By using both synthetic and field VSP data, we demonstrate that images derived from transmissions complement those from reflections. As an important application, we show that transmitted waves can illuminate zones above highly de-viated or horizontal wells, a region not imaged by reflection data. Because all of these benefits are obtained without extra data acquisition cost, we believe transmission imag-ing techniques will become widely adopted by the oil in-dustry.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. A13-A17 ◽  
Author(s):  
Valeri Korneev ◽  
Andrey Bakulin

The virtual source method (VSM) has been proposed as a practical approach to reduce distortions of seismic images caused by shallow, heterogeneous overburden. VSM is demanding at the acquisition stage because it requires placing downhole geophones below the most complex part of the heterogeneous overburden. Where such acquisition is possible, however, it pays off later at the processing stage because it does not require knowledge of the velocity model above the downhole receivers. This paper demonstrates that VSM can be viewed as an application of the Kirchhoff-Helmholtz integral (KHI) with an experimentally measured Green’s function. Direct measurement of the Green’s function ensures the effectiveness of the method in highly heterogeneous subsurface conditions.


Sign in / Sign up

Export Citation Format

Share Document