Waves and rays - II. Reflection and transmission. Seismic methods: Reading: Today: p Next Lecture: p Seismic rays obey Snell s Law

Similar documents
Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

INF-GEO Solutions, Geometrical Optics, Part 1

REFLECTION AND REFRACTION

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Types of Waves Transverse Shear. Waves. The Wave Equation

REFLECTION AND REFRACTION

CALCULUS BASIC SUMMER REVIEW

GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE. Hui QI, Yong SHI, Jingfu NAN

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 6 - Small Signal Stability

The Pendulum. Purpose

Repetition: Refractive Index

Lecture 11 Simple Linear Regression

Introduction to Signals and Systems, Part V: Lecture Summary

ADVANCED DIGITAL SIGNAL PROCESSING

Physics Supplement to my class. Kinetic Theory

567. Research of Dynamics of a Vibration Isolation Platform

The Scattering Matrix

OBJECTIVES. Chapter 1 INTRODUCTION TO INSTRUMENTATION FUNCTION AND ADVANTAGES INTRODUCTION. At the end of this chapter, students should be able to:

Free Space Optical Wireless Communications under Turbulence Channel Effect

Section 19. Dispersing Prisms

Section 19. Dispersing Prisms

1the 1it is said to be overdamped. When 1, the roots of

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

Cork Institute of Technology Bachelor of Science (Honours) in Applied Physics and Instrumentation-Award - (NFQ Level 8)

School of Mechanical Engineering Purdue University. ME375 Frequency Response - 1

Simple Linear Regression

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION

Principle Of Superposition

EE 505. Lecture 29. ADC Design. Oversampled

Scattering at an Interface:

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE

Lecture 3-7 Semiconductor Lasers.

EXPERIMENT OF SIMPLE VIBRATION

Numerical Methods in Fourier Series Applications

Chapter 7: The z-transform. Chih-Wei Liu

Assessment and Modeling of Forests. FR 4218 Spring Assignment 1 Solutions

PH 411/511 ECE B(k) Sin k (x) dk (1)

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5

A.1 Algebra Review: Polynomials/Rationals. Definitions:

Diversity Combining Techniques

Linear time invariant systems

CHAPTER NINE. Frequency Response Methods

EL 675 UHF Propagation for Modern Wireless Systems. Henry L. Bertoni Polytechnic University

Last time: Moments of the Poisson distribution from its generating function. Example: Using telescope to measure intensity of an object

ELEC1200: A System View of Communications: from Signals to Packets Lecture 3

II. Descriptive Statistics D. Linear Correlation and Regression. 1. Linear Correlation

Wind Energy Explained, 2 nd Edition Errata

Linear Regression Models

2C09 Design for seismic and climate changes

22. Theory of Multilayer Films

PH 411/511 ECE B(k) Sin k (x) dk (1)

THE MODELLING OF A REFLECTORLESS RANGE FINDER MAXIMUM RANGE WITH PHONG MODEL

Chapter 35 - Refraction. A PowerPoint Presentation by Paul E. Tippens, Professor of Physics Southern Polytechnic State University

Practical Spectral Anaysis (continue) (from Boaz Porat s book) Frequency Measurement

6.1 Analysis of frequency selective surfaces

Module 11: Applications : Linear prediction, Speech Analysis and Speech Enhancement Prof. Eliathamby Ambikairajah Dr. Tharmarajah Thiruvaran School

Parasitic Resistance L R W. Polysilicon gate. Drain. contact L D. V GS,eff R S R D. Drain

PHYS-3301 Lecture 10. Wave Packet Envelope Wave Properties of Matter and Quantum Mechanics I CHAPTER 5. Announcement. Sep.

DESCRIPTION OF THE SYSTEM

arxiv:physics/ v1 [physics.pop-ph] 29 Mar 2005

Voltage controlled oscillator (VCO)

COMM 602: Digital Signal Processing

Wave Phenomena Physics 15c

FIR Filter Design: Part II

Thin Film Interference

Chapter 2 Motion and Recombination of Electrons and Holes

Lecture 9: Diffusion, Electrostatics review, and Capacitors. Context

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and

Measurement uncertainty of the sound absorption

STP 226 ELEMENTARY STATISTICS

Global geophysics and wave propagation

Chapter 2 Systems and Signals

EEO 401 Digital Signal Processing Prof. Mark Fowler

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

Response Variable denoted by y it is the variable that is to be predicted measure of the outcome of an experiment also called the dependent variable

SPEC/4/PHYSI/SPM/ENG/TZ0/XX PHYSICS PAPER 1 SPECIMEN PAPER. 45 minutes INSTRUCTIONS TO CANDIDATES

Characterization and Simulation of HF Propagation in a realistic disturbed and horizontally inhomogeneous Ionosphere

Lecture Chapter 6: Convergence of Random Sequences

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

Chapter 2. Asymptotic Notation

UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA

Fast Monte Carlo Method in Stochastic Modelling of Charged Particle Multiplication

Two or more points can be used to describe a rigid body. This will eliminate the need to define rotational coordinates for the body!

Module 4. Signal Representation and Baseband Processing. Version 2 ECE IIT, Kharagpur

Synopsis of Euler s paper. E Memoire sur la plus grande equation des planetes. (Memoir on the Maximum value of an Equation of the Planets)

Dr. Seeler Department of Mechanical Engineering Fall 2009 Lafayette College ME 479: Control Systems and Mechatronics Design and Analysis

Chapter If n is odd, the median is the exact middle number If n is even, the median is the average of the two middle numbers

11 Correlation and Regression

Section 1 of Unit 03 (Pure Mathematics 3) Algebra

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare

Taylor Series (BC Only)

Chapter 2 Motion and Recombination of Electrons and Holes

