Fast computation of seabed spherical-wave reflection coefficients in geoacoustic inversion

Jorge E. Quijano, Stan E. Dosso, Jan Dettmer, Charles Holland

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

This paper develops a fast numerical approach to computing spherical-wave reflection coefficients (SWRCs) for layered seabeds, which provides substantial savings in computation time when used as the forward model for geoacoustic inversion of broadband seabed reflectivity data. The approach exploits the Sommerfeld-integral representation of SWRCs as the Hankel transform of a function proportional to the plane-wave reflection coefficient (PWRC), and applies Levin integration to the rapidly oscillating integrand cast as the product of a (pre-computed) media-independent matrix and a vector involving PWRCs at a sparse sampling of integration angles. Compared to conventional Simpson's rule integration for computation of the SWRC, the Levin integration yields speed-up factors of an order of magnitude or more. Further, it results in reduced memory requirements for storage of pre-computed quantities, a desirable property when a graphics processing unit (GPU) is used for parallel computation of SWRCs. The paper applies trans-dimensional Bayesian inversion to investigate the impact of forward modeling in terms of PWRCs and SWRCs on the estimation of geoacoustic parameters and uncertainties. Model comparisons are quantified in simulated- and measured-data inversions by comparing the estimated geoacoustic parameters to the true parameters or core measurements, respectively, and by calculating the deviance information criterion for model selection.

Original languageEnglish (US)
Pages (from-to)2106-2117
Number of pages12
JournalJournal of the Acoustical Society of America
Volume138
Issue number4
DOIs
StatePublished - Oct 1 2015

Fingerprint

spherical waves
wave reflection
inversions
reflectance
Inversion
Waves
casts
plane waves
sampling
broadband
requirements
products
matrices

All Science Journal Classification (ASJC) codes

  • Arts and Humanities (miscellaneous)
  • Acoustics and Ultrasonics

Cite this

@article{0db61d6407524347ab1a97ffcf3f5a2e,
title = "Fast computation of seabed spherical-wave reflection coefficients in geoacoustic inversion",
abstract = "This paper develops a fast numerical approach to computing spherical-wave reflection coefficients (SWRCs) for layered seabeds, which provides substantial savings in computation time when used as the forward model for geoacoustic inversion of broadband seabed reflectivity data. The approach exploits the Sommerfeld-integral representation of SWRCs as the Hankel transform of a function proportional to the plane-wave reflection coefficient (PWRC), and applies Levin integration to the rapidly oscillating integrand cast as the product of a (pre-computed) media-independent matrix and a vector involving PWRCs at a sparse sampling of integration angles. Compared to conventional Simpson's rule integration for computation of the SWRC, the Levin integration yields speed-up factors of an order of magnitude or more. Further, it results in reduced memory requirements for storage of pre-computed quantities, a desirable property when a graphics processing unit (GPU) is used for parallel computation of SWRCs. The paper applies trans-dimensional Bayesian inversion to investigate the impact of forward modeling in terms of PWRCs and SWRCs on the estimation of geoacoustic parameters and uncertainties. Model comparisons are quantified in simulated- and measured-data inversions by comparing the estimated geoacoustic parameters to the true parameters or core measurements, respectively, and by calculating the deviance information criterion for model selection.",
author = "Quijano, {Jorge E.} and Dosso, {Stan E.} and Jan Dettmer and Charles Holland",
year = "2015",
month = "10",
day = "1",
doi = "10.1121/1.4930186",
language = "English (US)",
volume = "138",
pages = "2106--2117",
journal = "Journal of the Acoustical Society of America",
issn = "0001-4966",
publisher = "Acoustical Society of America",
number = "4",

}

Fast computation of seabed spherical-wave reflection coefficients in geoacoustic inversion. / Quijano, Jorge E.; Dosso, Stan E.; Dettmer, Jan; Holland, Charles.

In: Journal of the Acoustical Society of America, Vol. 138, No. 4, 01.10.2015, p. 2106-2117.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Fast computation of seabed spherical-wave reflection coefficients in geoacoustic inversion

AU - Quijano, Jorge E.

AU - Dosso, Stan E.

AU - Dettmer, Jan

AU - Holland, Charles

PY - 2015/10/1

Y1 - 2015/10/1

N2 - This paper develops a fast numerical approach to computing spherical-wave reflection coefficients (SWRCs) for layered seabeds, which provides substantial savings in computation time when used as the forward model for geoacoustic inversion of broadband seabed reflectivity data. The approach exploits the Sommerfeld-integral representation of SWRCs as the Hankel transform of a function proportional to the plane-wave reflection coefficient (PWRC), and applies Levin integration to the rapidly oscillating integrand cast as the product of a (pre-computed) media-independent matrix and a vector involving PWRCs at a sparse sampling of integration angles. Compared to conventional Simpson's rule integration for computation of the SWRC, the Levin integration yields speed-up factors of an order of magnitude or more. Further, it results in reduced memory requirements for storage of pre-computed quantities, a desirable property when a graphics processing unit (GPU) is used for parallel computation of SWRCs. The paper applies trans-dimensional Bayesian inversion to investigate the impact of forward modeling in terms of PWRCs and SWRCs on the estimation of geoacoustic parameters and uncertainties. Model comparisons are quantified in simulated- and measured-data inversions by comparing the estimated geoacoustic parameters to the true parameters or core measurements, respectively, and by calculating the deviance information criterion for model selection.

AB - This paper develops a fast numerical approach to computing spherical-wave reflection coefficients (SWRCs) for layered seabeds, which provides substantial savings in computation time when used as the forward model for geoacoustic inversion of broadband seabed reflectivity data. The approach exploits the Sommerfeld-integral representation of SWRCs as the Hankel transform of a function proportional to the plane-wave reflection coefficient (PWRC), and applies Levin integration to the rapidly oscillating integrand cast as the product of a (pre-computed) media-independent matrix and a vector involving PWRCs at a sparse sampling of integration angles. Compared to conventional Simpson's rule integration for computation of the SWRC, the Levin integration yields speed-up factors of an order of magnitude or more. Further, it results in reduced memory requirements for storage of pre-computed quantities, a desirable property when a graphics processing unit (GPU) is used for parallel computation of SWRCs. The paper applies trans-dimensional Bayesian inversion to investigate the impact of forward modeling in terms of PWRCs and SWRCs on the estimation of geoacoustic parameters and uncertainties. Model comparisons are quantified in simulated- and measured-data inversions by comparing the estimated geoacoustic parameters to the true parameters or core measurements, respectively, and by calculating the deviance information criterion for model selection.

UR - http://www.scopus.com/inward/record.url?scp=84945232194&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84945232194&partnerID=8YFLogxK

U2 - 10.1121/1.4930186

DO - 10.1121/1.4930186

M3 - Article

AN - SCOPUS:84945232194

VL - 138

SP - 2106

EP - 2117

JO - Journal of the Acoustical Society of America

JF - Journal of the Acoustical Society of America

SN - 0001-4966

IS - 4

ER -