scholarly journals True-amplitude migration through regularized extended linearized waveform inversion

Geophysics ◽  
2021 ◽  
pp. 1-138
Author(s):  
Ettore Biondi ◽  
Mark A. Meadows ◽  
Biondo Biondi

The ability to create subsurface images whose amplitudes are proportional to the elastic wavefield variations recorded within seismic data as a function of reflection angle is fundamental for performing accurate amplitude-versus-offset (AVO) analysis and inversion. A process that generates such images is commonly referred to as true-amplitude migration. We demonstrate how the extended subsurface-offset image space is able to preserve the elastic behavior of the primary reflections when these events are acoustically migrated with a reverse-time-migration (RTM) approach performed in a least-squares fashion. Using a single-interface model, we show how the angle-domain image amplitude variations from an extended-offset acoustically migrated image closely follow the theoretical elastic Zoeppritz response even at the critical angle. Furthermore, we present a subsalt synthetic test in which single-component ocean-bottom-node (OBN) data are employed within a regularized linearized waveform inversion procedure. In this test, we highlight the ability of the acoustic extended-angle image domain to preserve the correct elastic amplitude variations of the reflected events from three subsalt sand lenses. The proposed method allows the accurate inversion of elastic-wave data for subsurface parameter variations that are critical for reservoir characterization in oil and gas exploration and production. We demonstrate its performance on an ocean-bottom-node (OBN) field dataset recorded in the Gulf of Mexico in which the AVO response of a potential gas-bearing prospect is correctly retrieved.

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. S261-S270 ◽  
Author(s):  
Daniel Rocha ◽  
Paul Sava ◽  
Antoine Guitton

We have developed a least-squares reverse time migration (LSRTM) method that uses an energy-based imaging condition to obtain faster convergence rates when compared with similar methods based on conventional imaging conditions. To achieve our goal, we also define a linearized modeling operator that is the proper adjoint of the energy migration operator. Our modeling and migration operators use spatial and temporal derivatives that attenuate imaging artifacts and deliver a better representation of the reflectivity and scattered wavefields. We applied the method to two Gulf of Mexico field data sets: a 2D towed-streamer benchmark data set and a 3D ocean-bottom node data set. We found LSRTM resolution improvement relative to RTM images, as well as the superior convergence rate obtained by the linearized modeling and migration operators based on the energy norm, coupled with inversion preconditioning using image-domain nonstationary matching filters.


Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. R63-R77 ◽  
Author(s):  
Denes Vigh ◽  
Kun Jiao ◽  
Dave Watts ◽  
Dong Sun

Recent computational improvements allowed us to simulate elastic wavefields in a 3D manner and undertake the challenge of executing elastic full-waveform inversion (EFWI). The 3D SEG/EAGE overthrust synthetic data were used to run the initial tests, which included using all three components for the simulation. The inversion targeted two regions: the channel system and the overthrusted zone, which proved the effectiveness of EFWI to delineate geology in terms of [Formula: see text] and [Formula: see text] velocity fields. For the field data experiment to demonstrate the technologies, we elected to use a Gulf of Mexico ocean bottom cable data set, which allowed us to take advantage of relatively large offsets along with the 4C acquisition. The input data were minimally processed mostly through noise removal, and the initial model was a Gaussian smoothed version of grid tomography output, which is done by a prestack migrated gather flattening process. During EWFI, a multiscale approach was followed to ensure convergence, and the early stages of the [Formula: see text]/[Formula: see text] ratio were constrained by the mud rock-line ratio. When the last sets of inversions were executed, this constraint was eliminated to ensure the simultaneous update of the [Formula: see text] and [Formula: see text] velocity fields. The density was kept constant to keep the inversion at a simple level, which allowed us to draw essential conclusions. The velocity fields were validated through an imaging algorithm of the elastic reverse time migration, and the imaging shows clear structural improvements when inputting the inverted velocities in conjunction with the measurements. If full-waveform inversion can provide multiple earth parameters, the user can use the process to detect gas zones along with sand and shale content of the subsurface, which will further assist the drilling decisions. We achieved this by simulating the earth more accurately with the elastic wave propagation in the algorithms.


2021 ◽  
Vol 40 (5) ◽  
pp. 348-356
Author(s):  
Cheryl Mifflin ◽  
Drew Eddy ◽  
Brad Wray ◽  
Lin Zheng ◽  
Nicolas Chazalnoel ◽  
...  

The story of seismic imaging over BHP's Shenzi Gulf of Mexico production field follows the history of offshore seismic imaging, from 2D to 3D narrow-azimuth streamer acquisition and to its leading the wide-azimuth movement with the Shenzi rich-azimuth (RAZ) survey. Each RAZ reprocessing project over the last 15 years applied the latest processing technology, culminating in hundreds of scenario tests to refine the salt model, but eventually the RAZ data reached a technical limit. A new ocean-bottom-node (OBN) survey acquired in 2020 has produced a step-change improvement over the legacy RAZ image. The uplift can be attributed to several factors. First, an OBN feasibility and survey design study demonstrated that a core of dense nodes combined with sparse nodes would improve the accuracy and resolution of the full-waveform inversion (FWI) solution. Second, the OBN data acquired following the survey design and employing FWI as the main model-building tool realized the predicted improvement. The result was a substantial change to the complex salt model, verified by a salt proximity survey as well as other salt markers, and improvement in imaging over the entire field. In addition to the improvement arising from a more accurate FWI velocity model, the steep-dip imaging also benefited from the new full-azimuth and long-offset data. However, the best steep-dip and fault imaging comes from the FWI image, a direct estimation of reflectivity from the FWI velocity. As the maximum frequency used by FWI moves toward the maximum frequency of the final reverse time migration (RTM), the FWI image approaches the resolution necessary to compete as the primary interpretation volume. Its subsalt illumination surpassed that of the RTM and even the least-squares RTM volumes. These imaging improvements are providing a new understanding of the faults and stratigraphic relationships of the field.


