The effects of topography on seismic ground motion of Wenchuan earthquake were simulated based upon elastodynamic equation for 3D wave propagation, the spectral-element method, ASTER DEM model and parallel computation environment and incorporated realistic 3D topography and attenuation of media. The simulated results show that topographic effects due to great rough terrain on two horizontal components of seismic wave are bigger than those on vertical component. The max PGA value in zone near epicentre is 671 cm/s2, while large PGA values are almost located near mountain tops and ridges. It can be concluded that the effects of 3D realistic topography on seismic ground motion are more complex than those of the existing 2D generalized model.