3D dynamic Green’s functions in a multilayered poroelastic half-space

The complete 3D dynamic Green’s functions in the multilayered poroelastic media are presented in this study. A method of potentials in cylindrical coordinate system is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to 3D wave propagation problems are obtained. After that, a three vector base and the propagator matrix method are introduced to treat 3D wave propagation problems in the stratified poroelastic half-space disturbed by buried sources. It is known that the original propagator algorithm has the loss-of-precision problem when the waves become evanescent. At present, an orthogonalization procedure is inserted into the matrix propagation loop to avoid the numerical difficulty of the original propagator algorithm. At last, the validity of the present approach for accurate and efficient calculating 3D dynamic Green’s functions of a multilayered poroelastic half-space is confirmed by comparing the numerical results with the known exact analytical solutions of a uniform poroelastic half-space.