Residual Moveout Analysis in a 3D dataset from Amplitude Variations with Offfset

Similar documents
Enhanced Imaging of Subsurface by Pre-Stack Merging of Multi-Vintage 3D Data Sets and its Prestack Time Migration A Case Study

Seamless Pre-Stack Integration of Streamer Datasets with Q-Marine Dataset A Case Study from Eastern Coast of India

Pre-stack Merging of 3D Vintages : A case History from Assam Arakan Fold Belt

High resolution residual velocity analysis and application

Selection of an optimised multiple attenuation scheme for a west coast of India data set

Common-angle processing using reflection angle computed by kinematic pre-stack time demigration

Enhanced imaging through 3D volumetric refraction static in Rajasthan area: a case study.

Comparison of Kirchhoff and Wave Equation PSDM : A Case Study from West Patan Area

Seismic Time Processing. The Basis for Modern Seismic Exploration

Seismic Reflection Method

Flexi Binning : A Solution for Irregular Acquisition Geometry

11th Biennial International Conference & Exposition

An implementable alternate survey geometry design for ideal Land 3D Seismic Acquisition

Th SRS3 07 A Global-scale AVO-based Pre-stack QC Workflow - An Ultra-dense Dataset in Tunisia

Using Similarity Attribute as a Quality Control Tool in 5D Interpolation

The Lesueur, SW Hub: Improving seismic response and attributes. Final Report

PSPRO AZIMUTHAL. Big Data Software for Interactive Analysis and Interpretation of Full-Azimuth Seismic Gathers. Power to Decide

Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction

Improved Imaging through Pre-stack Trace Interpolation for missing offsets of OBC data A case study from North Tapti area of West Coast, India

Seismic Modeling, Migration and Velocity Inversion

Velocity versus Offset (VVO) Estimation Using Local Event Correlation and Its Application in Seismic Processing and Analysis

5D leakage: measuring what 5D interpolation misses Peter Cary*and Mike Perz, Arcis Seismic Solutions

Practical implementation of SRME for land multiple attenuation

Tu N Maximising the Value of Seismic Data for Continued Mature Field Development in the East Shetland Basin

Machine learning assisted velocity auto-picking K.J. Smith, Shearwater GeoServices, OpenCPS

Summary. Introduction

A Short Narrative on the Scope of Work Involved in Data Conditioning and Seismic Reservoir Characterization

C-Wave processing of 3D seismic data of Cambay Basin-A case study

Angle Gathers for Gaussian Beam Depth Migration

Y015 Complementary Data-driven Methods for Interbed Demultiple of Land Data

Value Addition from Reprocessing of Seismic Data: An Analysis of Some of the Case Histories

Offset Scaling. Irfan Saputra. December 2008/CE8R3 SAMPLE IMAGE

3-D vertical cable processing using EOM

3D Seismic Data Merging A Case History in Indian Context S. Basu, S.N. Dalei and D.P. Sinha SPIC ONGC Mumbai,

Obstacles in the analysis of azimuth information from prestack seismic data Anat Canning* and Alex Malkin, Paradigm.

Th ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data

Foolproof AvO. Abstract

A stochastic approach of Residual Move Out Analysis in seismic data processing

VISTA. Desktop seismic data processing software

P071 Land Data Regularization and Interpolation Using Azimuth Moveout (AMO)

Amplitude-preserved low velocity noise suppression

We G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation

Ray-based tomography with limited picking

SeisSpace Software. SeisSpace enables the processor to be able focus on the science instead of being a glorified data manager.

Anisotropy-preserving 5D interpolation by hybrid Fourier transform

Azimuth Moveout (AMO) for data regularization and interpolation. Application to shallow resource plays in Western Canada

96 Alkhalifah & Biondi

Tu N Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction

Multiple suppression with land data: A case study

Internal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction

Anatomy of common scatterpoint (CSP) gathers formed during equivalent offset prestack migration (EOM)

Coherent partial stacking by offset continuation of 2-D prestack data

Interactive Velocity Analysis (aaspi_iva): Program: Precompute Velocity Semblance, Interactive Velocity Analysis, Velocity Interpolation

Analysis of different interpolation methods for uphole data using Surfer software

N. Kazemi* (University of Tehran) & H.R. SiahKoohi (University of Tehran)

Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco

We Solutions for Scattered Surface-wave Attenuation in the Western Desert of Egypt

G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry

Summary. Figure 1: Simplified CRS-based imaging workflow. This paper deals with the boxes highlighted in green.

It is well known among the seismic data processing

C014 Shot Based Pre-Processing Solutions for a WATS Survey An Example from a Field Trial in Green Canyon Gulf of Mexico

High definition tomography brings velocities to light Summary Introduction Figure 1:

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Equivalent offset migration: the implementation and application update

G034 Improvements in 4D Seismic Processing - Foinaven 4 Years On

