We studied the 20 December 2010 Rigan earthquake that occurred in southeastern Iran, using the ground‐motion simulation based on the stochastic finite‐fault model. This earthquake was recorded by 23 stations. The method requires region‐specific attenuation and generic site parameters. We first found an empirical function to describe the amplitude decay with distance of shear waves in the Rigan region of Iran. This function results from the analysis of spectral amplitudes of horizontal components from digital records from the Rigan earthquakes. The value of QS shows a dependence on frequency in the range of 1–24 Hz for the Rigan region. The average frequency‐dependent relationship estimated for the region is QS=(99±19)f(0.58±0.12). We then obtained a spectral decay parameter from the slope of smoothed amplitude of the Fourier spectra of acceleration at higher frequencies. The estimated zero‐distance Kappa value is 0.06. In the next step, we used these approximate values and the stochastic finite‐fault model in an iterative manner to optimally model the strong ground motion and rupture characteristics of the main event in terms of peak ground acceleration, peak ground velocity, peak ground displacement, Arias intensity, the amplitude Fourier spectra, and the response spectra of ground shaking. We compared our simulated results with recorded ones in both the frequency and time domain. Our near‐field estimates for the main event is in good conformity with the corresponding near‐field estimates reported by other studies.