h 2 XTV inversion Rayfract help Page 1 of 6

Size: px
Start display at page:

Download "h 2 XTV inversion Rayfract help Page 1 of 6"

Transcription

1 XTV inersion Rayfract help Page of 6 BXTV inersion Delta-t-V and XTV inersion assume that traces are gathered by Common Mid-Point CMP (Diebold and Stoffa 98). Also, these inersion methods rely on the suppression of layer dip effects on apparent elocity, when sorting traces in a CMP gather by unsigned offset (Diebold and Stoffa 98). This suppression fails at locations of strong refractor curature. Barton and Baker (3) try to obtain improed apparent elocities by regarding the shot direction. To enable the realistic interpretation of sudden apparent elocity increase in CMP sorted traeltime cures, (Winkelmann 998) proposes to extend the Delta-t-V gradient layer inersion (Gebrande 986), with Dix inersion and Intercept Time inersion. These two additional inersion methods allow the modeling of constant-elocity layers. As described by (Winkelmann 998) and (Gawlas ), the XTV inersion reconstructs the D elocity s. depth function (z) below a CMP based on XTV data triples, with alues X = (reduced) unsigned offset, T = (reduced) time and V = apparent elocity. These data triples sample a (reduced) CMP sorted traeltime cure. The inersion is based on the layer stripping principle. With offset we always mean unsigned offset in the following. XTV inersion uses three separate methods, for inersion of a data triple into a model layer : Modified Dix inersion Intercept Time inersion Gradient layer inersion (original Delta-t-V method) The inersion starts with the first XTV triple at the smallest offset X, as determined from the original CMP cure. Once the first layer has been determined by one of aboe methods, offset and time for all other triples are reduced to the bottom of this first layer. Then the XTV triple at the next smallest offset X is inerted into a second layer, and remaining triples are reduced to the bottom of this second layer. This triple inersion process is continued iteratiely, until all XTV triples hae been processed. We first describe each of these three inersion methods in detail. We then describe the oerall XTV algorithm, which decides what inersion method should be applied to the current XTV triple. All inersion methods are described for the two-layer case only. Since the XTV inersion is based on layer stripping, we only need to deal with two layers (current oerburden, current refractor) during each iteration and triple-to-layer inersion step. The modified Dix inersion assumes reflection of a ray. The layer thickness h is determined as follows : h Vt is the unsigned offset X between shot point and receier. t is the traeltime T between shot point and receier, separated by offset. V is the measured apparent elocity, at the bottom of the layer i.e. at offset. ()

2 XTV inersion Rayfract help Page of 6 The aerage layer elocity is obtained with V t () Velocities at the top and at the bottom of the modeled layer are both set to this aerage elocity. The Intercept Time inersion assumes critical refraction of a ray, with both oerburden and basement layer haing a constant elocity. Using intercept time as determined from the XTV triple with t V (3) the layer thickness h is determined with h ( ) ( ) V (4) is the direct wae elocity (for first XTV triple) or the elocity as modeled for the bottom of the preiously determined layer (for reduced XTV triples). Velocities at the top and at the bottom of the modeled layer are both set to. The Gradient layer inersion method or Delta-t-V method assumes a diing wae ray and a constant elocity-gradient equal to a. So the layer's elocity-depth function is ( z) az (5) Since the elocity gradient is assumed constant, the diing first-break ray follows a circular arc, between shot point and receier. Based on the two equations

3 XTV inersion Rayfract help Page 3 of 6 a ( V) V (6) and t( V ) a arch( V ) (7) the elocity at the top of the modeled layer is determined numerically. Then the layer thickness h is h V V (8) and the elocity at the bottom of the modeled layer is set to the measured elocity V. can be smaller than the elocity measured at the bottom of the oerlying layer. So the gradient layer inersion can recognize elocity inersions, at least in some situations. Also, the gradient layer inersion does not use the intercept time. Layer stripping is done using equations as described by Gibson et al. (979). XTV triples are sorted by unsigned offset X, for layer inersion as described here. Winkelmann (998) proposes the following XTV inersion algorithm : Use Intercept Time inersion if the apparent elocity V increases suddenly, between adjacent XTV triples. The Minimum elocity ratio required for application of the Intercept Time layer inersion can be adjusted by the user. Intercept Time inersion can be disabled with XTV parameter Enable Intercept Time layer inersion. Otherwise use Dix inersion if both the aerage elocity (determined with Dix inersion) and the elocity at the top of the modeled layer (determined with Gradient layer inersion) are smaller than the elocity modeled for the bottom of the preiously determined layer. You may disable Dix inersion with XTV parameter Enable Modified Dix layer inersion. Otherwise use Gradient layer inersion. A candidate XTV triple is optionally rejected if its apparent elocity V or intercept time are larger than aerage alues for the next three XTV triples. This triple filtering has the aim of suppressing reflections erroneously picked as first breaks.

