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

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

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

An investigation of the free surface effect

AVAZ and VVAZ practical analysis to estimate anisotropic properties

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

Phase-shift modelling for HTI media

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

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

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

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

Velocity and VTI anisotropy scanning in multicomponent seismic data

Exact elastic impedance in orthorhombic media

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

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

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

Acoustic Anisotropy Measurements and Interpretation in Deviated Wells

Constrained inversion of P-S seismic data

THE COMPOUND ANGLE IDENTITIES

Introduction to Seismology Spring 2008

Wave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces

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

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

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

An empirical method for estimation of anisotropic parameters in clastic rocks

Numerical Modeling for Different Types of Fractures

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

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

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

P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties

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

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.

Body-wave radiation patterns and AVO in transversely isotropic media

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

PEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity

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

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

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

Plane-wave migration in tilted coordinates

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

Shallow P and S velocity structure, Red Deer, Alberta

Elastic wavefield separation for VTI media

Borehole Geophysics. Acoustic logging measurements

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

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

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

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

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

Body-wave radiation patterns and AVO in transversely isotropic media

Seismology and Seismic Imaging

Radiation pattern in homogeneous and transversely isotropic attenuating media

3.2 Projectile Motion

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

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

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

Waveform inversion and time-reversal imaging in attenuative TI media

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

Phase and group velocity measurements from physically modeled transmission gathers

Integrals in cylindrical, spherical coordinates (Sect. 15.7)

Surface Waves and Free Oscillations. Surface Waves and Free Oscillations

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

Reflection and Transmission coefficient for VTI media

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

ANISOTROPIC PRESTACK DEPTH MIGRATION: AN OFFSHORE AFRICA CASE STUDY

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

Examples of prestack depth migration in TI media

Nonhyperbolic Reflection Moveout for Orthorhombic Media

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

Inversion for Anisotropic Velocity Parameter

SEG Houston 2009 International Exposition and Annual Meeting

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

Math Calculus II Homework # Due Date Solutions

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

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

Parametric Equations and Polar Coordinates

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.

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

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

Chapter 1. Introduction EARTH MODEL BUILDING

CHAPTER 4 Stress Transformation

ABSTRACT INTRODUCTION

NMO ellipse for a stratified medium with laterally varying velocity

Modelling of linearized Zoeppritz approximations

3-D description of normal moveout in anisotropic inhomogeneous media

( ) Trigonometric identities and equations, Mixed exercise 10

Reflection moveout and parameter estimation for horizontal transverse isotropy

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

Differential Equations: Homework 8

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates

Electromagnetic Waves Across Interfaces

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

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

On anelliptic approximations for qp velocities in VTI media

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

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

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

Edwin Soeryadjaya Problem Theoretical 3: Mirage

Math 120 Answers for Homework 14

Lecture 6, September 1, 2017

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

Pseudo-acoustic wavefield propagation for anisotropic media

Transcription:

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

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

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

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

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

[ ] ) ( 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

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

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

Angles and offsets included in the algorithm

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

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

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

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

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

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

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

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

δ=.5ε

δ=.5ε

δ=.75ε

δ= 1.ε

δ= 1.5ε

offset/depth δ= 1.5ε

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

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

For ε=.1 Displacement from NORSARD (m) Displacement from linear equations (m) Displacement from exact equations (m) δ=. 36.1 44.18 316.4 δ=.15 14.3 139.16 163.58 δ=.1 47 41.56 49.53 δ=.5-5 -56.5-49.39 δ=. -146-151.6-14.15 δ= -.5-44 -37-18.68 Table 1, NORSAR D experiments in VTI media, with ε=.1

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.

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.

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

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