Mapping the P-S Conversion Point in VTI Media. * Jianli Yang Don C. Lawton

Size: px
Start display at page:

Download "Mapping the P-S Conversion Point in VTI Media. * Jianli Yang Don C. Lawton"

Transcription

1 Mapping the P-S Conversion Point in VTI Media * Jianli Yang Don C. Lawton

2 Outline! Introduction! Theory! Numerical modeling methodology and results! NORSARD anisotropy ray mapping! Discussion and conclusions! Future work! Acknowledgement

3 Source MP Receiver P-wave S-wave The geometry of converted wave obeying Snell s law

4 Source MD Receiver P-wave S-wave P-S trajectory The conversion point traces a trajectory in the multilayered model

5 Source θ φ Elliptical wavefront Ray Spherical wavefront The definitions of the phase angle and ray angle

6 [ ] ) ( D θ εsin 1 α ) ( v * P θ θ + + = + = ) ( D β α θ εsin β α 1 β ) ( v * SV θ θ [ ] θ sin 1 β (θ) v SH γ + = = 1 θ sin ) α β (1 ε)ε α β 4(1 θ θcos sin ) α β (1 4δ 1 β α 1 1 ) ( D 1 4 * * θ [ ] ) dθ dv v tan θ (1 ) dθ dv v 1 (tan θ φ(θ) tan + = Thomsen s exact equations

7 v p v sv v sh () ( 4 θ = α 1 + δ sin θ cos θ + ε θ ) sin α cos β () θ = β 1 + ( ε δ ) sin θ θ () ( θ = β 1 + γ θ) sin tanφ = tanθ sinθ cosθ 1 dv v() θ dθ Thomsen s linear approximations

8 v ε = P ( π ) α α δ V = P (π 4) V P (π ) V P () V P () γ v = SH ( π ) β β Thomsen s definition of the anisotropy parameters

9 Angles and offsets included in the algorithm

10 Calculate the P- wave ray parameter for θ P Find the corresponding by Snell s law θ S VTI Calculate the and φ P V P Calculate the and φ S V S Isotropic Calculate XP Calculate XS XP + XS = offset X P (VTI) - X P (Isotropic) = displacement

11 δ=. δ=.1 δ=.5 ε=.1, exact equations

12 δ=. δ=. δ=.1 δ=.1 δ=.5 δ=.5 ε=.1, Thomsen s linear approximation

13 Source Receiver -1 - MP -3-4 Isotropic raypath m VTI raypath m ε=., δ=.5, offset/depth=1

14 Source Receiver MP Isotropic raypath VTI raypath ε=., δ=.1, offset/depth=1

15 Source Receiver MP VTI raypath Isotropic raypath ε=., δ=., offset/depth=1

16 Source Receiver MP Isotropic raypath VTI raypath ε=., δ=.15, offset/depth=1

17 Source Receiver VTI raypath MP Isotropic raypath ε=., δ=.5, offset/depth=1

18 δ=.5ε

19 δ=.5ε

20 δ=.75ε

21 δ= 1.ε

22 δ= 1.5ε

23 offset/depth δ= 1.5ε

24 P wave Isotropic case S wave Isotropic VTI The VTI model designed for NORSARD experiment

25 An example of the synthetic seismogram obtained from NORSARD anisotropy ray tracing on the model and displayed by PROMAX

26 For ε=.1 Displacement from NORSARD (m) Displacement from linear equations (m) Displacement from exact equations (m) δ= δ= δ= δ= δ= δ= Table 1, NORSAR D experiments in VTI media, with ε=.1

27 Discussion and Conclusions! The location of the conversion point in VTI media is different to that in the isotropic case.! The displacement of the conversion point is dependent on the offset/depth, velocity ratio, anisotropic parameters ε and δ.! When ε is greater than δ, the conversion point is displaced towards the source relative to its location in the isotropic case.

28 Discussion and Conclusions! When ε is less than δ, the conversion point moves towards the receiver compared to its location in isotropic case.! Results using linear approximations are similar to those obtained from NORSAR code.! Accurate placement of the conversion point is necessary for P-S survey design and data processing.

29 Future work! Further investigation of the relation between the displacement of the conversion point and Vp/Vs! Apply results of this work in the 3-C seismic survey design! Compare results using Thomsen s γ effective

30 Acknowledgements! We thank Dr. Larry Lines and Dr. Jim Brown for valuable suggestions! CREWES Sponsors financial support is also greatly appreciated

Mapping the conversion point in vertical transversely isotropic (VTI) media

Mapping the conversion point in vertical transversely isotropic (VTI) media Mapping the conversion point in vertical transversely isotropic (VTI) media Jianli Yang and Don. C. Lawton Conversion-point mapping ABSTRACT The important aspect of converted-wave (P-S) seismology is that

More information

Estimation of Thomsen s anisotropy parameters from compressional. and converted wave surface seismic traveltime data using NMO

Estimation of Thomsen s anisotropy parameters from compressional. and converted wave surface seismic traveltime data using NMO Important Notice This copy may be used only for the purposes of research and private study, and any use of the copy for a purpose other than research or private study may require the authorization of the

More information

An investigation of the free surface effect

An investigation of the free surface effect An investigation of the free surface effect Nasser S. Hamarbitan and Gary F. Margrave, An investigation of the free surface effect ABSTRACT When P and S seismic waves are incident on a solid-air interface

More information

AVAZ and VVAZ practical analysis to estimate anisotropic properties

AVAZ and VVAZ practical analysis to estimate anisotropic properties AVAZ and VVAZ practical analysis to estimate anisotropic properties Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Seismic anisotropic properties, such as orientation

More information

UNIVERSITY OF CALGARY. A comparison of different methods for estimating Thomsen's anisotropy parameters. Chunyan Xiao A THESIS

UNIVERSITY OF CALGARY. A comparison of different methods for estimating Thomsen's anisotropy parameters. Chunyan Xiao A THESIS Important Notice This copy may be used only for the purposes of research and private study, and any use of the copy for a purpose other than research or private study may require the authorization of the

More information

Phase-shift modelling for HTI media

Phase-shift modelling for HTI media HTI-Modelling Phase-shift modelling for HTI media R. K. Sharma and R. J. Ferguson ABSTRACT Fractures play an important role in hydrocarbon production as they determine the pathways and volume of crustal

More information

Estimation of stiffness coefficients of an orthorhombic physical model from group velocity measurements

Estimation of stiffness coefficients of an orthorhombic physical model from group velocity measurements Estimation of stiffness coefficients Estimation of stiffness coefficients of an orthorhombic physical model from group velocity measurements Faranak Mahmoudian, Gary Margrave, P.F. Daley, and Joe Wong

More information

AVAZ inversion for fracture orientation and intensity: a physical modeling study

AVAZ inversion for fracture orientation and intensity: a physical modeling study AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian*, Gary F. Margrave, and Joe Wong, University of Calgary. CREWES fmahmoud@ucalgary.ca Summary We present

More information

P-wave and S-wave near-surface characterization in NEBC. Liliana Zuleta and Don C. Lawton 1 st December, 2011

P-wave and S-wave near-surface characterization in NEBC. Liliana Zuleta and Don C. Lawton 1 st December, 2011 P-wave and S-wave near-surface characterization in NEBC Liliana Zuleta and Don C. Lawton 1 st December, 2011 Outline Objective heory and procedure Velocity and depth analysis SH data analysis / P-wave

More information

Far-field radiation from seismic sources in 2D attenuative anisotropic media

Far-field radiation from seismic sources in 2D attenuative anisotropic media CWP-535 Far-field radiation from seismic sources in 2D attenuative anisotropic media Yaping Zhu and Ilya Tsvankin Center for Wave Phenomena, Department of Geophysics, Colorado School of Mines, Golden,

More information

Velocity and VTI anisotropy scanning in multicomponent seismic data

Velocity and VTI anisotropy scanning in multicomponent seismic data PS multi-parameter scan Velocity and VTI anisotropy scanning in multicomponent seismic data Christopher O. Ogiesoba *, James E. Gaiser **, and Robert R. Stewart * ABSTRACT We present a prestack method

More information

Exact elastic impedance in orthorhombic media

Exact elastic impedance in orthorhombic media Exact elastic impedance in orthorhombic media F. Zhang (hina University of Petroleum), X.Y. Li (hina University of Petroleum, British Geological Survey) SUMMARY onventional elastic/ray impedance approximations

More information

Some comments on common-asymptotic-conversion-point (CACP) sorting of converted-wave data in isotropic, laterally inhomogeneous media

Some comments on common-asymptotic-conversion-point (CACP) sorting of converted-wave data in isotropic, laterally inhomogeneous media GEOPHYSICS, VOL. 70, NO. 3 (MAY-JUNE 2005); P. U29 U36, 5 FIGS. 0.90/.925750 Some comments on common-asymptotic-conversion-point (CACP) sorting of converted-wave data in isotropic, laterally inhomogeneous

More information

Pitfall in AVO Anisotropic Modeling: Plane Wave vs Spherical Wave. Qing Li May, 2003

Pitfall in AVO Anisotropic Modeling: Plane Wave vs Spherical Wave. Qing Li May, 2003 Pitfall in AVO Anisotropic Modeling: Plane Wave vs Spherical Wave Qing Li May, 2003 A Client s Model A client uses the following model to model anisotropic AVO effects and raised the question of discrepancy

More information

Exact Seismic Velocities for VTI and HTI Media and Extended Thomsen Formulas for Stronger Anisotropies. Abstract

Exact Seismic Velocities for VTI and HTI Media and Extended Thomsen Formulas for Stronger Anisotropies. Abstract Submitted to : Geophysics Exact Seismic Velocities for VTI and HTI Media and Extended Thomsen Formulas for Stronger Anisotropies James G. Berryman 1, 1 University of California, Lawrence Berkeley National

More information

Acoustic Anisotropy Measurements and Interpretation in Deviated Wells

Acoustic Anisotropy Measurements and Interpretation in Deviated Wells Acoustic Anisotropy Measurements and Interpretation in Deviated Wells X. M. Tang, and D. Patterson, Houston Technology Center, Baker Atlas, Houston, Texas, USA ABSTRACT Many acoustic anisotropy measurements

More information

Constrained inversion of P-S seismic data

Constrained inversion of P-S seismic data PS Inversion Constrained inversion of P-S seismic data Robert J. Ferguson, and Robert R. Stewart ABSTRACT A method to estimate S-wave interval velocity, using P-S seismic data is presented. The method

More information

THE COMPOUND ANGLE IDENTITIES

THE COMPOUND ANGLE IDENTITIES TRIGONOMETRY THE COMPOUND ANGLE IDENTITIES Question 1 Prove the validity of each of the following trigonometric identities. a) sin x + cos x 4 4 b) cos x + + 3 sin x + 2cos x 3 3 c) cos 2x + + cos 2x cos