4 XTV inersion Rayfract help Page 4 of 6 Intercept Time layer inersion for multiple adjacent XTV triples as required for Intercept Time layer inersion is not clearly defined, especially if the preious layer has been obtained with Intercept Time inersion as well. You may suppress application of our Intercept Time inersion to multiple adjacent XTV triples, with XTV parameter Allow adjacent Intercept time layer inersion. If the preious XTV triple was inerted with Intercept Time layer inersion as well, then the elocity for the bottom of the preious layer may be assumed to be the preious, or the apparent elocity V of the preious XTV triple, or any alue between these two elocities. Alternatiely, can be determined by interpolating between the preious and the apparent elocity V of the current XTV triple. You may specify how should be determined with XTV parameters Oerlying layer elocity step and Current layer elocity step. If the resulting exceeds the preious apparent elocity V, used for inersion of the current XTV triple is reset to V of the preious triple. The XTV inersion ensures that the sum of Oerlying layer elocity step and Current layer elocity step does not exceed %. If you want to use apparent elocities V (from preious and current XTV triples) exclusiely and disregard preious, for determination of the current based on step parameters Oerlying layer elocity step and Current layer elocity step, just enable XTV option Prefer measured layer top elocity oer inerted. If this option is enabled, the apparent elocity V as obtained for the preious XTV triple and layer is taken as an estimate for the elocity at the top of the current layer. In analogy, the elocity at the top of the preious layer is estimated with the apparent elocity V of the preious-preious XTV triple. Parameters for XTV inersion Minimum elocity ratio The minimum elocity ratio (between apparent elocity V of the current XTV triple and the preious triple) required for application of the Intercept time layer inersion. If the actual ratio is smaller, the Intercept Time layer inersion method will not be applied. Instead, Dix inersion or Gradient layer inersion will be applied to the current XTV triple. Valid ratio alues range from. to.5. This ratio is regarded if XTV parameter Enable XTV Intercept Time layer inersion is checked only.

5 XTV inersion Rayfract help Page 5 of 6 Enable Intercept Time layer inersion Check this XTV option if you want to enable XTV Intercept Time layer inersion Enable Modified Dix layer inersion Check this XTV option if you want to enable Dix layer inersion Allow adjacent Intercept time layer inersion Check this XTV option to enable application of our Intercept Time inersion method, for multiple adjacent XTV triples. Oerlying layer elocity step This XTV parameter is used if Allow adjacent Intercept time layer inersion is enabled only. This step parameter may ary between alues % and %. The elocity step determines how needed for Intercept Time inersion of the current XTV triple is obtained, from the preious XTV triple and by interpolation between the preious (step %) and the preious apparent elocity V (step %). Current layer elocity step This XTV parameter is used if Allow adjacent Intercept time layer inersion is enabled only. This step parameter may ary between alues % and 99%. The elocity step determines how needed for Intercept Time inersion of the current XTV triple is obtained, by interpolation between as obtained with step parameter Oerlying layer elocity step (current step %) and the current apparent elocity V (current step %). Prefer measured layer top elocity oer inerted Check this XTV option to use apparent elocities V (belonging to preious and current XTV triples) exclusiely and disregard the preious, for determination of the current, based on step parameters Oerlying layer elocity step and Current layer elocity step. Oerburden elocity is needed for our Intercept Time two-layer case inersion method. If this option is enabled, the apparent elocity V as obtained for the preious XTV triple and layer is taken as an estimate for the elocity at the top of the current layer. Suppress elocity artefacts Enable this option to suppress the generation of processing artefacts, i.e. unrealistic elocity ariations. Use best for medium and high coerage profiles. See (Winkelmann 998), top of page 36. If enabled, a candidate ray will be used for modeling of an incremental layer if the ray specific apparent elocity and intercept time (as modeled by local regression on CMP cure, at ray specific offset) both are lower than the mean of apparent elocity and intercept time, as estimated for the next three higher CMP offsets. This triple filtering has the aim of suppressing reflections erroneously picked as first breaks. If this setting is disabled, no

6 XTV inersion Rayfract help Page 6 of 6 candidate ray selection, i.e. filtering / enforcing of CMP traeltime cure continuity, based on apparent elocity and intercept time will occur. References J. B. Diebold and P. L. Stoffa 98. The traeltime equation, tau-p mapping, and inersion of common midpoint data. SEG Geophysics, olume 46, pp Roland A. Winkelmann 998 Ph. D. Thesis. Entwicklung und Anwendung eines Wellenfelderfahrens zur Auswertung on CMP-sortierten Refraktionseinsaetzen. Akademischer Verlag Muenchen, Munich. ISBN Gawlas, Peter Florian Ph. D. Thesis. Möglichkeiten eines DMO-Prozesses in der CMP- Refraktionsseismik. LMU Munich: Faculty of Geosciences. See H. Gebrande 986. CMP-Refraktionsseismik. Paper presented (in German) at Mintrop Seminar / Uni-Kontakt Ruhr-Uniersitaet Bochum, Expanded abstract "Seismik auf neuen Wegen", pp H. Gebrande and H. Miller 985. Refraktionsseismik (in German). In: F. Bender (Editor), Angewandte Geowissenschaften II. Ferdinand Enke, Stuttgart; pp ISBN Bruce S. Gibson, Mark E. Odegard and George H. Sutton 979. Nonlinear least-squares inersion of traeltime data for a linear elocity-depth relationship. Geophysics, olume 44, pp Barton P. and Barker N. 3. Velocity imaging by tau-p transformation of refracted seismic traeltimes. Geophysical Prospecting, olume 5, pp Copyright (c) Intelligent Resources Inc. All rights resered.

