HAL will be down for maintenance from Friday, June 10 at 4pm through Monday, June 13 at 9am. More information
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 : Monday, March 14, 2022 - 3:31:01 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

42

Files downloads

3