On forward and inverse modelling in seismology: raytracing in inhomogeneous media

Smith, Peter M. (2006) On forward and inverse modelling in seismology: raytracing in inhomogeneous media. Masters thesis, Memorial University of Newfoundland.

[img] [English] PDF (Migrated (PDF/A Conversion) from original format: (application/pdf)) - Accepted Version
Available under License - The author retains copyright ownership and moral rights in this thesis. Neither the thesis nor substantial extracts from it may be printed or otherwise reproduced without the author's permission.

Download (14Mb)
  • [img] [English] PDF - Accepted Version
    Available under License - The author retains copyright ownership and moral rights in this thesis. Neither the thesis nor substantial extracts from it may be printed or otherwise reproduced without the author's permission.
    (Original Version)

Abstract

The first part of this thesis deals with forward modelling. We present a raytracing method based on the concept of simulated annealing: a computational tool based on physical principles used for obtaining optimal solutions of problems of in areas ranging from combinatorics to condensed matter physics. Our method solves for rays that render signal traveltime stationary, in accordance with Fermat's principle of stationary traveltime. We test this method for two types of media: layered inhomogeneous media and linearly inhomogeneous media. We show that rays and traveltimes generated from this algorithm for these models quantitatively agree with predicted results. -- The second part of the thesis deals with inverse modelling. In this part, we introduce the generalized form of Radon's transform and its adjoint operator. We show that by treating traveltime as Radon's transform acting on the slowness function along a ray, we can use the adjoint operator to recover qualitative information about a medium from collected traveltimes. This method of backprojection is presented as an application of our raytracing method. We calculate rays and their associated traveltimes between sources and receivers on a square lattice for layered- and linearly-inhomogeneous media and use the back-projection method to construct slowness functions for each set of data. We show that although the backprojection method does not retain the quantitative properties of the original medium, results indicate that qualitative properties of the medium can be resolved by this method.

Item Type: Thesis (Masters)
URI: http://research.library.mun.ca/id/eprint/6669
Item ID: 6669
Additional Information: Bibliography: leaves 115-116. Link to CD-ROM materials included in the sidebar, at the end of the thesis. Pages [vi], [xii], [xvi], [xviii], [122] and [136] are blank and have been omitted from the digital reproduction.
Department(s): Science, Faculty of > Earth Sciences
Date: 2006
Date Type: Submission
Library of Congress Subject Heading: Inhomogeneous materials; Simulated annealing (Mathematics); Seismic waves--Mathematical models

Actions (login required)

View Item View Item

Downloads

Downloads per month over the past year

View more statistics