More information

Introduction to Seismology Spring 2008

Introduction to Seismology Spring 2008 MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to

More information

Wave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces

Wave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces Lecture 5: Crystal Optics Outline 1 Homogeneous, Anisotropic Media 2 Crystals 3 Plane Waves in Anisotropic Media 4 Wave Propagation in Uniaxial Media 5 Reflection and Transmission at Interfaces Christoph

More information

Omm Al-Qura University Dr. Abdulsalam Ai LECTURE OUTLINE CHAPTER 3. Vectors in Physics

Omm Al-Qura University Dr. Abdulsalam Ai LECTURE OUTLINE CHAPTER 3. Vectors in Physics LECTURE OUTLINE CHAPTER 3 Vectors in Physics 3-1 Scalars Versus Vectors Scalar a numerical value (number with units). May be positive or negative. Examples: temperature, speed, height, and mass. Vector

More information

A simple way to improve PP and PS AVO approximations. Chuck Ursenbach CREWES Sponsors Meeting Thursday, December 1, 2005

A simple way to improve PP and PS AVO approximations. Chuck Ursenbach CREWES Sponsors Meeting Thursday, December 1, 2005 A simple way to improve PP and PS AVO approximations Chuck Ursenbach CREWES Sponsors Meeting Thursday, December 1, 2005 Overview Notes on spherical-wave modeling Reflectivity Explorer observations Theoretical

