In this paper, we extend the velocity-stress acoustoporoelastic formulation from 2D to 3D, and Padé acoustoelasticity to Padé acoustoporoelasticity. Applications to experimental data with porous rocks differentiate the acoustoelasticity, acoustoporoelasticity, and Padé acoustoporoelasticity in accuracy. The SSG-FD numerical method was used to simulate the propagation of elastic waves in both 3D acoustoporoelastic and Padé acoustoporoelastic media. By comparing the theoretical and calculated wave velocities, the presented numerical scheme is verified using plane-wave theoretical solutions.Â