The DeltatV 1D method for seismic refraction inversion: Theory

The DeltatV 1D method for seismic refraction inversion: Theory The DeltatV 1D method for seismic refraction inversion: Theory Siegfried R. Rohdewald, Intelligent Resources Inc. Suite 142, 757 West Hastings Street, B.C. V6C 1A1, Canada. E-mail: info@rayfract.com. Abstract

More information

JOANNEUM OT0608 refraction line : Smooth inversion vs. 1.5D XTV inversion :

JOANNEUM OT0608 refraction line : Smooth inversion vs. 1.5D XTV inversion : JOANNEUM OT0608 refraction line : Smooth inversion vs. 1.5D XTV inversion : OT0608 20% wavepath width, 20 WET iterations, RMS error 1.8 %, 1D-Gradient smooth initial model, Version 3.21 4 4 3 3 2 1 0 200

More information

Refraction Seismology. Chapter :: 6

Refraction Seismology. Chapter :: 6 Refraction Seismology Chapter :: 6 Snell s Law & Critical Refraction Because seismic sources radiate waes in all directions Some ray must hit interface at exactly the critical angle, i c This critically

More information

Creating a new Rayfract profile, importing the SEG-2 trace files and the Gremix header data :

Creating a new Rayfract profile, importing the SEG-2 trace files and the Gremix header data : RAYFRACT Tutorial specific to IGT line 14 : The seven SEG-2 formatted shot sorted binary trace data files and the Interpex Gremix.GRM file used to import first breaks as described below are available on

More information

The azimuth-dependent offset-midpoint traveltime pyramid in 3D HTI media

The azimuth-dependent offset-midpoint traveltime pyramid in 3D HTI media The azimuth-dependent offset-midpoint traeltime pyramid in 3D HTI media Item Type Conference Paper Authors Hao Qi; Stoas Alexey; Alkhalifah Tariq Ali Eprint ersion Publisher's Version/PDF DOI.9/segam3-58.

More information

Th Broadband Processing of Conventional Streamer Data - Optimized De-Ghosting in the Tau-P Domain

Th Broadband Processing of Conventional Streamer Data - Optimized De-Ghosting in the Tau-P Domain Th 08 14 Broadband Processing of Conentional Streamer Data - Optimized De-Ghosting in the Tau-P Domain H. Masoomzadeh* (TGS) & N. Woodburn (TGS) SUMMARY Conentional marine seismic data is affected by the

More information

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

GEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011 GEOPHYS 242: Near Surface Geophysical Imaging Class 5: Refraction Migration Methods Wed, April 13, 2011 Migration versus tomography Refraction traveltime and wavefield migration The theory of interferometry

More information

SUMMARY INTRODUCTION WHICH WAY?

SUMMARY INTRODUCTION WHICH WAY? Apparent horizontal in timelapse seismic images Dae Hale 1, Barbara Cox 2 and Paul Hatchell 2 1 Center for Wae Phenomena, Colorado School of Mines 2 Shell International E&P SUMMARY In addition to ertical

More information

Extension of delay time analysis for 3-D seismic refraction statics

Extension of delay time analysis for 3-D seismic refraction statics Extension of delay time analysis for 3-D Extension of delay time analysis for 3-D seismic refraction statics Jocelyn Dufour and Don C. Lawton ABSTRACT This current research is on the extension of delay

More information

An Evaluation of Methods and Available Software for Seismic Refraction Tomography Analysis

An Evaluation of Methods and Available Software for Seismic Refraction Tomography Analysis 21 An Evaluation of Methods and Available Software for Seismic Refraction Tomography Analysis Jacob R. Sheehan 1, William E. Doll 1 and Wayne A. Mandell 2 1 Oak Ridge National Laboratory, Oak Ridge, TN

More information

Geology 228 Applied Geophysics Lecture 4. Seismic Refraction (Reynolds, Ch. 4-6)

Geology 228 Applied Geophysics Lecture 4. Seismic Refraction (Reynolds, Ch. 4-6) Geology 8 Applied Geophysics Lecture 4 Seismic Refraction (Reynolds, Ch. 4-6) Seismic Methods Some fundamentals of seismic waes One dimensional wae equation the solution of a plane wae in unbounded uniform

More information

Geogiga Seismic Pro 8.0 Release Notes

Geogiga Seismic Pro 8.0 Release Notes Geogiga Seismic Pro 8.0 Release Notes Copyright 2015, All rights reserved. Table of Contents Introduction...1 Part 1 General Enhancements...3 Trace Display...4 Color Section Display...6 Curve Plotting...8

More information

Geogiga Seismic Pro 8.1 Release Notes

Geogiga Seismic Pro 8.1 Release Notes Geogiga Seismic Pro 8.1 Release Notes Copyright 2016, All rights reserved. Table of Contents Introduction...1 Part 1 General Enhancements...3 Well Data...4 Data Format...6 Data Import...7 Trace Magnifier...8

More information

Geogiga Seismic Pro 9.0 Release Notes

Geogiga Seismic Pro 9.0 Release Notes Geogiga Seismic Pro 9.0 Release Notes Copyright 2018, All rights reserved. Table of Contents Introduction...1 Part 1 New Modules...3 Modeling2D...3 Surface RT...4 Part 2 General Enhancements...5 Part 3

More information

