Multipathing, reciprocal traveltime fields and raylets

N. Rawlinson, M. Sambridge, J. Hauser

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

A new theory for the calculation of later seismic arrivals in heterogeneous media is presented. We introduce the concept of a ‘raylet’, which is a segment of a later arriving ray path between source and receiver defined by joint properties of forward and reciprocal traveltime fields. We show that all rays between a single source and receiver in arbitrary heterogeneous media can be divided into a unique set of raylets, any one of which can be used to construct the complete two-point path. A particularly useful property of raylets is that they correspond to stationary curves in the summed (forward and reciprocal) traveltime fields, i.e. adding the first (or second, or even later) arriving traveltime field from some point A to the first (or second, or even later) arriving traveltime field from some point B. We show that many raylets, each corresponding to a later arriving phase, require only earlier arrival traveltime fields for their construction. The theory describing the properties of raylets is the primary result of the paper. One practical consequence is that many later arrivals between source and receiver in heterogeneous media can be found from just two first-arrival traveltime fields, one from the source and the other (the reciprocal field) from the receiver.

The theory and a method for constructing later arrivals is demonstrated though numerical experiments in 2-D but holds without change in 3-D. We use a simple grid-based eikonal solver to compute forward and reciprocal first-arrival traveltime fields, and validate our results with a ray-based wave front construction (WFC) technique. In one test involving a simple wave front triplication (or swallowtail), all three arrivals are retrieved. In another example featuring severe velocity heterogeneities, 16 out of a total of 37 arrivals are found. The theory shows that combining second and higher order traveltime fields from source and receiver yields all raylets and hence all later arrivals. In practice only first-arrival traveltime fields are usually available and for this case we describe a simple procedure to find any remaining arrivals using intermediate artificial sources together with their first-arrival traveltime fields. The properties of raylets and their manifestation in the joint traveltime field appears to be a previously unrecognized feature and provides a novel new approach to the calculation of later arrivals.
Original languageEnglish
Pages (from-to)1077-1092
Number of pages16
JournalGeophysical Journal International
Volume181
Issue number2
Early online date19 Mar 2010
DOIs
Publication statusPublished - May 2010

Fingerprint

arrivals
heterogeneous medium
receivers
rays
Experiments
wave fronts
experiment
calculation
grids

Keywords

  • numerical solutions
  • computational seismology
  • theoretical seismology
  • wave propagation

Cite this

Multipathing, reciprocal traveltime fields and raylets. / Rawlinson, N.; Sambridge, M.; Hauser, J.

In: Geophysical Journal International, Vol. 181, No. 2, 05.2010, p. 1077-1092.

Research output: Contribution to journalArticle

Rawlinson, N. ; Sambridge, M. ; Hauser, J. / Multipathing, reciprocal traveltime fields and raylets. In: Geophysical Journal International. 2010 ; Vol. 181, No. 2. pp. 1077-1092.
@article{eed08ec348c3447886a490b71d3bc2c1,
title = "Multipathing, reciprocal traveltime fields and raylets",
abstract = "A new theory for the calculation of later seismic arrivals in heterogeneous media is presented. We introduce the concept of a ‘raylet’, which is a segment of a later arriving ray path between source and receiver defined by joint properties of forward and reciprocal traveltime fields. We show that all rays between a single source and receiver in arbitrary heterogeneous media can be divided into a unique set of raylets, any one of which can be used to construct the complete two-point path. A particularly useful property of raylets is that they correspond to stationary curves in the summed (forward and reciprocal) traveltime fields, i.e. adding the first (or second, or even later) arriving traveltime field from some point A to the first (or second, or even later) arriving traveltime field from some point B. We show that many raylets, each corresponding to a later arriving phase, require only earlier arrival traveltime fields for their construction. The theory describing the properties of raylets is the primary result of the paper. One practical consequence is that many later arrivals between source and receiver in heterogeneous media can be found from just two first-arrival traveltime fields, one from the source and the other (the reciprocal field) from the receiver.The theory and a method for constructing later arrivals is demonstrated though numerical experiments in 2-D but holds without change in 3-D. We use a simple grid-based eikonal solver to compute forward and reciprocal first-arrival traveltime fields, and validate our results with a ray-based wave front construction (WFC) technique. In one test involving a simple wave front triplication (or swallowtail), all three arrivals are retrieved. In another example featuring severe velocity heterogeneities, 16 out of a total of 37 arrivals are found. The theory shows that combining second and higher order traveltime fields from source and receiver yields all raylets and hence all later arrivals. In practice only first-arrival traveltime fields are usually available and for this case we describe a simple procedure to find any remaining arrivals using intermediate artificial sources together with their first-arrival traveltime fields. The properties of raylets and their manifestation in the joint traveltime field appears to be a previously unrecognized feature and provides a novel new approach to the calculation of later arrivals.",
keywords = "numerical solutions, computational seismology, theoretical seismology, wave propagation",
author = "N. Rawlinson and M. Sambridge and J. Hauser",
year = "2010",
month = "5",
doi = "10.1111/j.1365-246X.2010.04558.x",
language = "English",
volume = "181",
pages = "1077--1092",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "2",

}