Common Reflection Angle Migration (CRAM) for improved input to reservoir description an example from Mumbai High Field

GeOlympus. GeoPACS GeoTR GeoStaR GeoWZ

Pre-stack deghosting for variable-depth streamer data. R. Soubaras* (CGGVeritas) Summary

Challenges of pre-salt imaging in Brazil s Santos Basin: A case study on a variable-depth streamer data set Summary

Th Towards Improved Time-lapse Seismic Repetition Accuracy by Use of Multimeasurement Streamer Reconstruction

GG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material!

GeoDepth. Velocity Determination, Modeling and Imaging

Seismic Data Enhancement While Interpretation A Challenge For Better Imagination of Structure and Reservoir Characterization With Confidence

Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS

Strategies to preserve the data fidelity at far offset Bo Zhang*, Tao Zhao, Jie Qi and Kurt J. Marfurt, The University of Oklahoma

Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic

Pre-Stack Seismic Data Analysis with 3D Visualization A Case Study*

Geogiga Seismic Pro 8.3 Release Notes

GEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011

PS wave AVO aspects on processing, inversion, and interpretation

Methodology of calculation and application of scalars to balance amplitude variations across merge volume

Geogiga Seismic Pro 8.0 Release Notes

Introduction. Surface and Interbed Multtple Elimination

Chapter 1. 2-D field tests INTRODUCTION AND SUMMARY

Zero offset VSP processing for coal reflections

Short Note. DMO velocity analysis with Jacubowicz s dip-decomposition method. David Kessler and Wai-Kin Chan*

PLEASE NOTE THAT THERE ARE QUESTIONS IN SOME OF THESE SECTIONS THAT ARE TO BE TURNED IN AS HOMEWORK.

The power of stacking, Fresnel zones, and prestack migration

Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers

3D predictive deconvolution for wide-azimuth gathers

Seismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy

D ATA P R O C E S S I N G S E R V I C E S. S e i s U P S O F T W A R E & D E V E L O P M E N T

Azimuthal anisotropy? The time and depth imaging points of view: an imaging case history.

Chapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design

A case study for salt model building using CFP full azimuth data W. Gao*, Z. Guo, M. Guo, Q. Zhang, S. Hightower, G. Cloudy Jr. and Z.

3-D VSP: Survey design and processing

Wave-equation migration from topography: Imaging Husky

An Integrated Pre-Stack Depth Migration Workflow Using Model- Based Velocity Estimation And Refinement

Seismic Full Waveform Inversion and Modeling

Static corrections in challenging cases

Transcription:

P-270 Residual Moveout Analysis in a 3D dataset from Amplitude Variations with Offfset Summary Subhendu Dutta, ONGC The amplitude variations with offset (AVO) is a very useful tool which compliments the conventional stacked section with added information concerning the converted waves (e.g. Ostrander 1984, Shuey 1985, Hilterman 1990). In 2001, H.W.Swan has shown that the sensitivity of amplitude variations with offset (AVO) to Normal moveout (NMO) velocity errors can be exploited to obtain more accurate velocities from a fully automated procedure, which significantly more sensitive to velocity variations than the semblance coefficient (as defined by Neideel and Taner in 1971). In this paper, we have discussed a case study to find out the residual moveout corrections on time migrated gathers from amplitude variations with offset in a large volume of 3D land dataset in Cauvery Basin using Swan s technique. This procedure uses AVO techniques in order to automatically perform detailed residual moveout analysis. Residual moveout operates as conventional normal moveout (2nd order) analysis as well as 4th order analysis to correct time migrated gathers at large offsets which exhibit 4th order moveout. This procedure enables to generate residual moveout or residual velocity volumes and thus can update RMS velocity volume & moveout corrected gathers (flat gathers in time migrated domain). In this paper, we have also discussed about the benefits and accuracy of this procedure over the conventional technique of picking residual moveout corrections through residual velocity analysis semblance plot using time migrated gathers as input & RMS velocity volume (migration velocity) as the background velocity. Introduction The pre-requisite for 3D Pre stack time migration (Kirchoff PSTM) for land data requires a conditioned decon cdp gather with applied residual statics and a near perfect RMS velocity model. The initial RMS velocity volume is usually obtained from DIX conversion of the of the stacking velocity in time domain for running of target oriented time migration (TOTM) in fixed line interval (inline or crossline). The product of TOTM is time migrated gathers (CRP gathers) and a PSTM stack. On these time migrated gathers, the actual RMS velocity is generally picked inversing the NMO on the time migrated gathers from the semblance plot and thus a near perfect RMS velocity volume is obtained after proper QC of the flatness of the TM gathers to run the final PSTM for full volume. After PSTM, again the flatness of the time migrated gathers are to be checked in fixed inline or cross line interval and it is seen that there is always a scope for residual velocity or residual moveout correction for fine tuning of the RMS velocity model as such. The procedure describe here provides an automatic calculation of residual moveout with respect to time migrated gathers after fixing parameters in the interactive mode on selected gathers and proper QC for flatness of gathers in different reflector level and thus a residual moveout volume can be prepared for the whole 3D volume. With the help of this residual moveout volume the earlier time migrated gathers can be corrected (PSTM flat gathers) and a final PSTM stack volume may be obtained. Study Area The area is situated to the South east of the Kovilkalappal field. It is bounded in the east and south by a NELP block operated by Niko Resources and in the North by Tiruttaraipundi high. Area covers a small oil field viz. Tulsapattinam. This campaign was acquired with the prime objective of delineation of Tulsapattinam pay (post rift Nannilam formation). A total of 9 exploratory wells have been drilled in and around the area. However one well A-1 falling in area A is the lone hydrocarbon bearing well in the area. The location map of the area is depicted in Fig-1. GPS, ONGC, Chennai, India * Email: sduttaind@yahoo.com

