A three-dimensional indirect boundary integral equation method for the scattering of seismic waves in a poroelastic layered half-space

A new indirect boundary integral equation method (IBIEM) is developed to solve the seismic wave scattering problem in a three-dimensional fluid-saturated layered half-space. Based on Biot's theory, the Green's functions of inclined circular loads in porous elastic layered half-space are firstly deduced. According to the single-layer potential theory, when establishing the scattered wave field, the uniform surface loads and fluid sources are located directly on the boundary surfaces of the scatterers, so as to avoid determining the optimal location of fictitious wave sources surface. At the same time, the radiation condition of wave in semi-infinite layered media can be accurately realized by using the dynamic Green's function of concentrated load, and the computational memory can be greatly reduced. The scattering of seismic waves by a canyon topography in saturated layered half-space is examined in detail. The numerical results indicate that with the increase of the porosity of the overlying soil, the amplitude of surface displacement can be amplified by 1∼6 times, and the amplification effect is more significant near the corner of the canyon, which can be attributed to the superposition of wave diffraction effect and resonance amplification effect of saturated soil layer.