Project 10.3 Vibrating Beams and Diving Boards

Time-Domain Representations of LTI Systems

Optics. n n. sin. 1. law of rectilinear propagation 2. law of reflection = 3. law of refraction

Analysis of MOS Capacitor Loaded Annular Ring MICROSTRIP Antenna

Transcription:

Seismic methods: Waves ad rays - II Readig: Today: p7-33 Net Lecture: p33-43 Reflectio ad trasmissio Seismic rays obey Sell s Law (just like i optics) The agle of icidece equals the agle of reflectio, ad the agle of trasmissio is related to the agle of icidece through the velocity ratio. si i si R si r

Reflectio ad trasmissio si i Seismic rays obey Sell s Law (just like i optics) The agle of icidece equals the agle of reflectio, ad the agle of trasmissio is related to the agle of icidece through the velocity ratio. But a coversio from to S or vice versa ca also occur. Still, the agles are determied by the velocity ratios. si R si r si R si r S S S S where p is the ray parameter ad is costat alog each ray. p Amplitudes reflected ad trasmitted The amplitude of the reflected, trasmitted ad coverted phases ca be calculated as a fuctio of the icidece agle usig Zoeppritz s equatios. Simple case: Normal icidece Reflectio coefficiet R C A A R i ρ ρ ρ + ρ Trasmissio coefficiet T C AT A i R C ρ ρ + ρ Reflectio ad trasmissio coefficiets for a specific impedace cotrast These coefficiets are determied by from the product of velocity ad desity the impedace of the material. R C usually small typically % of eergy is reflected.

Reflectio ad trasmissio You ca see: a direct wave, reflected ad trasmitted waves, plus multiples Normal move out (NMO) Reflectio from a sigle horizotal impedace cotrast: Arrival time or ( ) SR T h + T T0 + The arrival time curve is a hyperbola Note: a geophoe spread GG samples RR of the reflector. RR GG /

Normal move out (NMO) Arrival time curve is quadratic T T0 + So, if plot T vs. we ca determie the ad h from the slope ad itercept The importace of NMO Havig determied the layer velocity, we ca use the predicted quadratic shape to idetify reflectors The correct (shift traces) ad stack to ehace sigal to oise T NMO T T 0 T 0 Multiple layers Use Sell s Law to trace ray paths At each iterface si i si r p

NMO for layers Whe the offset is small w.r.t. reflector depth (<<h), the NMO curve is still a hyperbola T T0, + rms, where T 0, rms, k T 0, k k k T T k k k k h k k k k k3 Determie velocity structure oe layer at a time Diffractio From optics A sharp break i a reflector acts as a secodary source of a spherical wavefrot

Critical icidece si i si r whe >, r > i therefore, we ca icrease i util r 90 Whe r 90 i i C the critical agle si i C The critically refracted eergy travels alog the velocity iterface at cotiually refractig eergy back ito the upper medium at a agle i C a head wave Head wave Occurs due to a low to high velocity iterface Eergy travels alog the boudary at the higher velocity Eergy is cotiually refracted back ito the upper medium at a agle i C rovides costraits o the boudary depth e.g. Moho depth

Head wave You ca see: a head wave, trapped surface wave, divig body wave Factors affectig velocity Desity velocity typically icreases with desity (κ ad µ are depedat o ρ ad icrease more rapidly tha ρ) κ + 4 3 µ ρ S µ ρ orosity ad fluid saturatio Icreasig porosity reduces velocity. Fillig the porosity with fluid icreases the velocity. sat φ φ + F M oisso s ratio related to / S This is used to distiguish betwee rock/sedimet types. It is usually more sesitive tha just aloe. The sigificat variatios i sedimets are usually due to porosity variatios ad water saturatio. Water saturatio has o effect o S (for low porosities) but a sigificat effect o.

elocity ad desity Nafe-Drake curve sedimets ad sedimetary rocks igeous ad metamorphic rocks S This curve has bee approimated usig the epressio 4 ρ a (a is a costat: 670 whe ρ i km/m 3 ad i km/s) elocity ad desity Birch s Law A liear relatioship betwee velocity ad desity v aρ + b Three pressures Crust ad matle rock observatios

Typical rock velocity rages Usig velocity aloe to determie rock type is problematic to impossible. Seismic sources Rifles ad gus Cheap Repeatable fire ito water filled hole Shallow targets 0-50m Cosider Eergy iput Repeatability Cost Coveiece Sledge hammer Cheap Repeatable oce plate is stable (ad with traiig!) Targets 5-50m Weight drops Cheap Repeatable automated Targets > 50m

Seismic sources ibroseis No pulse, frequecy sweep Sigificat sigal with stackig/decovolutio Cosider Eergy iput Repeatability Cost Coveiece Eplosives arious sizes target depth Safety ad epese ca be a issue Air gus At sea ery repeatable Large array for big sigal Seismic receivers Geophoes Cylidrical coil suspeded i a magetic field The iertia of the coil causes motio relative to the maget geeratig a electrical sigal Geophoes are sesitive to velocity Istrumet respose The relatio betwee the iput groud motio ad the output electrical sigal Natural frequecy The frequecy which produces the maimum amplitude output Dampig Reduces the amplitude of the atural frequecy respose ad prevets ifiite oscillatios Wat a flat respose Hydrophoes Used at sea Use piezoelectric mierals to sese pressure variatios

Deploymet Importat cosideratios Need good couplig to the groud spike Mii-arrays to reduce surface wave oise Offset of geophoes Small offsets Near-vertical icidece retais -eergy High resolutio of subsurface reflectors Seismic reflectio aalysis Large offsets Improves velocity sesitivity rovides horizotal averages oly Seismic refractio aalysis