moveout (2nd order) analysis as well as 4 th order analysis to correct time migrated gathers at large offsets which exhibit 4 th order moveout. This procedure enables to generate residual moveout or velocity volumes and thus can update RMS velocity volume & moveout corrected gathers (flat gathers in time migrated domain). Auto-picking is computationally more efficient than it would be with semblance because only one correlation trace is calculated for each iteration, instead of one for every trial velocity. This efficiency makes velocity picking of every gather in a 3D survey feasible. Figure -1: Location Map of the area Reasons for 3D survey The proposed area is covered by earlier by 3D seismic survey in the year 1998-99 with (20m X 40m) Bin size, with the prime objective of delineation of Tulsapattinam pay (shallower play). But the quality of 3D seismic data is not sufficient to resolve deeper play or map the envisaged Stratigraphic features. Additionally, the frequency content in the available data is also low, as a result various seismic attribute attempted with Tulsapttinam 3D did not produce realistic output. As a result the high potential area remained under explored all these years. Thus this part of area is proposed for repeat acquisition with (20mX20m) Bin size. Theoretical Background The amplitude variations with offset (AVO) is a very useful tool which compliments the conventional stacked section with added information concerning the converted waves (e.g. Ostrander 1984, Shuey 1985, Hilterman 1990). In 2001, H.W.Swan has shown that the sensitivity of amplitude variations with offset (AVO) to Normal moveout (NMO) velocity errors can be exploited to obtain more accurate velocities from a fully automated procedure, which significantly more sensitive to velocity variations than the semblance coefficient (as defined by Neideel and Taner in 1971). This procedure uses AVO techniques in order to automatically perform detailed residual moveout analysis. Residual moveout operates as conventional normal This method consists of the following steps, AVO analysis, removal of NMO stretch errors, calculation of joint statistics between the analytic intercepts & gradients and picking an NMO velocity to null this indicator. Methodology & Procedure Input Data The input data for this procedure is time migrated gathers. These gathers have been obtained from PSTM (Kirchhoff). The conditioned decon gather & a RMS velocity model were the pre-requisite for first run of PSTM. There after RMS velocity was picked on time migrated gathers inversing NMO and a near perfect RMS velocity model was prepared. After that, PSTM for the full 3D volume was performed with this conditioned decon gather & RMS velocity. Fig-2 below illustrates the relationship for data with no AVO, but with a small error in stacking velocity. As per Swan s theory, a small stacking velocity error aligns data curve t (θ) and produces time misalignments ( t) and amplitude variations ( S). The time misalignment t is a function of the shot-receiver offset x, true rms velocity V, zero offset traveltime t₀, and velocity error V. 2

pstm gather in TM (time migrated) domain before and after pre conditioning, a calculated residual moveout section in one inline from the 3D volume and the residual section after smoothening. In this case only 2nd order nmo correction is applied. In this case the input gather is corrected with 2nd order correction only. Fig-5 shows the gathers after 4th order correction, which shows the flatness of the gathers in the far offset also. Figure-2 Processing Algorithm This procedure uses AVO techniques in order to automatically perform detailed residual moveout analysis. Automatic residual moveout operates as conventional normal moveout (2nd order) and as 4th order correction to correct gathers at large offsets which exhibit 4th order moveout. The basic workflow is mentioned in Fig-3. Figure- 4 Figure- 3 The input time migrated gathers requires some pre conditioning before AVO analysis and calculation for residual moveout in the gather if any. Fig-4 shows the input Figure- 5 Fig-6 shows the input gathers and the output flat pstm gathers after this process. 3

