Attenuation of free-surface multiples by up/down deconvolution for marine towed-streamer data

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. V129-V138 ◽  
Author(s):  
Mariusz Majdański ◽  
Clément Kostov ◽  
Ed Kragh ◽  
Ian Moore ◽  
Mark Thompson ◽  
...  

Free-surface-related multiples in marine seismic data are commonly attenuated using adaptive subtraction of the predicted multiple energy. An alternative method, based on deconvolution of the upgoing wavefield by the downgoing wavefield, was previously applied to ocean-bottom data. We apply the deconvolution method to towed-streamer data acquired in an over/under configuration. We also use direct arrival deconvolution that results in source wavelet designature only, as a benchmark to verify the full multiple deconvolution result. Detailed synthetic data analysis, including sensitivity tests, explains each data processing step and its effects on the final result. We then apply this verified preprocessing sequence to field data from the Kristin area of the North Sea, with a focus on the direct arrival prediction using the near-field hydrophone method. Prestack evaluation of the results shows that the method applied to the field data provides designature, source-side deghosting, and attenuation of multiples. We show comparable stacked results from our method and from 2D iterative surface-related multiple elimination. The workflow has the benefit that it does not require an adaptive subtraction step or iterative application. However, an accurate direct arrival prediction is essential for the successful application of the method. This prediction is obtained using near-field hydrophone measurements that can be recorded with some commercial acquisition systems.

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1303-1309 ◽  
Author(s):  
Ola Eiken ◽  
Geir Ultveit Haugen ◽  
Michel Schonewille ◽  
Adri Duijndam

Seismic reservoir monitoring has become an important tool in the management of many fields. Monitoring subtle changes in the seismic properties of a reservoir caused by production places strong demands on seismic repeatability. A lack of repeatability limits how frequently reservoir changes can be monitored or the applicability of seismic monitoring at all. In this paper we show that towing many streamers with narrow separation, combined with cross‐line interpolation of data onto predefined sail lines, can give highly repeatable marine seismic data. Results from two controlled zero time lag monitoring experiments in the North Sea demonstrate high sensitivity to changing water level and variations in lateral positions. After corrections by deterministic tidal time shifts and spatial interpolation of the irregularly sampled streamer data, relative rms difference amplitude levels are as low as 12% for a deep, structurally complex field and as low as 6% for a shallow, structurally simple field. Reducing the degree of nonrepeatability to as low as 6% to 12% allows monitoring of smaller reflectivity changes. In terms of reservoir management this has three important benefits: (1) reservoirs with small seismic changes resulting from production can be monitored, (2) reservoirs with large seismic changes can be monitored more frequently, and (3) monitoring data can be used more quantitatively.


Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 327-341 ◽  
Author(s):  
Lasse Amundsen

This paper presents a new, wave‐equation based method for eliminating the effect of the free surface from marine seismic data without destroying primary amplitudes and without any knowledge of the subsurface. Compared with previously published methods which require an estimate of the source wavelet, the present method has the following characteristics: it does not require any information about the marine source array and its signature, it does not rely on removal of the direct wave from the data, and it does not require any explicit deghosting. Moreover, the effect of the source signature is removed from the data in the multiple elimination process by deterministic signature deconvolution, replacing the original source signature radiated from the marine source array with any desired wavelet (within the data frequency‐band) radiated from a monopole point source. The fundamental constraint of the new method is that the vertical derivative of the pressure or the vertical component of the particle velocity is input to the free‐surface demultiple process along with pressure recordings. These additional data are routinely recorded in ocean‐bottom seismic surveys. The method can be applied to conventional towed streamer pressure data recorded in the water column at a depth which is greater than the depth of the source array only when the pressure derivative can be estimated, or even better, is measured. Since the direct wave and its source ghost is part of the free‐ surface demultiple, designature process, the direct arrival must be properly measured for the method to work successfully. In the case when the geology is close to horizontally layering, the free‐surface multiple elimination method greatly simplifies, reducing to a well‐known deterministic deconvolution process which can be applied to common shot gathers (or common receiver gathers or common midpoint gathers when source array variations are negligible) in the τ-p domain or frequency‐wavenumber domain.


2020 ◽  
Author(s):  
I. Gregory ◽  
Z. Dobo ◽  
F. Ebrahim ◽  
J. Sinden ◽  
P. McDonnell ◽  
...  

Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1293-1303 ◽  
Author(s):  
Luc T. Ikelle ◽  
Lasse Amundsen ◽  
Seung Yoo

