The solution of scattering of SH-wave was given by the method which to use an circular boundary of large radius to approximate the straight boundary of surface layer. With the theory of Helmholtz, the general solution of the Biot’s wave function was achieved. Utilizing the complex function method and the boundary conditions, we could transform the present problem into the problem in which we needed to solve the infinite linear algebraic equations with unknown coefficients. Numerical results of the ground motion were provided. The effects of incident wave number, incident angle, buried depth, stiffness of inclusion and layer on ground motion were analyzed qualitatively. The results in numerical examples show that greatly impact on ground motion can occur due to strong interaction between the surface layer, circular cavity and inclusion.