Abstract: When seismic waves pass through the underground scatterers, multi-order scattered waves will be generated, and the analysis of their seismic response characteristics is beneficial to infer the distribution and properties of the scatterers. Based on the two-dimensional scalar wave motion equation, combined with the seismic scattering theory and the Born approximation theory, the multi-order scattering wave equation is deduced. The numerical simulation of the two-point scatterer model and the complex scatterer model is carried out by using the high-order finite difference method, and the propagation law and wavefield characteristics of the multi-order scattered wave are analyzed. And then the single track records of multi-order scattering and each scattering are extracted and compared with the reference single track records, which verifies the accuracy of the scattering wave equation derived in this paper.