The inverse scattering multiple attenuation (ISMA) algorithm for ocean‐bottom seismic (OBS) data can be formulated in the form of a series expansion for each of the four components of OBS data. Besides the actual data, which constitute the first term of the series, each of the other terms is computed as a multidimensional convolution of OBS data with streamer data, and aims at removing one specific order of multiples. If the streamer data do not contain free‐surface multiples, we found that the computation of only the second term of the series is needed to predict and remove all orders of multiples, whatever the water depth. As the computation of the various terms of the series is the most expensive part of ISMA, this result can produce significant savings in computation time, even in data storage, as we no longer need to store the various terms of the series. For example, if the streamer data contained free‐surface multiples, OBS seismic data of 6‐s duration, corresponding to a geological model of the subsurface with 250‐m water depth, require the computation of five terms of the series for each of the four components of OBS data. With the new implementation, in which the streamer data do not contain free‐surface multiples, we need the computation of only one term of the series for each component of the OBS data. The saving in CPU time for this particular case is at least fourfold. The estimation of the inverse source signature, which is an essential part of ISMA, also benefits from the reduction of the number of terms needed for the demultiple to two because it becomes a linear inverse problem instead of a nonlinear one. Assuming that the removal of multiple events produces a significant reduction in the energy of the data, the optimization of this problem leads to a stable, noniterative analytic solution. We have also adapted these results to the implementation of ISMA for vertical‐cable (VC) data. This implementation is similar to that for OBS data. The key difference is that the basic model in VC imaging assumes that data consist of receiver ghosts of primaries instead of the primaries themselves. We have used the following property to achieve this goal. The combination of VC data with surface seismic data, which do not contain free‐surface multiples, allows us to predict free‐surface multiples and receiver ghosts as well as the receiver ghosts of primary reflections. However, if the direct wave arrivals are removed from the VC data, this combination will not predict the receiver ghosts of primary reflections. The difference between these two predictions produces data containing only receiver ghosts of primaries.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. V213-V225 ◽  
Author(s):  
Shaohuan Zu ◽  
Hui Zhou ◽  
Yangkang Chen ◽  
Shan Qu ◽  
Xiaofeng Zou ◽  
...  

We have designed a periodically varying code that can avoid the problem of the local coherency and make the interference distribute uniformly in a given range; hence, it was better at suppressing incoherent interference (blending noise) and preserving coherent useful signals compared with a random dithering code. We have also devised a new form of the iterative method to remove interference generated from the simultaneous source acquisition. In each iteration, we have estimated the interference using the blending operator following the proposed formula and then subtracted the interference from the pseudodeblended data. To further eliminate the incoherent interference and constrain the inversion, the data were then transformed to an auxiliary sparse domain for applying a thresholding operator. During the iterations, the threshold was decreased from the largest value to zero following an exponential function. The exponentially decreasing threshold aimed to gradually pass the deblended data to a more acceptable model subspace. Two numerically blended synthetic data sets and one numerically blended practical field data set from an ocean bottom cable were used to demonstrate the usefulness of our proposed method and the better performance of the periodically varying code over the traditional random dithering code.


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. V457-V471
Author(s):  
Thomas Andre Larsen Greiner ◽  
Volodya Hlebnikov ◽  
Jan Erik Lie ◽  
Odd Kolbjørnsen ◽  
Andreas Kjelsrud Evensen ◽  
...  

Seismic exploration in complex geologic settings and shallow geologic targets has led to a demand for higher spatial and temporal resolution in the final migrated image. Conventional marine seismic and wide-azimuth data acquisition lack near-offset coverage, which limits imaging in these settings. A new marine source-over-cable survey, with split-spread configuration, known as TopSeis, was introduced in 2017 to address the shallow-target problem. However, wavefield reconstruction in the near offsets is challenging in the shallow part of the seismic record due to the high temporal frequencies and coarse sampling that leads to severe spatial aliasing. We have investigated deep learning as a tool for the reconstruction problem, beyond spatial aliasing. Our method is based on a convolutional neural network (CNN) approach trained in the wavelet domain that is used to reconstruct the wavefield across the streamers. We determine the performance of the proposed method on broadband synthetic data and TopSeis field data from the Barents Sea. From our synthetic example, we find that the CNN can be learned in the inline direction and applied in the crossline direction, and that the approach preserves the characteristics of the geologic model in the migrated section. In addition, we compare our method to an industry-standard Fourier-based interpolation method, in which the CNN approach shows an improvement in the root-mean-square (rms) error close to a factor of two. In our field data example, we find that the approach reconstructs the wavefield across the streamers in the shot domain, and it displays promising characteristics of a reconstructed 3D wavefield.


Geophysics ◽  
1988 ◽  
Vol 53 (10) ◽  
pp. 1347-1354 ◽  
Author(s):  
Song Cao ◽  
Ian Lerche ◽  
Christian Hermanrud

Using nonlinear inverse techniques, we show that the change in borehole temperature with time, after mud circulation has stopped, can be used to provide very precise estimates of true formation temperatures and of mud temperature at the time circulation stopped. In addition, ruggedly stable estimates can also be made of the thermal invasion distance, the formation thermal conductivity, and the efficiency of heating the mud by the thermal recovery wave. It is well known that convective heat flow into (or out of) the formation influences the thermal recovery. We show that the flow rate at the borehole can also be obtained approximately from the borehole temperature measurements using the inverse method. Fidelity and reproducibility of the inverse procedure arc examined using synthetic data. Applications to field data from three wells in the Cooper basin of Australia and four wells in the North Sea confirm the accuracy of the method in satisfying the observed data, in determining true formation temperatures, and in assessing the magnitude of fluid penetration into the formation.


Sign in / Sign up

Export Citation Format

Share Document