Abstract: Due to the tectonic action, the interface of layered rock may be partially folded. The investigation of the influence of irregular interface on the scattering of cavities is of great significance for the seismic safety evaluation of the surface structure. Based on a substructure method, a complex scattering problem is transformed into a radiation problem and the evaluation of wave response of a layered half-space with regular boundary conditions (free field). The Fourier transform and dual variables are used to obtain the first-order ordinary differential wave equation. Besides, soil layers are merged by using the precise integration algorithm and load boundary conditions are applied to obtain the Green's function, and the dynamic stiffness is obtained. Based on a precise integration algorithm, the coefficient matrix of an improved transfer matrix method is proposed, and the wave motion of the layered half-space is calculated. The improved transfer matrix method has no restrictions on soil layer thickness and the number of layers. The accuracy and effectiveness of the proposed method are validated by comparing them with those of the results of the previous study. The scattering effects of horseshoe holes embedded in a complex layered half-space are investigated. The results show that the effect of local folds on the magnitude of surface displacement is related to factors such as incident wave type, incident wave frequency, and local fold geometry. The peak surface displacement is affected by the combined action of horseshoe-shaped cavity and local folds, and its characteristics are not obviously related to the type of incident wave.