The ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.

Similar documents
(1) Cov(, ) E[( E( ))( E( ))]

Common MidPoint (CMP) Records and Stacking

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

Chapters 2 Kinematics. Position, Distance, Displacement

Real-Time Systems. Example: scheduling using EDF. Feasibility analysis for EDF. Example: scheduling using EDF

Cyclone. Anti-cyclone

8. Queueing systems lect08.ppt S Introduction to Teletraffic Theory - Fall

Fresnel Equations cont.

Determination of Antoine Equation Parameters. December 4, 2012 PreFEED Corporation Yoshio Kumagae. Introduction

Interval Estimation. Consider a random variable X with a mean of X. Let X be distributed as X X

Midterm Exam. Tuesday, September hour, 15 minutes

4. THE DENSITY MATRIX

FORCED VIBRATION of MDOF SYSTEMS

PHYS 1443 Section 001 Lecture #4

Continuous Time Markov Chains

Fresnel Dragging Explained

For the plane motion of a rigid body, an additional equation is needed to specify the state of rotation of the body.

Quantum Mechanics II Lecture 11 Time-dependent perturbation theory. Time-dependent perturbation theory (degenerate or non-degenerate starting state)

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits.

Least squares and motion. Nuno Vasconcelos ECE Department, UCSD

Fundamentals of Speech Recognition Suggested Project The Hidden Markov Model

Lecture 3 Topic 2: Distributions, hypothesis testing, and sample size determination

The Poisson Process Properties of the Poisson Process

Modeling and Predicting Sequences: HMM and (may be) CRF. Amr Ahmed Feb 25

QR factorization. Let P 1, P 2, P n-1, be matrices such that Pn 1Pn 2... PPA

Chapter 8. Simple Linear Regression

The Properties of Probability of Normal Chain

θ = θ Π Π Parametric counting process models θ θ θ Log-likelihood: Consider counting processes: Score functions:

Chapter Lagrangian Interpolation

Displacement, Velocity, and Acceleration. (WHERE and WHEN?)

Square law expression is non linear between I D and V GS. Need to operate in appropriate region for linear behaviour. W L

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

The Linear Regression Of Weighted Segments

AML710 CAD LECTURE 12 CUBIC SPLINE CURVES. Cubic Splines Matrix formulation Normalised cubic splines Alternate end conditions Parabolic blending

NEWTON S SECOND LAW OF MOTION

Optimal Eye Movement Strategies in Visual Search (Supplement)

Fault Tolerant Computing. Fault Tolerant Computing CS 530 Probabilistic methods: overview

72 Calculus and Structures

C(p, ) 13 N. Nuclear reactions generate energy create new isotopes and elements. Notation for stellar rates: p 12

Physics Notes - Ch. 2 Motion in One Dimension

HYPOTHESIS TESTING. four steps

INSTANTANEOUS VELOCITY

Solution set Stat 471/Spring 06. Homework 2

Bianchi Type II Stiff Fluid Tilted Cosmological Model in General Relativity

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

Chapter 12: Velocity, acceleration, and forces

NUMERICAL EVALUATION of DYNAMIC RESPONSE

Unit 1 Test Review Physics Basics, Movement, and Vectors Chapters 1-3

Final Exam Applied Econometrics

SYRIAN SEISMIC CODE :

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

t = s D Overview of Tests Two-Sample t-test: Independent Samples Independent Samples t-test Difference between Means in a Two-sample Experiment

PubH 7440 Spring 2010 Midterm 2 April

TP A.14 The effects of cut angle, speed, and spin on object ball throw

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

WebAssign HW Due 11:59PM Tuesday Clicker Information

EMA5001 Lecture 3 Steady State & Nonsteady State Diffusion - Fick s 2 nd Law & Solutions

Mathematical Formulation

UNIT 1 ONE-DIMENSIONAL MOTION GRAPHING AND MATHEMATICAL MODELING. Objectives

Solution of Impulsive Differential Equations with Boundary Conditions in Terms of Integral Equations

Linear Regression Linear Regression with Shrinkage

Motion in Two Dimensions

Solutions from Chapter 9.1 and 9.2

SPECIAL CONSIDERATIONS FOR VOLUMETRIC Z-TEST FOR PROPORTIONS

