Skip to Main content Skip to Navigation
Journal articles

Discrete wavenumber solutions to numerical wave propagation in piecewise heterogeneous media - I. Theory of two-dimensional SH case

Abstract : A semi-analytical, semi-numerical method of seismogram synthesis is presented for piecewise heterogeneous media resulting from an arbitrary source. The method incorporates the discrete wavenumber Green's function representation into the boundary-volume integral equation numerical techniques. The presentation is restricted to 2-D antiplane motion (SH waves). To model different parts of the media to a necessary accuracy, the incident, boundary-scattering and volume-scattering waves are separately formulated in the discrete wavenumber domain and handled flexibly at various accuracies using approximation methods. These waves are accurately superposed through the generalized Lippmann-Schwinger integral (GLSI) equation. The full-waveform boundary method is used for the boundary-scattering wave to accurately simulate the reflection/transmission across strong-contrast boundaries. Meanwhile for volume heterogeneities, the following four flexible approaches have been developed in the numerical modelling scheme present here, with a great saving of computing time and memory: the solution implicitly for the volume-scattering wave with high accuracy to model subtle effects of volume heterogeneities; the solution semi-explicitly for the volume-scattering wave using the average Fresnel-radius approximation to volume integrations to reduce numerical burden by making the coefficient matrix sparser; the solution explicitly for the volume-scattering wave using the first-order Born approximation for smooth volume heterogeneities; and the solution explicitly for the volume-scattering wave using the second-order/high-order Born approximation for practical volume heterogeneities. These solutions are tested for dimensionless frequency responses to a heterogeneous alluvial valley where the velocity is perturbed randomly in the range of ca 5-20 per cent, which is not rare in most complex near-surface areas. Numerical experiments indicate that several times of site amplification can be expected as a result of heterogeneities introduced in a homogeneous valley. The test also confirms that the first-order Born approximation to the volume-scattering wave is strictly valid for velocity perturbation less than 10 per cent and approximately used for up to 15 per cent for general applications. The second-order Born approximation to the volume-scattering wave is strictly valid for velocity perturbation less than 15 per cent and approximately used for up to 20 per cent for general applications.
Document type :
Journal articles
Complete list of metadata

https://hal-insu.archives-ouvertes.fr/insu-03607131
Contributor : Nathalie POTHIER Connect in order to contact the contributor
Submitted on : Sunday, March 13, 2022 - 10:48:39 AM
Last modification on : Friday, September 30, 2022 - 4:07:24 AM

File

157-2-481.pdf
Publisher files allowed on an open archive

Licence


Distributed under a Creative Commons Attribution 4.0 International License

Identifiers

Collections

Citation

Li-Yun Fu, Michel Bouchon. Discrete wavenumber solutions to numerical wave propagation in piecewise heterogeneous media - I. Theory of two-dimensional SH case. Geophysical Journal International, Oxford University Press (OUP), 2004, 157, pp.481-498. ⟨10.1111/j.1365-246X.2004.02135.x⟩. ⟨insu-03607131⟩

Share

Metrics

Record views

45

Files downloads

6