Vector-acoustic reverse time migration of Volve ocean-bottom cable data set without up/down decomposed wavefields

Geophysics ◽  
2015 ◽  
Vol 80 (4) ◽  
pp. S137-S150 ◽  
Author(s):  
Matteo Ravasi ◽  
Ivan Vasconcelos ◽  
Andrew Curtis ◽  
Alexander Kritski
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 ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB175-WB182 ◽  
Author(s):  
Yan Huang ◽  
Bing Bai ◽  
Haiyong Quan ◽  
Tony Huang ◽  
Sheng Xu ◽  
...  

The availability of wide-azimuth data and the use of reverse time migration (RTM) have dramatically increased the capabilities of imaging complex subsalt geology. With these improvements, the current obstacle for creating accurate subsalt images now lies in the velocity model. One of the challenges is to generate common image gathers that take full advantage of the additional information provided by wide-azimuth data and the additional accuracy provided by RTM for velocity model updating. A solution is to generate 3D angle domain common image gathers from RTM, which are indexed by subsurface reflection angle and subsurface azimuth angle. We apply these 3D angle gathers to subsalt tomography with the result that there were improvements in velocity updating with a wide-azimuth data set in the Gulf of Mexico.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB27-WB39 ◽  
Author(s):  
Zheng-Zheng Zhou ◽  
Michael Howard ◽  
Cheryl Mifflin

Various reverse time migration (RTM) angle gather generation techniques have been developed to address poor subsalt data quality and multiarrival induced problems in gathers from Kirchhoff migration. But these techniques introduce new problems, such as inaccuracies in 2D subsurface angle gathers and edge diffraction artifacts in 3D subsurface angle gathers. The unique rich-azimuth data set acquired over the Shenzi field in the Gulf of Mexico enabled the generally artifact-free generation of 3D subsurface angle gathers. Using this data set, we carried out suprasalt tomography and salt model building steps and then produced 3D angle gathers to update the subsalt velocity. We used tilted transverse isotropy RTM with extended image condition to generate full 3D subsurface offset domain common image gathers, which were subsequently converted to 3D angle gathers. The angle gathers were substacked along the subsurface azimuth axis into azimuth sectors. Residual moveout analysis was carried out, and ray-based tomography was used to update velocities. The updated velocity model resulted in improved imaging of the subsalt section. We also applied residual moveout and selective stacking to 3D angle gathers from the final migration to produce an optimized stack image.


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.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. T167-T174 ◽  
Author(s):  
Dan Kosloff ◽  
Reynam C. Pestana ◽  
Hillel Tal-Ezer

A new scheme for the calculation of spatial derivatives has been developed. The technique is based on recursive derivative operators that are generated by an [Formula: see text] fit in the spectral domain. The use of recursive operators enables us to extend acoustic and elastic wave simulations to shorter wavelengths. The method is applied to the numerical solution of the 2D acoustic wave equation and to the solution of the equations of 2D dynamic elasticity in an isotropic medium. An example of reverse-time migration of a synthetic data set shows that the numerical dispersion can be significantly reduced with respect to schemes that are based on finite differences. The method is tested for the solutions of the equations of dynamic elasticity by comparing numerical and analytic solutions to Lamb’s problem.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. C295-C307 ◽  
Author(s):  
Pengfei Yu ◽  
Jianhua Geng ◽  
Chenlong Wang