TY - JOUR

T1 - Multipathing, reciprocal traveltime fields and raylets

AU - Rawlinson, N.

AU - Sambridge, M.

AU - Hauser, J.

PY - 2010/5

Y1 - 2010/5

N2 - A new theory for the calculation of later seismic arrivals in heterogeneous media is presented. We introduce the concept of a ‘raylet’, which is a segment of a later arriving ray path between source and receiver defined by joint properties of forward and reciprocal traveltime fields. We show that all rays between a single source and receiver in arbitrary heterogeneous media can be divided into a unique set of raylets, any one of which can be used to construct the complete two-point path. A particularly useful property of raylets is that they correspond to stationary curves in the summed (forward and reciprocal) traveltime fields, i.e. adding the first (or second, or even later) arriving traveltime field from some point A to the first (or second, or even later) arriving traveltime field from some point B. We show that many raylets, each corresponding to a later arriving phase, require only earlier arrival traveltime fields for their construction. The theory describing the properties of raylets is the primary result of the paper. One practical consequence is that many later arrivals between source and receiver in heterogeneous media can be found from just two first-arrival traveltime fields, one from the source and the other (the reciprocal field) from the receiver.The theory and a method for constructing later arrivals is demonstrated though numerical experiments in 2-D but holds without change in 3-D. We use a simple grid-based eikonal solver to compute forward and reciprocal first-arrival traveltime fields, and validate our results with a ray-based wave front construction (WFC) technique. In one test involving a simple wave front triplication (or swallowtail), all three arrivals are retrieved. In another example featuring severe velocity heterogeneities, 16 out of a total of 37 arrivals are found. The theory shows that combining second and higher order traveltime fields from source and receiver yields all raylets and hence all later arrivals. In practice only first-arrival traveltime fields are usually available and for this case we describe a simple procedure to find any remaining arrivals using intermediate artificial sources together with their first-arrival traveltime fields. The properties of raylets and their manifestation in the joint traveltime field appears to be a previously unrecognized feature and provides a novel new approach to the calculation of later arrivals.

AB - A new theory for the calculation of later seismic arrivals in heterogeneous media is presented. We introduce the concept of a ‘raylet’, which is a segment of a later arriving ray path between source and receiver defined by joint properties of forward and reciprocal traveltime fields. We show that all rays between a single source and receiver in arbitrary heterogeneous media can be divided into a unique set of raylets, any one of which can be used to construct the complete two-point path. A particularly useful property of raylets is that they correspond to stationary curves in the summed (forward and reciprocal) traveltime fields, i.e. adding the first (or second, or even later) arriving traveltime field from some point A to the first (or second, or even later) arriving traveltime field from some point B. We show that many raylets, each corresponding to a later arriving phase, require only earlier arrival traveltime fields for their construction. The theory describing the properties of raylets is the primary result of the paper. One practical consequence is that many later arrivals between source and receiver in heterogeneous media can be found from just two first-arrival traveltime fields, one from the source and the other (the reciprocal field) from the receiver.The theory and a method for constructing later arrivals is demonstrated though numerical experiments in 2-D but holds without change in 3-D. We use a simple grid-based eikonal solver to compute forward and reciprocal first-arrival traveltime fields, and validate our results with a ray-based wave front construction (WFC) technique. In one test involving a simple wave front triplication (or swallowtail), all three arrivals are retrieved. In another example featuring severe velocity heterogeneities, 16 out of a total of 37 arrivals are found. The theory shows that combining second and higher order traveltime fields from source and receiver yields all raylets and hence all later arrivals. In practice only first-arrival traveltime fields are usually available and for this case we describe a simple procedure to find any remaining arrivals using intermediate artificial sources together with their first-arrival traveltime fields. The properties of raylets and their manifestation in the joint traveltime field appears to be a previously unrecognized feature and provides a novel new approach to the calculation of later arrivals.

KW - numerical solutions

KW - computational seismology

KW - theoretical seismology

KW - wave propagation

U2 - 10.1111/j.1365-246X.2010.04558.x

DO - 10.1111/j.1365-246X.2010.04558.x

M3 - Article

VL - 181

SP - 1077

EP - 1092

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 2

ER -