Velocity Model Sampling and Interpolation: The Gradient Problem

Velocity Model Sampling and Interpolation: The Gradient Problem January 11, 2013 1 2 3 4 5 Kirchhoff migration uses raytracing equations for prestack depth migration. d x dt = v( x)2 p d p dt = v( x) v( x) These equations estimate the paths of seismic energy through

More information

Seismic Reflection Method

Seismic Reflection Method Seismic Reflection Method 1/GPH221L9 I. Introduction and General considerations Seismic reflection is the most widely used geophysical technique. It can be used to derive important details about the geometry

More information

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

Short Note. DMO velocity analysis with Jacubowicz s dip-decomposition method. David Kessler and Wai-Kin Chan* GEOPHYSICS, VOL. 58, NO. 10 (OCTOBER 1993); P. 1517-1524,9 FIGS. Short Note DMO velocity analysis with Jacubowicz s dip-decomposition method David Kessler and Wai-Kin Chan* INTRODUCTION Dip-moveout (DMO)

More information

Downloaded 09/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at A new parameterization for anisotropy update in full waeform inersion Nuno V. da Sila*, Andrew Ratcliffe, Graham Conroy, Vetle Vinje, Geoff Body, CGG Summary Full Waeform Inersion (FWI is now used regularly

More information

Seismic Imaging: Prestack

Seismic Imaging: Prestack Seismic Imaging: Prestack John C. Bancroft DEPARTMENT OF GEOLOGY & GEOPHYSICS, UNIVERSITY OF CALGARY, CALGARY, ALBERTA CSEG RECORDER NOV 2002 VOL. 26 NO. 09 Article Bios References Print Abstract The kinematics

More information

GeOlympus. GeoPACS GeoTR GeoStaR GeoWZ

GeOlympus. GeoPACS GeoTR GeoStaR GeoWZ GeOlympus GeoPACS GeoTR GeoStaR GeoWZ GeOlympus Workstation Work station Reflection time Refracted waves first arrivals Additional information Static and NMO corrections Layer by layer delay and statics

More information

Headwave Stacking in Terms of Partial Derivative Wavefield

Headwave Stacking in Terms of Partial Derivative Wavefield Geosystem Engineering, 7(1), 21-26 (March 2004) Headwave Stacking in Terms of Partial Derivative Wavefield Changsoo Shin School of Civil, Urban and Geosystem Engineering, Seoul National University, San

More information

Missing trace interpolation and its enhancement of seismic processes

Missing trace interpolation and its enhancement of seismic processes Missing trace interpolation Missing trace interpolation and its enhancement of seismic processes Wai-kin Chan and Robert R. Stewart ABSTRACT Many multi-channel seismic algorithms assume that the input

More information

Prestack Depth Migration with the Gabor Transform

Prestack Depth Migration with the Gabor Transform Prestac Depth Migration with the Gabor ransform Yongwang Ma* Uniersity of Calgary Calgary Alberta Canada yongma@ucalgary.ca and Gary Margrae Uniersity of Calgary Calgary Alberta Canada Abstract Summary

More information

Review of seismic imaging: Prestack

Review of seismic imaging: Prestack Review of seismic imaging: Prestack John C. Bancroft Published in the CSEG Recorder - November 001 Review of seismic migration: Prestack ABSTRACT The kinematics of prestack data considers an arbitrary

More information

SEG/New Orleans 2006 Annual Meeting

SEG/New Orleans 2006 Annual Meeting 3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in

More information

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

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary

More information

Seismic Modeling, Migration and Velocity Inversion

Seismic Modeling, Migration and Velocity Inversion Seismic Modeling, Migration and Velocity Inversion Full Waveform Inversion Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 18, 2014 Bee Bednar (Panorama Technologies)

More information

H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography

H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography M. Cavalca* (Schlumberger - Westerngeco) & R.P. Fletcher (Schlumberger - Westerngeco) SUMMARY Estimation of the

More information

It is widely considered that, in regions with significant

It is widely considered that, in regions with significant Multifocusing-based multiple attenuation Alex Berkovitch and Kostya Deev, Geomage Evgeny Landa, OPERA It is widely considered that, in regions with significant geologic complexity, methods which work directly

More information

Multiple suppression with land data: A case study

Multiple suppression with land data: A case study Multiple suppression with land data: A case study Gabriel Alvarez 1 ABSTRACT Some of the most important methods for multiple suppression are based on the moveout difference between the hyperbolas corresponding

More information

Genetic Algorithm for Seismic Velocity Picking

Genetic Algorithm for Seismic Velocity Picking Proceedings of International Joint Conference on Neural Networks, Dallas, Texas, USA, August 4-9, 2013 Genetic Algorithm for Seismic Velocity Picking Kou-Yuan Huang, Kai-Ju Chen, and Jia-Rong Yang Abstract

More information

Practical implementation of SRME for land multiple attenuation

Practical implementation of SRME for land multiple attenuation Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface

More information

Tunnel Excavation Disturbed Zone (EDZ) imaging with Rayfract version 3.35 :

Tunnel Excavation Disturbed Zone (EDZ) imaging with Rayfract version 3.35 : Tunnel Excavation Disturbed Zone (EDZ) imaging with Rayfract version 3.35 : Fig. 1 : left : Trace Shot gather, right : Refractor Shot breaks. Shows fit between picked times (solid colored curves) and modeled

More information

Reflection Seismic Processing

Reflection Seismic Processing Reflection Seismic Processing Objective Processing Systems General concept of CMP processing Processing Tools SEGY and similar file formats General features of ProMAX Reading: ProMAX manual (Introduction)

More information

On the Scattering Effect of Lateral Discontinuities on AVA Migration

On the Scattering Effect of Lateral Discontinuities on AVA Migration On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G

More information

Statics: the abusive power of trimming

Statics: the abusive power of trimming John C. Bancroft, Alan Richards, and Charles P. Ursenbach Statics: the abusive power of trimming ABSTRACT The application of trim statics directly to each seismic trace can be very dangerous as any seismic

More information

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

GG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material! GG450 April 6, 2010 Seismic Reflection I Today s material comes from p. 32-33 and 81-116 in the text book. Please read and understand all of this material! Back to seismic waves Last week we talked about

More information

Planning Land 3-D Seismic Surveys

Planning Land 3-D Seismic Surveys Planning Land 3-D Seismic Surveys Andreas Cordsen, Mike Galbraith, and John Peirce Edited by Bob A. Hardage Series Editor: Stephen J. Hill Geophysical Developments Series No. 9 Society of Exploration Geophysicists

More information

Downward propagation: using true near-surface velocity model, one-way wave- equation extrapolator:

Downward propagation: using true near-surface velocity model, one-way wave- equation extrapolator: Class 9: Beyond Seismic Statics Corrections Wed, Oct 7, 2009 Dynamic corrections: wave-equation and wavefield datuming methods Skipping the near-surface by interferometry (no model needed) Common-Focus

More information

Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples

Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint

More information

L 5 Seismic Method. Courtesy of ExxonMobil. Mitchum et al., 1977b

L 5 Seismic Method. Courtesy of ExxonMobil. Mitchum et al., 1977b Courtesy of ExxonMobil L 5 Seismic Method AAPG 1977 reprinted with permission of the AAPG whose permission is required for further use. Mitchum et al., 1977b Basic Exploration Workflow Identify Opportunities

More information

G009 Multi-dimensional Coherency Driven Denoising of Irregular Data

G009 Multi-dimensional Coherency Driven Denoising of Irregular Data G009 Multi-dimensional Coherency Driven Denoising of Irregular Data G. Poole* (CGGVeritas Services (UK) Ltd) SUMMARY Many land and ocean bottom datasets suffer from high levels of noise which make the

More information

Equivalent offset migration: the implementation and application update

Equivalent offset migration: the implementation and application update Equivalent offset migration: the implementation and application update Xinxiang Li, Yong Xu and John C. Bancroft INTRODUCTION We have some improvements about equivalent offset migration (EOM) method. In

More information

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

Anatomy of common scatterpoint (CSP) gathers formed during equivalent offset prestack migration (EOM) Anatomy of CSP gathers Anatomy of common scatterpoint (CSP) gathers formed during equivalent offset prestack migration (EOM) John C. Bancroft and Hugh D. Geiger SUMMARY The equivalent offset method of

More information

IXRefraX Tutorial. Version How to import and interpret seismic refraction data using IXRefraX 2005, 2006 Interpex Limited

IXRefraX Tutorial. Version How to import and interpret seismic refraction data using IXRefraX 2005, 2006 Interpex Limited IXRefraX Tutorial Version 1.02 How to import and interpret seismic refraction data using IXRefraX 2005, 2006 Interpex Limited Import & Interpretation Sequence Use File/Create Spread to import data. Use

More information

Geogiga Seismic Pro 8.3 Release Notes

Geogiga Seismic Pro 8.3 Release Notes Geogiga Seismic Pro 8.3 Release Notes Copyright 2017, All rights reserved. Table of Contents Introduction...1 Part 1 Utility Modules...2 Part 2 Reflection Modules...4 Updates in SF Imager...5 Updates in

More information

Anisotropy-preserving 5D interpolation by hybrid Fourier transform

Anisotropy-preserving 5D interpolation by hybrid Fourier transform Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,

More information

Target-oriented wavefield tomography: A field data example

Target-oriented wavefield tomography: A field data example Target-oriented wavefield tomography: A field data example Yaxun Tang and Biondo Biondi ABSTRACT We present a strategy for efficient migration velocity analysis in complex geological settings. The proposed

More information

Seismic refraction surveys

Seismic refraction surveys Seismic refraction surveys Seismic refraction surveys generate seismic waves that are refracted back to Earth s surface from velocity and density discontinuities at depth. Uses : Small-scale: geotechnical

More information

GEOPHYS 242: Near Surface Geophysical Imaging. Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011

GEOPHYS 242: Near Surface Geophysical Imaging. Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011 GEOPHYS 242: Near Surface Geophysical Imaging Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011 Wavefront tracing methods - speed versus accuracy Inversion algorithms - any magic approach?

More information

Applications of the tau-p Inversion Method. Shoba Maraj

Applications of the tau-p Inversion Method. Shoba Maraj Applications of the tau-p Inversion Method Shoba Maraj Final Report Theoretical Seismology Geop 523 Due 6th May 2011 Introduction Generally, velocity increases with depth in the Earth. They are related

More information

Master Thesis in Geophysics STEFAN JANSEN. Niels Bohr Institute, University of Copenhagen

Master Thesis in Geophysics STEFAN JANSEN. Niels Bohr Institute, University of Copenhagen PARAMETER INVESTIGATION FOR SUBSURFACE TOMOGRAPHY WITH REFRACTION SEISMIC DATA Master Thesis in Geophysics STEFAN JANSEN Niels Bohr Institute, University of Copenhagen Supervisors: Klaus Mosegaard, University

More information

WET inversion of reverse VSP survey, with Rayfract 3.25

WET inversion of reverse VSP survey, with Rayfract 3.25 WET inversion of reverse VSP survey, with Rayfract 3.25 Start up Rayfract Standard 3.25 via desktop icon. Select File New Profile. Set File name to ZIVKO13 and click Save button. Specify Station spacing

More information

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

Coherent partial stacking by offset continuation of 2-D prestack data Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced

More information

Th G Surface-offset RTM Gathers - What Happens When the Velocity is Wrong?

Th G Surface-offset RTM Gathers - What Happens When the Velocity is Wrong? Th G103 01 Surface-offset RTM Gathers - What Happens When the Velocity is Wrong? J.P. Montel* (CGG) SUMMARY Surface offset migrated common image gathers (SOCIGs) built by wave equation migration (WEM)

More information

High resolution residual velocity analysis and application

High resolution residual velocity analysis and application ABSTRACT Ye Zheng*, Graham Roberts and Andrew Ratcliffe Veritas DGC Inc. 2200-715 5th Avenue SW, Calgary, AB, T2P 5A2 Ye_zheng@veritasdgc.com Fidelity of velocity analysis is very important for seismic

More information

SUMMARY. denoise the original data at each iteration. This can be

SUMMARY. denoise the original data at each iteration. This can be A comparison of D reconstruction methods Aaron Stanton*, Nadia Kreimer, David Bonar, Mostafa Naghizadeh, and Mauricio Sacchi, Department of Physics, University of Alberta SUMMARY A comparison is made between

More information

y Input k y2 k y1 Introduction

y Input k y2 k y1 Introduction A regularization workflow for the processing of cross-spread COV data Gordon Poole 1, Philippe Herrmann 1, Erika Angerer 2, and Stephane Perrier 1 1 : CGG-Veritas, 2 : OMV Data regularization is critical

More information

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

Improved Imaging through Pre-stack Trace Interpolation for missing offsets of OBC data A case study from North Tapti area of West Coast, India P-256 Improved Imaging through Pre-stack Trace Interpolation for missing offsets of OBC data A case study from North Tapti area of West Coast, India M Lal, CPS Rana, Ramji Pathak, BN Bhatta, DP Sinha,

More information

When is anti-aliasing needed in Kirchhoff migration? a

When is anti-aliasing needed in Kirchhoff migration? a When is anti-aliasing needed in Kirchhoff migration? a a Published in SEP report, 80, 467-476 (1994) Dimitri Bevc and David E. Lumley ABSTRACT We present criteria to determine when numerical integration

More information

= 0) is the 2-D Fourier transform of the field (1),, z = 0). The phase ) is defined in the dispersion relation as where

