**7. References**


© 2012 Qian, licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use,

© 2012 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution,

distribution, and reproduction in any medium, provided the original work is properly cited.

and reproduction in any medium, provided the original work is properly cited.

**An Efficient Approach for** 

Additional information is available at the end of the chapter

Zheng-Hua Qian

**1. Introduction** 

http://dx.doi.org/10.5772/47806

**Seismoacoustic Scattering Simulation** 

**Based on Boundary Element Method** 

amplification decreases as the amplitude of the incident wave increases.

A seismogram synthesis is the theoretically calculated ground motion that would be recorded for a given crust structure and seismic source, which is an important tool that can be applied to the study of earthquake source from strong motion records. At the early stage of seismology, the reflectivity method [1] and the generalized ray method [2] were widely used for simulating seismic wave excitation and propagation, where a laterally homogeneous model is considered. The real Earth's crust, however, is a laterally inhomogeneous structure indeed. To study the site effect of strong ground motion, a stratified structure with irregular interfaces is an eligible model for simulating seismic wave propagation. For the study of seismic waves propagating through a sediment filled basin in the case of rigid grains, the work [3-4] is noteworthy, where a model was proposed to reproduct a nonlinear effect experimentally observed for real seismic waves: site

In the last three decades, the study of seismic wave excitation and propagation in laterally inhomogeneous media has become extensively important in seismology and geophysics, including both analytical solutions and non-analytical results. The analytical solutions, however, only exist for a few special cases [5], such as the semi-cylindrical canyon and alluvial valleys, semi-elliptical canyon, and hemispherical canyon. The methods based on high-frequency approximation [6](such as asymptotic ray theory, Gaussian beam method, and so on) can give a proper solution by including the contribution from the neighboring rays but are only appropriate for handling body wave propagation under relatively high frequency. The methods based on plane wave decomposition (such as discrete wave number method [7-10]) can provide a complete wave field for any finite frequency but usually involve the truncation of matrices and vectors having an infinite dimension. Nowadays, the
