A reasonable three-dimensional complex fault model was established by comprehensively considering the three-dimensional seismic structure model, aftershock distribution, and surface rupture investigation. Two possible rupture modes were adopted, and based on parallel non negative least squares method and multi time window technology, combined with far-eld, near-eld, GPS, and surface rupture data, the precise rupture process of the Wenchuan earthquake was inverted, and a reasonable rupture mode was provided. The results are as follows: (1) The occurrence of bilateral rupture at the intersection of the Beichuan fault and the Xiaoyudong fault is more in line with the actual situation for the Wenchuan earthquake. When the bilateral rupture occurs, dislocations that match the observed surface rupture are generated in the near surface area from Hongkou to Yingxiu. The synthetic records amplitude of the station in the near eld behind the rupture direction is closer to the observation record.(2) The joint inversion results show that the duration of the Wenchuan earthquake rupture reached 100s, and the released seismic moment was 1.058×10 21 N•m. There are 5 asperities on the fault plane, indicating that this earthquake is composed of at least 5 subevents. The slip is mainly distributed on the Beichuan Fault, indicating that the Beichuan Fault is the main rupture fault.