= 0) is the 2-D Fourier transform of the field (1),, z = 0). The phase ) is defined in the dispersion relation as where GEOPHYSICS, VOL. 61, NO. 5 (SEPTEMBER-OCTOBER 1996); P. 1412 1416, 4 FIGS. Short Note Prestack migration by split-step DSR. Alexander Mihai Popovici* INTRODUCTION The double-square-root (DSR) prestack

More information

PS User Guide Series Utility

PS User Guide Series Utility PS User Guide Series 2015 Utility Prepared By Choon B. Park, Ph.D. January 2015 Table of Contents Page 1. Overview 2 2. Format SEG-2 Files 5 3. Format TEXT data 7 4. Make Common Seismic Gathers 9 5. Make

More information

Tutorial on Multichannel Analysis of Surface Wave (MASW) in the RadExPro Software

Tutorial on Multichannel Analysis of Surface Wave (MASW) in the RadExPro Software Tutorial on Multichannel Analysis of Surface Wave (MASW) in the RadExPro Software (Revision of 23.11.2016) DECO Geophysical SC Moscow State University Science Park 1-77 Leninskie Gory Moscow 119992, Russia

More information

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic

More information

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

SeisSpace Software. SeisSpace enables the processor to be able focus on the science instead of being a glorified data manager. SeisSpace Software OVERVIEW Delivering large volumes of data quickly and accurately remains the hallmark of any seismic processing system. With intuitive analysis tools, state-of-the-art geophysical algorithms