More information

TTI Anisotropic Depth Migration: Which Tilt Estimate Should We Use?

TTI Anisotropic Depth Migration: Which Tilt Estimate Should We Use? TTI Anisotropic Depth Migration: Which Tilt Estimate Should We Use? Francois Audebert* CGG Americas, Calgary, Alberta, Canada faudebert@cgg.com and Volker Dirks CGG Americas, Calgary, Alberta, Canada Abstract

More information

An empirical method for estimation of anisotropic parameters in clastic rocks

An empirical method for estimation of anisotropic parameters in clastic rocks An empirical method for estimation of anisotropic parameters in clastic rocks YONGYI LI, Paradigm Geophysical, Calgary, Alberta, Canada Clastic sediments, particularly shale, exhibit transverse isotropic

More information

Numerical Modeling for Different Types of Fractures

Numerical Modeling for Different Types of Fractures umerical Modeling for Different Types of Fractures Xiaoqin Cui* CREWES Department of Geoscience University of Calgary Canada xicui@ucalgary.ca and Laurence R. Lines Edward S. Krebes Department of Geoscience

More information

P185 TTI Anisotropic Depth Migration - Which Tilt Estimate Should We Use?

P185 TTI Anisotropic Depth Migration - Which Tilt Estimate Should We Use? P18 TTI Anisotropic Depth Migration - Which Tilt Estimate Should We Use? F.S. Audebert* (CGG Americas Inc.), A. Pettenati (Ecole Supérieure d' Electricité) & V. Dirks (CGG Americas Inc) SUMMARY We perform

More information

Snell s law in transversely isotropic media using linearized group velocities and related quantities

Snell s law in transversely isotropic media using linearized group velocities and related quantities Snell's law using group angles and velocities Snell s law in transversely isotropic media using linearized group velocities and related quantities P.F. Daley ABSTRACT Using a linearized approximation for

More information

Linearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada

Linearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada Shahpoor Moradi,Kristopher. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada SUMMRY Study of linearized reflectivity is very important for amplitude versus offset VO) analysis.

