The authors present a series of rupture models of the Wenchuan M_W7.9 earthquake. A simulated annealing algorithm is used to determine the finite-fault model that minimizes the objective function described in terms of wavelet coefficients. With the approach, the authors can simultaneously invert for the slip amplitude, slip direction, rise time and rupture velocity. At the first step, experiments conducted on synthetic data are used to assess the ability to recover rupture slip details and teleseismic body waves. The authors investigate the resolvability of teleseismic body inversion on dip angle, rupture velocity and maximum rupture depth. Then 4 slip models are obtained for the Wenchuan earthquake by applying the method on 4 single fault planes which are specified by different dip angle. The authors analyze the resolvability of finite fault analysis by using telesesimic data and how to choose a reasonable dip angle for the Wenchuan earthquake. Two asperities are observed on different solutions, and likelihood of change dip is discussed.