Dynamic response of a multi-scale layered saturated porous half-space due to seismic dislocation source by using a revised dynamic stiffness matrix method

Abstract

The study investigating the dynamic response of seismic dislocation source in a layered saturated porous half-space, particularly the multi-scale characteristic of half-space, from crustal to geotechnical scale, is still lacking in the literature. particularly in will lead to the computation bottlenecks for multi-scale seismic wave propagation. Based on Biot's theory of wave propagation in a fluid-saturated porous solid, a dynamic stiffness matrix method is proposed to investigate the dynamic response of a multi-scale layered saturated porous half-space due to the seismic dislocation source. Assembling each layer stiffness matrix in the frequency-wavenumber domain, the global dynamic stiffness matrix is constructed to avoid errors accumulative compared to the reflection and transmission matrix method. The divergent index terms of the stiffness matrix with respect to thickness and wavenumber are extracted to revise the dynamic stiffness matrix adapted for multi-scale half-space. The seismic dislocation source action is equal to the external forces imposed on the total saturated porous half-space, which can obtain the dynamic response by using revised dynamic stiffness matrix method. The accuracy of the proposed method is verified via the comparison between the calculated results and of three published literatures. The multi-scale model with realistic near-surface fine structure is selected to investigate the effects of superficial fine structures on the dynamic responses. The numerical simulations results show that the presence of low-velocity saturated porous layers has a prominent impact on ground motion. The dynamic response seems to be sensitive to the porosity conditions, especially for high frequency, which cause the difference in variation of peak displacement value.