More information

P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties

P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties P.W. Wild* (Ikon Science Ltd), M. Kemper (Ikon Science Ltd), L. Lu (Ikon Science Ltd) & C.D. MacBeth (Heriot Watt University)

More information

Václav Bucha. Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague. SW3D meeting June 6-7, 2016 C OM S TR 3 D

Václav Bucha. Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague. SW3D meeting June 6-7, 2016 C OM S TR 3 D Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department

More information

Fermat s Principle. Fermat s Principle states that a ray of light in a medium will follow the path which takes the least amount of time.

Fermat s Principle. Fermat s Principle states that a ray of light in a medium will follow the path which takes the least amount of time. Homework Fermat s Principle Fermat s Principle states that a ray of light in a medium will follow the path which takes the least amount of time. Solution: The traversal time for the path is T = where ds

More information

Body-wave radiation patterns and AVO in transversely isotropic media

Body-wave radiation patterns and AVO in transversely isotropic media GEOPHYSICS, VOL. 60, NO. 5 (SEPTEMBER-OCTOBER 1995); P. 1409-1425, 11 FIGS. Downloaded 10/31/13 to 138.67.12.93. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/

More information

EXAM. Practice for Second Exam. Math , Fall Nov 4, 2003 ANSWERS

EXAM. Practice for Second Exam. Math , Fall Nov 4, 2003 ANSWERS EXAM Practice for Second Eam Math 135-006, Fall 003 Nov 4, 003 ANSWERS i Problem 1. In each part, find the integral. A. d (4 ) 3/ Make the substitution sin(θ). d cos(θ) dθ. We also have Then, we have d/dθ

More information

PEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity

PEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical

More information

Stanford Exploration Project, Report 123, October 31, 2005, pages 25 49

Stanford Exploration Project, Report 123, October 31, 2005, pages 25 49 Stanford Exploration Project, Report 123, October 31, 2005, pages 25 49 24 Stanford Exploration Project, Report 123, October 31, 2005, pages 25 49 Residual moveout in anisotropic angle-domain common image

More information

MATH 1080 Test 2 -Version A-SOLUTIONS Fall a. (8 pts) Find the exact length of the curve on the given interval.

MATH 1080 Test 2 -Version A-SOLUTIONS Fall a. (8 pts) Find the exact length of the curve on the given interval. MATH 8 Test -Version A-SOLUTIONS Fall 4. Consider the curve defined by y = ln( sec x), x. a. (8 pts) Find the exact length of the curve on the given interval. sec x tan x = = tan x sec x L = + tan x =

More information

We Challenges in shale-reservoir characterization by means of AVA and AVAZ

We Challenges in shale-reservoir characterization by means of AVA and AVAZ We-06-15 Challenges in shale-reservoir characterization by means of AVA and AVAZ N.C. Banik* (WesternGeco), M. Egan (WesternGeco), A. Koesoemadinata (WesternGeco) & A. Padhi (WesternGeco) SUMMARY In most

More information

Plane-wave migration in tilted coordinates

Plane-wave migration in tilted coordinates Stanford Exploration Project, Report 124, April 4, 2006, pages 1 16 Plane-wave migration in tilted coordinates Guojian Shan and Biondo Biondi ABSTRACT Plane-wave migration in tilted coordinates is powerful

More information

Azimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen

Azimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen Azimuthal AVO and Curvature Jesse Kolb* David Cho Kris Innanen Azimuthal AVO What is azimuthal AVO?: Analysis of incidence angle and azimuthal amplitude variations of reflection coefficients; Measures

More information

Shallow P and S velocity structure, Red Deer, Alberta

Shallow P and S velocity structure, Red Deer, Alberta Shallow P and S velocity structure, Red Deer, Alberta P & S velocity structure Don C. Lawton, Meredith A. McArthur, Rachel T. Newrick and Sarah E. Trend ABSTRACT A multioffset vertical seismic profile

More information

Elastic wavefield separation for VTI media

Elastic wavefield separation for VTI media CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically

More information

Borehole Geophysics. Acoustic logging measurements

