The feature extraction of a fracture surface is inaccurate due to the lack of local fragments and damage of the geometric features of ornamentation.Thus, in this paper, a method was proposed to solve this problem based on SURF feature descriptor and Jaccard distance. First, Canny operator was used to extract the contour lines of the debris edges and model surface, and a multi-scale space was constructed to extract the feature points of the fault surface. Second, to solve the high redundancy in constructing feature descriptors and high delay in Euclidean distance calculation, SURF feature descriptors with low redundancy was constructed, and then compared the similarity of feature points with Jaccard distance to determine the optimal adjacency relationship of fragments. Finally, the parameters of rigid body were calculated by ICP method, and the fragments were assembled accurately. The experimental results showed that the running time of the algorithm increased by 12%-16%, and the splicing error was at most 0.750 mm. Compared to the traditional method, the proposed method can effectively reduce the considerably large splicing gap and penetration caused by the damage of the fragment model,as well as the splicing error, and achieve fragment splicing efficiently.