Local inclusion topography has significant influence on seismic wave propagation, and the propagation characteristics of seismic waves in poroelastic soils are obviously different from those in single-phase media. Based on Biot’s theory, the scattering of plane P1 wave by inclusion in a three-dimensional poroelastic half-space is studied by using the indirect boundary element method (IBEM). The scattering field is constructed by introducing a virtual wave near the interface between inclusion and half-space and the surface of half-space, and the virtual wave density is obtained by establishing boundary integral equation based on the boundary conditions. The effects of the depth, geometric characteristics, boundary permeability, porosity, incident frequency, and incident angle of the inclusion on elastic wave scattering are systematically analyzed. The results show that due to the soil skeleton-pore water coupling effect, when the porosity is n = 0.3, the surface displacement amplitude of dry soil is larger than that of poroelastic soil. When the porosity is n = 0.36, the surface displacement amplitude of poroelastic soil is larger than that of dry soil. The surface displacement amplitude of poroelastic-drained condition is slightly larger than that of undrained condition. With the increase of inclusion depth, the scattering of elastic wave by inclusion decreases gradually. When P1 wave is incident, the surface displacement amplitude at the depth of H = 0.5 can be increased up to three times as much as that at the depth of H = 1.5. As the inclusion becomes narrower and flatter, the scattering of elastic waves by inclusion decreases gradually. When the ratio between height and length is S = 2/5, the surface displacement magnitude can reach up to 9.5.