Multiple prediction and subtraction from apparent slowness relations in 2D synthetic and field ocean-bottom cable data

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. V89-V99 ◽  
Author(s):  
Juanjuan Cao ◽  
George McMechan

A target-oriented algorithm is developed for the prediction of multiples recorded on ocean-bottom cables by utilizing apparent slowness relations in common-source and common-receiver gathers. It is based on combining offsets and times of direct waves and primary reflections to predict multiples by matching apparent slownesses at all source and receiver locations; all higher-order multiples can be predicted by matching apparent slownesses alternately in common-source and common-receiver gathers. No knowledge of the subsurface velocity is required. Traveltimes of the direct waves and primary reflections need to be picked from common-source gathers. The subtraction of multiples involves flattening the predicted times of the multiple events, subtracting a local spatial average trace from each trace, within a fixed time window containing the wavelet of the multiple, and then shifting the data back to its original times. Tests of synthetic and field data indicate that the proposed method predicts multiples very well and removes them from seismic data efficiently with negligible affect on the primary reflections, as long as the primary and multiple reflections do not overlap in time and slowness over substantial windows in the domain in which the removal is done.

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. B311-B324 ◽  
Author(s):  
Laura Gassner ◽  
Tobias Gerach ◽  
Thomas Hertweck ◽  
Thomas Bohlen

Evidence for gas-hydrate occurrence in the Western Black Sea is found from seismic measurements revealing bottom-simulating reflectors (BSRs) of varying distinctness. From an ocean-bottom seismic data set, low-resolution traveltime-tomography models of P-wave velocity [Formula: see text] are constructed. They serve as input for acoustic full-waveform inversion (FWI), which we apply to derive high-resolution parameter models aiding the interpretation of the seismic data for potential hydrate and gas deposits. Synthetic tests indicate the applicability of the FWI approach to robustly reconstruct [Formula: see text] models with a typical hydrate and gas signature. Models of S-wave velocity [Formula: see text] containing a hydrate signature can only be reconstructed when the parameter distribution of [Formula: see text] is already well-known. When we add noise to the modeled data to simulate field-data conditions, it prevents the reconstruction of [Formula: see text] completely, justifying the application of an acoustic approach. We invert for [Formula: see text] models from field data of two parallel profiles of 14 km length with a distance of 1 km. Results indicate a characteristic velocity trend for hydrate and gas occurrence at BSR depth in the first of the analyzed profiles. We find no indications for gas accumulations below the BSR on the second profile and only weak indications for hydrate. These differences in the [Formula: see text] signature are consistent with the reflectivity behavior of the migrated seismic streamer data of both profiles in which a zone of high-reflectivity amplitudes is coincident with the potential gas zone derived from the FWI result. Calculating saturation estimates for the potential hydrate and gas zones yields values of up to 30% and 1.2%, respectively.


Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. S439-S452 ◽  
Author(s):  
Matteo Ravasi

Marchenko redatuming is a revolutionary technique to estimate Green’s functions from virtual sources in the subsurface using only data measured at the earth’s surface, without having to place either sources or receivers in the subsurface. This goal is achieved by crafting special wavefields (so-called focusing functions) that can focus energy at a chosen point in the subsurface. Despite its great potential, strict requirements on the reflection response such as knowledge and accurate deconvolution of the source wavelet (including absolute scaling factor) and co-location of sources and receivers have so far challenged the application of Marchenko redatuming to real-world scenarios. I combine the coupled Marchenko equations with a one-way version of the Rayleigh integral representation to obtain a new redatuming scheme that handles internal as well as free-surface multiples using dual-sensor, band-limited seismic data (with an unknown source signature) from any acquisition system that presents arbitrarily located sources above a line of regularly sampled receivers—for example, ocean-bottom, source-over-spread streamer, and horizontal borehole seismic data. The redatuming scheme is validated using synthetic and field data, and the retrieved subsurface wavefields are used for improved structural imaging and taken as input for the computation of true-amplitude angle gathers, which can lead to more accurate amplitude-versus-angle interpretation and velocity analysis.


