A numerical method is proposed to study the scattering of seismic shear waves induced by the presence of underground cavities in homogeneous soils. The method is based on the superposition of two solutions: the solution of the free-wave propagation problem in a uniform half-space, easily determined analytically, and the solution of the wave scattering problem due to the cave presence, evaluated numerically by means of an ad hoc code implemented by using the ANSYS Parametric Design Language. In the two-dimensional setting, this technique is applied to the case of a single cave, placed at a certain depth from the ground level. The frequency spectrum of the seismic shear oscillation on the ground surface is determined for different dimensions and depths of the cave and compared with the spectrum registered without caves. The influence of the cave dimensions and depth on the spectrum amplification is analyzed and discussed.