Quasi-P (qP)-wavefield separation is a crucial step for elastic P-wave imaging in anisotropic media. It is, however, notoriously challenging to quickly and accurately obtain separated qP-wavefields. Based on the concepts of the trace of the stress tensor and the pressure fields defined in isotropic media, we have developed a new method to rapidly separate the qP-wave in a transversely isotropic medium with a vertical symmetry axis (VTI) by synthesized pressure from ocean-bottom seismic (OBS) data as a preprocessing step for elastic reverse time migration (ERTM). Another key aspect of OBS data elastic wave imaging is receiver-side 4C records back extrapolation. Recent studies have revealed that receiver-side tensorial extrapolation in isotropic media with ocean-bottom 4C records can sufficiently suppress nonphysical waves produced during receiver-side reverse time wavefield extrapolation. Similarly, the receiver-side 4C records tensorial extrapolation was extended to ERTM in VTI media in our studies. Combining a separated qP-wave by synthesizing pressure and receiver-side wavefield reverse time tensorial extrapolation with the crosscorrelation imaging condition, we have developed a robust, fast, flexible, and elastic imaging quality improved method in VTI media for OBS data.


Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. Q15-Q26 ◽  
Author(s):  
Giovanni Angelo Meles ◽  
Kees Wapenaar ◽  
Andrew Curtis

State-of-the-art methods to image the earth’s subsurface using active-source seismic reflection data involve reverse time migration. This and other standard seismic processing methods such as velocity analysis provide best results only when all waves in the data set are primaries (waves reflected only once). A variety of methods are therefore deployed as processing to predict and remove multiples (waves reflected several times); however, accurate removal of those predicted multiples from the recorded data using adaptive subtraction techniques proves challenging, even in cases in which they can be predicted with reasonable accuracy. We present a new, alternative strategy to construct a parallel data set consisting only of primaries, which is calculated directly from recorded data. This obviates the need for multiple prediction and removal methods. Primaries are constructed by using convolutional interferometry to combine the first-arriving events of upgoing and direct-wave downgoing Green’s functions to virtual receivers in the subsurface. The required upgoing wavefields to virtual receivers are constructed by Marchenko redatuming. Crucially, this is possible without detailed models of the earth’s subsurface reflectivity structure: Similar to the most migration techniques, the method only requires surface reflection data and estimates of direct (nonreflected) arrivals between the virtual subsurface sources and the acquisition surface. We evaluate the method on a stratified synclinal model. It is shown to be particularly robust against errors in the reference velocity model used and to improve the migrated images substantially.


2019 ◽  
Vol 16 (3) ◽  
pp. 327-337
Author(s):  
Ying-Ming Qu ◽  
Chong-Peng Huang ◽  
Chang Liu ◽  
Chang Zhou ◽  
Zhen-Chun Li ◽  
...  

Geophysics ◽  
2021 ◽  
pp. 1-67
Author(s):  
Yuzhu Liu ◽  
Weigang Liu ◽  
Zheng Wu ◽  
Jizhong Yang

Reverse time migration (RTM) has been widely used for imaging complex subsurface structures in oil and gas exploration. However, because only the adjoint of the forward Born modeling operator is applied to the seismic data in RTM, the output migration profile is biased in terms of the amplitude. To help partially balance the amplitude performance, the RTM image can be preconditioned with the inverse of the diagonal of the Hessian operator. Yet, existing preconditioning methods do not correctly consider the receiver-side effects, assuming that the receiver coverage is infinite or the velocity model is constant. We therefore provide a comparative study aiming to give a clearer understanding on the importance of incorporating the receiver-side effects by developing a frequency-domain scattering-integral reverse time migration (SI-RTM). In the proposed SI-RTM, the diagonal of the Hessian operator is explicitly computed in its exact formulation, and the source-side wavefield and receiver-side Green’s functions are obtained by solving the two-way wave equation. The computational cost is relatively affordable when compared with the more expensive least-squares RTM. In the comparative counterpart, the diagonal of the Hessian operator is approximated by the source-side illumination. We perform two synthetic numerical examples using an overthrust model and a complex reservoir model; the final migration images were significantly improved when the receiver-side effects were accurately considered. A third application of SI-RTM on one field data set acquired from the East China Sea further demonstrates the importance of incorporating the receiver-side effects in normalizing the RTM image. Findings of this study are expected to provide a theoretical basis for improving the ability of RTM imaging of subsurface structures, thereby critically advancing the application of geophysical techniques for imaging complex environments.


Sign in / Sign up

Export Citation Format

Share Document