2014 ◽  
Vol 32 (3) ◽  
pp. 395 ◽  
Author(s):  
Silmara L.R. Oliveira ◽  
Rosângela Corrêa Maciel ◽  
Michelângelo G. da Silva ◽  
Milton José Porsani

ABSTRACT. Short-period multiples attenuation is a difficult problem for shallow water marine seismic data processing. In the past few decades many filteringmethods have been developed to solve this problem and to improve the quality of seismic imaging. The Wiener-Levinson predictive deconvolution method is one of themost useful and well known filter methods used in the seismic data processing flow. It is a statistical approach to reduce redundancy along the time variable seismictrace, allowing us to both improve the time resolution and also attenuate multiple reflections of the seismic traces. One of the assumptions of the Wiener-Levinsonmethod is that the seismic wavelet is stationary along the entire seismic trace. However, this is not true for real seismic data and to bypass this limitation the methodis normally applied using fixed time windows, distributed along the seismic trace. The present study tested a new adaptive predictive deconvolution approach for theattenuation of short-period multiples. The new approach is based on a sliding window of fixed length that is shifted sample by sample along the entire seismic trace.At each position, a new filter is computed and applied. The implied systems of equations are solved by using a recursive Levinson-type algorithm. The main differencewith respect to the conventional Wiener-Levinson approach is that the filter is updated for each data sample along the trace and no assumption is imposed on the dataoutside the considered window. The new adaptive predictive deconvolution approach was tested using a seismic line of the Jequitinhonha Basin acquired by Petrobras.The results demonstrated that the new approach is very precise for the attenuation of short-period multiples, producing better results than the ones obtained fromthe conventional Wiener-Levinson predictive deconvolution approach. The results were obtained with filters of 25 coefficients, predictive distance of 5 samples andwindow length equal to 55 samples.Keywords: seismic processing, Jequitinhonha Basin, adaptive predictive deconvolution, multiple of attenuation, Wiener-Levinson deconvolution.RESUMO. A atenuação de reflexões múltiplas de curto período, presentes nos dados sísmicos adquiridos sobre lâmina d’água rasa, representa um grande problemado processamento de dados sísmicos marítimos. Nas últimas décadas, vários métodos de filtragem de dados sísmicos têm sido desenvolvidos com o propósito deatenuar reflexões múltiplas e melhorar a qualidade das seções sísmicas. O método de filtragem conhecido como deconvolução preditiva de Wiener-Levinson é bastante utilizado na indústria do petróleo. Ele permite melhorar a resolução temporal dos dados sísmicos e atenuar reflexões múltiplas, podendo ser visto como um método estatístico que remove a coerência temporal dos traços sísmicos. O método de Wiener-Levinson pressupõe que o pulso sísmico é estacionário, fato este que não ocorrenos dados sísmicos reais. Para contornar este problema, o método de Wiener-Levinson é normalmente aplicado utilizando-se janelas de tempo fixas, distribuídas ao longo do tempo de registro. No presente trabalho, empregamos um método de deconvolução preditiva adaptativa no qual as janelas de tempo deslizantes são deslocadas amostra a amostra ao longo de todo o traço sísmico. Os sistemas de equações são resolvidos com o algoritmo recursivo tipo-Levinson. Na deconvolução de Wiener-Levinson, com janelas de tempo fixa, os filtros são gerados e aplicados dentro de cada janela. Já na deconvolução preditiva adaptativa o algoritmo calcula um novo filtro a cada posição da janela deslizante. Para teste da nova abordagem utilizamos os dados sísmicos da Bacia de Jequitinhonha, cedidos pela Petrobras. Os melhores resultados foram obtidos com filtros de 25 coeficientes, distância de predição igual a 5 amostras e janela móvel de 55 amostras. Os resultados obtidos com a nova abordagem demonstram que a deconvolução preditiva adaptativa atua com eficácia na atenuação de múltiplas de curto período, apresentando resultados melhores queos gerados pelo método de deconvolução preditiva de Wiener-Levinson.Palavras-chave: processamento sísmico, Bacia do Jequitinhonha, deconvolução adaptativa, atenuação de múltiplas, deconvolução de Wiener-Levinson.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. S409-S418 ◽  
Author(s):  
M. Javad Khoshnavaz ◽  
Andrej Bóna ◽  
Aleksander Dzunic ◽  
Kevin Ung ◽  
Milovan Urosevic