2021 ◽  
Vol 9 ◽  
Author(s):  
Yanbao Zhang ◽  
Yike Liu ◽  
Jia Yi ◽  
Xuejian Liu

Nowadays the ocean bottom node (OBN) acquisition is widely used in oil and gas resource exploration and seismic monitoring. Conventional imaging algorithms of OBN data mainly focus on the processing of up-going primaries and down-going first-order multiples. Up-going multiples and higher-order down-going multiples are generally regarded as noise and should be eliminated or ignored in conventional migration methods. However, multiples carry abundant information about subsurface structures where primaries cannot achieve. To take full advantage of multiples, we propose a migration method using OBN down-going all-order multiples. And then the least-squares optimization algorithm is used to suppress crosstalks. Finally, a phase-encoding-based migration algorithm is developed to cut down the computational cost by blending several common receiver gathers together using random time delays and polarity reversals. Numerical experiments on the complex Marmousi model illustrate that the developed approach can enlarge the imaging area evidently, reduce the computational cost effectively, and enhance the image quality by suppressing crosstalks and improving the resolution.


2017 ◽  
Author(s):  
Bruno Pereira-Dias ◽  
André Bulcão ◽  
Djalma Soares Filho ◽  
Roberto Dias ◽  
Felipe Duarte ◽  
...  

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S333-S345 ◽  
Author(s):  
Pengfei Yu ◽  
Jianhua Geng ◽  
Xiaobo Li ◽  
Chenlong Wang

Conventionally, multicomponent geophones used to record the elastic wavefields in the solid seabed are necessary for ocean bottom seismic (OBS) data elastic reverse time migration (RTM). Particle velocity components are usually injected directly as boundary conditions in the elastic-wave equation in the receiver-side wavefield extrapolation step, which causes artifacts in the resulting elastic images. We have deduced a first-order acoustic-elastic coupled equation (AECE) by substituting pressure fields into the elastic velocity-stress equation (EVSE). AECE has three advantages for OBS data over EVSE when performing elastic RTM. First, the new equation unifies wave propagation in acoustic and elastic media. Second, the new equation separates P-waves directly during wavefield propagation. Third, three approaches are identified when using the receiver-side multicomponent particle velocity records and pressure records in elastic RTM processing: (1) particle velocity components are set as boundary conditions in receiver-side vectorial extrapolation with the AECE, which is equal to the elastic RTM using the conventional EVSE; (2) the pressure component may also be used for receiver-side scalar extrapolation with the AECE, and with which we can accomplish PP and PS images using only the pressure records and suppress most of the artifacts in the PP image with vectorial extrapolation; and (3) ocean-bottom 4C data can be simultaneously used for elastic images with receiver-side tensorial extrapolation using the AECE. Thus, the AECE may be used for conventional elastic RTM, but it also offers the flexibility to obtain PP and PS images using only pressure records.


2016 ◽  
Vol 4 (2) ◽  
pp. 28 ◽  
Author(s):  
Sunmonu Ayobami ◽  
Adabanija Adedapo ◽  
Adagunodo Aanuoluwa ◽  
Adeniji Ayokunnu

Hydrocarbon resources have become the most essential commodity contributing to any nation’s growth and development in the recent years. For the past decades now, the quest for hydrocarbon resources has been increasing in an arithmetic rate that its supply can no longer meets the demand for its consumption today. In petroleum industry, seismic and well log analyses play a vital role in oil and gas exploration and formation evaluation. This study is aimed to effectively characterize the reservoirs and analyze the by-passed pay in Philus Field, Niger-Delta, Nigeria in order to look into the economic viability and profitability of the volume of oil in the identified reservoir(s). The faults in the study area trend in NW-SE direction and dip towards the south. Seven reservoirs were mapped on Philus field. A discovery trap and a by-passed (new prospect) trap were mapped out on the field. The petrophysical analysis showed that porosity of Philus field was 0.24. The volumetric analysis showed that the Stock Tank Original Oil in Place of discovery trap (Philus field) ranged from 1.6 to 43.1 Mbbl while that of new prospect trap ranged from 18.1 to 211.3 Mbbl. It is recommended that the oil reserve of Philus field needs to be recalculated.


Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. R45-R55 ◽  
Author(s):  
Espen Birger Raknes ◽  
Wiktor Weibull

In reverse time migration (RTM) or full-waveform inversion (FWI), forward and reverse time propagating wavefields are crosscorrelated in time to form either the image condition in RTM or the misfit gradient in FWI. The crosscorrelation condition requires both fields to be available at the same time instants. For large-scale 3D problems, it is not possible, in practice, to store snapshots of the wavefields during forward modeling due to extreme storage requirements. We have developed an approximate wavefield reconstruction method that uses particle velocity field recordings on the boundaries to reconstruct the forward wavefields during the computation of the reverse time wavefields. The method is computationally effective and requires less storage than similar methods. We have compared the reconstruction method to a boundary reconstruction method that uses particle velocity and stress fields at the boundaries and to the optimal checkpointing method. We have tested the methods on a 2D vertical transversely isotropic model and a large-scale 3D elastic FWI problem. Our results revealed that there are small differences in the results for the three methods.


Sign in / Sign up

Export Citation Format

Share Document