More information

Wave-equation migration from topography: Imaging Husky

Wave-equation migration from topography: Imaging Husky Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought

More information

Angle-dependent reflectivity by wavefront synthesis imaging

Angle-dependent reflectivity by wavefront synthesis imaging Stanford Exploration Project, Report 80, May 15, 2001, pages 1 477 Angle-dependent reflectivity by wavefront synthesis imaging Jun Ji 1 ABSTRACT Elsewhere in this report, Ji and Palacharla (1994) show

More information

Writing Kirchhoff migration/modelling in a matrix form

Writing Kirchhoff migration/modelling in a matrix form Writing Kirchhoff migration/modelling in a matrix form Abdolnaser Yousefzadeh and John C. Bancroft Kirchhoff migration matrix ABSTRACT Kirchhoff prestack migration and modelling are linear operators. Therefore,

More information

CRS stack using global simultaneous multi-parameter optimization: a better velocity analysis for near-surface data.

CRS stack using global simultaneous multi-parameter optimization: a better velocity analysis for near-surface data. CRS stack using global simultaneous multi-parameter optimization: a better velocity analysis for near-surface data. Z. Heilmann*, G. Satta* and G.P. Deidda** *Center for Advanced Studies, Research, and

More information

Modelling, migration, and inversion using Linear Algebra

