In this study, the technique of stochastic finite- fault simulation method was used to simulate strong ground motion of Tehran earthquake. In this method, the respond spectrum on fault surface caused by each element is determined by the combination of an average spectrum of ω^2 and a unique Gowssy spectrum (Boor, 1983). In addition, to identify details of the source, path and the rupture spread of this earthquake, specified parameters were used to enter simulation software EXSIM. At first, acceleration coefficient of the main earthquake was calculated in different stations using GEOPSY software and spectral decay parameter was determined for this scenario. The fault model for this earthquake was reached through practical equations (Wells and Coopersmith, 1994) which shows a 127 Km causative rupture length and 29 Km causative rupture width, while strike and dip were estimated 263 and 75 degrees respectively. In the end, a comparison was made between this study and ground-motion prediction equation (GMPE) which was indicative of conformity in just some stations. This difference is justified by the tectonic conditions of equations and the limited nature of the data in this study.