PestoSeis’s documentation

The PestoSeis package is a set of functions to perform calculations of seismic traveltimes, rays and wave propagation in two dimensions (2D). PestoSeis covers the range from forward modeling wave propagation through a known model to reconstructing the structural map of an unknown medium from some observed (synthetic) data by providing algorithms to solve a 2D inverse problem by simple means.

Its main purpose is educational, it aims to be a simple tool to play around with some aspects of seismology. Due to the fact that many fields of research use the fundamentals of wave propagation, the functionalities within PestoSeis are not limited to seismology, but can also be applied for instance to problems in medical ultrasound (check out the tutorial on medical ultrasound). Due to the educational purpose of the code, no optimization for speed has been perfomed and it has not been tested extensively. Therefore, use of this code for research it is not encouraged.

The functions contained in this package aim at the following:

• Setting up a grid containing a velocity model;

• Computing traveltimes for given sources and receivers, optionally providing the entire 2D traveltime model;

• Tracing seismic rays given a 2D traveltime array and position of sources and receivers;

• Performing very simple traveltime tomography using rays;

• Computing acoustic and elastic wave propagation in 2D by finite differences;

• Performing basic seismic processing for exploration seismology.

Authors:

• Andrea Zunino

• Patrick Marty

• Ines E. Ulrich

Installation

See the Installation — loading the module for more details.

User guide

The User guide provides information about functionality and how to use the package with examples.

Traveltimes and rays in 2D

See Traveltimes and rays – using ttimerays.

Wave propagation in 2D

See Seismic wave propagation – using seismicwaves2d.

Tutorials

A set of examples in the form of Jupyter notebooks illustrating the capabilities of the code. See Tutorials.

For developers

See Developer zone.

Table of contents

• Installation — loading the module
• Traveltimes and rays – using ttimerays
  • Setup of grid parameters and models
  • Visualization
  • Computation of traveltimes
  • Rays
    • Trace rays in a 2D heterogeneous model
    • Trace straight rays
    • Trace rays in a horizontally layered medium
    • Ray tomography
• Seismic wave propagation – using seismicwaves2d
  • Acoustic wave equation in 2D
  • Elastic wave equation in 2D
  • Boundary conditions
• Tutorials
  • Tutorial 01 - Earthquake location
    • 1. Problem Setup
    • 2. Create list of trial source positions
    • 3. Calculate all traveltimes for all sources
    • 4. Traveltime Difference
    • 5. Plotting the Results
  • Tutorial 02 - Trace rays
    • 1. Problem Setup
    • 2. Compute the Traveltimes
    • 3. Trace the Rays
    • 4. Plot the Results
  • Tutorial 03 - Velocity inversion
    • 1. Problem Setup
    • 2. Trace Ray Paths
    • 3. Build the Tomography Matrix
    • 4. Perform the Least-Squares Inversion
    • 5. Plotting the Results
  • Tutorial 04 - Medical Ultrasound - Calculate travel times using the USCT setup
    • Problem Setup
    • Compute the Traveltimes
    • Trace the Rays
    • Plot the Results
    • References
  • Tutorial 05 - Acoustic Wave Propagation
    • Acoustic Forward Simulation
    • Preamble
    • Model Interpolation
    • Simulation Setup
    • Running the Simulation
    • Plotting the Results
    • References
  • Tutorial 06 - Elastic Wave Propagation
    • Elastic Forward Simulation
    • Preamble
    • Model Interpolation
    • Simulation Setup
    • Running the Simulation
    • Plotting the Results
    • References
  • Tutorial 07 - Processing seismic data I
    • Acoustic waves simulation
    • Create a movie
  • Tutorial 08 - Processing seismic data II
    • f-k Filtering
    • Preamble
    • Import Data
    • Transform to Frequency-Wavenumber Domain
    • Remove the Ground Roll
    • Convert back to Time-Space Domain
• API reference
  • Reflection seismology - reflectionseismo
    • agc()
    • calcreflectivity()
    • fwdconvolve()
    • geometrical_spreading()
    • imgshotgath()
    • nmocorrection()
    • wiggle()
  • Seismic waves - seismicwaves2d
    • animateacousticwaves()
    • animateelasticwaves()
    • gaussource()
    • rickersource()
    • solveacoustic2D()
    • solveelastic2D()
  • Traveltimes and rays - ttimerays
    • buildtomomat()
    • lininv()
    • plotgrid()
    • plotrays()
    • plotttimemod()
    • plotvelmod()
    • rollmod()
    • setupgrid()
    • traceall_straight_rays()
    • traceallrays()
    • tracerayhorlay()
    • traveltime()
    • unrollmod()
• Developer zone
  • All members including private ones

Indices and tables

• Index

• Module Index

• Search Page