摘要

Based on the integral expression of half-analytical solution for elastic waves in three-dimensional (3D) stratified porous media derived in the companion paper, we carry out numerical simulations by using the discrete wave-number method in this study. Comparing the analytical solutions of the elastic waves generated by a point source of single-force or explosion in homogeneous full-space porous media with our numerical solutions, we find they are identical. We also confirm the special reciprocity of the solid displacement by the numerical simulation on a two-layer half-space model. Thus, we verify the validity and reliability of our half-analytical solution, the numerical simulation method and the program code used by us. Then we carry out sensitivity study to discuss the impact of the media parameter on the reflection waves induced by an explosion point source on the interface. In the numerical simulation of a model of vertical seismic profile, we find that solid displacement, seepage displacement and pore-fluid pressure all can reflect the variation of the media properties. Meanwhile, we also discuss the dynamic compatibility phenomenon in porous media. It can be concluded that the developed numerical simulation method based on our half-analytical solution provides an alternative and effective instrument for the study of elastic waves' propagation in 3D stratified porous media.