Numerical simulation of three-dimensional wave equation requires large computer memory and high calculation speed.Parallel algorithm based on Message Passing Interface(MPI) was adopted to simulate Rayleigh wave by using staggered-grid finite difference method.The model region to be simulated was divided into several sub-regions.Then all processors collaborate each other and the simulation in sub-regions is completed by each processor in parallel.Accordingly,the purposes of accelerate the computing speed and expand the scale of the model are achieved.Acoustic-elastic boundary approach was adopted to implement the traction-free boundary in the simulation.Feasibility and correctness of the simulation method were certificated by comparing the single-channel seismic records from the heterogeneous isotropic medium model obtained by numerical simulation with the results of analytical solution and the snapshots.The method is further verified by the comparison of dispersion curves derived from wave records and analytical solution of a three-layer velocity incremented model.