Advanced

Broad band synthetic seismograms for a spherical inhomogeneity in a many-layered elastic half-space

Karlsson, Anders LU and Boström, Anders (1987) In Geophysical Journal International 89(2). p.527-547
Abstract
The propagation of elastic waves in a many-layered elastic half-space is considered. One of the layers contains a bounded inhomogeneity which in the numerical applications is taken as a sphere. The half-space is excited by an explosion, modelled as a sudden isotropic point source.



First the time-harmonic problem is solved using the null field approach (the T-matrix method). Starting from surface integral representations containing the free space Green's dyadic and field expansions in plane and spherical vector wave functions, the null field approach leads to a set of algebraic equations whose solution can be given in a form that naturally has a multiple scattering interpretation. The null field approach thus has a... (More)
The propagation of elastic waves in a many-layered elastic half-space is considered. One of the layers contains a bounded inhomogeneity which in the numerical applications is taken as a sphere. The half-space is excited by an explosion, modelled as a sudden isotropic point source.



First the time-harmonic problem is solved using the null field approach (the T-matrix method). Starting from surface integral representations containing the free space Green's dyadic and field expansions in plane and spherical vector wave functions, the null field approach leads to a set of algebraic equations whose solution can be given in a form that naturally has a multiple scattering interpretation. The null field approach thus has a building-block structure, i.e. the transition matrix of the inhomogeneity and the reflection and transmission coefficients of the interfaces are used as parts in the total solution. It is noted that in the absence of the inhomogeneity the null field approach yields a more decoupled set of equations than the classical methods. The main practical limitation of the null field approach, which is, in principle, an exact method, is that only low and intermediate frequencies can be treated. In the numerical applications the highest frequency corresponds to a sphere diameter of about 12 S wavelengths. The transformations from the frequency domain to the time domain is performed with an FFT algorithm and to limit the bandwidth appropriate filters are placed at the receivers.



The numerical examples show synthetic seismograms consisting of data from 10 observation points at increasing distances from the source. The layers have been chosen relatively thick so that the reflections from the different interfaces can be separated in time. As long as the reflections from the spherical inhomogeneity are also separated in time from other reflections they can mostly be well recognized. As the distance between the source and inhomogeneity increases, it becomes progressively more difficult to see any influence of the inhomogeneity in the seismograms. (Less)
Please use this url to cite or link to this publication:
author
publishing date
type
Contribution to journal
publication status
published
subject
in
Geophysical Journal International
volume
89
issue
2
pages
527 - 547
publisher
Wiley-Blackwell
external identifiers
  • scopus:0023480480
ISSN
0956-540X
DOI
10.1111/j.1365-246X.1987.tb05182.x
language
English
LU publication?
no
id
663b1f70-b8eb-4925-b81a-4bc9cc205927 (old id 1038697)
date added to LUP
2008-02-27 14:43:35
date last changed
2017-08-20 04:36:35
@article{663b1f70-b8eb-4925-b81a-4bc9cc205927,
  abstract     = {The propagation of elastic waves in a many-layered elastic half-space is considered. One of the layers contains a bounded inhomogeneity which in the numerical applications is taken as a sphere. The half-space is excited by an explosion, modelled as a sudden isotropic point source.<br/><br>
<br/><br>
First the time-harmonic problem is solved using the null field approach (the T-matrix method). Starting from surface integral representations containing the free space Green's dyadic and field expansions in plane and spherical vector wave functions, the null field approach leads to a set of algebraic equations whose solution can be given in a form that naturally has a multiple scattering interpretation. The null field approach thus has a building-block structure, i.e. the transition matrix of the inhomogeneity and the reflection and transmission coefficients of the interfaces are used as parts in the total solution. It is noted that in the absence of the inhomogeneity the null field approach yields a more decoupled set of equations than the classical methods. The main practical limitation of the null field approach, which is, in principle, an exact method, is that only low and intermediate frequencies can be treated. In the numerical applications the highest frequency corresponds to a sphere diameter of about 12 S wavelengths. The transformations from the frequency domain to the time domain is performed with an FFT algorithm and to limit the bandwidth appropriate filters are placed at the receivers.<br/><br>
<br/><br>
The numerical examples show synthetic seismograms consisting of data from 10 observation points at increasing distances from the source. The layers have been chosen relatively thick so that the reflections from the different interfaces can be separated in time. As long as the reflections from the spherical inhomogeneity are also separated in time from other reflections they can mostly be well recognized. As the distance between the source and inhomogeneity increases, it becomes progressively more difficult to see any influence of the inhomogeneity in the seismograms.},
  author       = {Karlsson, Anders and Boström, Anders},
  issn         = {0956-540X},
  language     = {eng},
  number       = {2},
  pages        = {527--547},
  publisher    = {Wiley-Blackwell},
  series       = {Geophysical Journal International},
  title        = {Broad band synthetic seismograms for a spherical inhomogeneity in a many-layered elastic half-space},
  url          = {http://dx.doi.org/10.1111/j.1365-246X.1987.tb05182.x},
  volume       = {89},
  year         = {1987},
}