Mechanics of Materials CIVL 3322 / MECH 3322

Density estimation III.

Integral Equations and their Relationship to Differential Equations with Initial Conditions

PHYSICS 151 Notes for Online Lecture #4

Scattering at an Interface: Oblique Incidence

Panel Data Regression Models

1. Consider an economy of identical individuals with preferences given by the utility function

Section 8 Convolution and Deconvolution

The textbook expresses the stock price as the present discounted value of the dividend paid and the price of the stock next period.

14. Poisson Processes

Q2.4 Average velocity equals instantaneous velocity when the speed is constant and motion is in a straight line.

Lecture 3. Diffusion. Reading: Chapter 3

Centroids Method of Composite Areas

Partial Molar Properties of solutions

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

Big O Notation for Time Complexity of Algorithms

CS623: Introduction to Computing with Neural Nets (lecture-10) Pushpak Bhattacharyya Computer Science and Engineering Department IIT Bombay

Key words: Fractional difference equation, oscillatory solutions,

Decomposition of Hadamard Matrices

Centroids & Moments of Inertia of Beam Sections

One-Dimensional Kinematics

Chapter 3 Kinematics in Two Dimensions

More Digital Logic. t p output. Low-to-high and high-to-low transitions could have different t p. V in (t)

Review Equations. Announcements 9/8/09. Table Tennis

Advanced time-series analysis (University of Lund, Economic History Department)

SUMMATION OF INFINITE SERIES REVISITED

Chapter 3: Vectors and Two-Dimensional Motion

s in boxe wers ans Put

Conservation of Momentum. The purpose of this experiment is to verify the conservation of momentum in two dimensions.

FALL HOMEWORK NO. 6 - SOLUTION Problem 1.: Use the Storage-Indication Method to route the Input hydrograph tabulated below.

Matlab and Python programming: how to get started

Physics 101: Lecture 03 Kinematics Today s lecture will cover Textbook Sections (and some Ch. 4)

Lecture VI Regression

Extremal graph theory II: K t and K t,t

Lecture 3 Probability review (cont d)

Transcription:

C. Trael me cures for mulple reflecors The ray pahs ad rael mes for mulple layers ca be compued usg ray-racg, as demosraed Lab. MATLAB scrp reflec_layers_.m performs smple ray racg. (m) ref(ms) ref(ms) ref(ms) 0.000 5. 9. 9.048 0.000 5.707 9.455 9.0 60.000 54.89 9.89 9.98 90.000 56.695 94.46 9.60 0.000 59.96 95.76 0.047 50.000 54.6 96.67 0.608 80.000 546.659 97.698.94 0.000 55.94 99.69.0 40.000 556.807 94.074.07 70.000 56.88 94.9 4.094 00.000 569.589 945.95 5.74 0.000 576.98 947.904 6.575 60.000 584.8 950.64 8.000 90.000 59. 95.609 9.545

The wo eamples boh show hree reflecos. Noe ha: (a) Rays are refraced a erfaces accordg o Sell s Law (b) The polary ad amplude of a refleco s due o he mpedace chage a he erface (c) Deeper reflecos are flaer (.e. hey ehb less ormal moeou) (d) The hgh amplude refleco from a shallow erface ca weake laer reflecos (from deeper srucure) (e)for a sho gaher, he ray pahs reflec (bouce) a dffere pos o he reflecor surface. (f)ths geomery, wh geophoes o each sde of he sho s called a spl-spread

I he real sho gaher oe ha: (a) The ercal scale s wo-way rael me secods. (b) Pose arrals hae he race shaded black. Ths helps he eye o correlae races whe ose s prese. Beer someme helps wh hs process, bu oo much has bee kow o hae he reerse effec. Calbrao s eeded o opmze hs ask. (c) Noe he drec wae, whch s he frs arral. (d) The groud roll s log lasg ad obscures reflecos a small offse (A). (e) Slower arrals plo a a seeper agle o he sesmogram (more me eeded o rael he same dsace) (f) Reflecos ca be recogzed by her parabolc shape ad are labelled B,C,D, ad E. Noe ha hey are eer he frs arral. Ths meas ha reflecos are recorded whe he groud s sll mog from earler arrals, ad reflecos ca be obscured by ose.