Borehole Geophysics. Acoustic logging measurements Acoustic logging measurements - Review of basic physics background - Concept of P- and S-wave measurements and logging tools - Tube waves - Seismic imaging - Synthetic seismograms - Field application examples

More information

Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC)

Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC) Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC) Summary We develop an elastic full waveform inversion

More information

Elastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave

Elastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous

More information

Detecting fractures using time-lapse 3C-3D seismic data

Detecting fractures using time-lapse 3C-3D seismic data data Zimin Zhang, Don C. Lawton and Robert R. Stewart ABSTRACT This report presents the interpretation of time-lapse 3C-3D seismic data for fracture detection in a Saskatchewan potash mine. Seismic interpretation

More information

Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story

Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear

More information

The exact eikonals (Hamiltonians) of the coupled quasi-compressional ( qp ) and quasi-shear ( qs

The exact eikonals (Hamiltonians) of the coupled quasi-compressional ( qp ) and quasi-shear ( qs Snell s law in transversely isotropic media P.F. aley Snell s law in transversely isotropic media ABSTRACT The problem of reflection and transmission of waves at a plane boundary separating two transversely

More information

Body-wave radiation patterns and AVO in transversely isotropic media

Body-wave radiation patterns and AVO in transversely isotropic media GEOPHYSICS, VOL. 60, NO. 5 (SEPTEMBER-OCTOBER 1995); P. 1409-1425, 11 FIGS. Body-wave radiation patterns and AVO in transversely isotropic media llya Tsvankin* ABSTRACT The angular dependence of reflection

More information

Seismology and Seismic Imaging

Seismology and Seismic Imaging Seismology and Seismic Imaging 4. Ray theory N. Rawlinson Research School of Earth Sciences, ANU Seismology lecture course p.1/23 The ray approximation Here, we consider the problem of how body waves (P

More information

Radiation pattern in homogeneous and transversely isotropic attenuating media

Radiation pattern in homogeneous and transversely isotropic attenuating media Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010

More information

3.2 Projectile Motion

3.2 Projectile Motion Motion in 2-D: Last class we were analyzing the distance in two-dimensional motion and revisited the concept of vectors, and unit-vector notation. We had our receiver run up the field then slant Northwest.

More information

Anisotropic Depth Migration and High-Resolution Tomography in Gulf of Mexico: A Case History

Anisotropic Depth Migration and High-Resolution Tomography in Gulf of Mexico: A Case History Anisotropic Depth Migration and High-Resolution Tomography in Gulf of Mexico: A Case History Gary Rodriguez, Sherry Yang, Diane Yang, Quincy Zhang, Steve Hightower, TGS Summary We present a case study

More information

Seismic modelling and monitoring of carbon storage in a shallow sandstone formation

Seismic modelling and monitoring of carbon storage in a shallow sandstone formation Seismic modelling and monitoring of carbon storage in a shallow sandstone formation Virginia C. Vera*, University of Calgary, Calgary, Alberta vcvera@ucalgary.ca and Don C. Lawton, University of Calgary,

More information

An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models

An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models Item Type Conference Paper Authors Waheed Umair bin; Alkhalifah Tariq Ali Eprint version Pre-print

More information

Waveform inversion and time-reversal imaging in attenuative TI media

Waveform inversion and time-reversal imaging in attenuative TI media Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas

More information

Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS

Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy

More information

Phase and group velocity measurements from physically modeled transmission gathers

Phase and group velocity measurements from physically modeled transmission gathers Phase and group velocity measurements Phase and group velocity measurements from physically modeled transmission gathers Faranak Mahmoudian, Gary Margrave, and Joe Wong ABSTRACT Physical model data have

More information

Integrals in cylindrical, spherical coordinates (Sect. 15.7)

Integrals in cylindrical, spherical coordinates (Sect. 15.7) Integrals in clindrical, spherical coordinates (Sect. 15.7 Integration in spherical coordinates. Review: Clindrical coordinates. Spherical coordinates in space. Triple integral in spherical coordinates.

More information

Surface Waves and Free Oscillations. Surface Waves and Free Oscillations

Surface Waves and Free Oscillations. Surface Waves and Free Oscillations Surface waves in in an an elastic half spaces: Rayleigh waves -Potentials - Free surface boundary conditions - Solutions propagating along the surface, decaying with depth - Lamb s problem Surface waves

More information

DMO processing for mode-converted waves in a medium with a linear increase in velocity with depth

DMO processing for mode-converted waves in a medium with a linear increase in velocity with depth P-SV DMO processinq for a linear velocity gradient DMO processing for mode-converted waves in a medium with a linear increase in velocity with depth Shaowu Wang, John C. Bancroft, and Don C. Lawton ABSTRACT

More information

Reflection and Transmission coefficient for VTI media

Reflection and Transmission coefficient for VTI media Reflection and Transmission coefficient for VTI media R. K. Sharma and R. J. Ferguson ABSTRACT VTI, R & T coefficients Presently, we obtain the reflection (R and transmission (T coefficients of plane waves

More information

Basic Ray Tracing. Rick Aster and Sue Bilek. October 3, 2003

Basic Ray Tracing. Rick Aster and Sue Bilek. October 3, 2003 Basic Ray Tracing Rick Aster and Sue Bilek October 3, 3 A key observation that we can make about a seismic signal is its arrival time. From systematic observations of arrival times, we can deduce useful

More information

ANISOTROPIC PRESTACK DEPTH MIGRATION: AN OFFSHORE AFRICA CASE STUDY

ANISOTROPIC PRESTACK DEPTH MIGRATION: AN OFFSHORE AFRICA CASE STUDY Copyright 000 by the Society of Exploration Geophysicists ANISOTROPIC PRESTACK DEPTH MIGRATION: AN OFFSHORE AFRICA CASE STUDY Philippe Berthet *, Paul Williamson *, Paul Sexton, Joachim Mispel * * Elf

More information

3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie Zhang, Wei Zhang, University of Science and Technology of China (USTC)

3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie Zhang, Wei Zhang, University of Science and Technology of China (USTC) Downloaded 01/03/14 to 16.01.198.34. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ 3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie

More information

Examples of prestack depth migration in TI media

Examples of prestack depth migration in TI media Eamples of prestack depth migration in TI media Robert J. Ferguson and Gary F. Margrave ABSTRACT Wave field etrapolation by nonstationary phase shift can be formulated to allow velocity variation with

More information

Nonhyperbolic Reflection Moveout for Orthorhombic Media

Nonhyperbolic Reflection Moveout for Orthorhombic Media Nonhyperbolic Reflection Moveout for Orthorhombic Media AbdulFattah Al-Dajani and M. Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary Sciences Massachusetts Institute of

More information

Interval anisotropic parameters estimation in a least squares sense Case histories from West Africa

Interval anisotropic parameters estimation in a least squares sense Case histories from West Africa P-263 Summary Interval anisotropic parameters estimation in a least squares sense Patrizia Cibin*, Maurizio Ferla Eni E&P Division (Milano, Italy), Emmanuel Spadavecchia - Politecnico di Milano (Milano,

More information

Inversion for Anisotropic Velocity Parameter

Inversion for Anisotropic Velocity Parameter Chapter 5 Inversion for Anisotropic Velocity Parameter Ben Aggarwala 1, Ellis Cumberbatch 2, Jeff Grossman 3, Michael Lamoureux 4, Vlad Shapiro 5, Mark Solomonovitch 6, Paul Webster 7 This report describes

More information

SEG Houston 2009 International Exposition and Annual Meeting

SEG Houston 2009 International Exposition and Annual Meeting TTI/VTI anisotropy parameters estimation by focusing analysis, Part I: theory Jun ai*, Yang He, Zhiming Li, in Wang, Manhong Guo TGS-Nopec Geophysical ompany, itywest lvd. Suite, Houston, TX 7742, US Summary

More information

C3 Exam Workshop 2. Workbook. 1. (a) Express 7 cos x 24 sin x in the form R cos (x + α) where R > 0 and 0 < α < 2

C3 Exam Workshop 2. Workbook. 1. (a) Express 7 cos x 24 sin x in the form R cos (x + α) where R > 0 and 0 < α < 2 C3 Exam Workshop 2 Workbook 1. (a) Express 7 cos x 24 sin x in the form R cos (x + α) where R > 0 and 0 < α < 2 π. Give the value of α to 3 decimal places. (b) Hence write down the minimum value of 7 cos

More information

Math Calculus II Homework # Due Date Solutions

Math Calculus II Homework # Due Date Solutions Math 35 - Calculus II Homework # - 007.08.3 Due Date - 007.09.07 Solutions Part : Problems from sections 7.3 and 7.4. Section 7.3: 9. + d We will use the substitution cot(θ, d csc (θ. This gives + + cot

More information

AVAZ inversion for fracture orientation and intensity: a physical modeling study

AVAZ inversion for fracture orientation and intensity: a physical modeling study AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian and Gary F Margrave ABSTRACT AVAZ inversion We present a pre-stack amplitude inversion of P-wave data

More information

Prevailing-frequency approximation of the coupling ray theory for S waves

Prevailing-frequency approximation of the coupling ray theory for S waves Prevailing-frequency approximation of the coupling ray theory for S waves Petr Bulant & Luděk Klimeš Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague S EI S MIC

More information

Parametric Equations and Polar Coordinates

Parametric Equations and Polar Coordinates Parametric Equations and Polar Coordinates Parametrizations of Plane Curves In previous chapters, we have studied curves as the graphs of functions or equations involving the two variables x and y. Another

More information

1. Evaluate the integrals. a. (9 pts) x e x/2 dx. Solution: Using integration by parts, let u = x du = dx and dv = e x/2 dx v = 2e x/2.

1. Evaluate the integrals. a. (9 pts) x e x/2 dx. Solution: Using integration by parts, let u = x du = dx and dv = e x/2 dx v = 2e x/2. MATH 8 Test -SOLUTIONS Spring 4. Evaluate the integrals. a. (9 pts) e / Solution: Using integration y parts, let u = du = and dv = e / v = e /. Then e / = e / e / e / = e / + e / = e / 4e / + c MATH 8

More information

P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients

P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),

More information

- 1 - θ 1. n 1. θ 2. mirror. object. image

- 1 - θ 1. n 1. θ 2. mirror. object. image TEST 5 (PHY 50) 1. a) How will the ray indicated in the figure on the following page be reflected by the mirror? (Be accurate!) b) Explain the symbols in the thin lens equation. c) Recall the laws governing

More information

Chapter 1. Introduction EARTH MODEL BUILDING

Chapter 1. Introduction EARTH MODEL BUILDING Chapter 1 Introduction Seismic anisotropy in complex earth subsurface has become increasingly important in seismic imaging due to the increasing offset and azimuth in modern seismic data. To account for

More information

CHAPTER 4 Stress Transformation

CHAPTER 4 Stress Transformation CHAPTER 4 Stress Transformation ANALYSIS OF STRESS For this topic, the stresses to be considered are not on the perpendicular and parallel planes only but also on other inclined planes. A P a a b b P z

More information

ABSTRACT INTRODUCTION

ABSTRACT INTRODUCTION Migration velocity analysis for P-S data Miqrationvelocity analysis for P-S data Wai-Kin Chan and Robert R. Stewart ABSTRACT Migration velocity analysis (MVA) in depth (Reshef, 1992) has been successfully

More information

NMO ellipse for a stratified medium with laterally varying velocity

NMO ellipse for a stratified medium with laterally varying velocity CWP-685 NMO ellipse for a stratified medium with laterally varying velocity Mamoru Takanashi 1, 2 & Ilya Tsvankin 1 1 Center for Wave Phenomena, Geophysics Department, Colorado School of Mines, Golden,

More information

Modelling of linearized Zoeppritz approximations

Modelling of linearized Zoeppritz approximations Modelling of linearized Zoeppritz approximations Arnim B. Haase Zoeppritz approximations ABSTRACT The Aki and Richards approximations to Zoeppritz s equations as well as approximations by Stewart, Smith

More information

3-D description of normal moveout in anisotropic inhomogeneous media

3-D description of normal moveout in anisotropic inhomogeneous media GEOPHYSICS, VOL. 63, NO. 3 (MAY-JUNE 1998); P. 1079-1092,10 FIGS. 3-D description of normal moveout in anisotropic inhomogeneous media Vladimir Grechka* and Ilya Tsvankin* ABSTRACT We present a new equation

More information

( ) Trigonometric identities and equations, Mixed exercise 10

( ) Trigonometric identities and equations, Mixed exercise 10 Trigonometric identities and equations, Mixed exercise 0 a is in the third quadrant, so cos is ve. The angle made with the horizontal is. So cos cos a cos 0 0 b sin sin ( 80 + 4) sin 4 b is in the fourth

More information

Reflection moveout and parameter estimation for horizontal transverse isotropy

Reflection moveout and parameter estimation for horizontal transverse isotropy GEOPHYSICS, VOL. 62, NO. 2 (MARCH-APRIL 1997); P. 614 629, 7 FIGS. Reflection moveout and parameter estimation for horizontal transverse isotropy Ilya Tsvankin ABSTRACT Transverse isotropy with a horizontal

More information

Calculus First Semester Review Name: Section: Evaluate the function: (g o f )( 2) f (x + h) f (x) h. m(x + h) m(x)

Calculus First Semester Review Name: Section: Evaluate the function: (g o f )( 2) f (x + h) f (x) h. m(x + h) m(x) Evaluate the function: c. (g o f )(x + 2) d. ( f ( f (x)) 1. f x = 4x! 2 a. f( 2) b. f(x 1) c. f (x + h) f (x) h 4. g x = 3x! + 1 Find g!! (x) 5. p x = 4x! + 2 Find p!! (x) 2. m x = 3x! + 2x 1 m(x + h)

More information

Differential Equations: Homework 8

Differential Equations: Homework 8 Differential Equations: Homework 8 Alvin Lin January 08 - May 08 Section.6 Exercise Find a general solution to the differential equation using the method of variation of parameters. y + y = tan(t) r +

More information

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates CBE 6333, R. Levicky 1 Orthogonal Curvilinear Coordinates Introduction. Rectangular Cartesian coordinates are convenient when solving problems in which the geometry of a problem is well described by the

More information

Electromagnetic Waves Across Interfaces

Electromagnetic Waves Across Interfaces Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph

More information

9.1. Click here for answers. Click here for solutions. PARAMETRIC CURVES

9.1. Click here for answers. Click here for solutions. PARAMETRIC CURVES SECTION 9. PARAMETRIC CURVES 9. PARAMETRIC CURVES A Click here for answers. S Click here for solutions. 5 (a) Sketch the curve b using the parametric equations to plot points. Indicate with an arrow the

More information

Synthetic seismograms, Synthetic sonic logs, Synthetic Core. Larry Lines and Mahbub Alam

Synthetic seismograms, Synthetic sonic logs, Synthetic Core. Larry Lines and Mahbub Alam Synthetic seismograms, Synthetic sonic logs, Synthetic Core Larry Lines and Mahbub Alam Synthetic Seismograms Synthetic seismograms range from 1-D model seismograms that are least general but most economical

More information

On anelliptic approximations for qp velocities in VTI media

On anelliptic approximations for qp velocities in VTI media Geophysical Prospecting, 2004, 52, 247 259 On anelliptic approximations for qp velocities in VTI media Sergey Fomel Bureau of Economic Geology, The University of Texas, Austin, University Station, Box

More information

AVAZ inversion for fracture orientation and intensity: A physical modeling study. Faranak Mahmoudian Gary Margrave

AVAZ inversion for fracture orientation and intensity: A physical modeling study. Faranak Mahmoudian Gary Margrave AVAZ inversion for fracture orientation and intensity: A physical modeling study Faranak Mahmoudian Gary Margrave Objective Fracture orientation: direction of fracture planes Fracture intensity: number

More information

ScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach

ScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 4305 4312 GHGT-12 Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including

More information

COMPARISON OF OPTICAL AND ELASTIC BREWSTER S ANGLES TO PROVIDE INVUITIVE INSIGHT INTO PROPAGATION OF P- AND S-WAVES. Robert H.

COMPARISON OF OPTICAL AND ELASTIC BREWSTER S ANGLES TO PROVIDE INVUITIVE INSIGHT INTO PROPAGATION OF P- AND S-WAVES. Robert H. COMPARISON OF OPTICAL AND ELASTIC BREWSTER S ANGLES TO PROVIDE INVUITIVE INSIGHT INTO PROPAGATION OF P- AND S-WAVES Robert H. Tatham Department of Geological Sciences The University of Texas at Austin

More information

Edwin Soeryadjaya Problem Theoretical 3: Mirage

Edwin Soeryadjaya Problem Theoretical 3: Mirage The refractive index of the air varies with temperature. Cold air is denser than warm air and has therefore a greater refractive index. Thus a temperature gradient in the atmosphere is always associated

More information

Math 120 Answers for Homework 14

Math 120 Answers for Homework 14 Math 0 Answrs for Homwork. Substitutions u = du = d d = du a d = du = du = u du = u + C = u = arctany du = +y dy dy = + y du b arctany arctany dy = + y du = + y + y arctany du = u du = u + C = arctan y

More information

Lecture 6, September 1, 2017

Lecture 6, September 1, 2017 Engineering Mathematics Fall 07 Lecture 6, September, 07 Escape Velocity Suppose we have a planet (or any large near to spherical heavenly body) of radius R and acceleration of gravity at the surface of

More information

LECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # & 2 #

LECTURE 5 - Wave Equation Hrvoje Tkalčić  2 # & 2 # LECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # "t = ( $ + 2µ ) & 2 # 2 % " 2 (& ' u r ) = µ "t 2 % & 2 (& ' u r ) *** N.B. The material presented in these lectures is from the principal textbooks, other

More information

Pseudo-acoustic wavefield propagation for anisotropic media

Pseudo-acoustic wavefield propagation for anisotropic media Anisotropic wavefield propagation Pseudo-acoustic wavefield propagation for anisotropic media Ben D. Wards, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Reverse-time migration (RTM) is a powerful

More information