According to the Biot wave theory and the complex function combined with multi coordinate method, seismic wave dynamic response of underground holes in the half space saturated soil is studied. Given P wave on two circular holes of the dynamic stress concentration factor of the numerical results, the distribution situations of the dynamic stress concentration factors are discussed with the variation of the hole spacing, incidence angle and the dimensionless wave number. The results show that, the effect of the underground hole spacing to the interaction effects of the holes is obvious. When the hole spacing reaches up to 6 times the radius of the hole, the peak value of the dynamic stress concentration factor decreases obviously, and the peripheral distribution of offset is weakened. With increasing of the incident angle, the hole peripheral stress concentration coefficient showed the trend of first increasing then decreasing, distribution also occurs at certain angle deflection. In low frequency waves, the hole stress peak value will influence mutually more, while in the high frequency input, interaction between holes will reduce.