Ray tracing in geophysics
|
|
- Kathlyn Maxwell
- 5 years ago
- Views:
Transcription
1 Ray tracing in geophysics Seminar Author: Klemen Kunstelj Adviser: dr. Andrej Gosar and dr. Jure Bajc
2 Abstract Alp 2002 is an international seismological project. The main issue is to construct the structure of the Lithosphere in the area of Suth-Eastern Alps. The experimental part took place in July 2002 and the analysis of these results is planned to be completed aproximately in two years. In the seminar some basic modelling methods that are used in the analysis of the Alp2002 data are presented. Figure 1: Network of geophones with seismograms in Slovenia during Alp2002
3 Figure 2: Seismograms of an explosion in Vojnik during Alp2002 gathered on the network of seismographs in Slovenia 1
4 CONTENTS 2 Contents 1 Introduction 3 2 Ray theory and seismic waves Wave equations for P- and S-waves Ray equation Ray propagation in media of constant velocity Ray tracing with RAYAMP 6 4 Ray tracing with JIVE3D Model parametrization The forward problem Two-point ray tracing Ray perturbation theory (RPT) Frechet derivatives The inverse problem Conclusion 14
5 1 Introduction 3 1 Introduction The main issue of Alp 2002 is to discover in detail the tectonics and geodynamics at the junction of Europian, Adriatic (fragment of African tectonic plate) and Tisza plates (contact zone between South-Eastern Alps, Dinarides and Panonian basin) by using combined seismic refraction/wide-angle reflection method. For this two techniques the maximal distance between the source (explosion) and the targets (geophone with seismograph) is around few hundred km. The velocity of seismic waves is in general increasing with depth. Waves are refracted at discontinuities and travel in deeper regions almost horizontally (refraction) or they are reflected under blunt angle (wide-angle reflection). Because seismic waves travel long distances, projects as Alp 2002 give us information about the lithosphere. The Earth Lithosphere is in general characterized by three discontinuities: Moho discontinuity between Earth s crust and mantle. It is named after Croatian seismologist Mohorovičič, who discovered it by analysing the earthquake in 1909 in the region of Kolpa. Moho is characterized by an increase of P-wave velocity (compresional waves or prima waves) from km/s in the crust to km/s in the mantle, and an increase of S-wave velocity (shear waves) from km/s to 4.8 km/s.moho is situated km below continents, 5 8 km below the ocean floor, and km below certain mountains regions. Conrad discontinuity between upper part (granite) and lower part (bazalt) of the crust, is located at depth km and is characterized by an increase of P-wave velocity from km/s to 6.5 km/s. It is not observed everywhere. Third discontinuity which is often seen on a seismic data is a junction between the sediments and magmathic base. In the framework of the Alp2002 project a network of 12 profiles of total length 4100 km was spreading over seven countries (Austria, Chech republic, Hungary, Croatia, Slovenia, Italy and Germany) portable seismographs were deployed and 31 strong (300 kg) explosions fired. In Slovenia 127 seismographs were placed along five profiles of total length 575 km and two explosions were fired; one near Vojnik and the other near Gradin. The collected data will allow construction of the three-dimensional model of the lithosphere.
6 2 Ray theory and seismic waves 4 2 Ray theory and seismic waves 2.1 Wave equations for P- and S-waves Displacements u(x, t) in Navier s equation ρ 2 u = (λ + µ) ( u) + µ 2 u t 2 form a vector field in an elastic medium. We use Helmholtz theorem in order to represent u in terms of a curl-free scalar potential φ and divergence-free vector potential ψ u = φ + ψ. (1) The Navier equation may than be re-written as [ (λ + 2µ) 2 φ ρ 2 φ t 2 ] + [ µ 2 ψ ρ 2 ψ t 2 ] = 0, (2) where ρ, λ and µ are taken to be costant troughout the medium and the change in the gravitational field strength are ignored. Equation (2) is setisfied, if 2 φ 1 2 φ λ + 2µ α 2 t = 0, α = (3) 2 ρ and 2 ψ 1 2 ψ µ β 2 t = 0, β = 2 ρ. (4) Equations (3) and (4) are two independent wave equations that describe compressional waves traveling at the P -wave velocity α and shear waves which travel at the S-wave velocity β. They represent a good aproximation to the equation of motion for heterogeneous media provided that the elastic propertis of the medium do not change significantly over a single wavelength. If the elastic properties do change rapidly, both wave equations become coupled and P - and S-waves may no longer be considered to propagate independently. In the layer-interface formalism both wave equations are normally treated as independent within each layer, but converted waves may be produced at intefaces, where there is a sudden change in the elastic properties of the medium. 2.2 Ray equation The procedure of execution the ray equation is similar in geophysics and optics. The solution of (3) for harmonic waves of frequency ω with constant amplitude and zero initial phase can be written (for P -waves) φ(x, t) = φ 0 (x)e iωt (5)
7 2.2 Ray equation 5 without the loss of generality. If φ 0 (x) is defined as φ 0 (x) = A(x)e ik 0S(x), (6) where k 0 = ω/α 0 and α 0 is a reference wave speed, equation (5) will form a solution to the wave equation provided that certain constraints on A(x) and S(x) (called the eikonal) describe the spatial variation of amplitude and phase within a general wave field. Let us now introduce the concept of a ray, which describes the path of a wave packet through the isotropic medium, being at all times perpendicular to the wavefront. It is clear that the unit vector ˆp = α α 0 S (7) is always perpendicular to the wavefront in the direction of the increasing phase. Like the travel time T (x), which is defined as a time taken for wavefront to travel from a reference point x 0 to the arbitrary point x, the eikonal is defined relative to the phase at the reference point S(x) = α 0 T (x) and the eikonal equation can be re-written in terms of the travel time and the wave speed v(x), ( S) 2 = α2 0 = ( T ) 2 = 1 α 2 v(x). (8) 2 The ray path may be described by the function x(s) where s is the curvilinear distance from the reference point along the ray path. Let us also define a second function p(s) for which p = T. This is called the slowness vector, because its magnitude at a point x(s) is equal to the reciprocal of the velocity at that point. Applying the condition that the ray path is orthogonal to the wavefronts yields dx = v T = vp. (9) ds By differentiating the eikonal equation (8) with respect to s and combining the result with (9) we obtain the ray equation ( ) d 1 dx = ds v ds ( 1 v ). (10) Once the ray trajectory has been obtained, travel times may be calculated by integrating ( ) 1 T = T ds = ds. (11) raypath raypath v(x(s)) So the travel time from A to B is equal to the travel time from B to A. This principle of reciprocity may be used to improve the efficiency of ray-tracing for certain problem types.
8 2.3 Ray propagation in media of constant velocity Ray propagation in media of constant velocity In this case we will take a layer of thickness H and velocity v over a halfspace of velocity v. This situation allows the study of ray trajectories and travel times due to the presence of a surface of contact between the layer and a half-space. The thickness of the layer is the parameter that scales a spatial dimension of the model. For this reason, the application of ray theory is valid, if the wave lengths are much smaller than the thickness of the layer. This example represents, in a simplified way, the situation of the Earth s crust over the upper mantle for small distances, for which the flat-earth approximation is valid. We consider the case v < v and the focus F at the Figure 3: Seismic waves surface (h = 0)(figure 3). There are three types of rays travelling from F to P, which is a distance x from the focus F: (a) direct rays; (b) rays reflected on the surface between the layer and the half-space (FCP); and (c) criticaly refracted rays or head waves, i.e. rays with a critical angle of incidence i c, propagate a certain distance horizontally through the half-space and come back to the free surface with the same angle of incidence i c (FBDP). In the next sections of this seminar I will present some possible approaches of modelling the seismic data. 3 Ray tracing with RAYAMP Rayamp is a ray tracing modelling program that is used to obtain 2D models without inversion stage. To define the velocity structure within a 2D model, there are two types of boudaries, model and divider boundaries. A model boundary is a straight line of arbitrary dip. It has assigned a costant velocity along its length and a nonzero velocity gradient normal to its length. A divider boundary is assigned a velocity zero and it separates two regions
9 3 Ray tracing with RAYAMP 7 Figure 4: Profil from Rijeka to Vojnik with different velocity and velocity gradient. Blocks may thus be defined, in which the velocity, magnitude and direction of velocity gradient are arbitrary. The ray path within a given block is a circular arc (because of constant velocity gradient), for which the travel time and the distance traveled may be calculated using very simple analytical expresions. The source may be located along any model boundary, targets are usually located on the surface of the model. If the incidence to a boundary is at the Figure 5: Profil from Vojnik to Ivanič grad critical angle, the head waves may occur. Beyond the point of intersection of
10 4 Ray tracing with JIVE3D 8 the critical ray with the boundary, the head waves are simulated by shooting critically refracted rays off the boundary at regular intervals along its length. If the incident angle is more than the critical, then the ray reflects from the boundary. For the precritical and multiple reflections, only the boundaries at which reflection is studied need to be specified. At all other boundaries the behaviour of a ray is controlled by the angle of incidence at the boundary. Thus, if no precritical or multiple reflections are studied, then a single specification of the range of take-off angles gives all wide-angle reflections, turning rays and head waves. The corresponding travel-time curve is divided into branches such that the distance along each branch increases or decreases monotonically with distance. The family of rays associated with each traveltime branch is labeled with a unique identification number, which is used in the synthetic seismogram routine for purposes of interpolation within a given ray family. 4 Ray tracing with JIVE3D Modelling with Jive3D can be divided into forward modelling and inversion stjpg. At each iteration of the algorithm, a set of synthetic travel-time data are produced from a working velocity model, and the Frechet derivatives (subsection 4.2.3), which link small changes in model parameters to small changes in travel-time data are calculated. These synthetic data and Frechet derivatives are then passed to the inversion stage, which compares the synthetic data with the given real data and calculates a new model based on a set of linear aproximations until the model converges to a point that optimizes the specified norm for smoothness and best fit. 4.1 Model parametrization An important feature of any tomographic inversion method is its approach to model parametrization. We require a parametrization that is able to describe models in the conventional layer-interface formalism. Each layer and interface is then parametrized by a grid of velocity and depth nodes, allowing the inclusion of reflected and refracted travel times. To define a 3D field of seismic velocity and an interface in a 3D velocity model, we will use qudratic B-splines in three-dimensions B2 3D (x, x ) = β 2 (x 1, x 1)β 2 (x 2, x 2)β 2 (x 3, x 3) and cubic B-splines in two-dimensions B3 2D (x, x ) = β 3 (x 1, x 1)β 3 (x 2, x 2). The β functions are basis functions for interpolation, also called B-spline functions.
11 4.2 The forward problem 9 Velocity field is then given by v(x) = 27 i=1 B 3D 2 (x, x i )v i, (12) where v i are the velocity B-spline coefficients that have influence at position x; x i are the node positions corresponding to those coefficients. Similarly, an interface in a 3D velocity model (modelled as a 2D surface), constructed using cubic B-splines is given by z(x) = 16 i=1 B 2D 3 (x, x i )z i, (13) where z i are depth B-spline coefficients that have influence at position x; x i are equivalent node positions. The interfaces represent discontinuities in seismic velocity at which reflections and refractions may occur. 4.2 The forward problem Two-point ray tracing In the section 2.2 it was stated that the ray equation may be used to determine the trajectory of a ray through a velocity model given the appropriate boundary conditions. These conditons normally take one of the following forms: 1. Source position and direction of propagation specified - find the position at which the ray returns to the surface (if at all) and the travel time; 2. Source position and direction specified - find the position within the model at which a given travel time is specified; 3. Start and end points specified - find all rays which join the two together and their travel times. In order to ensure that all two-point solutions may be found for any possible ray phase (including reflections and refractions through multiple layers), the shooting method is used in JIVE3D. This approach is relatively easy to implement for the models using the ray perturbation theory (RPT). It may be applied in a consistent manner to wide-angle reflection, wide-angle refraction and normal incidence reflection phases. All possible ray paths from source to receiver are being explored.
12 4.2 The forward problem Ray perturbation theory (RPT) RPT is a powerful set of tehniques allowing ray paths, travel times, amplitudes, paraxial rays and even waveforms to be approximated rapidly by applying perturbations to models, for which a previous solution has been calculated. Two mathematical formulations of RPT have emerged: Hamiltonian formulation and Lagrangian formulation. In both, a similar procedure is adopted. A ray is traced through a reference medium cell by cell for which an analytical solution exists and the quantities sought by the procedure (travel time, amplitude, final position and slowness vector) are calculated for reference ray (figure 6). One or more perturbations to the system are then applied and their effect on the ray is estimsted to first order, producing the perturbed ray. When ray-tracing through isotropic media, the Hamiltonian may be defined as follows H(x, p, τ) = 1 2 (p2 u 2 (x)), (14) where u(x) is slowness function, or the reciprocal of the seismic velocity; x and p are the position and the slowness vector perspectively, and the parameter τ is a measure of the extent of propagation along the ray, and may be defined by its differential relation to the travel time dt = u 2 dτ. Hamilton s equations yield a representation of the ray equations, dx dτ = ph, dp dτ = xh. (15) In order to apply RPT, the function u 2 is written in the form u 2 = u 2 ref + u 2 where u 2 ref is a reference function for which rays may be traced analytically and u 2 is perturbation that will be applied. This perturbation of u 2 results in a perturbation of the Hamiltonian H = 1 2 u2. The perturbations to the position and slowness now read x(τ) = x(τ 0 ) + (τ τ 0 ) p(τ 0 ) τ τ 0 (τ τ ) x ( u 2 (x))dτ, (16) p(τ) = p(τ 0 ) + 1 x ( u 2 (x))dτ. (17) 2 τ 0 The travel time as a function of τ is given by T (τ) = τ τ τ 0 u 2 (x)dτ (18) and the perturbation in travel time is Taylor expansion up to first order T (τ) = τ τ 0 u 2 (x(τ ))dτ + τ τ 0 Γ 0 x(τ )dτ, (19)
13 4.2 The forward problem 11 where Γ 0 is the gradient of the slowness squared function (Γ 0 = u 2 (x 0 )). These results are all expressed as polynomials in τ therefore the integrals in (16), (17),and (19) are straightforward. The equations above provide a mechanism for tracing a ray through the medium described in subsection 4.1 using quadratic B-splines Frechet derivatives Once the final position, final slowness and travel time have been determined for a ray within a single cell in the grid, the Frechet derivatives relating the travel time to each of the model parameters defining the velocity within this cell must be obtained (see figure 6). Frechet derivatives for ray solutions to Figure 6: A represantation of cell ray-tracing the two-point problem describe a change in ray path as well as travel time in response to small changes in the velocity or slowness squared parameter (see figure 7). Frechet derivative of the travel time with respect to the j-th model is given by T = s2 u(x) u(x)ds = ds (20) m j m j x 0 (s) s 1 m j where (s 1, s 2 ) is the range along the ray over which the parameter m j influences the travel time. Making use of the identities udτ = ds and d(u 2 ) = 2udu, we obtain T = 1 τ2 u 2 (x) dτ (21) m j 2 τ 1 m j
14 4.2 The forward problem 12 Figure 7: Ray geometry for Frechet derivatives calculated under (a) initial value boundary condition and (b) two-point boundary conditions. Shaded area is part of velocity field that is perturbed. The new ray path is shown as a dotted line. which can be rewritten, by substituting x(τ) = x ref (τ)+ x(τ), as a reference and a perturbation term, where and ( ) Tref m j ( T m j = 1 2 ) ( ) ( ) T Tref T = +, (22) m j m j m j τref τ 0 = 1 2 ( u 2 (x 0 ) + Γ ) 0 (x ref x 0 ) dτ (23) m j m j τfinal τ 0 ( ) Γ0 x + u2 (x 0 ) dτ. (24) m j m j The 3D quadratic B-spline used to produce the slowness squared field in a cell is u 2 (x) = Since Γ 0 = u 2 (x 0 ), we have 27 j=1 B 3D 2 (x, x j )m j. (25) u 2 (x 0 ) m j = B 3D 2 (x, x j ), Γ 0 m j = B 3D 2 (x, x j ). (26)
15 4.3 The inverse problem 13 So far, these expressions have been derived for the case, in which the model parameters describing the seismic velocity field contain B-spline coefficients in units of u 2 (represented as m u 2). When solving the inverse problem using regularized inversion (subsection 4.3), the model parameters must be converted to units of velocity using the relation m v = 1/ m u 2 and the Frechet derivatives T = 2 T m v (m v ) 2 m u 2 (27) Beginning with the differential expression for the change in travel time δt = pδx along a short segment of a ray δx, for which the slowness vector is p, a simple expression for the change in travel time for a ray crossing an interface due to perturbation in that interface may be obtained from Snell law; δt = δp 3 δz, where δp 3 is the change in the vertical component of the slowness vector at the reflection/refraction and δz is the change in the interface depth. The interface is defined by equation (13), so an expression for Frechet derivative may be obtained as T z = δp 3 = δp 3 B3 2D (x, x j ) (28) m j m j The results derived above allow a ray to be traced through a single cell in a model parametrized as described in subsection 4.1, producing final position and slowness vectors, travel times, and Frechet derivatives. 4.3 The inverse problem In the linearized inversion approach, the inverse problem may be formulated as: we have a model, described by the vector m, the components of which are the interface depths and velocity parameters included in the inversion. From this model, the forward modelling code provides a set of synthetic travel time data, contained in the vector t and a matrix A of Frechet derivatives, which measure the sensitivity of the model and synthetic travel times. We also have the set of real travel-time data, t real, the uncertainties in those data, σ, and information on the geometric arrangement of model parameters. In order to measure the fit of the current model to the data, we define the traveltime residual vector r = t real t, which is a subject of the optimization by least-squares formulation. For example I would like to present the evolution of salt dome inversion from the starting to the final model. We can notice from pictures bellow that χ 2 is decreasing step by step, so we are getting more and more realistic two dimensional model of investigating area.
16 5 Conclusion 14 5 Conclusion Figure 8: From starting to the final model We conclude by comparing both modelling techiques. Rayamp is only used for 2D forward modelling, which means, that it can not produce inversions to the starting model. If more realistic models are to be obtained with Rayamp, one would have to include some sort of alghoritm for automatic search of the best set of parameters, which are defining the model, such as; different model boundaries, velocity field and smoothing criteria. The procedure of Jive3D is more practical, it also includes inversion and is therefore able to constuct the final model as an evolution of the starting model.
17 REFERENCES 15 References [1] Agustin Udias: Principles of Seismology (Cambridge University Press 1999); [2] James William Douglas Hobro: Three-dimensional tomographic inversion of combined reflection and refraction seismic travel-time data (Chuchill College Cambridge 1999); [3] Andrej Gosar: Exploration of the South-Eastern Alps lithosphere with 3D refraction seismic (Project Alp2002) data acqusition in Slovenia (Geologija); [4] G. D. Spence, K. P. Whittall, R. M. Clowes: Practical Synthetic Seismograms for latterally varying media calculated by asymptotic ray theory (BSSA 1984; );
dp i = u ray c2 dt with t the travel time.
Three dimensional two point ray tracing using paraxial rays in cartessian coordinates Ronit Strahilevitz*, Paradigm Geophysical; Dan Kosloff, Department of Geophysics, Tel Aviv University, Israel and Paradigm
More informationGEOPHYS 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 informationVidale, 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 informationNon-linear surface wave phase velocity inversion based on ray theory
Geophys. J. Int. (22) 151, 583 596 Non-linear surface wave phase velocity inversion based on ray theory Marianne Bruneton, 1 Véronique Farra, 2 Helle Anette Pedersen 1 and the SVEKALAPKO Seismic Tomography
More informationThe geometry of reflection and refraction Wave conversion and reflection coefficient
7.2.2 Reflection and Refraction The geometry of reflection and refraction Wave conversion and reflection coefficient The geometry of reflection and refraction A wave incident on a boundary separating two
More informationGG450 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 informationComparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media
Raytracing in VTI media Comparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media Jianli Yang and Don C. Lawton ABSTRACT Determination of the conversion
More informationIntroduction to seismic body waves Tomography
Introduction to seismic body waves Tomography Presented by Han Zhang Prepared by group work of Justin Wilgus and Han Zhang Outlines 1. Velocity derived from wave equation 2. Two commonly used reference
More informationAngle-gather time migration a
Angle-gather time migration a a Published in SEP report, 1, 141-15 (1999) Sergey Fomel and Marie Prucha 1 ABSTRACT Angle-gather migration creates seismic images for different reflection angles at the reflector.
More informationComparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media
Comparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media Jianli Yang* and Don C. Lawton CREWES, University of Calgary, 5 UniversityDR.NW, Calgary, AB, TN
More informationProcessing 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 informationP063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA
1 P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA I. TRINKS 1, S. SINGH 2, C. CHAPMAN 3, P. BARTON 1, M. BOSCH 4, A. CHERRETT 5 Addresses 1 Bullard Laboratories,
More informationEffects 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 informationEARS2232 Exploration Seismics. Exploration. 2006/2007 Semester 1. Teaching Team. Dr. Sebastian Rost (Module Leader)
EARS3 Exploration Seismics Teaching Team Dr. Sebastian Rost (Module Leader) Dr. Graham Stuart (Seismic Interpretation) 006/007 Semester 1 Ben Dando (Demonstrator) Objectives On completion of this module
More informationInterpolation within ray tubes - state of the art
Interpolation within ray tubes - state of the art Petr Bulant Department of Geophysics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Praha 2, Czech Republic, http://sw3d.cz/staff/bulant.htm
More informationChapter 2 Methodology of Seismic Tomography
Chapter Methodology of Seismic Tomography Abstract In this chapter we introduce the tomographic methods which are widely used to study three-dimensional (3-D) seismic velocity attenuation and anisotropy
More informationAnalytical solutions for the ray perturbation in depth-varying media
Geophys. J. Inr. (1993) 114,281-288 Analytical solutions for the ray perturbation in depth-varying media Beverley J. Moore School of Marhematical Sciences, Uniuersiry of Technology, Sydney, PO Box 123,
More informationGEOPHYSICS. Time to depth conversion and seismic velocity estimation using time migration velocities
Time to depth conversion and seismic velocity estimation using time migration velocities Journal: Manuscript ID: Manuscript Type: Date Submitted by the Author: Geophysics draft Letters n/a Complete List
More informationAngle Gathers for Gaussian Beam Depth Migration
Angle Gathers for Gaussian Beam Depth Migration Samuel Gray* Veritas DGC Inc, Calgary, Alberta, Canada Sam Gray@veritasdgc.com Abstract Summary Migrated common-image-gathers (CIG s) are of central importance
More informationA comparison of exact, approximate, and linearized ray tracing methods in transversely isotropic media
Ray Tracing in TI media A comparison of exact, approximate, and linearized ray tracing methods in transversely isotropic media Patrick F. Daley, Edward S. Krebes, and Laurence R. Lines ABSTRACT The exact
More informationRay methods in rough models
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 86 Ray methods in rough models Thorbjørn Rekdal and Biondo Biondi 1 ABSTRACT First arrival traveltime methods have proven to fail to image
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationSuccessful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment.
Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment. Sergey Birdus 1, Dean Criddle 2, Alexey Artyomov 1, Li Li 1, Qingbing Tang 1,
More informationSeismic 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 informationGeometric Acoustics in High-Speed Boundary Layers
Accepted for presentation at the 9th International Symposium on Shock Waves. Madison, WI. July -9,. Paper #8 Geometric Acoustics in High-Speed Boundary Layers N. J. Parziale, J. E. Shepherd, and H. G.
More informationTravel Time Tomography using Neural Networks
Travel Time Tomography using Neural Networks Yoshiya Oda Tokyo Metropolitan University, Japan Tomohisa Ishiyama Tokyo Metropolitan University, Japan Shinya Yokono Tokyo Metropolitan University, Japan SUMMARY:
More informationFermat s principle and ray tracing in anisotropic layered media
Ray tracing in anisotropic layers Fermat s principle and ray tracing in anisotropic layered media Joe Wong ABSTRACT I consider the path followed by a seismic signal travelling through velocity models consisting
More informationPEAT SEISMOLOGY Lecture 16: Seismic Tomography I
PEAT8002 - SEISMOLOGY Lecture 16: Seismic Tomography I Nick Rawlinson Research School of Earth Sciences Australian National University Introduction Seismic data represent one of the most valuable resources
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationA new method for determination of a wave-ray trace based on tsunami isochrones
Bull. Nov. Comp. Center, Math. Model. in Geoph., 13 (2010), 93 101 c 2010 NCC Publisher A new method for determination of a wave-ray trace based on tsunami isochrones An.G. Marchuk Abstract. A new method
More informationGeometrical Optics INTRODUCTION. Wave Fronts and Rays
Geometrical Optics INTRODUCTION In this experiment, the optical characteristics of mirrors, lenses, and prisms will be studied based on using the following physics definitions and relationships plus simple
More informationProgress 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(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
More informationFinite-difference elastic modelling below a structured free surface
FD modelling below a structured free surface Finite-difference elastic modelling below a structured free surface Peter M. Manning ABSTRACT This paper shows experiments using a unique method of implementing
More informationSUMMARY. 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 informationAmplitude within the Fresnel zone for the zero-offset case
Amplitude within the Fresnel zone for the zero-offset case Shuang Sun and John C. Bancroft ABSTRACT Reflection energy from a linear reflector comes from an integrant over an aperture often described by
More informationKinematics of Reverse Time S-G Migration
Kinematics of Reverse Time S-G Migration William W. Symes TRIP Seminar Rice University September 2003 www.trip.caam.rice.edu 1 Agenda: explore prestack focussing properties of RT S-G migration, proper
More informationThe Effect of a Dipping Layer on a Reflected Seismic Signal.
The Effect of a Dipping Layer on a Reflected Seismic Signal. Edward O. Osagie, Ph.D. Lane College, Jackson, TN 38305 E-mail: eosagiee@yahoo.com ABSTRACT Modern seismic reflection prospecting is based on
More informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More informationThe Interface Fresnel Zone revisited
The Interface Fresnel Zone revisited Nathalie Favretto-Cristini, Paul Cristini and Eric de Bazelaire ABSTRACT We determine the part of reflectors which actually affects the reflected wavefield, which is
More informationHigh Frequency Wave Scattering
High Frequency Wave Scattering University of Reading March 21st, 2006 - Scattering theory What is the effect of obstacles or inhomogeneities on an incident wave? - Scattering theory What is the effect
More informationCoupled Wave Field Migration. Coupled Wave Field Migration
Coupled Wave Field Migration Coupled Wave Field Migration Vp=1650 m/s T=750 ms Vp=4000 m/s Surface seismic exploration uses reflected and refracted waves for imaging purposes. In this case, special conditions
More informationAPPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA.
APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. BY: ISAAC KUMA YEBOAH. Department of Engineering, Regent University
More informationP052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges
P052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges M. Ayzenberg* (Norwegian University of Science and Technology, A. Aizenberg (Institute of Petroleum Geology and Geophysics,
More informationA A14 A24 A34 A44 A45 A46
Computation of qp- and qs-wave Rays raveltimes Slowness Vector and Polarization in General Anisotropic Media Junwei Huang 1 * Juan M. Reyes-Montes 1 and R. Paul Young 2 1 Applied Seismology Consultants
More informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
More informationSEG/San Antonio 2007 Annual Meeting
Yaofeng He* and Ru-Shan Wu Modeling and Imaging Laboratory, Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA, 95064, USA Summary A one-way and one-return boundary
More informationSeismic modelling with the reflectivity method
Yongwang Ma, Luiz Loures, and Gary F. Margrave ABSTRACT Seismic modelling with the reflectivity method Numerical seismic modelling is a powerful tool in seismic imaging, interpretation and inversion. Wave
More informationModeling of diffractive optical elements for lens design.
Juan L. Rayces and Lan Lebich OA Applied Optics, 7421 Orangewood Ave., Garden Grove, A 92641. Abstract. The use of a standard aspheric profile to describe conventional optical elements in lens design programs
More informationPressure and shear wave separation of ocean bottom seismic data
Pressure and shear wave separation of ocean bottom seismic data Ohad Barak ABSTRACT PZ summation is a common method for separating the upgoing wavefield from the downgoing wavefield in data acquired by
More information2.710 Optics Spring 09 Solutions to Problem Set #1 Posted Wednesday, Feb. 18, 2009
MASSACHUSETTS INSTITUTE OF TECHNOLOGY.70 Optics Spring 09 Solutions to Problem Set # Posted Wednesday, Feb. 8, 009 Problem : Spherical waves and energy conservation In class we mentioned that the radiation
More informationFull waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)
of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many
More informationA Wavelet-Based Method for Simulation of Seismic Wave Propagation
AGU Fall Meeting, A Wavelet-Based Method for Simulation of Seismic Wave Propagation Tae-Kyung Hong & B.L.N. Kennett Research School of Earth Sciences The Australian National University Abstract Seismic
More informationStratigraphic coordinates, a coordinate system tailored to seismic interpretation a
Stratigraphic coordinates, a coordinate system tailored to seismic interpretation a a Published in Geophysical Prospecting, v. 63, 1246-1255, (2015) Parvaneh Karimi and Sergey Fomel ABSTRACT In certain
More information96 Alkhalifah & Biondi
Stanford Exploration Project, Report 97, July 8, 1998, pages 95 116 The azimuth moveout operator for vertically inhomogeneous media Tariq Alkhalifah and Biondo L. Biondi 1 keywords: AMO, DMO, dip moveout
More informationH003 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 informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationReflection seismic Method - 2D
Reflection seismic Method - 2D Acoustic Impedance Seismic events Wavelets Convolutional model Resolution Stacking and Signal/Noise Data orders Reading: Sheriff and Geldart, Chapters 6, 8 Acoustic Impedance
More information3-D traveltime computation using Huygens wavefront tracing
GEOPHYSICS, VOL. 66, NO. 3 MAY-JUNE 2001); P. 883 889, 12 FIGS., 1 TABLE. 3-D traveltime computation using Huygens wavefront tracing Paul Sava and Sergey Fomel ABSTRACT Traveltime computation is widely
More informationReflectivity modeling for stratified anelastic media
Reflectivity modeling for stratified anelastic media Peng Cheng and Gary F. Margrave Reflectivity modeling ABSTRACT The reflectivity method is widely used for the computation of synthetic seismograms for
More informationP H Y L A B 1 : G E O M E T R I C O P T I C S
P H Y 1 4 3 L A B 1 : G E O M E T R I C O P T I C S Introduction Optics is the study of the way light interacts with other objects. This behavior can be extremely complicated. However, if the objects in
More informationIMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF?
IMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF? Summary Samuel H. Gray* Veritas DGC Inc., 715 Fifth Ave. SW, Suite 2200, Calgary, AB Sam_gray@veritasdgc.com Carl Notfors Veritas
More informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
More informationu = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially.
Seismic Elastic Modelling Satinder Chopra, Arcis Corporation, Calgary, Canada Summary Seismic modeling experiments were performed to understand the seismic responses associated with a variety of subsurface
More informationGraphics and Interaction Transformation geometry and homogeneous coordinates
433-324 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationForce-Control for the Automated Footwear Testing System
Chapter 2 Force-Control for the Automated Footwear Testing System Problem presented by: Dr. Gerald K. Cole (Biomechanigg Research Inc.) Mentors: Greg Lewis (University of Ontario IT), Rex Westbrook (University
More informationTRAVELTIME TOMOGRAPHY (CONT)
30 March 005 MODE UNIQUENESS The forward model TRAVETIME TOMOGRAPHY (CONT di = A ik m k d = Am (1 Data vecto r Sensitivit y matrix Model vector states the linearized relationship between data and model
More informationCOMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates
COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationMathematics of Multidimensional Seismic Imaging, Migration, and Inversion
N. Bleistein J.K. Cohen J.W. Stockwell, Jr. Mathematics of Multidimensional Seismic Imaging, Migration, and Inversion With 71 Illustrations Springer Contents Preface List of Figures vii xxiii 1 Multidimensional
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More informationNUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS
- 465 - NUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS LIU, Z. Q. SUN, J. G. * SUN, H. * LIU, M. C. * GAO, Z. H. College for Geoexploration Science
More informationFast marching methods
1 Fast marching methods Lecture 3 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Metric discretization 2 Approach I:
More informationGeOlympus. 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 informationA NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE
Letter J. Phys. Earth, 34, 195-201, 1986 A NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE Ichiro NAKANISHI and Kaoru YAMAGUCHI
More informationSpectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications
Spectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications Lin Zheng*, Qi Zhao*, Qinya Liu, Bernd Milkereit, Giovanni Grasselli, University of Toronto, ON, Canada
More informationAn imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults
Earth Planets Space, 52, 3 11, 2000 An imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults Takumi Murakoshi 1, Hiroshi Takenaka
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationSeparation of diffracted waves in TI media
CWP-829 Separation of diffracted waves in TI media Yogesh Arora & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines Key words: Diffractions, Kirchhoff, Specularity, Anisotropy ABSTRACT
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationSeismic traveltime inversion for 2-D crustal velocity structure
Geophys. J. Int. (1992) 108, 16-34 Seismic traveltime inversion for 2-D crustal velocity structure C. A. Zelt* and R. B. Smith Department of Geology and Geophysics, University of Utah, Salt Lake City,
More informatione. \ 2-99 Calculation ofray Amplitudes: We seek a computationally tractable method to estimate ray
2-99 e. \ Michael C. Fehler, Steven T. Hildebrand, Lianjie Huang, and Douglas Alde MS 044; Los Alamos National Laboratiory; Los Alamos, NM 87545 USA Summary: Kirchhoff migration using ray tracing travel
More information3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines
3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose 3D angle decomposition methods from elastic reverse
More informationPHYS 202 Notes, Week 8
PHYS 202 Notes, Week 8 Greg Christian March 8 & 10, 2016 Last updated: 03/10/2016 at 12:30:44 This week we learn about electromagnetic waves and optics. Electromagnetic Waves So far, we ve learned about
More informationTravel Time Inversion in Seismic Tomography
Travel Time Inversion in Seismic Tomography S. Padina, D. Churchill, R. P. Bording March 2, 2006 1 Introduction The word tomography finds its origins in the two Greek words tomos, meaning section, and
More informationCHAPTER 6 Parametric Spline Curves
CHAPTER 6 Parametric Spline Curves When we introduced splines in Chapter 1 we focused on spline curves, or more precisely, vector valued spline functions. In Chapters 2 and 4 we then established the basic
More informationAbstract. Introduction. Eleventh International Congress of The Brazilian Geophysical Society
SS-traveltime parameters from PP- and PS-reflections Bjørn Ursin, Martin Tygel, and Einar Iversen Dept. of Petroleum Eng. and Appl. Geophysics, NTNU, S.P. Andersensvei 15A, NO-7491 Trondheim, Norway Dept.
More informationUT inspection of welding joins: issues and contribution of CIVA simulation tools
UT inspection of welding joins: issues and contribution of CIVA simulation tools GENERAL CONTEXT: ISSUES AND PROBLEMS OF WELD INSPECTION BY ULTRASONIC METHODS 67 th IIW Annual Assembly and International
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationAnisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS
Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Summary Anisotropic depth model building using surface seismic data alone is non-unique and
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationRotation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J. Fowler, WesternGeco.
otation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J Fowler, WesternGeco Summary Symmetry groups commonly used to describe seismic anisotropy
More informationA least-squares shot-profile application of time-lapse inverse scattering theory
A least-squares shot-profile application of time-lapse inverse scattering theory Mostafa Naghizadeh and Kris Innanen ABSTRACT The time-lapse imaging problem is addressed using least-squares shot-profile
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationDetermination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data
Determination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data Mohamed Amrouche and Hiroaki Yamanaka Dept. of Environmental Science and Technology, Tokyo Institute
More informationUNIVERSITY OF CALGARY. Traveltime Tomography in Isotropic and Transversely Isotropic Media. Marco Perez
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 informationReverse 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 informationSpline Curves. Spline Curves. Prof. Dr. Hans Hagen Algorithmic Geometry WS 2013/2014 1
Spline Curves Prof. Dr. Hans Hagen Algorithmic Geometry WS 2013/2014 1 Problem: In the previous chapter, we have seen that interpolating polynomials, especially those of high degree, tend to produce strong
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationLecture Ray Model of Light. Physics Help Q&A: tutor.leiacademy.org
Lecture 1201 Ray Model of Light Physics Help Q&A: tutor.leiacademy.org Reflection of Light A ray of light, the incident ray, travels in a medium. When it encounters a boundary with a second medium, part
More informationMODELING OF THREE-DIMENSIONAL PROPAGATION ON A COASTAL WEDGE WITH A SEDIMENT SUPPORTING SHEAR
MODELING OF THREE-DIMENSIONAL PROPAGATION ON A COASTAL WEDGE WITH A SEDIMENT SUPPORTING SHEAR Piotr Borejko Vienna University of Technology, Karlsplatz 13/E26/3, 14 Vienna, Austria Fax: + 43 1 588 1 21
More information