Technical Library

Utilizing a novel quantitative interpretation workflow to derisk shallow hydrocarbon prospects — a Barents Sea case study
Author: Laurent Olivier Feuilleaubois, Valerie Charoing, Andrea Maioli, Cyrille Reiser
First Break  1 March 2017The authors demonstrate that combining separated wavefield imaging with highresolution velocity computations improves the prestack amplitude analysis of shallow water hydrocarbon targets in the Barents Sea.

FWI for model updates in largecontrast media
Author: Sverre BrandsbergDahl, Nizar Chemingui, Alejandro Valenciano, Jaime RamosMartinez, Lingyun Qiu
The Leading Edge  1 January 2017We describe a new solution for recovering the longwavelength features of a velocity model in gradientbased fullwaveform inversion (FWI). The method uses reflected and transmitted wave modes to recover highresolution velocity models. The new FWI gradient enables reliable velocity updates deeper than the maximum penetration depth of diving waves and reduces the FWIdependency on recording ultralong offsets. We also discuss a new FWI regularization scheme that overcomes the limitations of the inversion in the presence of highcontrast geobodies and cycle skipping. The solution utilizes a priori information about the earth model in the regularization as an extra term in the objective function. The implementation makes use of the split Bregman method, making it efficient and accurate. Results from applying the new FWI gradient to field data show that we can combine both transmitted and reflected energy in a global FWI scheme to obtain highresolution velocity models without imprint of the reflectivity on the velocity updates. We illustrate the new regularization method's potential on the BP 2004 velocity benchmark model where our regularized FWI solution is capable of using a simple starting velocity model to deliver a highquality reconstruction of salt bodies.

New technique corrects distortions
Author: Volker Dirks, Senira Kattah
Hart E&P  1 November 2016WEI leads to improved amplitude fidelity and image resolution in the Jequitinhonha Basin offshore Brazil.

A robust FWI gradient for highresolution velocity model building
Author: Jaime RamosMartinez, Sean Crawley, Zuihong Zou, Alejandro Valenciano, Elena Klochikhina, Nizar Chemingui
SEG  1 October 2016We describe a robust method to produce longwavelength updates in gradientbased Full Waveform Inversion (FWI). The gradient is computed by applying dynamic weights in the velocity sensitivity kernel derived from impedance and velocity parameterization of the classical objective function. The new kernel implementation effectively eliminates the migration isochrones produced by the specular reflections and emphasizes the lowwavenumber components in the gradient in heterogeneous media. The new gradient is able to provide velocity updates beyond penetration depth of diving waves. We use a synthetic example to illustrate how this dynamically weighted FWI gradient successfully recovers the background velocity from precritical reflections. We apply the new approach to 2D and 3D dual sensor data from deepwater Gulf of Mexico. Results show how the dynamically weighted FWI gradient can combine both transmitted and reflected energy in a global FWI scheme and provide highresolution velocity models without migration imprint in the updates.

Quantification of the Resolution in the Inversion results from Towed Streamer EM Data
Author: Joel Skogman, Johan Mattsson,
SEG  1 October 2016The resolution of resistivity models from 3D inversion of Towed Streamer EM data is analyzed using the resolution matrix. In particular, a new resolution length derived from the resolution matrix is introduced. This measure is obtained by dividing the volume of an inversion cell with the corresponding diagonal element in the resolution matrix. Hence, the measure is independent on the actual cell discretization of the inversion domain. Resolution lengths are calculated both for a hypothetical case with noise free data and for resistivity sections obtained with regularized 3D inversion. The resolution of the ideal case with no noise constitutes the upper limit of the best possible resolving power of the Towed Streamer EM data. This limit is then compared against the results where smoothing regularization has been applied in the 3D inversion. A comparison is illustrated for inversion results of Towed Steamer EM data from a survey line in the Barents Sea.

Sparsity Promoting Morphological Decomposition for Coherent Noise Suppression: Application to Streamer Vibration Related Noise
Author: Pierre Turquais, Endrias Asgedom, Walter Söllner
SEG  1 October 2016We address the signal and noise separation problem where the noise is coherent. We use a dictionary learning method to learn a dictionary of unit vectors called atoms; each one representing an elementary waveform redundant in the noisy data. In such a learned dictionary, some atoms represent signal waveforms while others represent noise waveforms. Using a multivariate Gaussian classifier trained on a noise recording, the atoms representing noise waveforms are discriminated and separated from the atoms representing seismic waveforms and two subdictionaries are created; one describing the morphology of the signal, the other describing the morphology of the noise. Using these subdictionaries, a morphological component analysis problem is set to separate the seismic signal and the coherent noise. In contrast to fixing transforms for representing the noise and the signal, our method is entirely adapting to the morphology of the signal and the noise. We present an application for removing streamer vibration related noise and show successful denoising results on synthetic and field data examples.