Modelling, migration, and inversion using Linear Algebra Modelling, mid., and inv. using Linear Algebra Modelling, migration, and inversion using Linear Algebra John C. Bancroft and Rod Blais ABSRAC Modelling migration and inversion can all be accomplished using

More information

Ray based tomography using residual Stolt migration

Ray based tomography using residual Stolt migration Stanford Exploration Project, Report 11, November 11, 00, pages 1 15 Ray based tomography using residual Stolt migration Robert G. Clapp 1 ABSTRACT In complex areas, residual vertical movement is not an

More information

High-resolution Moveout Transform; a robust technique for modeling stackable seismic events Hassan Masoomzadeh* and Anthony Hardwick, TGS

High-resolution Moveout Transform; a robust technique for modeling stackable seismic events Hassan Masoomzadeh* and Anthony Hardwick, TGS High-resolution Moveout Transform; a robust technique for modeling stackable seismic events Hassan Masoomzadeh* and Anthony Hardwick, TGS Summary We propose a time-domain approach to transform a gather

More information

Building starting model for full waveform inversion from wide-aperture data by stereotomography

Building starting model for full waveform inversion from wide-aperture data by stereotomography Building starting model for full waveform inversion from wide-aperture data by stereotomography Vincent Prieux 1, G. Lambaré 2, S. Operto 1 and Jean Virieux 3 1 Géosciences Azur - CNRS - UNSA, France;

More information

Successes and challenges in 3D interpolation and deghosting of single-component marinestreamer

Successes and challenges in 3D interpolation and deghosting of single-component marinestreamer Successes and challenges in 3D interpolation and deghosting of single-component marinestreamer data James Rickett*, Schlumberger Gould Research Summary Combining deghosting with crossline interpolation

More information

IXRefraX Shootout Tutorial. Interpretation of Trimmed SAGEEP 2011 Shootout refraction tomography data using IXRefraX 2011 Interpex Limited

IXRefraX Shootout Tutorial. Interpretation of Trimmed SAGEEP 2011 Shootout refraction tomography data using IXRefraX 2011 Interpex Limited IXRefraX Shootout Tutorial Interpretation of Trimmed SAGEEP 2011 Shootout refraction tomography data using IXRefraX 2011 Interpex Limited Description There was a special Blind Test refraction tomography

More information

Multi-parameter Waveform Inversion with GPUs for the Cloud A Pipelined Implementation

Multi-parameter Waveform Inversion with GPUs for the Cloud A Pipelined Implementation Multi-parameter Waeform Inersion with GPUs for the Cloud A Pipelined Implementation Huy Le*, Stewart A. Lein, and Robert G. Clapp Geophysics Department, Stanford Uniersity March 28, 2018 Huy Le Multi-parameter

More information

Vidale, J., 1988, Finite-difference calculation of travel times, BSSA, Vo. 78, No. 6,

Vidale, J., 1988, Finite-difference calculation of travel times, BSSA, Vo. 78, No. 6, Class 4: First-Arrival Traveltime Tomography Mon, Sept 21, 2009 Wavefront tracing methods Inversion algorithms Model regularization Review of tomography case histories, presentation by Dr. Xianhuai Zhu

More information

VISTA. Desktop seismic data processing software

VISTA. Desktop seismic data processing software VISTA Desktop seismic data processing software VERSION 2017 VISTA desktop seismic data processing software VISTA Desktop seismic data processing software Comprehensive seismic processing and QC software

More information

SUMMARY. method to synthetic datasets is discussed in the present paper.

SUMMARY. method to synthetic datasets is discussed in the present paper. Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133

More information

Velocity-independent time-domain seismic imaging using local event slopes a

Velocity-independent time-domain seismic imaging using local event slopes a Velocity-independent time-domain seismic imaging using local event slopes a a Published in Geophysics, 72, no. 3, S139-S147, (2007) Sergey Fomel ABSTRACT I show that, by estimating local event slopes in

More information

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

Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave

More information

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

