Commun. Comput. Phys.,
Notice: Undefined index: year in /var/www/html/issue/abstract/readabs.php on line 20

Notice: Undefined index: ppage in /var/www/html/issue/abstract/readabs.php on line 21

Notice: Undefined index: issue in /var/www/html/issue/abstract/readabs.php on line 23
Volume 3.


Boundary Integral Modelling of Elastic Wave Propagation in Multi-Layered 2D Media with Irregular Interfaces

Enru Liu 1*, Zhongjie Zhang 2, Jianghua Yue 3, Andy Dobson 4

1 Department of Geophysics, China University of Mining and Technology, Xuzhou, 221008, China
2 Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China.
3 Department of Geophysics, China University of Mining and Technology, Xuzhou, 221008, China.
4 British Geological Survey, West Mains Road, Edinburgh EH9 3LA, Scotland, UK.

Received 7 August 2006; Accepted (in revised version) 6 June 2007
Available online 14 September 2007

Abstract

We present a semi-analytic method based on the propagation matrix formulation of indirect boundary element method to compute response of elastic (and acoustic) waves in multi-layered media with irregular interfaces. The method works recursively starting from the top-most free surface at which a stress-free boundary condition is applied, and the displacement-stress boundary conditions are then subsequently applied at each interface. The basic idea behind this method is the matrix formulation of the propagation matrix (PM) or more recently the reflectivity method as wide used in the geophysics community for the computation of synthetic seismograms in stratified media. The reflected and transmitted wave-fields between arbitrary shapes of layers can be computed using the indirect boundary element method (BEM, sometimes called IBEM). Like any standard BEM, the primary task of the BEM-based propagation matrix method (thereafter called PM-BEM) is the evaluation of element boundary integral of the Green's function, for which there are standard method that can be adapted. In addition, effective absorbing boundary conditions as used in the finite difference numerical method is adapted in our implementation to suppress the spurious arrivals from the artificial boundaries due to limited model space. To our knowledge, such implementation has not appeared in the literature. We present several examples in this paper to demonstrate the effectiveness of this proposed PM-BEM for modelling elastic waves in media with complex structure.

AMS subject classifications: 65N38, 86-08, 86A15

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Boundary element method, wave propagation, seismology, scattering.

*Corresponding author.
Email: enru.liu@exxonmobil.com (E. Liu), zhangzj@mail.iggcas.ac.cn (Z. Zhang), yuejh@cumt.edu.cn (J. Yue)
 

The Global Science Journal