Figure-8 Figure- 6 Fig-7 shows the second & fourth order residual move out curve. The scale of residual velocity or residual moveout is described in the diagram. The point towards right shows the 100% positive residual moveout whereas the leftside is negative moveout. The nature of the curve, describe the accuracy of actual RMS velocity pick from the time migrated gathers. Whereas in Fig-9, the iterations are 30 & 40 respectively. Sometimes more iterations are required for flatness of the gathers and more iterations require more computational time. Therefore the optimum iterations are to be chosen keeping in view the computational time & quality of the output gather & stack. The same exercise has been done manually to find out the residual moveout correction from residual velocity analysis semblance plot. Figure- 7 The figure-8 shows the flatness of the input pstm gathers after 10 iterations & 20 iterations respectively, thus customizing the parameters. Figure- 9 The fig-10 shows the residual velocity analysis or residual moveout analysis in the manual hand picked mode. The left side in the velocity panel is the time migrated section, in the middle residual time migrated move-out and in the right hand side of the semblance shows the pstm flat gathers after application of residual moveout analysis and application. 4

Figure 12 Figure- 10 The fig-11 shows the residual moveout analysis of the same time migrated gather from amplitude variation with technique. Figure 13 Fig-13 shows the residual move out analysis from amplitude variations with offset of the same gather after zooming it. Fig-12 & 13 shows the manual residual moveout correction and residual moveout correction from AVO studies with this technique are almost yields similar results. Figure -11 Fig-12 shows the residual move out analysis from hand picked of the same gather after zooming it. 5

Figure 14 Fig-14 above shows a panel where the two time migrated gathers are before and after pre conditioning with this technique and the stack panel shows stacks before application of residual move out (A), after application of residual move out which was hand picked (B), and the third one (C) is after application of residual move out which was calculated from amplitude variations with offset. The results of hand picked residuals and residual picked or computed from this technique is yielding same results. Figure 16 Fig-16 shows the comparisons from another inline section from the same 3D dataset. Figure 17 Fig-17 shows a full inline from the 3D dataset after computations and application of residual moveout correction from AVO studies. Conclusions Figure 15 The above fig-15 shows the stack comparison before application of RMO, after application of RMO (hand picked) and application of RMO from avo technique. The actual RMS velocity is generally picked manually from the semblance plot, inversing the NMO on the time migrated gathers derived from the first pass pre-stack migration. The accuracy of RMS velocity pick depends upon the flatness of the time migrated gathers. When picking the RMS velocities for a large volume of 3D 6

dataset with a very close grid velocity analysis say, (400mX400m), there is always a chance of left-out residual move-out correction in the dataset. It is legitimate to check the residual moveout correction in the time migrated gather before final stack of the data. In this paper, the residual move-out corrections have done both manually and using Swan s technique. This procedure uses AVO techniques in order to automatically perform detailed residual move-out analysis. At first, in the interactive mode, the parameters has to be designed, then the job can be performed through batch mode automatically, which reduces the manual labour. However, this process works better when signal to noise ratio in dataset is good. Auto-picking is computationally more efficient than it would be with semblance because only one correlation trace is calculated for each iteration, instead of one for every trial velocity. Moreover in turn it reduces the computational time. This efficiency makes velocity picking of every gather in a 3D survey feasible. This procedure enables to generate residual move-out or velocity volumes and thus can update RMS velocity volume & move-out corrected gathers (flat gathers in time migrated domain). The processed section reveals the sub-surface image better than earlier vintages and brings out an improvement in the seismic section compared to the earlier processed section. Automatic Residual Moveout Correction using Swan s Algorithm, J Castagna, M Backus, 1993. Offset dependent reflectivity theory and practice of AVO analysis, 345 p. H W Swan, 2001, Velocities from amplitude variation with offset, Geophysics v 66 p 1735-1743. M T Taner, F Koehler and R E Sheriff, 1979, Complex seismic trace analysis, Geophysics v 44 1041-1063. Dutta etal, 3D Pre Stack Merging A Case History in Cauvery Basin, P-100, SPG 2010, Hyderabad. Acknowledgements The authors sincerely thank Shri BSN Murthy, GM-HGS and Dr R C Iyer, GM(GP)-In-charge RCC for their help and constant encouragement to carry out this work and permitting to publish the paper. The authors also would like to thank Shri Kailash Prasad, DGM(GP) & Shri K.V.Krishnan, DGM(GP), RCC,Chennai for sharing his technical views on the subject. The authors are also grateful to the user group for their constant support through effective interactions in various stages of processing of the data. Views expressed in this paper are that of the author(s) only and may not necessarily be of ONGC. References Automatic Residual Moveout Correction using Swan s Algorithm, J Castagna, M Backus, 1993. Offset dependent reflectivity theory and practice of AVO analysis, 345 p. Dutta etal, Attenuation of High amplitude noise in pre stack data and preservation of relative amplitude, P-90, SPG 2008, Hyderabad. Efrain Melendeze etal., Merging of several large 3D s- Issues & Solutions,Kelman Technologies Inc. 7