G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry D.F. Halliday* (Schlumberger Cambridge Research), P.J. Bilsby (WesternGeco), J. Quigley (WesternGeco) & E. Kragh (Schlumberger

More information

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

Machine learning assisted velocity auto-picking K.J. Smith, Shearwater GeoServices, OpenCPS K.J. Smith, Shearwater GeoServices, OpenCPS Summary Velocity auto-picking can reduce time spent on processing large volumes of seismic data and increase the number of CMP gathers that are picked in a project.

More information

Seismic Trace Interpolation Using Adaptive Prediction Filters

Seismic Trace Interpolation Using Adaptive Prediction Filters Seismic Trace Interpolation Using Adaptive Prediction Filters Mostafa Naghizadeh and Mauricio Sacchi Signal Analysis and Imaging Group (SAIG) Department of Physics University of Alberta CSEG 2008 Calgary

More information

3D predictive deconvolution for wide-azimuth gathers

3D predictive deconvolution for wide-azimuth gathers 3D predictive deconvolution for wide-azimuth gathers P.Hugonnet(1), J.L.Boelle(2), P.Herrmann(1), F.Prat(1), S.Navion(1) (1) CGGVeritas, (2) Total E&P Summary Usual pre-stack predictive deconvolution solutions

More information

An analysis of seismic information obtained from crooked line seismic surveys in crystalline rocks, Australia

An analysis of seismic information obtained from crooked line seismic surveys in crystalline rocks, Australia An analysis of seismic information obtained from crooked line seismic surveys in crystalline rocks, Australia Milovan Urosevic* and Christopher Juhlin *Curtin University of Technology, Perth, Australia

More information

Tau-p: A Plane Wave Approach to the Analysis of Seismic Data

Tau-p: A Plane Wave Approach to the Analysis of Seismic Data Tau-p: A Plane Wave Approach to the Analysis of Seismic Data MODERN APPROACHES IN GEOPHYSICS formerly Seismology and Exploration Geophysics VOLUME 8 Managing Editor: G. NOLET, Department of Theoretical

More information

Refractor 8.1 User Guide

Refractor 8.1 User Guide Refractor 8.1 User Guide Copyright 2016, All rights reserved. Table of Contents Preface...1 Conventions Used in This Guide...1 Where to Find Information...1 Technical Support...2 Feedback...2 Chapter 1

More information

Berryhill s pure migration = Gardner s PSI

Berryhill s pure migration = Gardner s PSI Stanford Exploration Project, Report 82, May 11, 2001, pages 1 136 Berryhill s pure migration = Gardner s PSI Alexander M. Popovici 1 ABSTRACT Kinematically, Berryhill s pure migration or nonimaging shot-record

More information

Progress Report on: Interferometric Interpolation of 3D SSP Data

Progress Report on: Interferometric Interpolation of 3D SSP Data Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine

More information

CFP migration; practical aspects

CFP migration; practical aspects 13 CFP migration; practical aspects Jan Thorbecke 1 13.1 Introduction In the past year the Common Focus Point (CFP) technology has become an important paradigm in the DELPHI research program. The interpretation

More information

1D internal multiple prediction in a multidimensional world: errors and recommendations

1D internal multiple prediction in a multidimensional world: errors and recommendations 1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate

More information

IXRefraX Shootout Tutorial. Interpretation of SAGEEP 2011 Shootout refraction tomography data using IXRefraX 2011 Interpex Limited

IXRefraX Shootout Tutorial. Interpretation of SAGEEP 2011 Shootout refraction tomography data using IXRefraX 2011 Interpex Limited IXRefraX Shootout Tutorial Interpretation of SAGEEP 2011 Shootout refraction tomography data using IXRefraX 2011 Interpex Limited Description There was a special Blind Test refraction tomography on Monday

More information

DW Tomo 8.1 User Guide

DW Tomo 8.1 User Guide DW Tomo 8.1 User Guide Copyright 2016, All rights reserved. Table of Contents Preface...1 Conventions Used in This Guide...1 Where to Find Information...1 Technical Support...2 Feedback...2 Chapter 1 Introducing

More information

A MATLAB version of Equivalent Offset Migration (EOM)

A MATLAB version of Equivalent Offset Migration (EOM) A MATLAB version of Equivalent Offset Migration (EOM) John C. Bancroft ABSTRACT Equivalent offset migration continues to be a viable method of prestack migration that is used in anisotropic prestack depth

More information

Introduction to the interpretation of seismic refraction data within REFLEXW

Introduction to the interpretation of seismic refraction data within REFLEXW Introduction to the interpretation of seismic refraction data within REFLEXW In the following the complete interpretation of seismic refraction data is described including import of the seismic data, picking

More information

Experiment 6. Snell s Law. Use Snell s Law to determine the index of refraction of Lucite.

Experiment 6. Snell s Law. Use Snell s Law to determine the index of refraction of Lucite. Experiment 6 Snell s Law 6.1 Objectives Use Snell s Law to determine the index of refraction of Lucite. Observe total internal reflection and calculate the critical angle. Explain the basis of how optical

More information

Picking up the First Arrivals in VSP Based on Edge Detection

Picking up the First Arrivals in VSP Based on Edge Detection 2012 International Conference on Image, Vision and Computing (ICIVC 2012) IPCSIT vol. 50 (2012) (2012) IACSIT Press, Singapore DOI: 10.7763/IPCSIT.2012.V50.7 Picking up the First Arrivals in VSP Based

More information

Reverse time migration in midpoint-offset coordinates

Reverse time migration in midpoint-offset coordinates Stanford Exploration Project, Report 111, June 9, 00, pages 19 156 Short Note Reverse time migration in midpoint-offset coordinates Biondo Biondi 1 INTRODUCTION Reverse-time migration (Baysal et al., 198)

More information

Refraction Full-waveform Inversion in a Shallow Water Environment

Refraction Full-waveform Inversion in a Shallow Water Environment Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui

More information