Abstract: The randomness of artificial seismic vibration has important influences on the peak value of ground acceleration. Based on the results of soil layer shear wave velocity and soil body nonlinear statistics, three soil layer seismic response models are established. Taking the typical bedrock acceleration response spectrum of Heilongjiang region as the target one, 400 random seismic acceleration intervals are synthesized manually by 400 different random numbers for each target spectrum. Through the experience of soil layer seismic response, 10800 peak ground acceleration results are obtained from the three models. The distribution characteristics and differences of acceleration peaks under each target spectrum are analyzed statistically, and the reasons for this difference are explained preliminarily by comparing the Fourier spectrum with target response one. The results show that under the same profile and target response spectrum, there is a certain difference in the peak value of ground surface acceleration obtained by different random numbers, which shows a normal distribution. This difference is positively correlated with the response spectral value of the seismic targets corresponding to the Fourier spectrum peak periods of the surface earthquakes.