Please use this identifier to cite or link to this item: https://hdl.handle.net/10356/102531
Title: Frozen gaussian approximation for 3-D seismic wave propagation
Authors: Chai, Lihui
Tong, Ping
Yang, Xu
Keywords: DRNTU::Science::Mathematics
Fourier Analysis
Numerical Approximations And Analysis
Issue Date: 2016
Source: Chai, L., Tong, P., & Yang, X. (2017). Frozen Gaussian approximation for 3-D seismic wave propagation. Geophysical Journal International, 208(1), 59-74. doi:10.1093/gji/ggw368
Series/Report no.: Geophysical Journal International
Abstract: We present a systematic introduction on applying frozen Gaussian approximation (FGA) to compute synthetic seismograms in 3-D earth models. In this method, seismic wavefield is decomposed into frozen (fixed-width) Gaussian functions, which propagate along ray paths. Rather than the coherent state solution to the wave equation, this method is rigorously derived by asymptotic expansion on phase plane, with analysis of its accuracy determined by the ratio of short wavelength over large domain size. Similar to other ray-based beam methods (e.g. Gaussian beam methods), one can use relatively small number of Gaussians to get accurate approximations of high-frequency wavefield. The algorithm is embarrassingly parallel, which can drastically speed up the computation with a multicore-processor computer station. We illustrate the accuracy and efficiency of the method by comparing it to the spectral element method for a 3-D seismic wave propagation in homogeneous media, where one has the analytical solution as a benchmark. As another proof of methodology, simulations of high-frequency seismic wave propagation in heterogeneous media are performed for 3-D waveguide model and smoothed Marmousi model, respectively. The second contribution of this paper is that, we incorporate the Snell's law into the FGA formulation, and asymptotically derive reflection, transmission and free surface conditions for FGA to compute high-frequency seismic wave propagation in high contrast media. We numerically test these conditions by computing traveltime kernels of different phases in the 3-D crust-over-mantle model.
URI: https://hdl.handle.net/10356/102531
http://hdl.handle.net/10220/47267
ISSN: 0956-540X
DOI: http://dx.doi.org/10.1093/gji/ggw368
Rights: © 2016 The Authors. Published by Oxford University Press on behalf of The Royal Astronomical Society. This paper was published in Geophysical Journal International and is made available as an electronic reprint (preprint) with permission of The Authors 2016. Published by Oxford University Press on behalf of The Royal Astronomical Society. The published version is available at: [http://dx.doi.org/10.1093/gji/ggw368]. One print or electronic copy may be made for personal use only. Systematic or multiple reproduction, distribution to multiple locations via electronic or other means, duplication of any material in this paper for a fee or for commercial purposes, or modification of the content of the paper is prohibited and is subject to penalties under law.
Fulltext Permission: open
Fulltext Availability: With Fulltext
Appears in Collections:SPMS Journal Articles

Files in This Item:
File Description SizeFormat 
Frozen Gaussian approximation for 3-D seismic wave propagation.pdf5.06 MBAdobe PDFThumbnail
View/Open

Google ScholarTM

Check

Altmetric

Items in DR-NTU are protected by copyright, with all rights reserved, unless otherwise indicated.