Ieral elocy ad aerage elocy Cosder he mulplayer Earh model show aboe. The h layer has a hckess z ad a elocy. The ray raels hrough he layers o ad s refraced a each erface. I s he refleced by he h erface. The sesmc wae speds a me τ he h layer. The eral elocy s he acual elocy a specfc layer ad s defed as z / τ For he whole ray pah we ca defe a aerage elocy as V z τ τ τ Aoher way of aeragg he elocy s o use he roo-mea-square aerage. Ths s defed as: V rms, τ τ ad s eeded o compue eral eloces, as descrbed below. For he case of τ + τ V rms, τ + τ 4

Normal moe-ou for mulple layers ad he D equao To deerme he deph of he erfaces ad he eral eloces, a approach smlar o ha used C. s eeded. Noe ha we eed more ha jus he rael-me a ormal cdece, ad should use a plo of NMO s. offse (). Whe he sesmc sgals rael close o he ercal dreco, ca show ha he ormal moeou for he h refleco s: Δ Vrms, Ths s ery smlar o he equao C. for he ormal moeou aboe a sgle erface. The oly chage s ha he eral elocy ( ) has bee replaced by he r.m.s. elocy (.e. he aerage elocy dow o he h reflecor). I ca be show he roomea-square elocy, ges he correc aerage hs coe. Cosder he sesmogram wh hree reflecos ha s show aboe: For each refleco we ca measure he zero-offse rael mes (,, ). If ecessary a plo of ersus may be eeded (see C.) Ne he r.m.s. eloces are compued from he ormal moeou. Cosder he 4 h geophoe a dsace. The ormal moeou a hs alue of are Δ, Δ ad Δ. From hese NMO alues, he aerage (r.m.s.) eloces ca be compued. NMO Δ( ) ( ) V rms, Re-arragg hs equao ges V rms, Δ ( ) 5

For he oher reflecos s obous ha V rms, ( ad Δ ) V rms, Δ ( ) The e sage of he daa aalyss s o compue he eral (rue) eloces (,, ) from he aerage (r.m.s.) eloces (V rms,, V rms,, V rms,) The s acheed wh he D equao ha saes: rms, rms, Whle hs appears complcaed, s really a smple recurso equao. Ths meas ha f we kow he elocy oe layer, he equao ells us wha he elocy wll be he layer below. To sar hs process, we beg a he surface. For he frs refleco, he sesmc sgal raels oly Layer. Thus V rms, Now V rms, represes a aerage elocy of layers ad. Thus o fd we ca use he D equao wh rms, rms, Ths ca he be repeaed as may mes as ecessary. For he e layer dow we ca wre rms, rms, Eample : Look a he daa he able aboe for he hree refleco sesmogram Eample. The rael mes for he frs refleco are he same as for he eample C. The soluo o ha problem gae: 500 m s - ; z 400 m ad 0 0.5 s Noe ha V rms, 500 m s - 6

Cosder he secod refleco Choose a geophoe, Zero offse rael me, Normal moeou Δ r.m.s. elocy for Layers ad V rms, Use he D equao o compue he eral elocy for layer : rms, rms, Does hs agree wh he rue model? Now compue he hckess of he secod layer. To do hs oe ha he refleco speds a me - he secod layer. I hs me, raels a dsace z. Thus we ca wre z Re-arragg he equao ges z ( ) Repea hs for he hrd refleco Choose a geophoe, Zero offse rael me, Normal moeou Δ V rms, rms, rms, z ( ) 7

Mulple reflecos (a) The secod refleco s a mulple bouce from he erface a 0.4 km deph. (b) The rael me of he mulple a zero-offse s eacly double ha of he sgle refleco. Ths ca help defy a mulple. (c) Noe ha whe he sesmc eergy bouces a he Earh-Ar erface, he refleco coeffce ca be large. Remember ha R ρ ρ ρ + ρ Cosder he upward raellg wae as approaches he surface. ρ s ow he desy of he ar whch s esseally zero. Ths ges R - ad almos all eergy s refleced dowwards. (d) For a eample of mulples sgle chael sesmc daa, see Kearey Fgure 4.56 ad he dscusso C. 8