Adaptive multiple subtraction by statistical curvelet matching
Author: Truong Nguyen, Richard Dyer
SEG  1 October 2016In this paper, we present a novel algorithm that performs adaptive subtraction of a multiple model in the curvelet domain. The algorithm is based on the observation that the predicted model is an imperfect estimate of the actual multiples, containing two types of error: 1. a systematic error that manifests as an approximately constant phase and amplitude error within each subband and 2. a localized error that potentially varies from coefficient to coefficient. Adaptation of the model is automatically controlled by parameters provided by a statistical modelling step. Results show that the algorithm works well with different types of multiples and levels of noise.

Imaging the total wavefields by reflectivity inversion using amplitudenormalized wavefield decomposition: field data example
Author: Alba Ordoñez, University of Oslo, Walter Söllner, Tilman Klüver
SEG  1 October 2016Based on an example acquired with dualsensor towedstreamers and time and depth distributed sources, we image the total upand downgoing wavefields using primaries and multiples. The imaging framework is based on computing the subsurface impulse response (i.e., reflectivity). At every depth level, the latter can be obtained by inverting the matrix form of an integral equation defined in terms of the amplitudenormalized upgoing pressure and downgoing vertical velocity wavefields. This procedure gives the reflectivity matrix. The total upgoing wavefield used in the imaging scheme is composed of the scattered energy from primaries and multiples. The primary reflected wavefield is generated by a direct downgoing source wavefield, which is mostly passing the acquisition surface at offsets smaller than the nearest data channel. Hence, the most relevant part of the direct wavefield is not measured in the studied example. From the nearfield pressure measurements, we predict the missing direct arrivals needed to image the primary reflected wavefield; and we then synthesize the total downgoing wavefield by adding the downgoing scattered energy. By downward extrapolating the total up and downgoing wavefields, the information of the subsurface is extracted from the reflectivity matrix in the spatial and angular domains.

Postmigration image optimization: a Gulf of Mexico case study
Author: Sergey Frolov, Samuel Brown, Sriram Arasanipalai, Nizar Chemingui
SEG  1 October 2016There are many factors which affect the quality of a subsurface image: acquisition geometry, data quality, accuracy of the subsurface model, imaging algorithm, image postprocessing steps, etc. Despite improvements in seismic acquisition, such as full azimuth (FAZ) surveys, and advanced imaging algorithms such as reverse time migration (RTM), focusing energy beneath a complex salt overburden remains a challenging problem. We propose a postimaging workflow to mitigate migration artifacts, limited illumination problems, and velocity model inaccuracies that cause image distortion. RTM images derived with an inverse scattering imaging condition are decomposed into angle/azimuth domain gathers. The angles and azimuths which constructively interfere at each image point to produce the best image are chosen so as to remove both coherent and incoherent noise and produce an optimal stack. Finally, structureconformable filtering is applied. We apply this workflow to an RTM image from Garden Banks in the Gulf of Mexico .

3D inversion of the towed streamer EM data using the seismic constraints
Author: Michael S. Zhdanov, Masashi Endo, David Sunwall, Martin Čuma, JennyAnn Malmberg, Allan McKay, Tashi Tshering, Jonathan Midgley
SEG  1 October 2016In this paper we introduce a largescale 3D inversion technique for towed streamer electromagnetic (EM) data, which incorporates seismic constraints. The inversion algorithm is based on the integral equation (IE) forward modeling and utilizes a reweighted regularized conjugate gradient method with adaptive regularization to minimize the objective functional. We have also incorporated in the inversion the moving sensitivity domain approach in order to invert the entire largescale towed streamer EM survey data while keeping the accuracy and reducing the time and memory/storage of the computation. The developed algorithm and software can take into account the constraints based on seismic and welllog data, and provide the inversion guided by these constraints. Application of the developed method to the interpretation of the largescale towed streamer EM survey data acquired in the Barents Sea demonstrates its practical effectiveness.