Seismic imaging techniques often require an input velocity model. Velocity analysis is one of the most critical stages in seismic data processing. Standard ways to find the velocity model from seismic data in the time domain are constant velocity stack and semblance velocity analysis that may be time consuming and labor intensive. Oriented/velocity-less imaging using local event slopes is an alternative to the conventional imaging techniques. In some previous oriented techniques, seismic data must be sorted in two different domains, whereas seismic data are not always available in both domains and the use of interpolation is inevitable in such cases. Other methods are developed in terms of the higher order derivatives of traveltime with respect to offset, whereas estimation of the higher order derivatives is difficult to achieve with the required accuracy. We addressed the limitations by developing an oriented local slope based prestack time migration technique in only one domain: the common-source domain. The migration technique is developed for reflectors with small curvature. In the proposed approach, the need for the estimation of higher order derivatives is replaced by a point-to-point mapping of seismic data using the predictive painting technique. The theoretical contents of the proposed technique are tested on a simple synthetic data example and applied to a field data set.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. E25-E33 ◽  
Author(s):  
Anatoly Baumstein ◽  
Mohamed T. Hadidi

The wide success of 2D surface-related multiple elimination (SRME) in attenuating complex multiples in many cases has spurred efforts to apply the method in three dimensions. However, application of 3D SRME to conventional marine data is often impeded by severe crossline aliasing characteristic of marine acquisition geometries. We propose to overcome this limitation using a dip-moveout (DMO)-based procedure consisting of the following steps: resorting the data into common offsets to improve crossline sampling, performing DMO to eliminate azimuth variations in the common-offset domain, and efficiently implementing inverse shot-record DMO to reconstruct densely sampled shot records required for 3D SRME to predict multiples correctly. We use a field data example to demonstrate that the proposed shot reconstruction procedure leads to kinematically accurate reconstruction of primaries but may not be able to simultaneously position multiples correctly. The mispositioning of multiples becomes a problem when second- and higher-order multiples must be predicted. We propose to resolve this difficulty by using a layer-stripping approach to multiple prediction. Alternatively, an approximate algorithm that relies on adaptive subtraction to compensate for inaccurate positioning of predicted multiples can be used. Application of the latter approach is illustrated with a field data example, and its performance is evaluated quantitatively through a measurement of S/N ratio improvement. We demonstrate that a DMO-based implementation of 3D SRME outperforms conventional 2D SRME and can accurately predict and attenuate complex 3D multiples.


Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. S65-S70 ◽  
Author(s):  
Lele Zhang ◽  
Evert Slob

Internal multiple reflections have been widely considered as coherent noise in measured seismic data, and many approaches have been developed for their attenuation. The Marchenko multiple elimination (MME) scheme eliminates internal multiple reflections without model information or adaptive subtraction. This scheme was originally derived from coupled Marchenko equations, but it was modified to make it model independent. It filters primary reflections with their two-way traveltimes and physical amplitudes from measured seismic data. The MME scheme is applied to a deepwater field data set from the Norwegian North Sea to evaluate its success in removing internal multiple reflections. The result indicates that most internal multiple reflections are successfully removed and primary reflections masked by overlapping internal multiple reflections are recovered.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. V115-V130 ◽  
Author(s):  
Jun Matsushima ◽  
Makoto Suzuki ◽  
Ippei Matsugi ◽  
Yoshibumi Kato ◽  
Shuichi Rokugawa

