Abstract: On May 28, 2012, an M4.8 earthquake occurred at the border between Guye District and Luanxian County in Tangshan city. With the available three-dimensional velocity model and topography data, this paper applies curve grid finite difference method which considering the terrain effect to the earthquake simulation of ground motion, get a few of snapshots of the image speed V_X. It can be seen that the earthquake spread outward in the form of approximate circular, energy of the earthquake is uniformly distributed, that is in accordance with point-source seismic characteristics. By comparing calculated waveform with station recorded waveform, it is found that the calculated waveform and the bedrock area stations recorded waveform can be well compared, the beginning of P wave and S wave, phase and maximum amplitude of 4.8 earthquake waveform have good consistency, which proves the validity of the simulation results. The result of calculated waveform comparing with the recorded waveform which stations on the Quaternary strata coverage area and the transition region is not well, the beginning of P wave and S wave can correspond, but maximum amplitude and phase have larger deviation. After the arrival of S waves, the energy attenuation of the calculated wave is faster than that of stations recorded waveform, which may be caused by lack of Quaternary thickness data in P wave velocity structure model of the simulation calculation.