The earthquake slip distribution self-similarity is investigated in this study. We complied finite fault slip models for earthquakes in the Taiwan orogenic belt and global earthquakes to determine the slip distribution self-similarity. Forty-one earthquakes (19 Taiwan earthquakes and 22 global earthquakes) in the M_w = 4.6 - 8.9 magnitude range were examined. The fault slip exhibited self-similar scaling between the rupture slip and area. The average area ratio (R_s) and slip ratio (R_d) follows a scaling of R_s = 10^{a - n(Rd)}. Slip self-similarity implies that a fault rupture exhibits fractal behavior. The scaling exponent can be considered as a measure for the roughness degree of the slip distribution on the fault surface. This study suggests that the slip distribution for large earthquakes (M_w > 7.0) tends to have a more homogeneous slip. Scaling exponents can provide insight into earthquake rupture mechanics and the scaling of heterogeneous slips on the fault surface provides a basis for ground mo­tion simulation for a finite fault for an earthquake scenario, particularly for near-fault motion.