It is important to obtain reliable attenuation results from experimental data to elucidate the physical mechanism responsible for ultrasonic wave attenuation. For attenuation estimation, a time window is often used to compute the frequencies of the direct-arrival waveforms. However, the effect of windowing distorts the spectral distribution due to a spectral leakage effect, degrading the attenuation estimates. We propose a method that enables accurate measurement of ultrasonic attenuation using sweep signals under the assumptions that velocity dispersion can be ignored and the quality factor [Formula: see text] is not dependent on frequency. We obtained the spectral amplitude of the sweep signal in the frequency-time domain using the continuous wavelet transform and estimated attenuation in the time-scale spectrum domain using the spectral-ratio method. This method is independent of the effect of windowing, whereas the windowing effect underestimates the attenuation results. In the absence of noise, the estimated attenuation results using sweep signals are in perfect agreement with the given input values, and the accuracy of the estimated attenuation results from windowed pulse waveforms depends on the extraction window length. However, our numerical experiments demonstrated that the proposed method is largely influenced by the existence of overlapping sweep events such as multiple reflections between the source and receiver transducer. Thus, applicability of the proposed method is limited to highly attenuative media, in which overlapping events are much smaller than direct sweep signals because these multiple reflected events are largely attenuated. Application of the proposed method to laboratory experimental data yielded similar underestimation of the attenuation results due to the windowing effect in the case of highly attenuative media. We also evaluated the usefulness of observing compressed pulse waveforms with shorter duration from the crosscorrelation of sweep waveforms than the case of pulse generation.


2010 ◽  
Vol 7 (2) ◽  
pp. 149-157 ◽  
Author(s):  
Xiang-Chun Wang ◽  
Chang-Liang Xia ◽  
Xue-Wei Liu

2021 ◽  
Author(s):  
Rick Schrynemeeckers

Abstract Current offshore hydrocarbon detection methods employ vessels to collect cores along transects over structures defined by seismic imaging which are then analyzed by standard geochemical methods. Due to the cost of core collection, the sample density over these structures is often insufficient to map hydrocarbon accumulation boundaries. Traditional offshore geochemical methods cannot define reservoir sweet spots (i.e. areas of enhanced porosity, pressure, or net pay thickness) or measure light oil or gas condensate in the C7 – C15 carbon range. Thus, conventional geochemical methods are limited in their ability to help optimize offshore field development production. The capability to attach ultrasensitive geochemical modules to Ocean Bottom Seismic (OBS) nodes provides a new capability to the industry which allows these modules to be deployed in very dense grid patterns that provide extensive coverage both on structure and off structure. Thus, both high resolution seismic data and high-resolution hydrocarbon data can be captured simultaneously. Field trials were performed in offshore Ghana. The trial was not intended to duplicate normal field operations, but rather provide a pilot study to assess the viability of passive hydrocarbon modules to function properly in real world conditions in deep waters at elevated pressures. Water depth for the pilot survey ranged from 1500 – 1700 meters. Positive thermogenic signatures were detected in the Gabon samples. A baseline (i.e. non-thermogenic) signature was also detected. The results indicated the positive signatures were thermogenic and could easily be differentiated from baseline or non-thermogenic signatures. The ability to deploy geochemical modules with OBS nodes for reoccurring surveys in repetitive locations provides the ability to map the movement of hydrocarbons over time as well as discern depletion affects (i.e. time lapse geochemistry). The combined technologies will also be able to: Identify compartmentalization, maximize production and profitability by mapping reservoir sweet spots (i.e. areas of higher porosity, pressure, & hydrocarbon richness), rank prospects, reduce risk by identifying poor prospectivity areas, accurately map hydrocarbon charge in pre-salt sequences, augment seismic data in highly thrusted and faulted areas.


Sign in / Sign up

Export Citation Format

Share Document